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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-4109.6830	8259.3660	0.3276	0.3276
HKYG	16	-4114.0641	8260.1283	0.2238	0.5514
GTRI	20	-4110.4579	8260.9157	0.1509	0.7023
GTRIG	21	-4109.6642	8261.3285	0.1228	0.8251
HKYI	16	-4114.9250	8261.8501	0.0946	0.9197
HKYIG	17	-4114.0896	8262.1792	0.0803	1.0000
HKY	15	-4145.0802	8320.1604	0.0000	1.0000
GTR	19	-4141.4390	8320.8780	0.0000	1.0000
K2PG	13	-4165.4072	8356.8144	0.0000	1.0000
K2PI	13	-4166.1543	8358.3087	0.0000	1.0000
K2PIG	14	-4165.4402	8358.8804	0.0000	1.0000
SYMG	17	-4162.8006	8359.6013	0.0000	1.0000
SYMI	17	-4163.6683	8361.3367	0.0000	1.0000
SYMIG	18	-4162.8452	8361.6903	0.0000	1.0000
K2P	12	-4193.7223	8411.4446	0.0000	1.0000
SYM	16	-4190.8228	8413.6456	0.0000	1.0000
F81G	15	-4193.6700	8417.3400	0.0000	1.0000
F81I	15	-4194.1828	8418.3656	0.0000	1.0000
F81IG	16	-4193.3218	8418.6436	0.0000	1.0000
F81	14	-4222.1866	8472.3732	0.0000	1.0000
JC69G	12	-4233.8777	8491.7554	0.0000	1.0000
JC69I	12	-4234.1863	8492.3725	0.0000	1.0000
JC69IG	13	-4233.5058	8493.0116	0.0000	1.0000
JC69	11	-4260.3317	8542.6633	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4114.0641	8341.0822	0.6890	0.6890
HKYI	16	-4114.9250	8342.8039	0.2913	0.9803
HKYIG	17	-4114.0896	8348.1927	0.0197	0.9999
GTRG	20	-4109.6830	8360.5584	0.0000	1.0000
GTRI	20	-4110.4579	8362.1081	0.0000	1.0000
GTRIG	21	-4109.6642	8367.5805	0.0000	1.0000
HKY	15	-4145.0802	8396.0547	0.0000	1.0000
GTR	19	-4141.4390	8417.0108	0.0000	1.0000
K2PG	13	-4165.4072	8422.5894	0.0000	1.0000
K2PI	13	-4166.1543	8424.0837	0.0000	1.0000
K2PIG	14	-4165.4402	8429.7151	0.0000	1.0000
SYMG	17	-4162.8006	8445.6148	0.0000	1.0000
SYMI	17	-4163.6683	8447.3502	0.0000	1.0000
SYMIG	18	-4162.8452	8452.7635	0.0000	1.0000
K2P	12	-4193.7223	8472.1600	0.0000	1.0000
F81G	15	-4193.6700	8493.2343	0.0000	1.0000
F81I	15	-4194.1828	8494.2599	0.0000	1.0000
SYM	16	-4190.8228	8494.5995	0.0000	1.0000
F81IG	16	-4193.3218	8499.5975	0.0000	1.0000
F81	14	-4222.1866	8543.2079	0.0000	1.0000
JC69G	12	-4233.8777	8552.4708	0.0000	1.0000
JC69I	12	-4234.1863	8553.0879	0.0000	1.0000
JC69IG	13	-4233.5058	8558.7866	0.0000	1.0000
JC69	11	-4260.3317	8598.3191	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/chr1_2946.phylip = 1 - 1164;
 Partition Dummy = 1:591-tomodel-phylip/chr1_2946.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/chr1_2946.phylip = 1 - 1164;
 Partition Dummy = 1:591-tomodel-phylip/chr1_2946.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 (HKYG)]

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

