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

Input data from file "591-tomodel-phylip/chr2_24864.phylip" (7 taxa, 645 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	-1991.4052	4015.6766	0.4564	0.4564
HKYI	16	-1991.7401	4016.3465	0.3265	0.7828
HKYIG	17	-1991.3832	4017.7425	0.1624	0.9453
GTRG	20	-1989.9700	4021.2861	0.0276	0.9729
GTRI	20	-1990.4383	4022.2228	0.0173	0.9902
GTRIG	21	-1989.9377	4023.3585	0.0098	1.0000
HKY	15	-2010.0025	4050.7681	0.0000	1.0000
GTR	19	-2008.2939	4055.8037	0.0000	1.0000
K2PG	13	-2031.9668	4090.5104	0.0000	1.0000
K2PI	13	-2032.2622	4091.1013	0.0000	1.0000
K2PIG	14	-2031.9661	4092.5989	0.0000	1.0000
SYMG	17	-2030.1930	4095.3622	0.0000	1.0000
SYMI	17	-2030.6538	4096.2837	0.0000	1.0000
SYMIG	18	-2030.1822	4097.4570	0.0000	1.0000
K2P	12	-2053.8645	4132.2226	0.0000	1.0000
F81G	15	-2051.0742	4132.9115	0.0000	1.0000
F81I	15	-2051.8459	4134.4550	0.0000	1.0000
F81IG	16	-2051.0663	4134.9988	0.0000	1.0000
SYM	16	-2053.1000	4139.0663	0.0000	1.0000
F81	14	-2069.7918	4168.2503	0.0000	1.0000
JC69G	12	-2085.2368	4194.9674	0.0000	1.0000
JC69I	12	-2085.9788	4196.4513	0.0000	1.0000
JC69IG	13	-2085.2282	4197.0333	0.0000	1.0000
JC69	11	-2105.4406	4233.2983	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1991.4052	4014.8104	0.4453	0.4453
HKYI	16	-1991.7401	4015.4802	0.3186	0.7638
HKYIG	17	-1991.3832	4016.7664	0.1674	0.9313
GTRG	20	-1989.9700	4019.9400	0.0343	0.9655
GTRI	20	-1990.4383	4020.8766	0.0214	0.9870
GTRIG	21	-1989.9377	4021.8754	0.0130	1.0000
HKY	15	-2010.0025	4050.0050	0.0000	1.0000
GTR	19	-2008.2939	4054.5877	0.0000	1.0000
K2PG	13	-2031.9668	4089.9335	0.0000	1.0000
K2PI	13	-2032.2622	4090.5244	0.0000	1.0000
K2PIG	14	-2031.9661	4091.9323	0.0000	1.0000
SYMG	17	-2030.1930	4094.3861	0.0000	1.0000
SYMI	17	-2030.6538	4095.3076	0.0000	1.0000
SYMIG	18	-2030.1822	4096.3643	0.0000	1.0000
K2P	12	-2053.8645	4131.7290	0.0000	1.0000
F81G	15	-2051.0742	4132.1483	0.0000	1.0000
F81I	15	-2051.8459	4133.6919	0.0000	1.0000
F81IG	16	-2051.0663	4134.1325	0.0000	1.0000
SYM	16	-2053.1000	4138.2000	0.0000	1.0000
F81	14	-2069.7918	4167.5836	0.0000	1.0000
JC69G	12	-2085.2368	4194.4737	0.0000	1.0000
JC69I	12	-2085.9788	4195.9576	0.0000	1.0000
JC69IG	13	-2085.2282	4196.4564	0.0000	1.0000
JC69	11	-2105.4406	4232.8813	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1991.4052	4086.3184	0.5696	0.5696
HKYI	16	-1991.7401	4086.9882	0.4075	0.9771
HKYIG	17	-1991.3832	4092.7437	0.0229	1.0000
GTRG	20	-1989.9700	4109.3250	0.0000	1.0000
GTRI	20	-1990.4383	4110.2616	0.0000	1.0000
GTRIG	21	-1989.9377	4115.7296	0.0000	1.0000
HKY	15	-2010.0025	4117.0437	0.0000	1.0000
GTR	19	-2008.2939	4139.5035	0.0000	1.0000
K2PG	13	-2031.9668	4148.0338	0.0000	1.0000
K2PI	13	-2032.2622	4148.6247	0.0000	1.0000
K2PIG	14	-2031.9661	4154.5018	0.0000	1.0000
SYMG	17	-2030.1930	4170.3633	0.0000	1.0000
SYMI	17	-2030.6538	4171.2849	0.0000	1.0000
SYMIG	18	-2030.1822	4176.8108	0.0000	1.0000
K2P	12	-2053.8645	4185.3600	0.0000	1.0000
F81G	15	-2051.0742	4199.1871	0.0000	1.0000
F81I	15	-2051.8459	4200.7306	0.0000	1.0000
F81IG	16	-2051.0663	4205.6405	0.0000	1.0000
SYM	16	-2053.1000	4209.7080	0.0000	1.0000
F81	14	-2069.7918	4230.1531	0.0000	1.0000
JC69G	12	-2085.2368	4248.1047	0.0000	1.0000
JC69I	12	-2085.9788	4249.5886	0.0000	1.0000
JC69IG	13	-2085.2282	4254.5567	0.0000	1.0000
JC69	11	-2105.4406	4282.0430	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/chr2_24864.phylip = 1 - 645;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24864.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/chr2_24864.phylip = 1 - 645;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24864.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/chr2_24864.phylip = 1 - 645;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24864.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

