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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1779.5180	3599.0360	0.5674	0.5674
GTRIG	21	-1778.9886	3599.9772	0.3544	0.9218
HKYI	16	-1786.2960	3604.5920	0.0353	0.9571
HKYIG	17	-1785.5953	3605.1906	0.0261	0.9832
GTRG	20	-1783.0714	3606.1429	0.0162	0.9995
HKYG	16	-1790.5117	3613.0234	0.0005	1.0000
GTR	19	-1797.5205	3633.0411	0.0000	1.0000
HKY	15	-1805.0208	3640.0415	0.0000	1.0000
F81I	15	-1812.6776	3655.3552	0.0000	1.0000
F81IG	16	-1812.0002	3656.0004	0.0000	1.0000
F81G	15	-1815.7334	3661.4668	0.0000	1.0000
SYMI	17	-1819.4868	3672.9736	0.0000	1.0000
SYMIG	18	-1818.9320	3673.8639	0.0000	1.0000
K2PI	13	-1824.1701	3674.3401	0.0000	1.0000
K2PIG	14	-1823.3886	3674.7772	0.0000	1.0000
SYMG	17	-1822.8086	3679.6172	0.0000	1.0000
K2PG	13	-1828.0799	3682.1597	0.0000	1.0000
F81	14	-1829.4040	3686.8080	0.0000	1.0000
SYM	16	-1836.7907	3705.5814	0.0000	1.0000
K2P	12	-1842.3933	3708.7866	0.0000	1.0000
JC69I	12	-1847.9154	3719.8308	0.0000	1.0000
JC69IG	13	-1847.2092	3720.4184	0.0000	1.0000
JC69G	12	-1850.9979	3725.9958	0.0000	1.0000
JC69	11	-1864.6689	3751.3377	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1786.2960	3683.2118	0.9266	0.9266
HKYIG	17	-1785.5953	3688.7241	0.0589	0.9855
HKYG	16	-1790.5117	3691.6432	0.0137	0.9991
GTRI	20	-1779.5180	3697.3107	0.0008	0.9999
GTRIG	21	-1778.9886	3703.1657	0.0000	1.0000
GTRG	20	-1783.0714	3704.4176	0.0000	1.0000
HKY	15	-1805.0208	3713.7476	0.0000	1.0000
GTR	19	-1797.5205	3726.4021	0.0000	1.0000
F81I	15	-1812.6776	3729.0613	0.0000	1.0000
F81IG	16	-1812.0002	3734.6202	0.0000	1.0000
F81G	15	-1815.7334	3735.1729	0.0000	1.0000
K2PI	13	-1824.1701	3738.2187	0.0000	1.0000
K2PIG	14	-1823.3886	3743.5695	0.0000	1.0000
K2PG	13	-1828.0799	3746.0383	0.0000	1.0000
F81	14	-1829.4040	3755.6003	0.0000	1.0000
SYMI	17	-1819.4868	3756.5072	0.0000	1.0000
SYMIG	18	-1818.9320	3762.3112	0.0000	1.0000
SYMG	17	-1822.8086	3763.1508	0.0000	1.0000
K2P	12	-1842.3933	3767.7515	0.0000	1.0000
JC69I	12	-1847.9154	3778.7957	0.0000	1.0000
SYM	16	-1836.7907	3784.2012	0.0000	1.0000
JC69IG	13	-1847.2092	3784.2970	0.0000	1.0000
JC69G	12	-1850.9979	3784.9607	0.0000	1.0000
JC69	11	-1864.6689	3805.3889	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/chrZ_11372.phylip = 1 - 1006;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_11372.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/chrZ_11372.phylip = 1 - 1006;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_11372.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/chrZ_11372.phylip = 1 - 1006;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_11372.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

