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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-2008.6025	4050.2200	0.6784	0.6784
HKYIG	17	-2008.6316	4052.4071	0.2273	0.9057
GTRI	20	-2006.9053	4055.3895	0.0512	0.9569
HKYG	16	-2011.9834	4056.9818	0.0231	0.9800
GTRIG	21	-2006.8634	4057.4670	0.0181	0.9981
GTRG	20	-2010.1890	4061.9570	0.0019	1.0000
F81I	15	-2030.6708	4092.2356	0.0000	1.0000
F81IG	16	-2030.6927	4094.4003	0.0000	1.0000
F81G	15	-2033.3077	4097.5092	0.0000	1.0000
GTR	19	-2039.9245	4119.2748	0.0000	1.0000
HKY	15	-2044.3589	4119.6117	0.0000	1.0000
F81	14	-2065.2306	4159.2419	0.0000	1.0000
SYMG	17	-2083.5216	4202.1871	0.0000	1.0000
SYMI	17	-2083.6154	4202.3748	0.0000	1.0000
SYMIG	18	-2083.2535	4203.7879	0.0000	1.0000
K2PI	13	-2096.0934	4218.8621	0.0000	1.0000
K2PIG	14	-2096.0993	4220.9792	0.0000	1.0000
K2PG	13	-2097.6654	4222.0061	0.0000	1.0000
JC69I	12	-2109.1523	4242.8824	0.0000	1.0000
JC69IG	13	-2109.1541	4244.9835	0.0000	1.0000
JC69G	12	-2110.5487	4245.6751	0.0000	1.0000
SYM	16	-2119.4871	4271.9891	0.0000	1.0000
K2P	12	-2131.3774	4287.3326	0.0000	1.0000
JC69	11	-2143.2975	4309.0830	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2008.6025	4049.2051	0.6521	0.6521
HKYIG	17	-2008.6316	4051.2632	0.2330	0.8852
GTRI	20	-2006.9053	4053.8106	0.0652	0.9504
GTRIG	21	-2006.8634	4055.7269	0.0250	0.9754
HKYG	16	-2011.9834	4055.9669	0.0222	0.9976
GTRG	20	-2010.1890	4060.3780	0.0024	1.0000
F81I	15	-2030.6708	4091.3417	0.0000	1.0000
F81IG	16	-2030.6927	4093.3854	0.0000	1.0000
F81G	15	-2033.3077	4096.6154	0.0000	1.0000
GTR	19	-2039.9245	4117.8489	0.0000	1.0000
HKY	15	-2044.3589	4118.7178	0.0000	1.0000
F81	14	-2065.2306	4158.4613	0.0000	1.0000
SYMG	17	-2083.5216	4201.0432	0.0000	1.0000
SYMI	17	-2083.6154	4201.2308	0.0000	1.0000
SYMIG	18	-2083.2535	4202.5070	0.0000	1.0000
K2PI	13	-2096.0934	4218.1868	0.0000	1.0000
K2PIG	14	-2096.0993	4220.1986	0.0000	1.0000
K2PG	13	-2097.6654	4221.3308	0.0000	1.0000
JC69I	12	-2109.1523	4242.3046	0.0000	1.0000
JC69IG	13	-2109.1541	4244.3082	0.0000	1.0000
JC69G	12	-2110.5487	4245.0973	0.0000	1.0000
SYM	16	-2119.4871	4270.9742	0.0000	1.0000
K2P	12	-2131.3774	4286.7549	0.0000	1.0000
JC69	11	-2143.2975	4308.5950	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2008.6025	4118.2508	0.9299	0.9299
HKYIG	17	-2008.6316	4124.6243	0.0384	0.9683
HKYG	16	-2011.9834	4125.0126	0.0316	1.0000
GTRI	20	-2006.9053	4140.1177	0.0000	1.0000
GTRIG	21	-2006.8634	4146.3494	0.0000	1.0000
GTRG	20	-2010.1890	4146.6852	0.0000	1.0000
F81I	15	-2030.6708	4156.0721	0.0000	1.0000
F81G	15	-2033.3077	4161.3457	0.0000	1.0000
F81IG	16	-2030.6927	4162.4311	0.0000	1.0000
HKY	15	-2044.3589	4183.4482	0.0000	1.0000
GTR	19	-2039.9245	4199.8407	0.0000	1.0000
F81	14	-2065.2306	4218.8763	0.0000	1.0000
K2PI	13	-2096.0934	4274.2865	0.0000	1.0000
SYMG	17	-2083.5216	4274.4042	0.0000	1.0000
SYMI	17	-2083.6154	4274.5919	0.0000	1.0000
K2PG	13	-2097.6654	4277.4304	0.0000	1.0000
SYMIG	18	-2083.2535	4280.1835	0.0000	1.0000
K2PIG	14	-2096.0993	4280.6136	0.0000	1.0000
JC69I	12	-2109.1523	4294.0889	0.0000	1.0000
JC69G	12	-2110.5487	4296.8816	0.0000	1.0000
JC69IG	13	-2109.1541	4300.4079	0.0000	1.0000
K2P	12	-2131.3774	4338.5392	0.0000	1.0000
SYM	16	-2119.4871	4340.0199	0.0000	1.0000
JC69	11	-2143.2975	4356.0639	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_27741.phylip = 1 - 553;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27741.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_27741.phylip = 1 - 553;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27741.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_27741.phylip = 1 - 553;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27741.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

