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

Input data from file "591-tomodel-phylip/chr13_696.phylip" (7 taxa, 941 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	-2744.9147	5529.8293	0.5362	0.5362
GTRIG	21	-2744.9175	5531.8349	0.1967	0.7329
GTRI	20	-2746.3365	5532.6730	0.1294	0.8623
HKYG	16	-2750.7400	5533.4801	0.0864	0.9487
HKYIG	17	-2750.7481	5535.4962	0.0315	0.9803
HKYI	16	-2752.2171	5536.4342	0.0197	1.0000
GTR	19	-2758.4552	5554.9103	0.0000	1.0000
HKY	15	-2763.9025	5557.8051	0.0000	1.0000
K2PG	13	-2775.0930	5576.1860	0.0000	1.0000
SYMG	17	-2771.2272	5576.4544	0.0000	1.0000
K2PIG	14	-2775.0928	5578.1856	0.0000	1.0000
SYMIG	18	-2771.2294	5578.4588	0.0000	1.0000
K2PI	13	-2776.5707	5579.1413	0.0000	1.0000
SYMI	17	-2772.6512	5579.3023	0.0000	1.0000
K2P	12	-2788.4587	5600.9174	0.0000	1.0000
SYM	16	-2784.6218	5601.2437	0.0000	1.0000
F81G	15	-2806.7025	5643.4051	0.0000	1.0000
F81IG	16	-2806.7042	5645.4085	0.0000	1.0000
F81I	15	-2807.9894	5645.9787	0.0000	1.0000
F81	14	-2818.2321	5664.4641	0.0000	1.0000
JC69G	12	-2827.3808	5678.7617	0.0000	1.0000
JC69IG	13	-2827.3868	5680.7737	0.0000	1.0000
JC69I	12	-2828.6624	5681.3249	0.0000	1.0000
JC69	11	-2839.0639	5700.1278	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2750.7400	5611.0312	0.7929	0.7929
HKYI	16	-2752.2171	5613.9853	0.1810	0.9739
HKYIG	17	-2750.7481	5617.8943	0.0256	0.9996
GTRG	20	-2744.9147	5626.7682	0.0003	0.9999
GTRI	20	-2746.3365	5629.6119	0.0001	0.9999
HKY	15	-2763.9025	5630.5092	0.0000	1.0000
GTRIG	21	-2744.9175	5633.6207	0.0000	1.0000
K2PG	13	-2775.0930	5639.1962	0.0000	1.0000
K2PI	13	-2776.5707	5642.1516	0.0000	1.0000
K2PIG	14	-2775.0928	5646.0428	0.0000	1.0000
GTR	19	-2758.4552	5647.0023	0.0000	1.0000
SYMG	17	-2771.2272	5658.8525	0.0000	1.0000
K2P	12	-2788.4587	5659.0807	0.0000	1.0000
SYMI	17	-2772.6512	5661.7003	0.0000	1.0000
SYMIG	18	-2771.2294	5665.7038	0.0000	1.0000
SYM	16	-2784.6218	5678.7948	0.0000	1.0000
F81G	15	-2806.7025	5716.1092	0.0000	1.0000
F81I	15	-2807.9894	5718.6829	0.0000	1.0000
F81IG	16	-2806.7042	5722.9596	0.0000	1.0000
F81	14	-2818.2321	5732.3213	0.0000	1.0000
JC69G	12	-2827.3808	5736.9250	0.0000	1.0000
JC69I	12	-2828.6624	5739.4882	0.0000	1.0000
JC69IG	13	-2827.3868	5743.7840	0.0000	1.0000
JC69	11	-2839.0639	5753.4442	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/chr13_696.phylip = 1 - 941;
 Partition Dummy = 1:591-tomodel-phylip/chr13_696.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/chr13_696.phylip = 1 - 941;
 Partition Dummy = 1:591-tomodel-phylip/chr13_696.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/chr13_696.phylip = 1 - 941;
 Partition Dummy = 1:591-tomodel-phylip/chr13_696.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

