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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2118.7596	4269.5191	0.5272	0.5272
GTRI	20	-2115.8277	4271.6555	0.1812	0.7083
GTRG	20	-2116.0738	4272.1476	0.1416	0.8500
GTRIG	21	-2115.7519	4273.5037	0.0719	0.9219
HKYI	16	-2121.0923	4274.1845	0.0512	0.9730
HKYIG	17	-2120.7528	4275.5057	0.0264	0.9995
GTR	19	-2123.0636	4284.1271	0.0004	0.9998
HKY	15	-2127.7365	4285.4730	0.0002	1.0000
F81G	15	-2140.2725	4310.5451	0.0000	1.0000
F81I	15	-2140.6804	4311.3608	0.0000	1.0000
F81IG	16	-2140.0206	4312.0412	0.0000	1.0000
F81	14	-2147.7255	4323.4510	0.0000	1.0000
K2PG	13	-2149.8899	4325.7799	0.0000	1.0000
K2PI	13	-2151.3010	4328.6019	0.0000	1.0000
SYMI	17	-2147.7926	4329.5853	0.0000	1.0000
SYMG	17	-2147.9509	4329.9019	0.0000	1.0000
K2PIG	14	-2150.9986	4329.9973	0.0000	1.0000
SYMIG	18	-2147.1628	4330.3256	0.0000	1.0000
K2P	12	-2159.4988	4342.9975	0.0000	1.0000
SYM	16	-2156.7197	4345.4394	0.0000	1.0000
JC69G	12	-2169.7225	4363.4449	0.0000	1.0000
JC69I	12	-2169.8564	4363.7127	0.0000	1.0000
JC69IG	13	-2169.3235	4364.6469	0.0000	1.0000
JC69	11	-2177.5476	4377.0952	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2118.7596	4349.9703	0.9047	0.9047
HKYI	16	-2121.0923	4354.6358	0.0878	0.9925
HKY	15	-2127.7365	4360.8960	0.0038	0.9963
HKYIG	17	-2120.7528	4360.9851	0.0037	1.0000
GTRI	20	-2115.8277	4372.2195	0.0000	1.0000
GTRG	20	-2116.0738	4372.7116	0.0000	1.0000
GTRIG	21	-2115.7519	4379.0960	0.0000	1.0000
GTR	19	-2123.0636	4379.6629	0.0000	1.0000
F81G	15	-2140.2725	4385.9681	0.0000	1.0000
F81I	15	-2140.6804	4386.7838	0.0000	1.0000
K2PG	13	-2149.8899	4391.1465	0.0000	1.0000
F81IG	16	-2140.0206	4392.4924	0.0000	1.0000
F81	14	-2147.7255	4393.8458	0.0000	1.0000
K2PI	13	-2151.3010	4393.9685	0.0000	1.0000
K2PIG	14	-2150.9986	4400.3921	0.0000	1.0000
K2P	12	-2159.4988	4403.3359	0.0000	1.0000
SYMI	17	-2147.7926	4415.0647	0.0000	1.0000
SYMG	17	-2147.9509	4415.3813	0.0000	1.0000
SYMIG	18	-2147.1628	4420.8332	0.0000	1.0000
JC69G	12	-2169.7225	4423.7833	0.0000	1.0000
JC69I	12	-2169.8564	4424.0511	0.0000	1.0000
SYM	16	-2156.7197	4425.8906	0.0000	1.0000
JC69IG	13	-2169.3235	4430.0135	0.0000	1.0000
JC69	11	-2177.5476	4432.4054	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_14674.phylip = 1 - 1128;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_14674.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_14674.phylip = 1 - 1128;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_14674.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_14674.phylip = 1 - 1128;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_14674.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

