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

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

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-2624.2804	5289.6632	0.7317	0.7317
GTRIG	21	-2624.2498	5291.7137	0.2625	0.9942
HKYG	16	-2633.9894	5300.6890	0.0030	0.9972
HKYI	16	-2634.8953	5302.5008	0.0012	0.9984
HKYIG	17	-2634.0042	5302.8084	0.0010	0.9994
GTRI	20	-2631.3666	5303.8357	0.0006	1.0000
SYMG	17	-2655.9589	5346.7178	0.0000	1.0000
SYMIG	18	-2655.8793	5348.6539	0.0000	1.0000
K2PG	13	-2661.7581	5349.9895	0.0000	1.0000
K2PIG	14	-2661.6969	5351.9407	0.0000	1.0000
SYMI	17	-2663.5676	5361.9353	0.0000	1.0000
K2PI	13	-2668.3670	5363.2074	0.0000	1.0000
F81G	15	-2674.7231	5380.0719	0.0000	1.0000
F81I	15	-2675.6065	5381.8388	0.0000	1.0000
F81IG	16	-2674.7312	5382.1727	0.0000	1.0000
GTR	19	-2674.0888	5387.1736	0.0000	1.0000
HKY	15	-2679.1229	5388.8716	0.0000	1.0000
JC69G	12	-2707.1091	5438.6235	0.0000	1.0000
JC69IG	13	-2707.1221	5440.7175	0.0000	1.0000
JC69I	12	-2708.2712	5440.9476	0.0000	1.0000
SYM	16	-2705.8792	5444.4687	0.0000	1.0000
K2P	12	-2712.7663	5449.9378	0.0000	1.0000
F81	14	-2719.0753	5466.6974	0.0000	1.0000
JC69	11	-2750.6007	5523.5439	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2624.2804	5288.5609	0.7215	0.7215
GTRIG	21	-2624.2498	5290.4995	0.2737	0.9952
HKYG	16	-2633.9894	5299.9788	0.0024	0.9976
HKYI	16	-2634.8953	5301.7906	0.0010	0.9985
HKYIG	17	-2634.0042	5302.0084	0.0009	0.9994
GTRI	20	-2631.3666	5302.7333	0.0006	1.0000
SYMG	17	-2655.9589	5345.9178	0.0000	1.0000
SYMIG	18	-2655.8793	5347.7586	0.0000	1.0000
K2PG	13	-2661.7581	5349.5162	0.0000	1.0000
K2PIG	14	-2661.6969	5351.3938	0.0000	1.0000
SYMI	17	-2663.5676	5361.1353	0.0000	1.0000
K2PI	13	-2668.3670	5362.7340	0.0000	1.0000
F81G	15	-2674.7231	5379.4461	0.0000	1.0000
F81I	15	-2675.6065	5381.2130	0.0000	1.0000
F81IG	16	-2674.7312	5381.4625	0.0000	1.0000
GTR	19	-2674.0888	5386.1775	0.0000	1.0000
HKY	15	-2679.1229	5388.2458	0.0000	1.0000
JC69G	12	-2707.1091	5438.2183	0.0000	1.0000
JC69IG	13	-2707.1221	5440.2441	0.0000	1.0000
JC69I	12	-2708.2712	5440.5424	0.0000	1.0000
SYM	16	-2705.8792	5443.7585	0.0000	1.0000
K2P	12	-2712.7663	5449.5326	0.0000	1.0000
F81	14	-2719.0753	5466.1505	0.0000	1.0000
JC69	11	-2750.6007	5523.2014	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2633.9894	5374.5889	0.6815	0.6815
HKYI	16	-2634.8953	5376.4007	0.2755	0.9570
HKYIG	17	-2634.0042	5381.2817	0.0240	0.9810
GTRG	20	-2624.2804	5381.8235	0.0183	0.9993
GTRIG	21	-2624.2498	5388.4253	0.0007	1.0000
GTRI	20	-2631.3666	5395.9960	0.0000	1.0000
K2PG	13	-2661.7581	5410.1369	0.0000	1.0000
K2PIG	14	-2661.6969	5416.6776	0.0000	1.0000
K2PI	13	-2668.3670	5423.3547	0.0000	1.0000
SYMG	17	-2655.9589	5425.1911	0.0000	1.0000
SYMIG	18	-2655.8793	5431.6950	0.0000	1.0000
SYMI	17	-2663.5676	5440.4086	0.0000	1.0000
F81G	15	-2674.7231	5449.3931	0.0000	1.0000
F81I	15	-2675.6065	5451.1600	0.0000	1.0000
F81IG	16	-2674.7312	5456.0726	0.0000	1.0000
HKY	15	-2679.1229	5458.1928	0.0000	1.0000
GTR	19	-2674.0888	5474.7771	0.0000	1.0000
JC69G	12	-2707.1091	5494.1759	0.0000	1.0000
JC69I	12	-2708.2712	5496.5000	0.0000	1.0000
JC69IG	13	-2707.1221	5500.8649	0.0000	1.0000
K2P	12	-2712.7663	5505.4902	0.0000	1.0000
SYM	16	-2705.8792	5518.3686	0.0000	1.0000
F81	14	-2719.0753	5531.4344	0.0000	1.0000
JC69	11	-2750.6007	5574.4959	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/chrZ_7997.phylip = 1 - 783;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7997.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/chrZ_7997.phylip = 1 - 783;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7997.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/chrZ_7997.phylip = 1 - 783;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7997.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

