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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1291.8565	2625.0832	0.3119	0.3119
GTRI	20	-1291.8592	2625.0888	0.3110	0.6228
HKYI	16	-1297.0922	2627.0661	0.1157	0.7386
HKYG	16	-1297.0923	2627.0663	0.1157	0.8543
GTRIG	21	-1291.8645	2627.2388	0.1061	0.9604
HKYIG	17	-1297.1109	2629.2153	0.0395	0.9999
GTR	19	-1301.2370	2641.7118	0.0001	1.0000
HKY	15	-1306.8386	2644.4538	0.0000	1.0000
F81I	15	-1309.9804	2650.7375	0.0000	1.0000
F81G	15	-1309.9814	2650.7396	0.0000	1.0000
F81IG	16	-1309.9814	2652.8446	0.0000	1.0000
K2PI	13	-1316.1907	2658.9686	0.0000	1.0000
K2PG	13	-1316.1972	2658.9816	0.0000	1.0000
SYMI	17	-1312.4306	2659.8547	0.0000	1.0000
SYMG	17	-1312.4307	2659.8549	0.0000	1.0000
K2PIG	14	-1316.1908	2661.0601	0.0000	1.0000
SYMIG	18	-1312.4307	2661.9735	0.0000	1.0000
F81	14	-1319.6754	2668.0292	0.0000	1.0000
SYM	16	-1322.3121	2677.5059	0.0000	1.0000
K2P	12	-1326.5020	2677.5064	0.0000	1.0000
JC69G	12	-1328.4799	2681.4622	0.0000	1.0000
JC69I	12	-1328.4801	2681.4626	0.0000	1.0000
JC69IG	13	-1328.4799	2683.5469	0.0000	1.0000
JC69	11	-1338.4064	2699.2373	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1291.8565	2623.7129	0.3280	0.3280
GTRI	20	-1291.8592	2623.7185	0.3271	0.6551
GTRIG	21	-1291.8645	2625.7290	0.1197	0.7748
HKYI	16	-1297.0922	2626.1844	0.0953	0.8702
HKYG	16	-1297.0923	2626.1846	0.0953	0.9655
HKYIG	17	-1297.1109	2628.2218	0.0344	0.9999
GTR	19	-1301.2370	2640.4740	0.0001	1.0000
HKY	15	-1306.8386	2643.6771	0.0000	1.0000
F81I	15	-1309.9804	2649.9608	0.0000	1.0000
F81G	15	-1309.9814	2649.9629	0.0000	1.0000
F81IG	16	-1309.9814	2651.9629	0.0000	1.0000
K2PI	13	-1316.1907	2658.3815	0.0000	1.0000
K2PG	13	-1316.1972	2658.3945	0.0000	1.0000
SYMI	17	-1312.4306	2658.8612	0.0000	1.0000
SYMG	17	-1312.4307	2658.8614	0.0000	1.0000
K2PIG	14	-1316.1908	2660.3816	0.0000	1.0000
SYMIG	18	-1312.4307	2660.8613	0.0000	1.0000
F81	14	-1319.6754	2667.3507	0.0000	1.0000
SYM	16	-1322.3121	2676.6242	0.0000	1.0000
K2P	12	-1326.5020	2677.0040	0.0000	1.0000
JC69G	12	-1328.4799	2680.9598	0.0000	1.0000
JC69I	12	-1328.4801	2680.9602	0.0000	1.0000
JC69IG	13	-1328.4799	2682.9598	0.0000	1.0000
JC69	11	-1338.4064	2698.8129	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1297.0922	2697.4172	0.4898	0.4898
HKYG	16	-1297.0923	2697.4174	0.4898	0.9796
HKYIG	17	-1297.1109	2703.9066	0.0191	0.9987
HKY	15	-1306.8386	2710.4579	0.0007	0.9994
GTRG	20	-1291.8565	2712.7539	0.0002	0.9996
GTRI	20	-1291.8592	2712.7594	0.0002	0.9998
K2PI	13	-1316.1907	2716.2581	0.0000	0.9999
K2PG	13	-1316.1972	2716.2711	0.0000	0.9999
F81I	15	-1309.9804	2716.7415	0.0000	1.0000
F81G	15	-1309.9814	2716.7436	0.0000	1.0000
GTRIG	21	-1291.8645	2719.2220	0.0000	1.0000
K2PIG	14	-1316.1908	2722.7103	0.0000	1.0000
F81IG	16	-1309.9814	2723.1957	0.0000	1.0000
GTR	19	-1301.2370	2725.0630	0.0000	1.0000
F81	14	-1319.6754	2729.6794	0.0000	1.0000
K2P	12	-1326.5020	2730.4285	0.0000	1.0000
JC69G	12	-1328.4799	2734.3844	0.0000	1.0000
JC69I	12	-1328.4801	2734.3847	0.0000	1.0000
SYMI	17	-1312.4306	2734.5460	0.0000	1.0000
SYMG	17	-1312.4307	2734.5463	0.0000	1.0000
JC69IG	13	-1328.4799	2740.8365	0.0000	1.0000
SYMIG	18	-1312.4307	2740.9982	0.0000	1.0000
JC69	11	-1338.4064	2747.7854	0.0000	1.0000
SYM	16	-1322.3121	2747.8570	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_2502.phylip = 1 - 634;
 Partition Dummy = 1:591-tomodel-phylip/chr9_2502.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_2502.phylip = 1 - 634;
 Partition Dummy = 1:591-tomodel-phylip/chr9_2502.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_2502.phylip = 1 - 634;
 Partition Dummy = 1:591-tomodel-phylip/chr9_2502.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

