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

Input data from file "683-tomodel-phylip/chr19_2703.phylip" (7 taxa, 1133 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	-2192.6347	4417.2693	0.3692	0.3692
HKYI	16	-2192.7886	4417.5771	0.3165	0.6856
HKYIG	17	-2192.9232	4419.8463	0.1018	0.7874
GTRI	20	-2190.0603	4420.1207	0.0887	0.8761
GTRG	20	-2190.2230	4420.4460	0.0754	0.9515
GTRIG	21	-2190.2551	4422.5101	0.0269	0.9784
HKY	15	-2196.6944	4423.3888	0.0173	0.9957
GTR	19	-2194.0908	4426.1816	0.0043	1.0000
F81G	15	-2220.6229	4471.2458	0.0000	1.0000
F81I	15	-2220.6812	4471.3624	0.0000	1.0000
F81IG	16	-2220.8039	4473.6079	0.0000	1.0000
F81	14	-2224.4502	4476.9005	0.0000	1.0000
K2PG	13	-2247.3100	4520.6200	0.0000	1.0000
K2PI	13	-2247.5849	4521.1697	0.0000	1.0000
SYMI	17	-2244.3139	4522.6277	0.0000	1.0000
SYMG	17	-2244.5606	4523.1213	0.0000	1.0000
K2PIG	14	-2247.6879	4523.3758	0.0000	1.0000
SYMIG	18	-2244.4426	4524.8851	0.0000	1.0000
K2P	12	-2251.7626	4527.5252	0.0000	1.0000
SYM	16	-2248.7874	4529.5749	0.0000	1.0000
JC69G	12	-2272.5739	4569.1479	0.0000	1.0000
JC69I	12	-2272.7870	4569.5740	0.0000	1.0000
JC69IG	13	-2272.8919	4571.7837	0.0000	1.0000
JC69	11	-2276.3053	4574.6105	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2192.6347	4497.7913	0.4064	0.4064
HKYI	16	-2192.7886	4498.0991	0.3485	0.7549
HKY	15	-2196.6944	4498.8782	0.2360	0.9909
HKYIG	17	-2192.9232	4505.4010	0.0090	1.0000
GTRI	20	-2190.0603	4520.7732	0.0000	1.0000
GTRG	20	-2190.2230	4521.0985	0.0000	1.0000
GTR	19	-2194.0908	4521.8015	0.0000	1.0000
GTRIG	21	-2190.2551	4528.1952	0.0000	1.0000
F81G	15	-2220.6229	4546.7352	0.0000	1.0000
F81I	15	-2220.6812	4546.8518	0.0000	1.0000
F81	14	-2224.4502	4547.3572	0.0000	1.0000
F81IG	16	-2220.8039	4554.1299	0.0000	1.0000
K2PG	13	-2247.3100	4586.0441	0.0000	1.0000
K2PI	13	-2247.5849	4586.5939	0.0000	1.0000
K2P	12	-2251.7626	4587.9167	0.0000	1.0000
K2PIG	14	-2247.6879	4593.8326	0.0000	1.0000
SYMI	17	-2244.3139	4608.1824	0.0000	1.0000
SYMG	17	-2244.5606	4608.6759	0.0000	1.0000
SYM	16	-2248.7874	4610.0969	0.0000	1.0000
SYMIG	18	-2244.4426	4615.4724	0.0000	1.0000
JC69G	12	-2272.5739	4629.5394	0.0000	1.0000
JC69I	12	-2272.7870	4629.9655	0.0000	1.0000
JC69	11	-2276.3053	4629.9694	0.0000	1.0000
JC69IG	13	-2272.8919	4637.2079	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/chr19_2703.phylip = 1 - 1133;
 Partition Dummy = 1:683-tomodel-phylip/chr19_2703.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/chr19_2703.phylip = 1 - 1133;
 Partition Dummy = 1:683-tomodel-phylip/chr19_2703.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/chr19_2703.phylip = 1 - 1133;
 Partition Dummy = 1:683-tomodel-phylip/chr19_2703.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

