
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:52 08/01/2011

Input data from file "261-aggressive-phylip/7611_hg18.phylip" (29 taxa, 442 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


(Output sorted by AICc. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	64	-1362.2723	2874.6137	0.6667	0.6667
GTRIG	65	-1362.1726	2877.1642	0.1863	0.8530
GTRI	64	-1363.8574	2877.7838	0.1366	0.9896
HKYG	60	-1372.3132	2883.8391	0.0066	0.9962
HKYIG	61	-1372.2024	2886.3100	0.0019	0.9982
HKYI	60	-1373.6890	2886.5906	0.0017	0.9998
GTR	63	-1371.9715	2891.2763	0.0002	1.0000
HKY	59	-1382.7189	2901.9719	0.0000	1.0000
SYMG	61	-1391.4600	2924.8252	0.0000	1.0000
K2PG	57	-1397.4189	2926.0565	0.0000	1.0000
SYMIG	62	-1391.3147	2927.2415	0.0000	1.0000
SYMI	61	-1393.1865	2928.2782	0.0000	1.0000
K2PIG	58	-1397.2629	2928.3953	0.0000	1.0000
F81G	59	-1396.3014	2929.1368	0.0000	1.0000
K2PI	57	-1399.1484	2929.5156	0.0000	1.0000
F81IG	60	-1396.1884	2931.5895	0.0000	1.0000
F81I	59	-1397.9801	2932.4943	0.0000	1.0000
SYM	60	-1402.2541	2943.7209	0.0000	1.0000
K2P	56	-1409.1556	2946.8931	0.0000	1.0000
F81	58	-1406.9699	2947.8094	0.0000	1.0000
JC69G	56	-1418.8801	2966.3421	0.0000	1.0000
JC69IG	57	-1418.7408	2968.7003	0.0000	1.0000
JC69I	56	-1420.5901	2969.7620	0.0000	1.0000
JC69	55	-1430.3652	2986.6889	0.0000	1.0000



(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1362.2723	2852.5447	0.6190	0.6190
GTRIG	65	-1362.1726	2854.3451	0.2516	0.8706
GTRI	64	-1363.8574	2855.7148	0.1269	0.9974
HKYG	60	-1372.3132	2864.6265	0.0015	0.9989
HKYIG	61	-1372.2024	2866.4047	0.0006	0.9995
HKYI	60	-1373.6890	2867.3780	0.0004	0.9999
GTR	63	-1371.9715	2869.9430	0.0001	1.0000
HKY	59	-1382.7189	2883.4379	0.0000	1.0000
SYMG	61	-1391.4600	2904.9200	0.0000	1.0000
SYMIG	62	-1391.3147	2906.6293	0.0000	1.0000
SYMI	61	-1393.1865	2908.3730	0.0000	1.0000
K2PG	57	-1397.4189	2908.8377	0.0000	1.0000
K2PIG	58	-1397.2629	2910.5259	0.0000	1.0000
F81G	59	-1396.3014	2910.6028	0.0000	1.0000
K2PI	57	-1399.1484	2912.2969	0.0000	1.0000
F81IG	60	-1396.1884	2912.3769	0.0000	1.0000
F81I	59	-1397.9801	2913.9602	0.0000	1.0000
SYM	60	-1402.2541	2924.5083	0.0000	1.0000
F81	58	-1406.9699	2929.9399	0.0000	1.0000
K2P	56	-1409.1556	2930.3113	0.0000	1.0000
JC69G	56	-1418.8801	2949.7603	0.0000	1.0000
JC69IG	57	-1418.7408	2951.4815	0.0000	1.0000
JC69I	56	-1420.5901	2953.1801	0.0000	1.0000
JC69	55	-1430.3652	2970.7304	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1372.3132	3110.1051	0.6876	0.6876
HKYI	60	-1373.6890	3112.8566	0.1737	0.8614
GTRG	64	-1362.2723	3114.3885	0.0808	0.9421
HKYIG	61	-1372.2024	3115.9746	0.0365	0.9787
GTRI	64	-1363.8574	3117.5586	0.0166	0.9952
GTRIG	65	-1362.1726	3120.2802	0.0042	0.9995
HKY	59	-1382.7189	3124.8252	0.0004	0.9999
GTR	63	-1371.9715	3127.6955	0.0001	1.0000
K2PG	57	-1397.4189	3142.0424	0.0000	1.0000
K2PI	57	-1399.1484	3145.5015	0.0000	1.0000
K2PIG	58	-1397.2629	3147.8218	0.0000	1.0000
F81G	59	-1396.3014	3151.9901	0.0000	1.0000
SYMG	61	-1391.4600	3154.4899	0.0000	1.0000
F81I	59	-1397.9801	3155.3475	0.0000	1.0000
F81IG	60	-1396.1884	3157.8555	0.0000	1.0000
SYMI	61	-1393.1865	3157.9429	0.0000	1.0000
K2P	56	-1409.1556	3159.4246	0.0000	1.0000
SYMIG	62	-1391.3147	3160.2906	0.0000	1.0000
F81	58	-1406.9699	3167.2359	0.0000	1.0000
SYM	60	-1402.2541	3169.9869	0.0000	1.0000
JC69G	56	-1418.8801	3178.8736	0.0000	1.0000
JC69I	56	-1420.5901	3182.2935	0.0000	1.0000
JC69IG	57	-1418.7408	3184.6862	0.0000	1.0000
JC69	55	-1430.3652	3195.7524	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7611_hg18.phylip = 1 - 442;
 Partition Dummy = 1:261-aggressive-phylip/7611_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7611_hg18.phylip = 1 - 442;
 Partition Dummy = 1:261-aggressive-phylip/7611_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7611_hg18.phylip = 1 - 442;
 Partition Dummy = 1:261-aggressive-phylip/7611_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

