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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	16	-930.6693	1895.0074	0.4441	0.4441
HKYG	16	-930.9569	1895.5825	0.3332	0.7773
HKYIG	17	-930.6762	1897.2355	0.1458	0.9231
HKY	15	-934.5651	1900.5980	0.0271	0.9502
GTRI	20	-929.0874	1900.7834	0.0247	0.9750
GTRG	20	-929.4017	1901.4122	0.0181	0.9930
GTRIG	21	-929.3922	1903.6629	0.0059	0.9989
GTR	19	-933.3101	1906.9732	0.0011	1.0000
K2PI	13	-947.7363	1922.5790	0.0000	1.0000
K2PG	13	-947.9966	1923.0996	0.0000	1.0000
K2PIG	14	-947.7426	1924.7656	0.0000	1.0000
SYMI	17	-946.6224	1929.1279	0.0000	1.0000
K2P	12	-952.1146	1929.1747	0.0000	1.0000
SYMG	17	-946.7975	1929.4781	0.0000	1.0000
SYMIG	18	-946.6256	1931.3623	0.0000	1.0000
SYM	16	-951.1234	1935.9156	0.0000	1.0000
F81I	15	-957.0798	1945.6275	0.0000	1.0000
F81G	15	-957.3264	1946.1206	0.0000	1.0000
F81IG	16	-957.0868	1947.8423	0.0000	1.0000
F81	14	-960.2811	1949.8426	0.0000	1.0000
JC69I	12	-972.8627	1970.6708	0.0000	1.0000
JC69G	12	-973.0701	1971.0857	0.0000	1.0000
JC69IG	13	-972.8677	1972.8417	0.0000	1.0000
JC69	11	-976.1569	1975.1113	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-930.6693	1893.3387	0.4251	0.4251
HKYG	16	-930.9569	1893.9137	0.3188	0.7439
HKYIG	17	-930.6762	1895.3524	0.1553	0.8992
GTRI	20	-929.0874	1898.1747	0.0379	0.9371
GTRG	20	-929.4017	1898.8035	0.0277	0.9647
HKY	15	-934.5651	1899.1301	0.0235	0.9882
GTRIG	21	-929.3922	1900.7844	0.0103	0.9985
GTR	19	-933.3101	1904.6202	0.0015	1.0000
K2PI	13	-947.7363	1921.4726	0.0000	1.0000
K2PG	13	-947.9966	1921.9932	0.0000	1.0000
K2PIG	14	-947.7426	1923.4851	0.0000	1.0000
SYMI	17	-946.6224	1927.2448	0.0000	1.0000
SYMG	17	-946.7975	1927.5950	0.0000	1.0000
K2P	12	-952.1146	1928.2292	0.0000	1.0000
SYMIG	18	-946.6256	1929.2512	0.0000	1.0000
SYM	16	-951.1234	1934.2468	0.0000	1.0000
F81I	15	-957.0798	1944.1596	0.0000	1.0000
F81G	15	-957.3264	1944.6527	0.0000	1.0000
F81IG	16	-957.0868	1946.1735	0.0000	1.0000
F81	14	-960.2811	1948.5621	0.0000	1.0000
JC69I	12	-972.8627	1969.7253	0.0000	1.0000
JC69G	12	-973.0701	1970.1402	0.0000	1.0000
JC69IG	13	-972.8677	1971.7353	0.0000	1.0000
JC69	11	-976.1569	1974.3137	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-930.6693	1954.7423	0.4585	0.4585
HKYG	16	-930.9569	1955.3174	0.3440	0.8025
HKY	15	-934.5651	1956.6961	0.1726	0.9751
HKYIG	17	-930.6762	1960.5938	0.0246	0.9997
K2PI	13	-947.7363	1971.3631	0.0001	0.9998
K2PG	13	-947.9966	1971.8837	0.0001	0.9999
K2P	12	-952.1146	1974.2820	0.0000	1.0000
GTRI	20	-929.0874	1974.9293	0.0000	1.0000
GTRG	20	-929.4017	1975.5581	0.0000	1.0000
K2PIG	14	-947.7426	1977.2133	0.0000	1.0000
GTR	19	-933.3101	1977.5371	0.0000	1.0000
GTRIG	21	-929.3922	1981.3767	0.0000	1.0000
SYMI	17	-946.6224	1992.4863	0.0000	1.0000
SYMG	17	-946.7975	1992.8364	0.0000	1.0000
SYM	16	-951.1234	1995.6505	0.0000	1.0000
SYMIG	18	-946.6256	1998.3303	0.0000	1.0000
F81I	15	-957.0798	2001.7256	0.0000	1.0000
F81G	15	-957.3264	2002.2187	0.0000	1.0000
F81	14	-960.2811	2002.2903	0.0000	1.0000
F81IG	16	-957.0868	2007.5772	0.0000	1.0000
JC69I	12	-972.8627	2015.7781	0.0000	1.0000
JC69G	12	-973.0701	2016.1930	0.0000	1.0000
JC69	11	-976.1569	2016.5288	0.0000	1.0000
JC69IG	13	-972.8677	2021.6258	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/chr6_3198.phylip = 1 - 343;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3198.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/chr6_3198.phylip = 1 - 343;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3198.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 591-tomodel-phylip/chr6_3198.phylip = 1 - 343;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3198.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

