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

Input data from file "591-tomodel-phylip/chr1_28829.phylip" (7 taxa, 908 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	-2585.6540	5211.3080	0.4298	0.4298
GTRI	20	-2585.6958	5211.3916	0.4122	0.8419
GTRIG	21	-2585.6540	5213.3081	0.1581	1.0000
HKYI	16	-2601.2833	5234.5666	0.0000	1.0000
HKYG	16	-2601.4122	5234.8244	0.0000	1.0000
HKYIG	17	-2601.4123	5236.8246	0.0000	1.0000
SYMI	17	-2615.7056	5265.4113	0.0000	1.0000
SYMG	17	-2615.8606	5265.7212	0.0000	1.0000
SYMIG	18	-2615.8607	5267.7213	0.0000	1.0000
K2PG	13	-2622.3780	5270.7559	0.0000	1.0000
K2PI	13	-2622.4727	5270.9454	0.0000	1.0000
K2PIG	14	-2622.3780	5272.7560	0.0000	1.0000
GTR	19	-2640.8764	5319.7528	0.0000	1.0000
F81I	15	-2650.6227	5331.2455	0.0000	1.0000
F81G	15	-2650.6528	5331.3056	0.0000	1.0000
F81IG	16	-2650.6528	5333.3057	0.0000	1.0000
HKY	15	-2657.8779	5345.7558	0.0000	1.0000
JC69G	12	-2668.8875	5361.7749	0.0000	1.0000
JC69I	12	-2668.9965	5361.9929	0.0000	1.0000
JC69IG	13	-2668.8875	5363.7750	0.0000	1.0000
SYM	16	-2673.2194	5378.4388	0.0000	1.0000
K2P	12	-2681.6548	5387.3097	0.0000	1.0000
F81	14	-2705.5539	5439.1077	0.0000	1.0000
JC69	11	-2725.3925	5472.7849	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2585.6540	5307.5328	0.4447	0.4447
GTRI	20	-2585.6958	5307.6165	0.4265	0.8712
HKYI	16	-2601.2833	5311.5465	0.0598	0.9310
HKYG	16	-2601.4122	5311.8043	0.0525	0.9835
GTRIG	21	-2585.6540	5314.3442	0.0148	0.9983
HKYIG	17	-2601.4123	5318.6157	0.0017	1.0000
K2PG	13	-2622.3780	5333.3021	0.0000	1.0000
K2PI	13	-2622.4727	5333.4916	0.0000	1.0000
K2PIG	14	-2622.3780	5340.1134	0.0000	1.0000
SYMI	17	-2615.7056	5347.2025	0.0000	1.0000
SYMG	17	-2615.8606	5347.5124	0.0000	1.0000
SYMIG	18	-2615.8607	5354.3237	0.0000	1.0000
F81I	15	-2650.6227	5403.4142	0.0000	1.0000
F81G	15	-2650.6528	5403.4742	0.0000	1.0000
F81IG	16	-2650.6528	5410.2856	0.0000	1.0000
GTR	19	-2640.8764	5411.1665	0.0000	1.0000
HKY	15	-2657.8779	5417.9244	0.0000	1.0000
JC69G	12	-2668.8875	5419.5099	0.0000	1.0000
JC69I	12	-2668.9965	5419.7279	0.0000	1.0000
JC69IG	13	-2668.8875	5426.3212	0.0000	1.0000
K2P	12	-2681.6548	5445.0446	0.0000	1.0000
SYM	16	-2673.2194	5455.4187	0.0000	1.0000
F81	14	-2705.5539	5506.4651	0.0000	1.0000
JC69	11	-2725.3925	5525.7086	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_28829.phylip = 1 - 908;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28829.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_28829.phylip = 1 - 908;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28829.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/chr1_28829.phylip = 1 - 908;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28829.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

