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

Input data from file "591-tomodel-phylip/chr6_5749.phylip" (7 taxa, 1132 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	-2979.1446	5990.2891	0.5432	0.5432
HKYI	16	-2980.0034	5992.0068	0.2301	0.7734
HKYIG	17	-2979.1469	5992.2937	0.1994	0.9727
GTRG	20	-2978.7119	5997.4239	0.0153	0.9881
GTRI	20	-2979.6012	5999.2024	0.0063	0.9944
GTRIG	21	-2978.7141	5999.4283	0.0056	1.0000
HKY	15	-3006.2319	6042.4637	0.0000	1.0000
GTR	19	-3005.8638	6049.7275	0.0000	1.0000
F81G	15	-3047.1589	6124.3177	0.0000	1.0000
F81I	15	-3047.7155	6125.4311	0.0000	1.0000
F81IG	16	-3047.1671	6126.3343	0.0000	1.0000
F81	14	-3074.3704	6176.7408	0.0000	1.0000
K2PG	13	-3083.9175	6193.8349	0.0000	1.0000
K2PI	13	-3084.6860	6195.3721	0.0000	1.0000
K2PIG	14	-3083.8931	6195.7862	0.0000	1.0000
SYMG	17	-3081.0285	6196.0570	0.0000	1.0000
SYMIG	18	-3080.9462	6197.8925	0.0000	1.0000
SYMI	17	-3081.9535	6197.9070	0.0000	1.0000
K2P	12	-3114.1076	6252.2151	0.0000	1.0000
SYM	16	-3112.4768	6256.9536	0.0000	1.0000
JC69G	12	-3141.4418	6306.8837	0.0000	1.0000
JC69I	12	-3142.1452	6308.2903	0.0000	1.0000
JC69IG	13	-3141.4275	6308.8549	0.0000	1.0000
JC69	11	-3170.5334	6363.0667	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2979.1446	6070.7970	0.6881	0.6881
HKYI	16	-2980.0034	6072.5146	0.2915	0.9796
HKYIG	17	-2979.1469	6077.8333	0.0204	1.0000
GTRG	20	-2978.7119	6098.0587	0.0000	1.0000
GTRI	20	-2979.6012	6099.8373	0.0000	1.0000
GTRIG	21	-2978.7141	6105.0949	0.0000	1.0000
HKY	15	-3006.2319	6117.9398	0.0000	1.0000
GTR	19	-3005.8638	6145.3306	0.0000	1.0000
F81G	15	-3047.1589	6199.7938	0.0000	1.0000
F81I	15	-3047.7155	6200.9072	0.0000	1.0000
F81IG	16	-3047.1671	6206.8421	0.0000	1.0000
F81	14	-3074.3704	6247.1852	0.0000	1.0000
K2PG	13	-3083.9175	6259.2476	0.0000	1.0000
K2PI	13	-3084.6860	6260.7847	0.0000	1.0000
K2PIG	14	-3083.8931	6266.2306	0.0000	1.0000
SYMG	17	-3081.0285	6281.5966	0.0000	1.0000
SYMI	17	-3081.9535	6283.4466	0.0000	1.0000
SYMIG	18	-3080.9462	6288.4638	0.0000	1.0000
K2P	12	-3114.1076	6312.5960	0.0000	1.0000
SYM	16	-3112.4768	6337.4614	0.0000	1.0000
JC69G	12	-3141.4418	6367.2646	0.0000	1.0000
JC69I	12	-3142.1452	6368.6712	0.0000	1.0000
JC69IG	13	-3141.4275	6374.2676	0.0000	1.0000
JC69	11	-3170.5334	6418.4159	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 591-tomodel-phylip/chr6_5749.phylip = 1 - 1132;
 Partition Dummy = 1:591-tomodel-phylip/chr6_5749.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 591-tomodel-phylip/chr6_5749.phylip = 1 - 1132;
 Partition Dummy = 1:591-tomodel-phylip/chr6_5749.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 591-tomodel-phylip/chr6_5749.phylip = 1 - 1132;
 Partition Dummy = 1:591-tomodel-phylip/chr6_5749.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

