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

Input data from file "591-tomodel-phylip/chr1_28903.phylip" (7 taxa, 515 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	-1200.7298	2434.5519	0.4447	0.4447
HKYG	16	-1200.8867	2434.8659	0.3801	0.8248
HKYIG	17	-1200.8867	2437.0049	0.1304	0.9552
GTRI	20	-1199.4997	2440.6998	0.0206	0.9757
GTRG	20	-1199.6260	2440.9523	0.0181	0.9939
GTRIG	21	-1199.6259	2443.1261	0.0061	1.0000
HKY	15	-1212.0601	2455.0822	0.0000	1.0000
GTR	19	-1211.0983	2461.7319	0.0000	1.0000
K2PI	13	-1220.9992	2468.7249	0.0000	1.0000
K2PG	13	-1221.2085	2469.1436	0.0000	1.0000
K2PIG	14	-1221.2085	2471.2570	0.0000	1.0000
F81I	15	-1220.4807	2471.9233	0.0000	1.0000
F81G	15	-1220.6534	2472.2686	0.0000	1.0000
F81IG	16	-1220.6533	2474.3990	0.0000	1.0000
SYMI	17	-1220.2817	2475.7947	0.0000	1.0000
SYMG	17	-1220.5050	2476.2414	0.0000	1.0000
SYMIG	18	-1220.5051	2478.3892	0.0000	1.0000
K2P	12	-1232.2009	2489.0233	0.0000	1.0000
F81	14	-1231.0262	2490.8923	0.0000	1.0000
SYM	16	-1231.5137	2496.1198	0.0000	1.0000
JC69I	12	-1238.6459	2501.9133	0.0000	1.0000
JC69G	12	-1238.8533	2502.3282	0.0000	1.0000
JC69IG	13	-1238.8533	2504.4332	0.0000	1.0000
JC69	11	-1249.2481	2521.0210	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1200.7298	2433.4596	0.4334	0.4334
HKYG	16	-1200.8867	2433.7735	0.3704	0.8038
HKYIG	17	-1200.8867	2435.7735	0.1363	0.9401
GTRI	20	-1199.4997	2438.9994	0.0272	0.9672
GTRG	20	-1199.6260	2439.2519	0.0239	0.9912
GTRIG	21	-1199.6259	2441.2519	0.0088	1.0000
HKY	15	-1212.0601	2454.1203	0.0000	1.0000
GTR	19	-1211.0983	2460.1966	0.0000	1.0000
K2PI	13	-1220.9992	2467.9984	0.0000	1.0000
K2PG	13	-1221.2085	2468.4170	0.0000	1.0000
K2PIG	14	-1221.2085	2470.4170	0.0000	1.0000
F81I	15	-1220.4807	2470.9614	0.0000	1.0000
F81G	15	-1220.6534	2471.3067	0.0000	1.0000
F81IG	16	-1220.6533	2473.3067	0.0000	1.0000
SYMI	17	-1220.2817	2474.5633	0.0000	1.0000
SYMG	17	-1220.5050	2475.0101	0.0000	1.0000
SYMIG	18	-1220.5051	2477.0101	0.0000	1.0000
K2P	12	-1232.2009	2488.4017	0.0000	1.0000
F81	14	-1231.0262	2490.0523	0.0000	1.0000
SYM	16	-1231.5137	2495.0275	0.0000	1.0000
JC69I	12	-1238.6459	2501.2918	0.0000	1.0000
JC69G	12	-1238.8533	2501.7067	0.0000	1.0000
JC69IG	13	-1238.8533	2503.7066	0.0000	1.0000
JC69	11	-1249.2481	2520.4961	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1200.7298	2501.3663	0.5283	0.5283
HKYG	16	-1200.8867	2501.6802	0.4516	0.9799
HKYIG	17	-1200.8867	2507.9243	0.0199	0.9998
HKY	15	-1212.0601	2517.7828	0.0001	1.0000
K2PI	13	-1220.9992	2523.1726	0.0000	1.0000
K2PG	13	-1221.2085	2523.5912	0.0000	1.0000
GTRI	20	-1199.4997	2523.8827	0.0000	1.0000
GTRG	20	-1199.6260	2524.1352	0.0000	1.0000
K2PIG	14	-1221.2085	2529.8353	0.0000	1.0000
GTRIG	21	-1199.6259	2530.3794	0.0000	1.0000
F81I	15	-1220.4807	2534.6239	0.0000	1.0000
F81G	15	-1220.6534	2534.9692	0.0000	1.0000
K2P	12	-1232.2009	2539.3317	0.0000	1.0000
GTR	19	-1211.0983	2540.8358	0.0000	1.0000
F81IG	16	-1220.6533	2541.2134	0.0000	1.0000
SYMI	17	-1220.2817	2546.7142	0.0000	1.0000
SYMG	17	-1220.5050	2547.1609	0.0000	1.0000
F81	14	-1231.0262	2549.4707	0.0000	1.0000
JC69I	12	-1238.6459	2552.2218	0.0000	1.0000
JC69G	12	-1238.8533	2552.6367	0.0000	1.0000
SYMIG	18	-1220.5051	2553.4051	0.0000	1.0000
JC69IG	13	-1238.8533	2558.8808	0.0000	1.0000
SYM	16	-1231.5137	2562.9341	0.0000	1.0000
JC69	11	-1249.2481	2567.1820	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/chr1_28903.phylip = 1 - 515;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28903.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/chr1_28903.phylip = 1 - 515;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28903.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/chr1_28903.phylip = 1 - 515;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28903.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

