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

Input data from file "591-tomodel-phylip/chr2_24859.phylip" (7 taxa, 881 characters)

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2218.9824	4477.9648	0.2519	0.2519
HKYG	16	-2223.0367	4478.0733	0.2386	0.4905
GTRI	20	-2219.3631	4478.7261	0.1722	0.6627
HKYI	16	-2223.4888	4478.9776	0.1518	0.8145
GTRIG	21	-2218.9663	4479.9327	0.0942	0.9086
HKYIG	17	-2222.9968	4479.9935	0.0913	1.0000
GTR	19	-2232.7093	4503.4187	0.0000	1.0000
HKY	15	-2236.7467	4503.4934	0.0000	1.0000
K2PG	13	-2241.0767	4508.1534	0.0000	1.0000
K2PI	13	-2241.6780	4509.3561	0.0000	1.0000
K2PIG	14	-2241.0096	4510.0193	0.0000	1.0000
SYMG	17	-2238.1384	4510.2768	0.0000	1.0000
SYMI	17	-2238.7496	4511.4992	0.0000	1.0000
SYMIG	18	-2238.0805	4512.1611	0.0000	1.0000
K2P	12	-2254.3940	4532.7881	0.0000	1.0000
SYM	16	-2251.1943	4534.3886	0.0000	1.0000
F81G	15	-2281.2014	4592.4028	0.0000	1.0000
F81I	15	-2281.6685	4593.3371	0.0000	1.0000
F81IG	16	-2281.1755	4594.3510	0.0000	1.0000
F81	14	-2293.3015	4614.6029	0.0000	1.0000
JC69G	12	-2296.4160	4616.8320	0.0000	1.0000
JC69I	12	-2296.9498	4617.8997	0.0000	1.0000
JC69IG	13	-2296.3793	4618.7586	0.0000	1.0000
JC69	11	-2308.4901	4638.9802	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2223.0367	4554.5702	0.5980	0.5980
HKYI	16	-2223.4888	4555.4745	0.3805	0.9786
HKYIG	17	-2222.9968	4561.2715	0.0210	0.9995
K2PG	13	-2241.0767	4570.3071	0.0002	0.9998
K2PI	13	-2241.6780	4571.5098	0.0001	0.9999
GTRG	20	-2218.9824	4573.5860	0.0000	0.9999
GTRI	20	-2219.3631	4574.3473	0.0000	1.0000
HKY	15	-2236.7467	4575.2093	0.0000	1.0000
K2PIG	14	-2241.0096	4576.9541	0.0000	1.0000
GTRIG	21	-2218.9663	4580.3349	0.0000	1.0000
K2P	12	-2254.3940	4590.1608	0.0000	1.0000
SYMG	17	-2238.1384	4591.5548	0.0000	1.0000
SYMI	17	-2238.7496	4592.7771	0.0000	1.0000
GTR	19	-2232.7093	4594.2588	0.0000	1.0000
SYMIG	18	-2238.0805	4598.2201	0.0000	1.0000
SYM	16	-2251.1943	4610.8856	0.0000	1.0000
F81G	15	-2281.2014	4664.1186	0.0000	1.0000
F81I	15	-2281.6685	4665.0529	0.0000	1.0000
F81IG	16	-2281.1755	4670.8479	0.0000	1.0000
JC69G	12	-2296.4160	4674.2047	0.0000	1.0000
JC69I	12	-2296.9498	4675.2724	0.0000	1.0000
JC69IG	13	-2296.3793	4680.9124	0.0000	1.0000
F81	14	-2293.3015	4681.5377	0.0000	1.0000
JC69	11	-2308.4901	4691.5719	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_24859.phylip = 1 - 881;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24859.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24859.phylip = 1 - 881;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24859.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/chr2_24859.phylip = 1 - 881;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24859.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

