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

Input data from file "591-tomodel-phylip/chrZ_11569.phylip" (7 taxa, 818 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	-2135.1569	4302.9929	0.4363	0.4363
HKYI	16	-2135.4825	4303.6441	0.3151	0.7514
HKYIG	17	-2134.9186	4304.6022	0.1952	0.9466
GTRG	20	-2133.8014	4308.6567	0.0257	0.9723
GTRI	20	-2134.1375	4309.3290	0.0184	0.9906
GTRIG	21	-2133.7579	4310.6765	0.0094	1.0000
K2PG	13	-2151.8995	4330.2518	0.0000	1.0000
K2PI	13	-2152.2221	4330.8969	0.0000	1.0000
K2PIG	14	-2151.8273	4332.1776	0.0000	1.0000
SYMG	17	-2151.3217	4337.4084	0.0000	1.0000
SYMI	17	-2151.6132	4337.9914	0.0000	1.0000
SYMIG	18	-2151.2576	4339.3713	0.0000	1.0000
F81G	15	-2160.6502	4351.8989	0.0000	1.0000
F81I	15	-2161.1576	4352.9137	0.0000	1.0000
F81IG	16	-2160.4481	4353.5753	0.0000	1.0000
HKY	15	-2162.6807	4355.9600	0.0000	1.0000
GTR	19	-2160.6035	4360.1593	0.0000	1.0000
JC69G	12	-2176.6142	4377.6160	0.0000	1.0000
JC69I	12	-2177.0500	4378.4876	0.0000	1.0000
JC69IG	13	-2176.4044	4379.2615	0.0000	1.0000
K2P	12	-2180.8852	4386.1579	0.0000	1.0000
SYM	16	-2180.3081	4393.2954	0.0000	1.0000
F81	14	-2187.3525	4403.2280	0.0000	1.0000
JC69	11	-2203.8798	4430.0871	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2135.1569	4302.3138	0.4277	0.4277
HKYI	16	-2135.4825	4302.9649	0.3088	0.7365
HKYIG	17	-2134.9186	4303.8372	0.1997	0.9362
GTRG	20	-2133.8014	4307.6028	0.0304	0.9666
GTRI	20	-2134.1375	4308.2750	0.0217	0.9883
GTRIG	21	-2133.7579	4309.5157	0.0117	1.0000
K2PG	13	-2151.8995	4329.7990	0.0000	1.0000
K2PI	13	-2152.2221	4330.4442	0.0000	1.0000
K2PIG	14	-2151.8273	4331.6546	0.0000	1.0000
SYMG	17	-2151.3217	4336.6434	0.0000	1.0000
SYMI	17	-2151.6132	4337.2264	0.0000	1.0000
SYMIG	18	-2151.2576	4338.5152	0.0000	1.0000
F81G	15	-2160.6502	4351.3004	0.0000	1.0000
F81I	15	-2161.1576	4352.3152	0.0000	1.0000
F81IG	16	-2160.4481	4352.8961	0.0000	1.0000
HKY	15	-2162.6807	4355.3615	0.0000	1.0000
GTR	19	-2160.6035	4359.2070	0.0000	1.0000
JC69G	12	-2176.6142	4377.2284	0.0000	1.0000
JC69I	12	-2177.0500	4378.1000	0.0000	1.0000
JC69IG	13	-2176.4044	4378.8087	0.0000	1.0000
K2P	12	-2180.8852	4385.7704	0.0000	1.0000
SYM	16	-2180.3081	4392.6162	0.0000	1.0000
F81	14	-2187.3525	4402.7049	0.0000	1.0000
JC69	11	-2203.8798	4429.7596	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2135.1569	4377.6236	0.5654	0.5654
HKYI	16	-2135.4825	4378.2747	0.4083	0.9737
HKYIG	17	-2134.9186	4383.8538	0.0251	0.9987
K2PG	13	-2151.8995	4390.9882	0.0007	0.9995
K2PI	13	-2152.2221	4391.6334	0.0005	1.0000
K2PIG	14	-2151.8273	4397.5506	0.0000	1.0000
GTRG	20	-2133.8014	4401.7400	0.0000	1.0000
GTRI	20	-2134.1375	4402.4122	0.0000	1.0000
GTRIG	21	-2133.7579	4408.3598	0.0000	1.0000
SYMG	17	-2151.3217	4416.6601	0.0000	1.0000
SYMI	17	-2151.6132	4417.2430	0.0000	1.0000
F81G	15	-2160.6502	4421.9034	0.0000	1.0000
F81I	15	-2161.1576	4422.9181	0.0000	1.0000
SYMIG	18	-2151.2576	4423.2387	0.0000	1.0000
HKY	15	-2162.6807	4425.9644	0.0000	1.0000
F81IG	16	-2160.4481	4428.2059	0.0000	1.0000
JC69G	12	-2176.6142	4433.7108	0.0000	1.0000
JC69I	12	-2177.0500	4434.5824	0.0000	1.0000
JC69IG	13	-2176.4044	4439.9980	0.0000	1.0000
K2P	12	-2180.8852	4442.2527	0.0000	1.0000
GTR	19	-2160.6035	4448.6373	0.0000	1.0000
SYM	16	-2180.3081	4467.9260	0.0000	1.0000
F81	14	-2187.3525	4468.6010	0.0000	1.0000
JC69	11	-2203.8798	4481.5351	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_11569.phylip = 1 - 818;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11569.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_11569.phylip = 1 - 818;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11569.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_11569.phylip = 1 - 818;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11569.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

