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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2187.8084	4415.6167	0.4764	0.4764
GTRI	20	-2188.1453	4416.2905	0.3401	0.8165
GTRIG	21	-2187.7625	4417.5251	0.1835	1.0000
HKYG	16	-2203.2306	4438.4612	0.0000	1.0000
HKYI	16	-2203.7375	4439.4750	0.0000	1.0000
HKYIG	17	-2202.9549	4439.9098	0.0000	1.0000
GTR	19	-2214.8496	4467.6992	0.0000	1.0000
HKY	15	-2233.0384	4496.0767	0.0000	1.0000
SYMG	17	-2233.3315	4500.6629	0.0000	1.0000
SYMI	17	-2233.5553	4501.1105	0.0000	1.0000
SYMIG	18	-2233.2825	4502.5651	0.0000	1.0000
K2PG	13	-2242.2768	4510.5537	0.0000	1.0000
K2PIG	14	-2241.6087	4511.2174	0.0000	1.0000
K2PI	13	-2242.6736	4511.3472	0.0000	1.0000
F81G	15	-2247.8093	4525.6187	0.0000	1.0000
F81I	15	-2247.9146	4525.8293	0.0000	1.0000
F81IG	16	-2247.8094	4527.6187	0.0000	1.0000
SYM	16	-2260.3988	4552.7976	0.0000	1.0000
K2P	12	-2272.0658	4568.1316	0.0000	1.0000
F81	14	-2277.1293	4582.2585	0.0000	1.0000
JC69G	12	-2281.4249	4586.8498	0.0000	1.0000
JC69I	12	-2281.5312	4587.0624	0.0000	1.0000
JC69IG	13	-2281.4360	4588.8720	0.0000	1.0000
JC69	11	-2310.7198	4643.4395	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2187.8084	4511.3736	0.4978	0.4978
GTRI	20	-2188.1453	4512.0474	0.3554	0.8532
HKYG	16	-2203.2306	4515.0667	0.0785	0.9317
HKYI	16	-2203.7375	4516.0805	0.0473	0.9790
GTRIG	21	-2187.7625	4518.0698	0.0175	0.9965
HKYIG	17	-2202.9549	4521.3032	0.0035	1.0000
GTR	19	-2214.8496	4558.6683	0.0000	1.0000
HKY	15	-2233.0384	4567.8944	0.0000	1.0000
K2PG	13	-2242.2768	4572.7956	0.0000	1.0000
K2PI	13	-2242.6736	4573.5892	0.0000	1.0000
K2PIG	14	-2241.6087	4578.2472	0.0000	1.0000
SYMG	17	-2233.3315	4582.0563	0.0000	1.0000
SYMI	17	-2233.5553	4582.5039	0.0000	1.0000
SYMIG	18	-2233.2825	4588.7463	0.0000	1.0000
F81G	15	-2247.8093	4597.4364	0.0000	1.0000
F81I	15	-2247.9146	4597.6469	0.0000	1.0000
F81IG	16	-2247.8094	4604.2243	0.0000	1.0000
K2P	12	-2272.0658	4625.5858	0.0000	1.0000
SYM	16	-2260.3988	4629.4031	0.0000	1.0000
JC69G	12	-2281.4249	4644.3039	0.0000	1.0000
JC69I	12	-2281.5312	4644.5166	0.0000	1.0000
F81	14	-2277.1293	4649.2883	0.0000	1.0000
JC69IG	13	-2281.4360	4651.1140	0.0000	1.0000
JC69	11	-2310.7198	4696.1058	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/chr12_5983.phylip = 1 - 887;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5983.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/chr12_5983.phylip = 1 - 887;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5983.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 (GTRG)]

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


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

