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

Input data from file "591-tomodel-phylip/chr1_31933.phylip" (7 taxa, 429 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	-1011.0036	2055.3276	0.4964	0.4964
HKYI	16	-1011.5900	2056.5004	0.2762	0.7726
HKYIG	17	-1010.8626	2057.2143	0.1933	0.9658
GTRG	20	-1010.0298	2062.1185	0.0166	0.9825
GTRI	20	-1010.6247	2063.3083	0.0092	0.9917
GTRIG	21	-1009.8876	2064.0454	0.0064	0.9980
HKY	15	-1017.6502	2066.4625	0.0019	0.9999
GTR	19	-1016.4476	2072.7534	0.0001	1.0000
F81G	15	-1023.2273	2077.6168	0.0000	1.0000
F81I	15	-1023.5333	2078.2288	0.0000	1.0000
F81IG	16	-1023.1740	2079.6683	0.0000	1.0000
F81	14	-1029.3734	2087.7613	0.0000	1.0000
SYMG	17	-1051.2917	2138.0725	0.0000	1.0000
SYMI	17	-1051.6364	2138.7619	0.0000	1.0000
SYMIG	18	-1051.2006	2140.0695	0.0000	1.0000
K2PG	13	-1058.4493	2143.7758	0.0000	1.0000
K2PI	13	-1058.7786	2144.4343	0.0000	1.0000
K2PIG	14	-1058.3794	2145.7734	0.0000	1.0000
SYM	16	-1057.8906	2149.1017	0.0000	1.0000
K2P	12	-1064.7519	2154.2538	0.0000	1.0000
JC69G	12	-1067.4507	2159.6514	0.0000	1.0000
JC69I	12	-1067.6745	2160.0991	0.0000	1.0000
JC69IG	13	-1067.4169	2161.7109	0.0000	1.0000
JC69	11	-1073.4209	2169.4749	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1011.0036	2054.0073	0.4809	0.4809
HKYI	16	-1011.5900	2055.1801	0.2675	0.7484
HKYIG	17	-1010.8626	2055.7253	0.2037	0.9521
GTRG	20	-1010.0298	2060.0597	0.0233	0.9754
GTRI	20	-1010.6247	2061.2494	0.0129	0.9883
GTRIG	21	-1009.8876	2061.7752	0.0099	0.9982
HKY	15	-1017.6502	2065.3003	0.0017	0.9999
GTR	19	-1016.4476	2070.8952	0.0001	1.0000
F81G	15	-1023.2273	2076.4545	0.0000	1.0000
F81I	15	-1023.5333	2077.0665	0.0000	1.0000
F81IG	16	-1023.1740	2078.3479	0.0000	1.0000
F81	14	-1029.3734	2086.7468	0.0000	1.0000
SYMG	17	-1051.2917	2136.5834	0.0000	1.0000
SYMI	17	-1051.6364	2137.2728	0.0000	1.0000
SYMIG	18	-1051.2006	2138.4012	0.0000	1.0000
K2PG	13	-1058.4493	2142.8987	0.0000	1.0000
K2PI	13	-1058.7786	2143.5572	0.0000	1.0000
K2PIG	14	-1058.3794	2144.7589	0.0000	1.0000
SYM	16	-1057.8906	2147.7813	0.0000	1.0000
K2P	12	-1064.7519	2153.5038	0.0000	1.0000
JC69G	12	-1067.4507	2158.9014	0.0000	1.0000
JC69I	12	-1067.6745	2159.3491	0.0000	1.0000
JC69IG	13	-1067.4169	2160.8338	0.0000	1.0000
JC69	11	-1073.4209	2168.8418	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1011.0036	2118.9906	0.6101	0.6101
HKYI	16	-1011.5900	2120.1634	0.3394	0.9495
HKYIG	17	-1010.8626	2124.7700	0.0339	0.9835
HKY	15	-1017.6502	2126.2222	0.0164	0.9999
F81G	15	-1023.2273	2137.3764	0.0001	0.9999
F81I	15	-1023.5333	2137.9884	0.0000	1.0000
GTRG	20	-1010.0298	2141.2888	0.0000	1.0000
GTRI	20	-1010.6247	2142.4786	0.0000	1.0000
F81IG	16	-1023.1740	2143.3313	0.0000	1.0000
F81	14	-1029.3734	2143.6072	0.0000	1.0000
GTRIG	21	-1009.8876	2147.0658	0.0000	1.0000
GTR	19	-1016.4476	2148.0629	0.0000	1.0000
K2PG	13	-1058.4493	2195.6976	0.0000	1.0000
K2PI	13	-1058.7786	2196.3562	0.0000	1.0000
K2PIG	14	-1058.3794	2201.6193	0.0000	1.0000
K2P	12	-1064.7519	2202.2413	0.0000	1.0000
SYMG	17	-1051.2917	2205.6282	0.0000	1.0000
SYMI	17	-1051.6364	2206.3176	0.0000	1.0000
JC69G	12	-1067.4507	2207.6389	0.0000	1.0000
JC69I	12	-1067.6745	2208.0866	0.0000	1.0000
SYMIG	18	-1051.2006	2211.5074	0.0000	1.0000
SYM	16	-1057.8906	2212.7646	0.0000	1.0000
JC69	11	-1073.4209	2213.5178	0.0000	1.0000
JC69IG	13	-1067.4169	2213.6328	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/chr1_31933.phylip = 1 - 429;
 Partition Dummy = 1:591-tomodel-phylip/chr1_31933.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/chr1_31933.phylip = 1 - 429;
 Partition Dummy = 1:591-tomodel-phylip/chr1_31933.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/chr1_31933.phylip = 1 - 429;
 Partition Dummy = 1:591-tomodel-phylip/chr1_31933.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

