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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: SYMG


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	-2039.1229	4119.3689	0.4878	0.4878
GTRI	20	-2039.6076	4120.3382	0.3005	0.7883
GTRIG	21	-2038.9581	4121.1531	0.1999	0.9882
SYMG	17	-2046.7951	4128.4051	0.0053	0.9936
SYMI	17	-2046.9615	4128.7380	0.0045	0.9981
SYMIG	18	-2046.7549	4130.4218	0.0019	1.0000
HKYI	16	-2061.4458	4155.6150	0.0000	1.0000
HKYG	16	-2061.5580	4155.8395	0.0000	1.0000
K2PI	13	-2065.6643	4157.8107	0.0000	1.0000
HKYIG	17	-2061.5580	4157.9310	0.0000	1.0000
K2PG	13	-2065.8247	4158.1315	0.0000	1.0000
GTR	19	-2060.2928	4159.6002	0.0000	1.0000
K2PIG	14	-2065.8247	4160.2064	0.0000	1.0000
SYM	16	-2066.5698	4165.8629	0.0000	1.0000
F81I	15	-2077.5720	4185.7814	0.0000	1.0000
F81G	15	-2077.5744	4185.7863	0.0000	1.0000
JC69I	12	-2081.2256	4186.8639	0.0000	1.0000
JC69G	12	-2081.2679	4186.9485	0.0000	1.0000
F81IG	16	-2077.6026	4187.9286	0.0000	1.0000
HKY	15	-2079.1867	4189.0109	0.0000	1.0000
JC69IG	13	-2081.3079	4189.0979	0.0000	1.0000
K2P	12	-2083.2168	4190.8463	0.0000	1.0000
F81	14	-2093.8841	4216.3253	0.0000	1.0000
JC69	11	-2097.4723	4217.2933	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2039.1229	4118.2459	0.4830	0.4830
GTRI	20	-2039.6076	4119.2152	0.2975	0.7804
GTRIG	21	-2038.9581	4119.9161	0.2095	0.9899
SYMG	17	-2046.7951	4127.5901	0.0045	0.9944
SYMI	17	-2046.9615	4127.9230	0.0038	0.9983
SYMIG	18	-2046.7549	4129.5098	0.0017	1.0000
HKYI	16	-2061.4458	4154.8916	0.0000	1.0000
HKYG	16	-2061.5580	4155.1161	0.0000	1.0000
HKYIG	17	-2061.5580	4157.1160	0.0000	1.0000
K2PI	13	-2065.6643	4157.3286	0.0000	1.0000
K2PG	13	-2065.8247	4157.6493	0.0000	1.0000
GTR	19	-2060.2928	4158.5855	0.0000	1.0000
K2PIG	14	-2065.8247	4159.6494	0.0000	1.0000
SYM	16	-2066.5698	4165.1395	0.0000	1.0000
F81I	15	-2077.5720	4185.1440	0.0000	1.0000
F81G	15	-2077.5744	4185.1488	0.0000	1.0000
JC69I	12	-2081.2256	4186.4512	0.0000	1.0000
JC69G	12	-2081.2679	4186.5358	0.0000	1.0000
F81IG	16	-2077.6026	4187.2052	0.0000	1.0000
HKY	15	-2079.1867	4188.3734	0.0000	1.0000
JC69IG	13	-2081.3079	4188.6158	0.0000	1.0000
K2P	12	-2083.2168	4190.4336	0.0000	1.0000
F81	14	-2093.8841	4215.7683	0.0000	1.0000
JC69	11	-2097.4723	4216.9445	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	17	-2046.7951	4206.5567	0.4858	0.4858
SYMI	17	-2046.9615	4206.8896	0.4113	0.8972
GTRG	20	-2039.1229	4211.1477	0.0489	0.9461
GTRI	20	-2039.6076	4212.1171	0.0301	0.9762
SYMIG	18	-2046.7549	4213.1214	0.0182	0.9945
GTRIG	21	-2038.9581	4217.4630	0.0021	0.9965
K2PI	13	-2065.6643	4217.7148	0.0018	0.9984
K2PG	13	-2065.8247	4218.0355	0.0016	0.9999
K2PIG	14	-2065.8247	4224.6806	0.0001	1.0000
HKYI	16	-2061.4458	4229.2131	0.0000	1.0000
HKYG	16	-2061.5580	4229.4375	0.0000	1.0000
HKYIG	17	-2061.5580	4236.0826	0.0000	1.0000
SYM	16	-2066.5698	4239.4610	0.0000	1.0000
JC69I	12	-2081.2256	4242.1923	0.0000	1.0000
JC69G	12	-2081.2679	4242.2769	0.0000	1.0000
K2P	12	-2083.2168	4246.1747	0.0000	1.0000
GTR	19	-2060.2928	4246.8422	0.0000	1.0000
JC69IG	13	-2081.3079	4249.0020	0.0000	1.0000
F81I	15	-2077.5720	4254.8203	0.0000	1.0000
F81G	15	-2077.5744	4254.8252	0.0000	1.0000
HKY	15	-2079.1867	4258.0498	0.0000	1.0000
F81IG	16	-2077.6026	4261.5267	0.0000	1.0000
JC69	11	-2097.4723	4268.0405	0.0000	1.0000
F81	14	-2093.8841	4280.7996	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/chr6_9783.phylip = 1 - 769;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9783.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/chr6_9783.phylip = 1 - 769;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9783.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 (SYMG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_9783.phylip = 1 - 769;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9783.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=Fixed(Equal);
END;


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

