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

Input data from file "591-tomodel-phylip/chr7_11510.phylip" (7 taxa, 987 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	-3155.9444	6351.8888	0.3728	0.3728
HKYG	16	-3160.0381	6352.0762	0.3395	0.7123
GTRIG	21	-3155.9445	6353.8890	0.1371	0.8495
HKYIG	17	-3160.0382	6354.0764	0.1249	0.9743
HKYI	16	-3163.1936	6358.3872	0.0145	0.9888
GTRI	20	-3159.4495	6358.8989	0.0112	1.0000
HKY	15	-3175.7152	6381.4304	0.0000	1.0000
GTR	19	-3172.1691	6382.3382	0.0000	1.0000
K2PG	13	-3186.7643	6399.5287	0.0000	1.0000
K2PIG	14	-3186.7644	6401.5288	0.0000	1.0000
SYMG	17	-3185.2705	6404.5411	0.0000	1.0000
K2PI	13	-3189.7741	6405.5481	0.0000	1.0000
SYMIG	18	-3185.2705	6406.5411	0.0000	1.0000
SYMI	17	-3188.4718	6410.9436	0.0000	1.0000
K2P	12	-3202.3672	6428.7345	0.0000	1.0000
SYM	16	-3201.0052	6434.0104	0.0000	1.0000
F81G	15	-3248.8506	6527.7012	0.0000	1.0000
F81IG	16	-3248.8507	6529.7014	0.0000	1.0000
F81I	15	-3251.3601	6532.7202	0.0000	1.0000
F81	14	-3262.7144	6553.4289	0.0000	1.0000
JC69G	12	-3271.4837	6566.9674	0.0000	1.0000
JC69IG	13	-3271.4838	6568.9676	0.0000	1.0000
JC69I	12	-3273.8730	6571.7460	0.0000	1.0000
JC69	11	-3285.3592	6592.7183	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3160.0381	6430.3909	0.9307	0.9307
HKYI	16	-3163.1936	6436.7019	0.0397	0.9703
HKYIG	17	-3160.0382	6437.2858	0.0296	0.9999
GTRG	20	-3155.9444	6449.7822	0.0001	1.0000
HKY	15	-3175.7152	6454.8505	0.0000	1.0000
GTRIG	21	-3155.9445	6456.6771	0.0000	1.0000
GTRI	20	-3159.4495	6456.7923	0.0000	1.0000
K2PG	13	-3186.7643	6463.1594	0.0000	1.0000
K2PI	13	-3189.7741	6469.1789	0.0000	1.0000
K2PIG	14	-3186.7644	6470.0541	0.0000	1.0000
GTR	19	-3172.1691	6475.3370	0.0000	1.0000
K2P	12	-3202.3672	6487.4705	0.0000	1.0000
SYMG	17	-3185.2705	6487.7505	0.0000	1.0000
SYMI	17	-3188.4718	6494.1530	0.0000	1.0000
SYMIG	18	-3185.2705	6494.6451	0.0000	1.0000
SYM	16	-3201.0052	6512.3251	0.0000	1.0000
F81G	15	-3248.8506	6601.1213	0.0000	1.0000
F81I	15	-3251.3601	6606.1402	0.0000	1.0000
F81IG	16	-3248.8507	6608.0161	0.0000	1.0000
F81	14	-3262.7144	6621.9543	0.0000	1.0000
JC69G	12	-3271.4837	6625.7035	0.0000	1.0000
JC69I	12	-3273.8730	6630.4820	0.0000	1.0000
JC69IG	13	-3271.4838	6632.5984	0.0000	1.0000
JC69	11	-3285.3592	6646.5597	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/chr7_11510.phylip = 1 - 987;
 Partition Dummy = 1:591-tomodel-phylip/chr7_11510.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/chr7_11510.phylip = 1 - 987;
 Partition Dummy = 1:591-tomodel-phylip/chr7_11510.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/chr7_11510.phylip = 1 - 987;
 Partition Dummy = 1:591-tomodel-phylip/chr7_11510.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

