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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1972.0351	3985.1189	0.3118	0.3118
GTRG	20	-1972.0562	3985.1610	0.3053	0.6171
HKYI	16	-1977.1744	3987.0247	0.1202	0.7373
HKYG	16	-1977.2099	3987.0957	0.1160	0.8533
GTRIG	21	-1972.0629	3987.2808	0.1058	0.9591
HKYIG	17	-1977.2100	3989.1811	0.0409	1.0000
K2PG	13	-1999.6563	4025.7630	0.0000	1.0000
K2PI	13	-1999.6745	4025.7995	0.0000	1.0000
SYMG	17	-1995.8316	4026.4244	0.0000	1.0000
GTR	19	-1993.7445	4026.4366	0.0000	1.0000
SYMI	17	-1995.8691	4026.4995	0.0000	1.0000
K2PIG	14	-1999.6604	4027.8413	0.0000	1.0000
SYMIG	18	-1995.8316	4028.5149	0.0000	1.0000
HKY	15	-1999.4405	4029.4764	0.0000	1.0000
K2P	12	-2024.1994	4072.7844	0.0000	1.0000
SYM	16	-2020.0873	4072.8505	0.0000	1.0000
F81I	15	-2042.3822	4115.3600	0.0000	1.0000
F81G	15	-2042.4615	4115.5186	0.0000	1.0000
F81IG	16	-2042.4507	4117.5773	0.0000	1.0000
JC69G	12	-2062.6006	4149.5868	0.0000	1.0000
JC69I	12	-2062.6614	4149.7085	0.0000	1.0000
JC69IG	13	-2062.5498	4151.5500	0.0000	1.0000
F81	14	-2062.9062	4154.3328	0.0000	1.0000
JC69	11	-2084.2737	4190.8733	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1972.0351	3984.0703	0.3247	0.3247
GTRG	20	-1972.0562	3984.1123	0.3180	0.6427
GTRIG	21	-1972.0629	3986.1258	0.1162	0.7589
HKYI	16	-1977.1744	3986.3489	0.1039	0.8628
HKYG	16	-1977.2099	3986.4199	0.1003	0.9631
HKYIG	17	-1977.2100	3988.4199	0.0369	1.0000
K2PG	13	-1999.6563	4025.3125	0.0000	1.0000
K2PI	13	-1999.6745	4025.3490	0.0000	1.0000
GTR	19	-1993.7445	4025.4890	0.0000	1.0000
SYMG	17	-1995.8316	4025.6632	0.0000	1.0000
SYMI	17	-1995.8691	4025.7383	0.0000	1.0000
K2PIG	14	-1999.6604	4027.3208	0.0000	1.0000
SYMIG	18	-1995.8316	4027.6631	0.0000	1.0000
HKY	15	-1999.4405	4028.8809	0.0000	1.0000
SYM	16	-2020.0873	4072.1747	0.0000	1.0000
K2P	12	-2024.1994	4072.3987	0.0000	1.0000
F81I	15	-2042.3822	4114.7644	0.0000	1.0000
F81G	15	-2042.4615	4114.9231	0.0000	1.0000
F81IG	16	-2042.4507	4116.9015	0.0000	1.0000
JC69G	12	-2062.6006	4149.2011	0.0000	1.0000
JC69I	12	-2062.6614	4149.3229	0.0000	1.0000
JC69IG	13	-2062.5498	4151.0995	0.0000	1.0000
F81	14	-2062.9062	4153.8123	0.0000	1.0000
JC69	11	-2084.2737	4190.5474	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1977.1744	4061.7367	0.5002	0.5002
HKYG	16	-1977.2099	4061.8077	0.4827	0.9829
HKYIG	17	-1977.2100	4068.5195	0.0168	0.9997
GTRI	20	-1972.0351	4078.3051	0.0001	0.9999
GTRG	20	-1972.0562	4078.3471	0.0001	1.0000
GTRIG	21	-1972.0629	4085.0724	0.0000	1.0000
K2PG	13	-1999.6563	4086.5651	0.0000	1.0000
K2PI	13	-1999.6745	4086.6016	0.0000	1.0000
K2PIG	14	-1999.6604	4093.2852	0.0000	1.0000
HKY	15	-1999.4405	4099.5570	0.0000	1.0000
SYMG	17	-1995.8316	4105.7628	0.0000	1.0000
SYMI	17	-1995.8691	4105.8379	0.0000	1.0000
SYMIG	18	-1995.8316	4112.4744	0.0000	1.0000
GTR	19	-1993.7445	4115.0121	0.0000	1.0000
K2P	12	-2024.1994	4128.9396	0.0000	1.0000
SYM	16	-2020.0873	4147.5625	0.0000	1.0000
F81I	15	-2042.3822	4185.4405	0.0000	1.0000
F81G	15	-2042.4615	4185.5992	0.0000	1.0000
F81IG	16	-2042.4507	4192.2893	0.0000	1.0000
JC69G	12	-2062.6006	4205.7420	0.0000	1.0000
JC69I	12	-2062.6614	4205.8637	0.0000	1.0000
JC69IG	13	-2062.5498	4212.3522	0.0000	1.0000
F81	14	-2062.9062	4219.7767	0.0000	1.0000
JC69	11	-2084.2737	4242.3765	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/chr6_5979.phylip = 1 - 822;
 Partition Dummy = 1:591-tomodel-phylip/chr6_5979.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/chr6_5979.phylip = 1 - 822;
 Partition Dummy = 1:591-tomodel-phylip/chr6_5979.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/chr6_5979.phylip = 1 - 822;
 Partition Dummy = 1:591-tomodel-phylip/chr6_5979.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

