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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: SYMG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-4408.9661	8857.9321	0.7284	0.7284
GTRIG	21	-4408.9660	8859.9321	0.2680	0.9964
SYMG	17	-4417.8622	8869.7244	0.0020	0.9984
GTRI	20	-4415.6553	8871.3107	0.0009	0.9993
SYMIG	18	-4417.8622	8871.7243	0.0007	1.0000
SYMI	17	-4424.4942	8882.9885	0.0000	1.0000
HKYG	16	-4436.4527	8904.9054	0.0000	1.0000
HKYIG	17	-4436.4527	8906.9055	0.0000	1.0000
HKYI	16	-4443.1411	8918.2822	0.0000	1.0000
K2PG	13	-4448.5727	8923.1454	0.0000	1.0000
K2PIG	14	-4448.5726	8925.1452	0.0000	1.0000
K2PI	13	-4455.5052	8937.0104	0.0000	1.0000
GTR	19	-4449.5741	8937.1482	0.0000	1.0000
SYM	16	-4461.3725	8954.7451	0.0000	1.0000
HKY	15	-4475.9999	8981.9998	0.0000	1.0000
K2P	12	-4489.7456	9003.4912	0.0000	1.0000
F81G	15	-4545.5825	9121.1649	0.0000	1.0000
F81IG	16	-4545.5827	9123.1655	0.0000	1.0000
F81I	15	-4550.7072	9131.4145	0.0000	1.0000
JC69G	12	-4563.7116	9151.4232	0.0000	1.0000
JC69IG	13	-4563.7118	9153.4235	0.0000	1.0000
JC69I	12	-4568.8374	9161.6749	0.0000	1.0000
F81	14	-4577.3046	9182.6092	0.0000	1.0000
JC69	11	-4596.6296	9215.2591	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	17	-4417.8622	8956.0419	0.8232	0.8232
GTRG	20	-4408.9661	8959.4821	0.1474	0.9705
SYMIG	18	-4417.8622	8963.1193	0.0239	0.9945
GTRIG	21	-4408.9660	8966.5595	0.0043	0.9987
SYMI	17	-4424.4942	8969.3059	0.0011	0.9998
GTRI	20	-4415.6553	8972.8607	0.0002	1.0000
HKYG	16	-4436.4527	8986.1454	0.0000	1.0000
K2PG	13	-4448.5727	8989.1529	0.0000	1.0000
HKYIG	17	-4436.4527	8993.2229	0.0000	1.0000
K2PIG	14	-4448.5726	8996.2302	0.0000	1.0000
HKYI	16	-4443.1411	8999.5222	0.0000	1.0000
K2PI	13	-4455.5052	9003.0179	0.0000	1.0000
GTR	19	-4449.5741	9033.6206	0.0000	1.0000
SYM	16	-4461.3725	9035.9850	0.0000	1.0000
HKY	15	-4475.9999	9058.1623	0.0000	1.0000
K2P	12	-4489.7456	9064.4212	0.0000	1.0000
F81G	15	-4545.5825	9197.3274	0.0000	1.0000
F81IG	16	-4545.5827	9204.4054	0.0000	1.0000
F81I	15	-4550.7072	9207.5769	0.0000	1.0000
JC69G	12	-4563.7116	9212.3532	0.0000	1.0000
JC69IG	13	-4563.7118	9219.4310	0.0000	1.0000
JC69I	12	-4568.8374	9222.6049	0.0000	1.0000
F81	14	-4577.3046	9253.6942	0.0000	1.0000
JC69	11	-4596.6296	9271.1116	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_7959.phylip = 1 - 1185;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7959.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_7959.phylip = 1 - 1185;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7959.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 BIC model (SYMG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_7959.phylip = 1 - 1185;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7959.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=Fixed(Equal);
END;


-------------------------------------------------------------
End of Output

