
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:09 04/04/2011

Input data from file "683-tomodel-phylip/chr1_14856.phylip" (7 taxa, 802 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1507.0513	3055.1781	0.3837	0.3837
GTRG	20	-1507.0968	3055.2692	0.3666	0.7504
GTRIG	21	-1507.0968	3057.3783	0.1277	0.8781
HKYI	16	-1513.4238	3059.5406	0.0433	0.9214
HKYG	16	-1513.4303	3059.5536	0.0430	0.9645
GTR	19	-1511.0805	3061.1330	0.0195	0.9840
HKYIG	17	-1513.4373	3061.6552	0.0150	0.9991
HKY	15	-1518.3001	3067.2108	0.0009	1.0000
F81I	15	-1535.1939	3100.9986	0.0000	1.0000
F81G	15	-1535.2164	3101.0436	0.0000	1.0000
F81IG	16	-1535.2165	3103.1259	0.0000	1.0000
F81	14	-1540.6322	3109.7981	0.0000	1.0000
SYMI	17	-1550.4918	3135.7642	0.0000	1.0000
SYMG	17	-1550.5110	3135.8025	0.0000	1.0000
K2PI	13	-1555.5304	3137.5227	0.0000	1.0000
K2PG	13	-1555.5306	3137.5230	0.0000	1.0000
SYMIG	18	-1550.5110	3137.8955	0.0000	1.0000
K2PIG	14	-1555.5455	3139.6248	0.0000	1.0000
SYM	16	-1556.0560	3144.8051	0.0000	1.0000
K2P	12	-1561.9689	3148.3332	0.0000	1.0000
JC69I	12	-1575.6695	3175.7345	0.0000	1.0000
JC69G	12	-1575.6808	3175.7571	0.0000	1.0000
JC69IG	13	-1575.6808	3177.8236	0.0000	1.0000
JC69	11	-1581.7787	3185.8916	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1507.0513	3054.1025	0.3881	0.3881
GTRG	20	-1507.0968	3054.1936	0.3708	0.7588
GTRIG	21	-1507.0968	3056.1937	0.1364	0.8952
HKYI	16	-1513.4238	3058.8476	0.0362	0.9314
HKYG	16	-1513.4303	3058.8606	0.0359	0.9674
GTR	19	-1511.0805	3060.1611	0.0188	0.9861
HKYIG	17	-1513.4373	3060.8746	0.0131	0.9992
HKY	15	-1518.3001	3066.6001	0.0008	1.0000
F81I	15	-1535.1939	3100.3879	0.0000	1.0000
F81G	15	-1535.2164	3100.4329	0.0000	1.0000
F81IG	16	-1535.2165	3102.4329	0.0000	1.0000
F81	14	-1540.6322	3109.2644	0.0000	1.0000
SYMI	17	-1550.4918	3134.9836	0.0000	1.0000
SYMG	17	-1550.5110	3135.0219	0.0000	1.0000
SYMIG	18	-1550.5110	3137.0219	0.0000	1.0000
K2PI	13	-1555.5304	3137.0608	0.0000	1.0000
K2PG	13	-1555.5306	3137.0611	0.0000	1.0000
K2PIG	14	-1555.5455	3139.0911	0.0000	1.0000
SYM	16	-1556.0560	3144.1121	0.0000	1.0000
K2P	12	-1561.9689	3147.9378	0.0000	1.0000
JC69I	12	-1575.6695	3175.3390	0.0000	1.0000
JC69G	12	-1575.6808	3175.3617	0.0000	1.0000
JC69IG	13	-1575.6808	3177.3617	0.0000	1.0000
JC69	11	-1581.7787	3185.5574	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1513.4238	3133.8413	0.4451	0.4451
HKYG	16	-1513.4303	3133.8544	0.4422	0.8874
HKY	15	-1518.3001	3136.9068	0.0961	0.9835
HKYIG	17	-1513.4373	3140.5555	0.0155	0.9990
GTRI	20	-1507.0513	3147.8447	0.0004	0.9994
GTRG	20	-1507.0968	3147.9358	0.0004	0.9998
GTR	19	-1511.0805	3149.2162	0.0002	1.0000
GTRIG	21	-1507.0968	3154.6230	0.0000	1.0000
F81I	15	-1535.1939	3170.6945	0.0000	1.0000
F81G	15	-1535.2164	3170.7395	0.0000	1.0000
F81	14	-1540.6322	3174.8839	0.0000	1.0000
F81IG	16	-1535.2165	3177.4267	0.0000	1.0000
K2PI	13	-1555.5304	3197.9932	0.0000	1.0000
K2PG	13	-1555.5306	3197.9935	0.0000	1.0000
K2P	12	-1561.9689	3204.1831	0.0000	1.0000
K2PIG	14	-1555.5455	3204.7106	0.0000	1.0000
SYMI	17	-1550.4918	3214.6645	0.0000	1.0000
SYMG	17	-1550.5110	3214.7028	0.0000	1.0000
SYM	16	-1556.0560	3219.1058	0.0000	1.0000
SYMIG	18	-1550.5110	3221.3899	0.0000	1.0000
JC69I	12	-1575.6695	3231.5843	0.0000	1.0000
JC69G	12	-1575.6808	3231.6070	0.0000	1.0000
JC69	11	-1581.7787	3237.1156	0.0000	1.0000
JC69IG	13	-1575.6808	3238.2941	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_14856.phylip = 1 - 802;
 Partition Dummy = 1:683-tomodel-phylip/chr1_14856.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_14856.phylip = 1 - 802;
 Partition Dummy = 1:683-tomodel-phylip/chr1_14856.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_14856.phylip = 1 - 802;
 Partition Dummy = 1:683-tomodel-phylip/chr1_14856.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

