
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:41 04/04/2011

Input data from file "591-tomodel-phylip/chr5_11657.phylip" (7 taxa, 1156 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2812.0121	5656.0243	0.5266	0.5266
HKYI	16	-2812.7743	5657.5486	0.2458	0.7724
HKYIG	17	-2811.9477	5657.8954	0.2066	0.9790
GTRG	20	-2811.8521	5663.7042	0.0113	0.9904
GTRI	20	-2812.6296	5665.2593	0.0052	0.9956
GTRIG	21	-2811.7879	5665.5757	0.0044	1.0000
HKY	15	-2826.7010	5683.4020	0.0000	1.0000
GTR	19	-2826.4544	5690.9089	0.0000	1.0000
K2PG	13	-2849.5997	5725.1994	0.0000	1.0000
K2PI	13	-2850.4182	5726.8365	0.0000	1.0000
K2PIG	14	-2849.5211	5727.0422	0.0000	1.0000
SYMG	17	-2848.2758	5730.5516	0.0000	1.0000
SYMI	17	-2849.0615	5732.1229	0.0000	1.0000
SYMIG	18	-2848.1967	5732.3934	0.0000	1.0000
K2P	12	-2864.4440	5752.8880	0.0000	1.0000
SYM	16	-2863.2475	5758.4951	0.0000	1.0000
F81G	15	-2867.7012	5765.4025	0.0000	1.0000
F81I	15	-2868.3784	5766.7568	0.0000	1.0000
F81IG	16	-2867.6482	5767.2965	0.0000	1.0000
F81	14	-2881.1876	5790.3751	0.0000	1.0000
JC69G	12	-2900.6320	5825.2640	0.0000	1.0000
JC69I	12	-2901.3311	5826.6621	0.0000	1.0000
JC69IG	13	-2900.5712	5827.1424	0.0000	1.0000
JC69	11	-2914.2491	5850.4981	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2812.0121	5736.8678	0.6675	0.6675
HKYI	16	-2812.7743	5738.3922	0.3115	0.9790
HKYIG	17	-2811.9477	5743.7917	0.0209	1.0000
HKY	15	-2826.7010	5759.1928	0.0000	1.0000
GTRG	20	-2811.8521	5764.7587	0.0000	1.0000
GTRI	20	-2812.6296	5766.3137	0.0000	1.0000
GTRIG	21	-2811.7879	5771.6829	0.0000	1.0000
GTR	19	-2826.4544	5786.9106	0.0000	1.0000
K2PG	13	-2849.5997	5790.8848	0.0000	1.0000
K2PI	13	-2850.4182	5792.5219	0.0000	1.0000
K2PIG	14	-2849.5211	5797.7803	0.0000	1.0000
K2P	12	-2864.4440	5813.5207	0.0000	1.0000
SYMG	17	-2848.2758	5816.4479	0.0000	1.0000
SYMI	17	-2849.0615	5818.0192	0.0000	1.0000
SYMIG	18	-2848.1967	5823.3424	0.0000	1.0000
SYM	16	-2863.2475	5839.3386	0.0000	1.0000
F81G	15	-2867.7012	5841.1933	0.0000	1.0000
F81I	15	-2868.3784	5842.5477	0.0000	1.0000
F81IG	16	-2867.6482	5848.1400	0.0000	1.0000
F81	14	-2881.1876	5861.1132	0.0000	1.0000
JC69G	12	-2900.6320	5885.8967	0.0000	1.0000
JC69I	12	-2901.3311	5887.2948	0.0000	1.0000
JC69IG	13	-2900.5712	5892.8278	0.0000	1.0000
JC69	11	-2914.2491	5906.0781	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11657.phylip = 1 - 1156;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11657.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11657.phylip = 1 - 1156;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11657.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11657.phylip = 1 - 1156;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11657.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

