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

Input data from file "683-tomodel-phylip/chrZ_7997.phylip" (7 taxa, 979 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	-2731.5491	5495.0981	0.5928	0.5928
HKYIG	17	-2731.4500	5496.9000	0.2408	0.8335
HKYI	16	-2732.8572	5497.7143	0.1602	0.9938
GTRI	20	-2732.3006	5504.6012	0.0051	0.9989
GTRG	20	-2734.1782	5508.3564	0.0008	0.9997
GTRIG	21	-2734.0807	5510.1613	0.0003	1.0000
HKY	15	-2756.6042	5543.2083	0.0000	1.0000
GTR	19	-2761.2763	5560.5526	0.0000	1.0000
K2PG	13	-2771.9703	5569.9406	0.0000	1.0000
K2PIG	14	-2771.8776	5571.7552	0.0000	1.0000
K2PI	13	-2773.1763	5572.3526	0.0000	1.0000
SYMG	17	-2772.5970	5579.1939	0.0000	1.0000
SYMIG	18	-2772.4973	5580.9947	0.0000	1.0000
SYMI	17	-2773.7739	5581.5478	0.0000	1.0000
F81G	15	-2788.2452	5606.4903	0.0000	1.0000
F81IG	16	-2788.1757	5608.3514	0.0000	1.0000
F81I	15	-2789.3268	5608.6535	0.0000	1.0000
K2P	12	-2797.8605	5619.7210	0.0000	1.0000
SYM	16	-2800.4664	5632.9327	0.0000	1.0000
F81	14	-2814.1014	5656.2028	0.0000	1.0000
JC69G	12	-2826.0871	5676.1741	0.0000	1.0000
JC69IG	13	-2826.0139	5678.0277	0.0000	1.0000
JC69I	12	-2827.1520	5678.3040	0.0000	1.0000
JC69	11	-2851.6922	5725.3844	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2731.5491	5573.2826	0.7659	0.7659
HKYI	16	-2732.8572	5575.8988	0.2070	0.9730
HKYIG	17	-2731.4500	5579.9710	0.0270	1.0000
GTRI	20	-2732.3006	5602.3319	0.0000	1.0000
GTRG	20	-2734.1782	5606.0870	0.0000	1.0000
GTRIG	21	-2734.0807	5612.7785	0.0000	1.0000
HKY	15	-2756.6042	5616.5063	0.0000	1.0000
K2PG	13	-2771.9703	5633.4656	0.0000	1.0000
K2PI	13	-2773.1763	5635.8775	0.0000	1.0000
K2PIG	14	-2771.8776	5640.1666	0.0000	1.0000
GTR	19	-2761.2763	5653.3967	0.0000	1.0000
SYMG	17	-2772.5970	5662.2649	0.0000	1.0000
SYMI	17	-2773.7739	5664.6189	0.0000	1.0000
SYMIG	18	-2772.4973	5668.9522	0.0000	1.0000
K2P	12	-2797.8605	5678.3594	0.0000	1.0000
F81G	15	-2788.2452	5679.7883	0.0000	1.0000
F81I	15	-2789.3268	5681.9515	0.0000	1.0000
F81IG	16	-2788.1757	5686.5359	0.0000	1.0000
SYM	16	-2800.4664	5711.1172	0.0000	1.0000
F81	14	-2814.1014	5724.6143	0.0000	1.0000
JC69G	12	-2826.0871	5734.8125	0.0000	1.0000
JC69I	12	-2827.1520	5736.9424	0.0000	1.0000
JC69IG	13	-2826.0139	5741.5526	0.0000	1.0000
JC69	11	-2851.6922	5779.1362	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_7997.phylip = 1 - 979;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_7997.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_7997.phylip = 1 - 979;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_7997.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_7997.phylip = 1 - 979;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_7997.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

