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

Input data from file "591-tomodel-phylip/chr2_17818.phylip" (7 taxa, 1037 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	-3221.7827	6483.5653	0.4379	0.4379
GTRI	20	-3222.7553	6485.5106	0.1656	0.6034
GTRIG	21	-3221.7613	6485.5225	0.1646	0.7680
HKYG	16	-3227.0080	6486.0159	0.1286	0.8966
HKYI	16	-3227.8596	6487.7192	0.0549	0.9514
HKYIG	17	-3226.9819	6487.9638	0.0486	1.0000
K2PG	13	-3257.0118	6540.0236	0.0000	1.0000
K2PIG	14	-3256.9587	6541.9173	0.0000	1.0000
K2PI	13	-3257.9880	6541.9761	0.0000	1.0000
SYMG	17	-3255.0931	6544.1861	0.0000	1.0000
SYMIG	18	-3255.0475	6546.0951	0.0000	1.0000
SYMI	17	-3256.3867	6546.7734	0.0000	1.0000
GTR	19	-3254.6164	6547.2328	0.0000	1.0000
HKY	15	-3259.5323	6549.0645	0.0000	1.0000
K2P	12	-3290.5203	6605.0407	0.0000	1.0000
SYM	16	-3288.7145	6609.4290	0.0000	1.0000
F81G	15	-3292.4533	6614.9066	0.0000	1.0000
F81I	15	-3292.9730	6615.9460	0.0000	1.0000
F81IG	16	-3292.3534	6616.7067	0.0000	1.0000
JC69G	12	-3320.0754	6664.1507	0.0000	1.0000
JC69I	12	-3320.7576	6665.5152	0.0000	1.0000
JC69IG	13	-3319.9941	6665.9883	0.0000	1.0000
F81	14	-3322.6515	6673.3030	0.0000	1.0000
JC69	11	-3350.7909	6723.5818	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3227.0080	6565.1213	0.6855	0.6855
HKYI	16	-3227.8596	6566.8246	0.2925	0.9780
HKYIG	17	-3226.9819	6572.0133	0.0218	0.9998
GTRG	20	-3221.7827	6582.4470	0.0001	1.0000
GTRI	20	-3222.7553	6584.3923	0.0000	1.0000
GTRIG	21	-3221.7613	6589.3483	0.0000	1.0000
K2PG	13	-3257.0118	6604.2967	0.0000	1.0000
K2PI	13	-3257.9880	6606.2492	0.0000	1.0000
K2PIG	14	-3256.9587	6611.1345	0.0000	1.0000
HKY	15	-3259.5323	6623.2258	0.0000	1.0000
SYMG	17	-3255.0931	6628.2356	0.0000	1.0000
SYMI	17	-3256.3867	6630.8229	0.0000	1.0000
SYMIG	18	-3255.0475	6635.0886	0.0000	1.0000
GTR	19	-3254.6164	6641.1705	0.0000	1.0000
K2P	12	-3290.5203	6664.3697	0.0000	1.0000
SYM	16	-3288.7145	6688.5344	0.0000	1.0000
F81G	15	-3292.4533	6689.0679	0.0000	1.0000
F81I	15	-3292.9730	6690.1073	0.0000	1.0000
F81IG	16	-3292.3534	6695.8121	0.0000	1.0000
JC69G	12	-3320.0754	6723.4798	0.0000	1.0000
JC69I	12	-3320.7576	6724.8442	0.0000	1.0000
JC69IG	13	-3319.9941	6730.2614	0.0000	1.0000
F81	14	-3322.6515	6742.5203	0.0000	1.0000
JC69	11	-3350.7909	6777.9668	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/chr2_17818.phylip = 1 - 1037;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17818.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/chr2_17818.phylip = 1 - 1037;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17818.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/chr2_17818.phylip = 1 - 1037;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17818.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

