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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: K2PI


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	16	-1743.3597	3519.7400	0.5447	0.5447
HKYG	16	-1744.9995	3523.0197	0.1057	0.6503
GTRI	20	-1740.9741	3523.5361	0.0816	0.7320
K2PI	13	-1748.5436	3523.7662	0.0728	0.8047
GTRG	20	-1741.1269	3523.8417	0.0701	0.8748
HKYIG	17	-1744.4329	3524.0162	0.0642	0.9390
GTRIG	21	-1741.0882	3525.9265	0.0247	0.9637
K2PG	13	-1750.2242	3527.1274	0.0136	0.9772
K2PIG	14	-1749.7225	3528.2301	0.0078	0.9850
SYMI	17	-1746.6518	3528.4539	0.0070	0.9920
SYMG	17	-1746.7332	3528.6168	0.0064	0.9984
SYMIG	18	-1747.0783	3531.4448	0.0016	1.0000
HKY	15	-1768.3813	3567.6615	0.0000	1.0000
GTR	19	-1764.1458	3567.7255	0.0000	1.0000
K2P	12	-1773.7456	3572.0721	0.0000	1.0000
SYM	16	-1770.0042	3573.0291	0.0000	1.0000
F81G	15	-1791.8223	3614.5434	0.0000	1.0000
F81I	15	-1792.6697	3616.2382	0.0000	1.0000
F81IG	16	-1791.7739	3616.5685	0.0000	1.0000
JC69G	12	-1796.5411	3617.6632	0.0000	1.0000
JC69IG	13	-1796.4832	3619.6455	0.0000	1.0000
JC69I	12	-1797.5687	3619.7184	0.0000	1.0000
F81	14	-1811.7046	3652.1943	0.0000	1.0000
JC69	11	-1816.5180	3655.5267	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1743.3597	3518.7193	0.5179	0.5179
GTRI	20	-1740.9741	3521.9482	0.1031	0.6210
HKYG	16	-1744.9995	3521.9990	0.1005	0.7215
GTRG	20	-1741.1269	3522.2538	0.0885	0.8099
HKYIG	17	-1744.4329	3522.8658	0.0651	0.8751
K2PI	13	-1748.5436	3523.0871	0.0583	0.9334
GTRIG	21	-1741.0882	3524.1765	0.0338	0.9672
K2PG	13	-1750.2242	3526.4483	0.0109	0.9781
SYMI	17	-1746.6518	3527.3035	0.0071	0.9852
K2PIG	14	-1749.7225	3527.4451	0.0066	0.9918
SYMG	17	-1746.7332	3527.4664	0.0065	0.9983
SYMIG	18	-1747.0783	3530.1567	0.0017	1.0000
GTR	19	-1764.1458	3566.2916	0.0000	1.0000
HKY	15	-1768.3813	3566.7626	0.0000	1.0000
K2P	12	-1773.7456	3571.4911	0.0000	1.0000
SYM	16	-1770.0042	3572.0084	0.0000	1.0000
F81G	15	-1791.8223	3613.6445	0.0000	1.0000
F81I	15	-1792.6697	3615.3393	0.0000	1.0000
F81IG	16	-1791.7739	3615.5479	0.0000	1.0000
JC69G	12	-1796.5411	3617.0822	0.0000	1.0000
JC69IG	13	-1796.4832	3618.9664	0.0000	1.0000
JC69I	12	-1797.5687	3619.1374	0.0000	1.0000
F81	14	-1811.7046	3651.4093	0.0000	1.0000
JC69	11	-1816.5180	3655.0360	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1748.5436	3579.1160	0.8223	0.8223
K2PG	13	-1750.2242	3582.4773	0.1532	0.9754
HKYI	16	-1743.3597	3587.6780	0.0114	0.9868
K2PIG	14	-1749.7225	3587.7840	0.0108	0.9976
HKYG	16	-1744.9995	3590.9577	0.0022	0.9998
HKYIG	17	-1744.4329	3596.1345	0.0002	1.0000
SYMI	17	-1746.6518	3600.5722	0.0000	1.0000
SYMG	17	-1746.7332	3600.7350	0.0000	1.0000
SYMIG	18	-1747.0783	3607.7352	0.0000	1.0000
GTRI	20	-1740.9741	3608.1465	0.0000	1.0000
GTRG	20	-1741.1269	3608.4522	0.0000	1.0000
GTRIG	21	-1741.0882	3614.6848	0.0000	1.0000
K2P	12	-1773.7456	3623.2101	0.0000	1.0000
HKY	15	-1768.3813	3631.4114	0.0000	1.0000
SYM	16	-1770.0042	3640.9671	0.0000	1.0000
GTR	19	-1764.1458	3648.1800	0.0000	1.0000
JC69G	12	-1796.5411	3668.8012	0.0000	1.0000
JC69I	12	-1797.5687	3670.8564	0.0000	1.0000
JC69IG	13	-1796.4832	3674.9953	0.0000	1.0000
F81G	15	-1791.8223	3678.2933	0.0000	1.0000
F81I	15	-1792.6697	3679.9881	0.0000	1.0000
F81IG	16	-1791.7739	3684.5066	0.0000	1.0000
JC69	11	-1816.5180	3702.4451	0.0000	1.0000
F81	14	-1811.7046	3711.7481	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_38361.phylip = 1 - 550;
 Partition Dummy = 1:591-tomodel-phylip/chr1_38361.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 AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_38361.phylip = 1 - 550;
 Partition Dummy = 1:591-tomodel-phylip/chr1_38361.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 (K2PI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_38361.phylip = 1 - 550;
 Partition Dummy = 1:591-tomodel-phylip/chr1_38361.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

