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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-4227.2761	8486.5522	0.5563	0.5563
HKYIG	17	-4227.2767	8488.5533	0.2045	0.7608
HKYG	16	-4228.6874	8489.3748	0.1356	0.8964
GTRI	20	-4225.4096	8490.8192	0.0659	0.9623
GTRIG	21	-4225.4153	8492.8306	0.0241	0.9864
GTRG	20	-4226.9866	8493.9732	0.0136	1.0000
HKY	15	-4259.7532	8549.5065	0.0000	1.0000
GTR	19	-4257.6191	8553.2381	0.0000	1.0000
SYMI	17	-4266.7324	8567.4647	0.0000	1.0000
K2PI	13	-4270.7532	8567.5065	0.0000	1.0000
SYMG	17	-4267.4152	8568.8304	0.0000	1.0000
SYMIG	18	-4266.6725	8569.3450	0.0000	1.0000
K2PIG	14	-4270.7598	8569.5197	0.0000	1.0000
K2PG	13	-4272.3567	8570.7134	0.0000	1.0000
SYM	16	-4300.1628	8632.3255	0.0000	1.0000
K2P	12	-4304.9571	8633.9142	0.0000	1.0000
F81I	15	-4302.2010	8634.4021	0.0000	1.0000
F81IG	16	-4302.2211	8636.4421	0.0000	1.0000
F81G	15	-4304.2906	8638.5813	0.0000	1.0000
F81	14	-4332.9734	8693.9469	0.0000	1.0000
JC69I	12	-4357.5767	8739.1534	0.0000	1.0000
JC69IG	13	-4357.5970	8741.1940	0.0000	1.0000
JC69G	12	-4359.4446	8742.8892	0.0000	1.0000
JC69	11	-4387.0833	8796.1667	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-4227.2761	8567.8998	0.7857	0.7857
HKYG	16	-4228.6874	8570.7225	0.1916	0.9773
HKYIG	17	-4227.2767	8574.9852	0.0227	1.0000
GTRI	20	-4225.4096	8592.5037	0.0000	1.0000
GTRG	20	-4226.9866	8595.6577	0.0000	1.0000
GTRIG	21	-4225.4153	8599.5994	0.0000	1.0000
HKY	15	-4259.7532	8625.7699	0.0000	1.0000
K2PI	13	-4270.7532	8633.6014	0.0000	1.0000
K2PG	13	-4272.3567	8636.8084	0.0000	1.0000
K2PIG	14	-4270.7598	8640.6988	0.0000	1.0000
GTR	19	-4257.6191	8649.8384	0.0000	1.0000
SYMI	17	-4266.7324	8653.8966	0.0000	1.0000
SYMG	17	-4267.4152	8655.2622	0.0000	1.0000
SYMIG	18	-4266.6725	8660.8611	0.0000	1.0000
K2P	12	-4304.9571	8694.9249	0.0000	1.0000
F81I	15	-4302.2010	8710.6655	0.0000	1.0000
SYM	16	-4300.1628	8713.6732	0.0000	1.0000
F81G	15	-4304.2906	8714.8447	0.0000	1.0000
F81IG	16	-4302.2211	8717.7898	0.0000	1.0000
F81	14	-4332.9734	8765.1260	0.0000	1.0000
JC69I	12	-4357.5767	8800.1641	0.0000	1.0000
JC69G	12	-4359.4446	8803.8999	0.0000	1.0000
JC69IG	13	-4357.5970	8807.2889	0.0000	1.0000
JC69	11	-4387.0833	8852.0932	0.0000	1.0000



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

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

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


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

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


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

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

