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

Input data from file "683-tomodel-phylip/chr1_15540.phylip" (7 taxa, 863 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	-2217.1454	4466.2908	0.4808	0.4808
HKYIG	17	-2217.1488	4468.2976	0.1763	0.6571
HKYI	16	-2218.2538	4468.5076	0.1587	0.8158
GTRG	20	-2214.6422	4469.2843	0.1076	0.9234
GTRIG	21	-2214.6448	4471.2896	0.0395	0.9629
GTRI	20	-2215.7808	4471.5615	0.0345	0.9974
HKY	15	-2223.6465	4477.2929	0.0020	0.9994
GTR	19	-2220.7522	4479.5043	0.0006	1.0000
F81G	15	-2260.5593	4551.1186	0.0000	1.0000
F81I	15	-2261.3037	4552.6074	0.0000	1.0000
F81IG	16	-2260.5655	4553.1309	0.0000	1.0000
F81	14	-2265.4728	4558.9457	0.0000	1.0000
SYMG	17	-2298.7123	4631.4246	0.0000	1.0000
SYMI	17	-2299.4019	4632.8038	0.0000	1.0000
SYMIG	18	-2298.7218	4633.4436	0.0000	1.0000
K2PG	13	-2304.9412	4635.8825	0.0000	1.0000
K2PIG	14	-2304.9405	4637.8809	0.0000	1.0000
K2PI	13	-2306.2797	4638.5595	0.0000	1.0000
SYM	16	-2304.2822	4640.5644	0.0000	1.0000
K2P	12	-2308.8904	4641.7808	0.0000	1.0000
JC69G	12	-2335.2681	4694.5362	0.0000	1.0000
JC69I	12	-2335.8367	4695.6734	0.0000	1.0000
JC69IG	13	-2335.2699	4696.5398	0.0000	1.0000
JC69	11	-2340.2846	4702.5691	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2217.1454	4542.4574	0.7102	0.7102
HKYI	16	-2218.2538	4544.6742	0.2344	0.9446
HKY	15	-2223.6465	4548.6991	0.0313	0.9759
HKYIG	17	-2217.1488	4549.2247	0.0241	1.0000
GTRG	20	-2214.6422	4564.4926	0.0000	1.0000
GTRI	20	-2215.7808	4566.7698	0.0000	1.0000
GTR	19	-2220.7522	4569.9522	0.0000	1.0000
GTRIG	21	-2214.6448	4571.2583	0.0000	1.0000
F81G	15	-2260.5593	4622.5248	0.0000	1.0000
F81I	15	-2261.3037	4624.0136	0.0000	1.0000
F81	14	-2265.4728	4625.5915	0.0000	1.0000
F81IG	16	-2260.5655	4629.2976	0.0000	1.0000
K2PG	13	-2304.9412	4697.7679	0.0000	1.0000
K2P	12	-2308.8904	4698.9057	0.0000	1.0000
K2PI	13	-2306.2797	4700.4449	0.0000	1.0000
K2PIG	14	-2304.9405	4704.5267	0.0000	1.0000
SYMG	17	-2298.7123	4712.3516	0.0000	1.0000
SYMI	17	-2299.4019	4713.7308	0.0000	1.0000
SYM	16	-2304.2822	4716.7310	0.0000	1.0000
SYMIG	18	-2298.7218	4719.1311	0.0000	1.0000
JC69G	12	-2335.2681	4751.6612	0.0000	1.0000
JC69I	12	-2335.8367	4752.7983	0.0000	1.0000
JC69	11	-2340.2846	4754.9337	0.0000	1.0000
JC69IG	13	-2335.2699	4758.4252	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/chr1_15540.phylip = 1 - 863;
 Partition Dummy = 1:683-tomodel-phylip/chr1_15540.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/chr1_15540.phylip = 1 - 863;
 Partition Dummy = 1:683-tomodel-phylip/chr1_15540.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/chr1_15540.phylip = 1 - 863;
 Partition Dummy = 1:683-tomodel-phylip/chr1_15540.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

