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

Input data from file "591-tomodel-phylip/chr3_2735.phylip" (7 taxa, 849 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	-2106.8034	4245.6069	0.3919	0.3919
HKYIG	17	-2106.2350	4246.4699	0.2545	0.6464
HKYI	16	-2107.9186	4247.8371	0.1285	0.7749
GTRG	20	-2104.0332	4248.0664	0.1146	0.8895
GTRIG	21	-2103.5130	4249.0260	0.0709	0.9604
GTRI	20	-2105.0964	4250.1928	0.0396	1.0000
HKY	15	-2135.3506	4300.7012	0.0000	1.0000
GTR	19	-2132.1850	4302.3700	0.0000	1.0000
F81G	15	-2154.6858	4339.3716	0.0000	1.0000
F81IG	16	-2154.0169	4340.0338	0.0000	1.0000
F81I	15	-2155.6621	4341.3241	0.0000	1.0000
K2PG	13	-2174.9587	4375.9174	0.0000	1.0000
K2PIG	14	-2174.1548	4376.3095	0.0000	1.0000
K2PI	13	-2175.9367	4377.8733	0.0000	1.0000
SYMG	17	-2174.5840	4383.1680	0.0000	1.0000
SYMIG	18	-2173.7581	4383.5162	0.0000	1.0000
SYMI	17	-2175.5834	4385.1668	0.0000	1.0000
F81	14	-2184.1362	4396.2723	0.0000	1.0000
K2P	12	-2205.1489	4434.2979	0.0000	1.0000
SYM	16	-2204.2733	4440.5465	0.0000	1.0000
JC69G	12	-2215.1658	4454.3317	0.0000	1.0000
JC69IG	13	-2214.3828	4454.7657	0.0000	1.0000
JC69I	12	-2216.1252	4456.2504	0.0000	1.0000
JC69	11	-2244.6416	4511.2831	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2106.8034	4321.5118	0.7202	0.7202
HKYI	16	-2107.9186	4323.7421	0.2361	0.9563
HKYIG	17	-2106.2350	4327.1189	0.0436	1.0000
GTRG	20	-2104.0332	4342.9476	0.0000	1.0000
GTRI	20	-2105.0964	4345.0739	0.0000	1.0000
GTRIG	21	-2103.5130	4348.6512	0.0000	1.0000
HKY	15	-2135.3506	4371.8621	0.0000	1.0000
GTR	19	-2132.1850	4392.5071	0.0000	1.0000
F81G	15	-2154.6858	4410.5324	0.0000	1.0000
F81I	15	-2155.6621	4412.4850	0.0000	1.0000
F81IG	16	-2154.0169	4415.9387	0.0000	1.0000
K2PG	13	-2174.9587	4437.5902	0.0000	1.0000
K2PI	13	-2175.9367	4439.5461	0.0000	1.0000
K2PIG	14	-2174.1548	4442.7264	0.0000	1.0000
F81	14	-2184.1362	4462.6891	0.0000	1.0000
SYMG	17	-2174.5840	4463.8170	0.0000	1.0000
SYMI	17	-2175.5834	4465.8158	0.0000	1.0000
SYMIG	18	-2173.7581	4468.9093	0.0000	1.0000
K2P	12	-2205.1489	4491.2266	0.0000	1.0000
JC69G	12	-2215.1658	4511.2604	0.0000	1.0000
JC69I	12	-2216.1252	4513.1791	0.0000	1.0000
JC69IG	13	-2214.3828	4516.4385	0.0000	1.0000
SYM	16	-2204.2733	4516.4515	0.0000	1.0000
JC69	11	-2244.6416	4563.4678	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/chr3_2735.phylip = 1 - 849;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2735.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/chr3_2735.phylip = 1 - 849;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2735.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/chr3_2735.phylip = 1 - 849;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2735.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

