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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1881.4255	3794.8510	0.3719	0.3719
HKYG	16	-1881.4765	3794.9530	0.3534	0.7252
HKYIG	17	-1881.4765	3796.9530	0.1300	0.8552
GTRI	20	-1879.2162	3798.4325	0.0620	0.9173
GTRG	20	-1879.2741	3798.5483	0.0586	0.9758
GTRIG	21	-1879.2741	3800.5482	0.0215	0.9974
F81I	15	-1888.3365	3806.6730	0.0010	0.9984
F81G	15	-1888.4063	3806.8126	0.0009	0.9993
F81IG	16	-1888.4062	3808.8124	0.0003	0.9997
HKY	15	-1889.5586	3809.1172	0.0003	1.0000
GTR	19	-1887.7437	3813.4874	0.0000	1.0000
F81	14	-1896.0234	3820.0469	0.0000	1.0000
K2PI	13	-1907.6849	3841.3699	0.0000	1.0000
K2PG	13	-1907.7265	3841.4529	0.0000	1.0000
K2PIG	14	-1907.7971	3843.5942	0.0000	1.0000
SYMI	17	-1906.6212	3847.2424	0.0000	1.0000
SYMG	17	-1906.6763	3847.3526	0.0000	1.0000
SYMIG	18	-1906.6763	3849.3525	0.0000	1.0000
JC69I	12	-1915.6648	3855.3297	0.0000	1.0000
JC69G	12	-1915.7362	3855.4725	0.0000	1.0000
K2P	12	-1915.9709	3855.9419	0.0000	1.0000
JC69IG	13	-1915.7362	3857.4725	0.0000	1.0000
SYM	16	-1914.9754	3861.9509	0.0000	1.0000
JC69	11	-1923.2725	3868.5450	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1881.4255	3872.2997	0.4877	0.4877
HKYG	16	-1881.4765	3872.4018	0.4634	0.9511
HKYIG	17	-1881.4765	3879.2423	0.0152	0.9662
F81I	15	-1888.3365	3879.2812	0.0149	0.9811
F81G	15	-1888.4063	3879.4208	0.0139	0.9949
HKY	15	-1889.5586	3881.7254	0.0044	0.9993
F81IG	16	-1888.4062	3886.2611	0.0005	0.9998
F81	14	-1896.0234	3887.8146	0.0002	1.0000
GTRI	20	-1879.2162	3895.2434	0.0000	1.0000
GTRG	20	-1879.2741	3895.3592	0.0000	1.0000
GTRIG	21	-1879.2741	3902.1997	0.0000	1.0000
K2PI	13	-1907.6849	3904.2970	0.0000	1.0000
K2PG	13	-1907.7265	3904.3800	0.0000	1.0000
GTR	19	-1887.7437	3905.4577	0.0000	1.0000
K2PIG	14	-1907.7971	3911.3619	0.0000	1.0000
JC69I	12	-1915.6648	3913.4163	0.0000	1.0000
JC69G	12	-1915.7362	3913.5591	0.0000	1.0000
K2P	12	-1915.9709	3914.0284	0.0000	1.0000
JC69IG	13	-1915.7362	3920.3996	0.0000	1.0000
JC69	11	-1923.2725	3921.7911	0.0000	1.0000
SYMI	17	-1906.6212	3929.5317	0.0000	1.0000
SYMG	17	-1906.6763	3929.6419	0.0000	1.0000
SYMIG	18	-1906.6763	3936.4824	0.0000	1.0000
SYM	16	-1914.9754	3939.3996	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 683-tomodel-phylip/chr8_4030.phylip = 1 - 935;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4030.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 683-tomodel-phylip/chr8_4030.phylip = 1 - 935;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4030.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 (HKYI)]

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

