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

Input data from file "683-tomodel-phylip/chr5_11657.phylip" (7 taxa, 1112 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	-2723.9006	5479.8011	0.3391	0.3391
HKYI	16	-2724.5100	5481.0199	0.1844	0.5235
GTRG	20	-2720.6136	5481.2273	0.1662	0.6897
HKYIG	17	-2723.7130	5481.4260	0.1505	0.8402
GTRI	20	-2721.2438	5482.4876	0.0885	0.9287
GTRIG	21	-2720.4601	5482.9201	0.0713	1.0000
HKY	15	-2744.5226	5519.0453	0.0000	1.0000
GTR	19	-2741.5969	5521.1939	0.0000	1.0000
K2PG	13	-2761.2208	5548.4416	0.0000	1.0000
K2PI	13	-2761.7571	5549.5142	0.0000	1.0000
K2PIG	14	-2761.0559	5550.1118	0.0000	1.0000
SYMG	17	-2760.4665	5554.9330	0.0000	1.0000
SYMI	17	-2760.9976	5555.9952	0.0000	1.0000
SYMIG	18	-2760.3068	5556.6137	0.0000	1.0000
F81G	15	-2768.8142	5567.6284	0.0000	1.0000
F81I	15	-2769.2691	5568.5381	0.0000	1.0000
F81IG	16	-2768.7024	5569.4048	0.0000	1.0000
K2P	12	-2781.8008	5587.6017	0.0000	1.0000
SYM	16	-2781.0514	5594.1028	0.0000	1.0000
F81	14	-2787.8588	5603.7175	0.0000	1.0000
JC69G	12	-2801.6494	5627.2988	0.0000	1.0000
JC69I	12	-2801.9333	5627.8666	0.0000	1.0000
JC69IG	13	-2801.4160	5628.8319	0.0000	1.0000
JC69	11	-2820.8099	5663.6198	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2723.9006	5560.0238	0.6330	0.6330
HKYI	16	-2724.5100	5561.2426	0.3441	0.9771
HKYIG	17	-2723.7130	5566.6626	0.0229	1.0000
GTRG	20	-2720.6136	5581.5056	0.0000	1.0000
GTRI	20	-2721.2438	5582.7659	0.0000	1.0000
GTRIG	21	-2720.4601	5588.2123	0.0000	1.0000
HKY	15	-2744.5226	5594.2540	0.0000	1.0000
K2PG	13	-2761.2208	5613.6225	0.0000	1.0000
K2PI	13	-2761.7571	5614.6951	0.0000	1.0000
GTR	19	-2741.5969	5616.4583	0.0000	1.0000
K2PIG	14	-2761.0559	5620.3067	0.0000	1.0000
SYMG	17	-2760.4665	5640.1695	0.0000	1.0000
SYMI	17	-2760.9976	5641.2318	0.0000	1.0000
F81G	15	-2768.8142	5642.8372	0.0000	1.0000
F81I	15	-2769.2691	5643.7469	0.0000	1.0000
SYMIG	18	-2760.3068	5646.8641	0.0000	1.0000
K2P	12	-2781.8008	5647.7686	0.0000	1.0000
F81IG	16	-2768.7024	5649.6274	0.0000	1.0000
F81	14	-2787.8588	5673.9124	0.0000	1.0000
SYM	16	-2781.0514	5674.3254	0.0000	1.0000
JC69G	12	-2801.6494	5687.4658	0.0000	1.0000
JC69I	12	-2801.9333	5688.0336	0.0000	1.0000
JC69IG	13	-2801.4160	5694.0128	0.0000	1.0000
JC69	11	-2820.8099	5718.7729	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/chr5_11657.phylip = 1 - 1112;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11657.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/chr5_11657.phylip = 1 - 1112;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11657.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/chr5_11657.phylip = 1 - 1112;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11657.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

