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

Input data from file "591-tomodel-phylip/chr7_4698.phylip" (7 taxa, 886 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	-1791.7568	3623.5137	0.4962	0.4962
GTRG	20	-1792.5293	3625.0587	0.2292	0.7254
GTRIG	21	-1791.7568	3625.5135	0.1826	0.9080
HKYI	16	-1797.9804	3627.9608	0.0537	0.9617
HKYIG	17	-1797.9800	3629.9600	0.0198	0.9814
HKYG	16	-1799.0432	3630.0864	0.0186	1.0000
SYMI	17	-1813.5495	3661.0991	0.0000	1.0000
SYMIG	18	-1812.8545	3661.7090	0.0000	1.0000
SYMG	17	-1814.4658	3662.9316	0.0000	1.0000
K2PI	13	-1818.5987	3663.1974	0.0000	1.0000
GTR	19	-1812.8480	3663.6959	0.0000	1.0000
K2PIG	14	-1818.5986	3665.1971	0.0000	1.0000
K2PG	13	-1820.5793	3667.1587	0.0000	1.0000
HKY	15	-1819.7844	3669.5688	0.0000	1.0000
F81I	15	-1820.4393	3670.8785	0.0000	1.0000
F81IG	16	-1820.4392	3672.8784	0.0000	1.0000
F81G	15	-1822.3025	3674.6050	0.0000	1.0000
SYM	16	-1834.8154	3701.6308	0.0000	1.0000
JC69I	12	-1839.9719	3703.9437	0.0000	1.0000
K2P	12	-1840.7108	3705.4216	0.0000	1.0000
JC69IG	13	-1839.9716	3705.9432	0.0000	1.0000
JC69G	12	-1841.0022	3706.0044	0.0000	1.0000
F81	14	-1841.8373	3711.6746	0.0000	1.0000
JC69	11	-1860.6078	3743.2156	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1797.9804	3704.5482	0.7246	0.7246
HKYG	16	-1799.0432	3706.6739	0.2503	0.9749
HKYIG	17	-1797.9800	3711.3342	0.0244	0.9993
GTRI	20	-1791.7568	3719.2480	0.0005	0.9997
GTRG	20	-1792.5293	3720.7930	0.0002	1.0000
K2PI	13	-1818.5987	3725.4247	0.0000	1.0000
GTRIG	21	-1791.7568	3726.0346	0.0000	1.0000
K2PG	13	-1820.5793	3729.3860	0.0000	1.0000
K2PIG	14	-1818.5986	3732.2112	0.0000	1.0000
HKY	15	-1819.7844	3741.3696	0.0000	1.0000
SYMI	17	-1813.5495	3742.4732	0.0000	1.0000
F81I	15	-1820.4393	3742.6793	0.0000	1.0000
SYMG	17	-1814.4658	3744.3058	0.0000	1.0000
F81G	15	-1822.3025	3746.4057	0.0000	1.0000
SYMIG	18	-1812.8545	3747.8699	0.0000	1.0000
F81IG	16	-1820.4392	3749.4659	0.0000	1.0000
GTR	19	-1812.8480	3754.6435	0.0000	1.0000
JC69I	12	-1839.9719	3761.3843	0.0000	1.0000
K2P	12	-1840.7108	3762.8622	0.0000	1.0000
JC69G	12	-1841.0022	3763.4450	0.0000	1.0000
JC69IG	13	-1839.9716	3768.1706	0.0000	1.0000
SYM	16	-1834.8154	3778.2183	0.0000	1.0000
F81	14	-1841.8373	3778.6887	0.0000	1.0000
JC69	11	-1860.6078	3795.8695	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 591-tomodel-phylip/chr7_4698.phylip = 1 - 886;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4698.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 591-tomodel-phylip/chr7_4698.phylip = 1 - 886;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4698.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 591-tomodel-phylip/chr7_4698.phylip = 1 - 886;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4698.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

