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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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
GTRG	20	-1944.0870	3929.3723	0.4785	0.4785
GTRI	20	-1944.8124	3930.8231	0.2316	0.7101
GTRIG	21	-1943.8693	3931.0586	0.2059	0.9160
HKYG	16	-1950.7200	3934.2116	0.0426	0.9586
HKYI	16	-1951.2364	3935.2445	0.0254	0.9840
HKYIG	17	-1950.6477	3936.1647	0.0160	1.0000
K2PG	13	-1978.5830	3983.6802	0.0000	1.0000
K2PI	13	-1978.8211	3984.1564	0.0000	1.0000
SYMG	17	-1975.4229	3985.7152	0.0000	1.0000
K2PIG	14	-1978.5732	3985.7405	0.0000	1.0000
SYMI	17	-1975.6960	3986.2613	0.0000	1.0000
SYMIG	18	-1975.3575	3987.6880	0.0000	1.0000
GTR	19	-1979.6985	3998.4796	0.0000	1.0000
HKY	15	-1986.3933	4003.4665	0.0000	1.0000
F81I	15	-1991.1186	4012.9171	0.0000	1.0000
F81IG	16	-1990.7039	4014.1795	0.0000	1.0000
F81G	15	-1992.7915	4016.2630	0.0000	1.0000
JC69G	12	-2014.7297	4053.8994	0.0000	1.0000
K2P	12	-2014.7554	4053.9509	0.0000	1.0000
JC69I	12	-2014.9063	4054.2527	0.0000	1.0000
JC69IG	13	-2014.7136	4055.9413	0.0000	1.0000
SYM	16	-2011.6098	4055.9911	0.0000	1.0000
F81	14	-2022.2463	4073.0866	0.0000	1.0000
JC69	11	-2046.6948	4115.7614	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1944.0870	3928.1741	0.4797	0.4797
GTRI	20	-1944.8124	3929.6248	0.2322	0.7119
GTRIG	21	-1943.8693	3929.7386	0.2194	0.9313
HKYG	16	-1950.7200	3933.4400	0.0345	0.9658
HKYI	16	-1951.2364	3934.4729	0.0206	0.9864
HKYIG	17	-1950.6477	3935.2953	0.0136	1.0000
K2PG	13	-1978.5830	3983.1661	0.0000	1.0000
K2PI	13	-1978.8211	3983.6423	0.0000	1.0000
SYMG	17	-1975.4229	3984.8459	0.0000	1.0000
K2PIG	14	-1978.5732	3985.1464	0.0000	1.0000
SYMI	17	-1975.6960	3985.3920	0.0000	1.0000
SYMIG	18	-1975.3575	3986.7150	0.0000	1.0000
GTR	19	-1979.6985	3997.3969	0.0000	1.0000
HKY	15	-1986.3933	4002.7867	0.0000	1.0000
F81I	15	-1991.1186	4012.2372	0.0000	1.0000
F81IG	16	-1990.7039	4013.4078	0.0000	1.0000
F81G	15	-1992.7915	4015.5831	0.0000	1.0000
JC69G	12	-2014.7297	4053.4593	0.0000	1.0000
K2P	12	-2014.7554	4053.5108	0.0000	1.0000
JC69I	12	-2014.9063	4053.8126	0.0000	1.0000
SYM	16	-2011.6098	4055.2195	0.0000	1.0000
JC69IG	13	-2014.7136	4055.4272	0.0000	1.0000
F81	14	-2022.2463	4072.4926	0.0000	1.0000
JC69	11	-2046.6948	4115.3896	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1950.7200	4006.7524	0.6102	0.6102
HKYI	16	-1951.2364	4007.7853	0.3641	0.9742
HKYIG	17	-1950.6477	4013.1898	0.0244	0.9986
GTRG	20	-1944.0870	4019.8146	0.0009	0.9995
GTRI	20	-1944.8124	4021.2653	0.0004	1.0000
GTRIG	21	-1943.8693	4025.9611	0.0000	1.0000
K2PG	13	-1978.5830	4042.7324	0.0000	1.0000
K2PI	13	-1978.8211	4043.2086	0.0000	1.0000
K2PIG	14	-1978.5732	4049.2948	0.0000	1.0000
SYMG	17	-1975.4229	4062.7403	0.0000	1.0000
SYMI	17	-1975.6960	4063.2864	0.0000	1.0000
SYMIG	18	-1975.3575	4069.1915	0.0000	1.0000
HKY	15	-1986.3933	4071.5170	0.0000	1.0000
F81I	15	-1991.1186	4080.9676	0.0000	1.0000
F81G	15	-1992.7915	4084.3135	0.0000	1.0000
GTR	19	-1979.6985	4084.4554	0.0000	1.0000
F81IG	16	-1990.7039	4086.7202	0.0000	1.0000
JC69G	12	-2014.7297	4108.4436	0.0000	1.0000
K2P	12	-2014.7554	4108.4951	0.0000	1.0000
JC69I	12	-2014.9063	4108.7969	0.0000	1.0000
JC69IG	13	-2014.7136	4114.9935	0.0000	1.0000
SYM	16	-2011.6098	4128.5319	0.0000	1.0000
F81	14	-2022.2463	4136.6409	0.0000	1.0000
JC69	11	-2046.6948	4165.7919	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_8735.phylip = 1 - 722;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8735.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_8735.phylip = 1 - 722;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8735.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_8735.phylip = 1 - 722;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8735.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

