
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:46 02/17/2011

Input data from file "183-phylip/chr2_8655.phylip" (29 taxa, 491 characters)

Minimum AIC  model: HKYIG
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	60	-1174.7714	2486.5660	0.3987	0.3987
HKYIG	61	-1173.7022	2487.0360	0.3152	0.7140
HKYG	60	-1175.2526	2487.5284	0.2464	0.9604
GTRI	64	-1172.7043	2492.9390	0.0165	0.9769
GTRIG	65	-1171.6681	2493.5244	0.0123	0.9892
GTRG	64	-1173.1251	2493.7807	0.0108	1.0000
K2PI	57	-1190.6694	2510.6090	0.0000	1.0000
K2PIG	58	-1189.5691	2510.9807	0.0000	1.0000
K2PG	57	-1191.2571	2511.7844	0.0000	1.0000
SYMI	61	-1189.1562	2517.9440	0.0000	1.0000
SYMIG	62	-1188.0706	2518.3936	0.0000	1.0000
SYMG	61	-1189.6376	2518.9070	0.0000	1.0000
F81I	59	-1195.1666	2524.7602	0.0000	1.0000
F81IG	60	-1194.1622	2525.3476	0.0000	1.0000
F81G	59	-1195.6772	2525.7813	0.0000	1.0000
HKY	59	-1196.4696	2527.3660	0.0000	1.0000
GTR	63	-1194.6203	2534.1259	0.0000	1.0000
JC69I	56	-1210.4330	2547.5756	0.0000	1.0000
JC69IG	57	-1209.4058	2548.0818	0.0000	1.0000
JC69G	56	-1210.9690	2548.6477	0.0000	1.0000
K2P	56	-1212.8967	2552.5031	0.0000	1.0000
SYM	60	-1211.0829	2559.1891	0.0000	1.0000
F81	58	-1216.8651	2565.5728	0.0000	1.0000
JC69	55	-1232.3970	2588.9548	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1173.7022	2469.4043	0.3480	0.3480
HKYI	60	-1174.7714	2469.5428	0.3247	0.6727
HKYG	60	-1175.2526	2470.5051	0.2007	0.8734
GTRIG	65	-1171.6681	2473.3362	0.0487	0.9221
GTRI	64	-1172.7043	2473.4085	0.0470	0.9691
GTRG	64	-1173.1251	2474.2501	0.0309	1.0000
K2PIG	58	-1189.5691	2495.1381	0.0000	1.0000
K2PI	57	-1190.6694	2495.3388	0.0000	1.0000
K2PG	57	-1191.2571	2496.5142	0.0000	1.0000
SYMIG	62	-1188.0706	2500.1413	0.0000	1.0000
SYMI	61	-1189.1562	2500.3123	0.0000	1.0000
SYMG	61	-1189.6376	2501.2753	0.0000	1.0000
F81IG	60	-1194.1622	2508.3244	0.0000	1.0000
F81I	59	-1195.1666	2508.3333	0.0000	1.0000
F81G	59	-1195.6772	2509.3543	0.0000	1.0000
HKY	59	-1196.4696	2510.9391	0.0000	1.0000
GTR	63	-1194.6203	2515.2407	0.0000	1.0000
JC69IG	57	-1209.4058	2532.8116	0.0000	1.0000
JC69I	56	-1210.4330	2532.8659	0.0000	1.0000
JC69G	56	-1210.9690	2533.9380	0.0000	1.0000
K2P	56	-1212.8967	2537.7934	0.0000	1.0000
SYM	60	-1211.0829	2542.1659	0.0000	1.0000
F81	58	-1216.8651	2549.7302	0.0000	1.0000
JC69	55	-1232.3970	2574.7939	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1174.7714	2721.3294	0.5708	0.5708
HKYG	60	-1175.2526	2722.2918	0.3528	0.9236
HKYIG	61	-1173.7022	2725.3874	0.0750	0.9987
K2PI	57	-1190.6694	2734.5361	0.0008	0.9994
K2PG	57	-1191.2571	2735.7116	0.0004	0.9999
K2PIG	58	-1189.5691	2738.5319	0.0001	1.0000
GTRI	64	-1172.7043	2741.9809	0.0000	1.0000
GTRG	64	-1173.1251	2742.8226	0.0000	1.0000
GTRIG	65	-1171.6681	2746.1051	0.0000	1.0000
F81I	59	-1195.1666	2755.9235	0.0000	1.0000
SYMI	61	-1189.1562	2756.2954	0.0000	1.0000
F81G	59	-1195.6772	2756.9445	0.0000	1.0000
SYMG	61	-1189.6376	2757.2584	0.0000	1.0000
HKY	59	-1196.4696	2758.5293	0.0000	1.0000
F81IG	60	-1194.1622	2760.1110	0.0000	1.0000
SYMIG	62	-1188.0706	2760.3208	0.0000	1.0000
JC69I	56	-1210.4330	2767.8668	0.0000	1.0000
JC69G	56	-1210.9690	2768.9389	0.0000	1.0000
JC69IG	57	-1209.4058	2772.0089	0.0000	1.0000
K2P	56	-1212.8967	2772.7943	0.0000	1.0000
GTR	63	-1194.6203	2779.6167	0.0000	1.0000
F81	58	-1216.8651	2793.1240	0.0000	1.0000
SYM	60	-1211.0829	2793.9525	0.0000	1.0000
JC69	55	-1232.3970	2805.5983	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_8655.phylip = 1 - 491;
 Partition Dummy = 1:183-phylip/chr2_8655.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_8655.phylip = 1 - 491;
 Partition Dummy = 1:183-phylip/chr2_8655.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_8655.phylip = 1 - 491;
 Partition Dummy = 1:183-phylip/chr2_8655.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

