
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:17 04/04/2011

Input data from file "683-tomodel-phylip/chr6_3227.phylip" (7 taxa, 1242 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	-4189.8123	8419.6246	0.6459	0.6459
GTRIG	21	-4189.8179	8421.6357	0.2363	0.8822
GTRI	20	-4191.5277	8423.0554	0.1162	0.9984
HKYG	16	-4200.2392	8432.4785	0.0010	0.9995
HKYIG	17	-4200.2442	8434.4884	0.0004	0.9998
HKYI	16	-4202.0884	8436.1768	0.0002	1.0000
GTR	19	-4205.8517	8449.7034	0.0000	1.0000
HKY	15	-4215.9741	8461.9482	0.0000	1.0000
K2PG	13	-4260.4850	8546.9700	0.0000	1.0000
K2PIG	14	-4260.3404	8548.6808	0.0000	1.0000
SYMG	17	-4257.6284	8549.2567	0.0000	1.0000
K2PI	13	-4262.4302	8550.8603	0.0000	1.0000
SYMIG	18	-4257.6337	8551.2674	0.0000	1.0000
SYMI	17	-4259.5030	8553.0059	0.0000	1.0000
K2P	12	-4277.4308	8578.8617	0.0000	1.0000
SYM	16	-4273.9599	8579.9197	0.0000	1.0000
F81G	15	-4290.3373	8610.6746	0.0000	1.0000
F81IG	16	-4290.3434	8612.6867	0.0000	1.0000
F81I	15	-4291.8013	8613.6026	0.0000	1.0000
F81	14	-4304.3196	8636.6392	0.0000	1.0000
JC69G	12	-4339.9852	8703.9704	0.0000	1.0000
JC69IG	13	-4339.9852	8705.9704	0.0000	1.0000
JC69I	12	-4341.6846	8707.3691	0.0000	1.0000
JC69	11	-4355.2679	8732.5357	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4200.2392	8514.4701	0.8250	0.8250
HKYI	16	-4202.0884	8518.1684	0.1298	0.9549
HKYIG	17	-4200.2442	8521.6045	0.0233	0.9782
GTRG	20	-4189.8123	8522.1141	0.0181	0.9962
GTRI	20	-4191.5277	8525.5450	0.0032	0.9995
GTRIG	21	-4189.8179	8529.2497	0.0005	1.0000
HKY	15	-4215.9741	8538.8154	0.0000	1.0000
GTR	19	-4205.8517	8547.0685	0.0000	1.0000
K2PG	13	-4260.4850	8613.5882	0.0000	1.0000
K2PI	13	-4262.4302	8617.4785	0.0000	1.0000
K2PIG	14	-4260.3404	8620.4235	0.0000	1.0000
SYMG	17	-4257.6284	8636.3729	0.0000	1.0000
SYMI	17	-4259.5030	8640.1220	0.0000	1.0000
K2P	12	-4277.4308	8640.3554	0.0000	1.0000
SYMIG	18	-4257.6337	8643.5080	0.0000	1.0000
SYM	16	-4273.9599	8661.9114	0.0000	1.0000
F81G	15	-4290.3373	8687.5418	0.0000	1.0000
F81I	15	-4291.8013	8690.4698	0.0000	1.0000
F81IG	16	-4290.3434	8694.6784	0.0000	1.0000
F81	14	-4304.3196	8708.3819	0.0000	1.0000
JC69G	12	-4339.9852	8765.4641	0.0000	1.0000
JC69I	12	-4341.6846	8768.8629	0.0000	1.0000
JC69IG	13	-4339.9852	8772.5886	0.0000	1.0000
JC69	11	-4355.2679	8788.9050	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 683-tomodel-phylip/chr6_3227.phylip = 1 - 1242;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3227.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 683-tomodel-phylip/chr6_3227.phylip = 1 - 1242;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3227.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 683-tomodel-phylip/chr6_3227.phylip = 1 - 1242;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3227.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

