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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2177.4878	4387.7187	0.5140	0.5140
HKYIG	17	-2177.1561	4389.1493	0.2513	0.7653
HKYI	16	-2179.0273	4390.7978	0.1102	0.8755
GTRIG	21	-2174.2779	4391.8267	0.0659	0.9414
GTRI	20	-2176.1320	4393.4179	0.0297	0.9712
GTRG	20	-2176.1633	4393.4804	0.0288	1.0000
K2PG	13	-2208.8471	4444.1894	0.0000	1.0000
K2PIG	14	-2208.4177	4445.4076	0.0000	1.0000
K2PI	13	-2210.4447	4447.3846	0.0000	1.0000
HKY	15	-2209.1467	4448.9483	0.0000	1.0000
SYMG	17	-2207.7322	4450.3016	0.0000	1.0000
GTR	19	-2205.7395	4450.5216	0.0000	1.0000
SYMIG	18	-2207.3359	4451.6088	0.0000	1.0000
SYMI	17	-2209.3159	4453.4690	0.0000	1.0000
F81G	15	-2225.1599	4480.9745	0.0000	1.0000
F81IG	16	-2224.9714	4482.6859	0.0000	1.0000
F81I	15	-2226.3975	4483.4499	0.0000	1.0000
K2P	12	-2242.8542	4510.1323	0.0000	1.0000
SYM	16	-2241.8309	4516.4049	0.0000	1.0000
JC69G	12	-2252.6169	4529.6576	0.0000	1.0000
JC69IG	13	-2252.3507	4531.1967	0.0000	1.0000
JC69I	12	-2253.9839	4532.3918	0.0000	1.0000
F81	14	-2255.5067	4539.5855	0.0000	1.0000
JC69	11	-2284.6177	4591.5935	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2177.4878	4386.9756	0.4917	0.4917
HKYIG	17	-2177.1561	4388.3121	0.2520	0.7437
HKYI	16	-2179.0273	4390.0546	0.1055	0.8491
GTRIG	21	-2174.2779	4390.5557	0.0821	0.9312
GTRI	20	-2176.1320	4392.2641	0.0349	0.9661
GTRG	20	-2176.1633	4392.3265	0.0339	1.0000
K2PG	13	-2208.8471	4443.6941	0.0000	1.0000
K2PIG	14	-2208.4177	4444.8354	0.0000	1.0000
K2PI	13	-2210.4447	4446.8894	0.0000	1.0000
HKY	15	-2209.1467	4448.2935	0.0000	1.0000
SYMG	17	-2207.7322	4449.4644	0.0000	1.0000
GTR	19	-2205.7395	4449.4791	0.0000	1.0000
SYMIG	18	-2207.3359	4450.6718	0.0000	1.0000
SYMI	17	-2209.3159	4452.6318	0.0000	1.0000
F81G	15	-2225.1599	4480.3197	0.0000	1.0000
F81IG	16	-2224.9714	4481.9427	0.0000	1.0000
F81I	15	-2226.3975	4482.7951	0.0000	1.0000
K2P	12	-2242.8542	4509.7084	0.0000	1.0000
SYM	16	-2241.8309	4515.6617	0.0000	1.0000
JC69G	12	-2252.6169	4529.2337	0.0000	1.0000
JC69IG	13	-2252.3507	4530.7014	0.0000	1.0000
JC69I	12	-2253.9839	4531.9679	0.0000	1.0000
F81	14	-2255.5067	4539.0133	0.0000	1.0000
JC69	11	-2284.6177	4591.2353	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2177.4878	4460.8754	0.7903	0.7903
HKYI	16	-2179.0273	4463.9544	0.1695	0.9598
HKYIG	17	-2177.1561	4466.8307	0.0402	1.0000
GTRI	20	-2176.1320	4484.6389	0.0000	1.0000
GTRG	20	-2176.1633	4484.7013	0.0000	1.0000
GTRIG	21	-2174.2779	4487.5492	0.0000	1.0000
K2PG	13	-2208.8471	4503.7377	0.0000	1.0000
K2PI	13	-2210.4447	4506.9330	0.0000	1.0000
K2PIG	14	-2208.4177	4509.4978	0.0000	1.0000
HKY	15	-2209.1467	4517.5746	0.0000	1.0000
SYMG	17	-2207.7322	4527.9830	0.0000	1.0000
SYMI	17	-2209.3159	4531.1503	0.0000	1.0000
SYMIG	18	-2207.3359	4533.8091	0.0000	1.0000
GTR	19	-2205.7395	4537.2351	0.0000	1.0000
F81G	15	-2225.1599	4549.6008	0.0000	1.0000
F81I	15	-2226.3975	4552.0762	0.0000	1.0000
F81IG	16	-2224.9714	4555.8425	0.0000	1.0000
K2P	12	-2242.8542	4565.1333	0.0000	1.0000
JC69G	12	-2252.6169	4584.6586	0.0000	1.0000
JC69I	12	-2253.9839	4587.3928	0.0000	1.0000
SYM	16	-2241.8309	4589.5616	0.0000	1.0000
JC69IG	13	-2252.3507	4590.7450	0.0000	1.0000
F81	14	-2255.5067	4603.6757	0.0000	1.0000
JC69	11	-2284.6177	4642.0414	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYG)]

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


[Mrbayes block for the best AICc model (HKYG)]

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


[Mrbayes block for the best BIC model (HKYG)]

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

