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

Input data from file "591-tomodel-phylip/chrZ_11555.phylip" (7 taxa, 859 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	-2276.6251	4593.2503	0.5655	0.5655
GTRIG	21	-2276.5435	4595.0871	0.2257	0.7913
GTRI	20	-2277.8893	4595.7786	0.1597	0.9510
HKYG	16	-2283.6644	4599.3288	0.0271	0.9781
HKYI	16	-2284.5259	4601.0518	0.0114	0.9895
HKYIG	17	-2283.6150	4601.2300	0.0105	1.0000
GTR	19	-2293.7325	4625.4649	0.0000	1.0000
HKY	15	-2299.8809	4629.7618	0.0000	1.0000
K2PG	13	-2303.6399	4633.2797	0.0000	1.0000
K2PI	13	-2304.2967	4634.5935	0.0000	1.0000
K2PIG	14	-2303.5975	4635.1951	0.0000	1.0000
SYMG	17	-2301.2708	4636.5417	0.0000	1.0000
SYMIG	18	-2301.1892	4638.3784	0.0000	1.0000
SYMI	17	-2302.2996	4638.5992	0.0000	1.0000
K2P	12	-2320.5177	4665.0353	0.0000	1.0000
SYM	16	-2318.3382	4668.6765	0.0000	1.0000
F81G	15	-2349.3524	4728.7048	0.0000	1.0000
F81I	15	-2350.1264	4730.2527	0.0000	1.0000
F81IG	16	-2349.3203	4730.6405	0.0000	1.0000
F81	14	-2364.5323	4757.0646	0.0000	1.0000
JC69G	12	-2367.6055	4759.2111	0.0000	1.0000
JC69I	12	-2368.3269	4760.6538	0.0000	1.0000
JC69IG	13	-2367.5746	4761.1491	0.0000	1.0000
JC69	11	-2383.1874	4788.3748	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2283.6644	4675.4211	0.6847	0.6847
HKYI	16	-2284.5259	4677.1441	0.2893	0.9740
HKYIG	17	-2283.6150	4682.0781	0.0245	0.9985
GTRG	20	-2276.6251	4688.3657	0.0011	0.9996
GTRI	20	-2277.8893	4690.8940	0.0003	0.9999
GTRIG	21	-2276.5435	4694.9582	0.0000	0.9999
K2PG	13	-2303.6399	4695.1047	0.0000	1.0000
K2PI	13	-2304.2967	4696.4185	0.0000	1.0000
HKY	15	-2299.8809	4701.0984	0.0000	1.0000
K2PIG	14	-2303.5975	4701.7758	0.0000	1.0000
GTR	19	-2293.7325	4715.8245	0.0000	1.0000
SYMG	17	-2301.2708	4717.3897	0.0000	1.0000
SYMI	17	-2302.2996	4719.4473	0.0000	1.0000
K2P	12	-2320.5177	4722.1045	0.0000	1.0000
SYMIG	18	-2301.1892	4723.9823	0.0000	1.0000
SYM	16	-2318.3382	4744.7688	0.0000	1.0000
F81G	15	-2349.3524	4800.0414	0.0000	1.0000
F81I	15	-2350.1264	4801.5893	0.0000	1.0000
F81IG	16	-2349.3203	4806.7328	0.0000	1.0000
JC69G	12	-2367.6055	4816.2803	0.0000	1.0000
JC69I	12	-2368.3269	4817.7231	0.0000	1.0000
JC69IG	13	-2367.5746	4822.9741	0.0000	1.0000
F81	14	-2364.5323	4823.6453	0.0000	1.0000
JC69	11	-2383.1874	4840.6883	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/chrZ_11555.phylip = 1 - 859;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11555.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/chrZ_11555.phylip = 1 - 859;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11555.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 591-tomodel-phylip/chrZ_11555.phylip = 1 - 859;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11555.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

