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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2116.6236	4273.2473	0.3801	0.3801
GTRIG	21	-2115.7728	4273.5457	0.3274	0.7075
GTRI	20	-2116.9306	4273.8613	0.2796	0.9870
HKYG	16	-2125.0299	4282.0598	0.0046	0.9917
HKYIG	17	-2124.0886	4282.1772	0.0044	0.9961
HKYI	16	-2125.1912	4282.3823	0.0039	1.0000
GTR	19	-2142.4016	4322.8033	0.0000	1.0000
SYMIG	18	-2149.7993	4335.5987	0.0000	1.0000
SYMG	17	-2150.8255	4335.6511	0.0000	1.0000
SYMI	17	-2150.9000	4335.8000	0.0000	1.0000
HKY	15	-2152.9700	4335.9400	0.0000	1.0000
K2PIG	14	-2156.2643	4340.5285	0.0000	1.0000
K2PI	13	-2157.3585	4340.7171	0.0000	1.0000
K2PG	13	-2157.3680	4340.7359	0.0000	1.0000
F81G	15	-2156.9989	4343.9978	0.0000	1.0000
F81I	15	-2157.1254	4344.2508	0.0000	1.0000
F81IG	16	-2158.1730	4348.3460	0.0000	1.0000
SYM	16	-2178.9741	4389.9481	0.0000	1.0000
JC69G	12	-2187.0233	4398.0466	0.0000	1.0000
JC69I	12	-2187.0710	4398.1420	0.0000	1.0000
K2P	12	-2187.3307	4398.6613	0.0000	1.0000
F81	14	-2186.9378	4401.8756	0.0000	1.0000
JC69IG	13	-2188.1084	4402.2168	0.0000	1.0000
JC69	11	-2218.0675	4458.1350	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2125.0299	4359.6279	0.5145	0.5145
HKYI	16	-2125.1912	4359.9504	0.4379	0.9523
HKYIG	17	-2124.0886	4364.5933	0.0430	0.9953
GTRG	20	-2116.6236	4370.2074	0.0026	0.9979
GTRI	20	-2116.9306	4370.8214	0.0019	0.9998
GTRIG	21	-2115.7728	4375.3538	0.0002	1.0000
K2PI	13	-2157.3585	4403.7411	0.0000	1.0000
K2PG	13	-2157.3680	4403.7600	0.0000	1.0000
K2PIG	14	-2156.2643	4408.4006	0.0000	1.0000
HKY	15	-2152.9700	4408.6600	0.0000	1.0000
GTR	19	-2142.4016	4414.9154	0.0000	1.0000
F81G	15	-2156.9989	4416.7179	0.0000	1.0000
F81I	15	-2157.1254	4416.9708	0.0000	1.0000
SYMG	17	-2150.8255	4418.0672	0.0000	1.0000
SYMI	17	-2150.9000	4418.2161	0.0000	1.0000
SYMIG	18	-2149.7993	4422.8628	0.0000	1.0000
F81IG	16	-2158.1730	4425.9141	0.0000	1.0000
JC69G	12	-2187.0233	4456.2227	0.0000	1.0000
JC69I	12	-2187.0710	4456.3181	0.0000	1.0000
K2P	12	-2187.3307	4456.8374	0.0000	1.0000
JC69IG	13	-2188.1084	4465.2408	0.0000	1.0000
SYM	16	-2178.9741	4467.5162	0.0000	1.0000
F81	14	-2186.9378	4469.7476	0.0000	1.0000
JC69	11	-2218.0675	4511.4631	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 683-tomodel-phylip/chrZ_8113.phylip = 1 - 942;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_8113.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 683-tomodel-phylip/chrZ_8113.phylip = 1 - 942;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_8113.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 683-tomodel-phylip/chrZ_8113.phylip = 1 - 942;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_8113.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

