
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:34 02/17/2011

Input data from file "183-phylip/chr1_15480.phylip" (29 taxa, 294 characters)

Minimum AIC  model: GTRG
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	60	-1231.6010	2614.6184	0.6599	0.6599
HKYIG	61	-1231.0405	2616.6844	0.2349	0.8949
GTRG	64	-1227.5274	2619.3868	0.0608	0.9557
GTRIG	65	-1226.9786	2621.5888	0.0202	0.9759
HKYI	60	-1235.1490	2621.7143	0.0190	0.9949
K2PG	57	-1241.7306	2625.4781	0.0029	0.9978
K2PIG	58	-1241.1663	2627.4560	0.0011	0.9989
GTRI	64	-1231.7899	2627.9116	0.0009	0.9997
K2PI	57	-1244.5988	2631.2145	0.0002	0.9999
SYMG	61	-1238.9928	2632.5891	0.0001	1.0000
SYMIG	62	-1238.5210	2634.8601	0.0000	1.0000
SYMI	61	-1241.9566	2638.5166	0.0000	1.0000
F81G	59	-1262.4989	2673.2542	0.0000	1.0000
F81IG	60	-1261.9625	2675.3413	0.0000	1.0000
JC69G	56	-1269.3569	2677.6505	0.0000	1.0000
JC69IG	57	-1268.8676	2679.7521	0.0000	1.0000
F81I	59	-1265.7732	2679.8029	0.0000	1.0000
JC69I	56	-1272.4082	2683.7531	0.0000	1.0000
HKY	59	-1267.7860	2683.8284	0.0000	1.0000
GTR	63	-1263.8426	2688.7460	0.0000	1.0000
K2P	56	-1277.3873	2693.7114	0.0000	1.0000
SYM	60	-1274.4525	2700.3213	0.0000	1.0000
F81	58	-1296.8157	2738.7547	0.0000	1.0000
JC69	55	-1303.8759	2743.6341	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1227.5274	2583.0549	0.3118	0.3118
HKYG	60	-1231.6010	2583.2021	0.2897	0.6015
GTRIG	65	-1226.9786	2583.9572	0.1986	0.8001
HKYIG	61	-1231.0405	2584.0809	0.1867	0.9868
HKYI	60	-1235.1490	2590.2980	0.0083	0.9951
GTRI	64	-1231.7899	2591.5797	0.0044	0.9995
K2PG	57	-1241.7306	2597.4611	0.0002	0.9997
K2PIG	58	-1241.1663	2598.3326	0.0002	0.9999
SYMG	61	-1238.9928	2599.9856	0.0001	0.9999
SYMIG	62	-1238.5210	2601.0419	0.0000	1.0000
K2PI	57	-1244.5988	2603.1975	0.0000	1.0000
SYMI	61	-1241.9566	2605.9131	0.0000	1.0000
F81G	59	-1262.4989	2642.9978	0.0000	1.0000
F81IG	60	-1261.9625	2643.9250	0.0000	1.0000
F81I	59	-1265.7732	2649.5465	0.0000	1.0000
JC69G	56	-1269.3569	2650.7138	0.0000	1.0000
JC69IG	57	-1268.8676	2651.7352	0.0000	1.0000
HKY	59	-1267.7860	2653.5720	0.0000	1.0000
GTR	63	-1263.8426	2653.6851	0.0000	1.0000
JC69I	56	-1272.4082	2656.8164	0.0000	1.0000
K2P	56	-1277.3873	2666.7747	0.0000	1.0000
SYM	60	-1274.4525	2668.9050	0.0000	1.0000
F81	58	-1296.8157	2709.6313	0.0000	1.0000
JC69	55	-1303.8759	2717.7518	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1231.6010	2804.2169	0.7327	0.7327
K2PG	57	-1241.7306	2807.4252	0.1473	0.8800
HKYIG	61	-1231.0405	2808.7793	0.0749	0.9549
HKYI	60	-1235.1490	2811.3128	0.0211	0.9759
K2PIG	58	-1241.1663	2811.9802	0.0151	0.9910
K2PI	57	-1244.5988	2813.1616	0.0084	0.9994
GTRG	64	-1227.5274	2818.8040	0.0005	0.9999
GTRIG	65	-1226.9786	2823.3899	0.0001	1.0000
SYMG	61	-1238.9928	2824.6840	0.0000	1.0000
GTRI	64	-1231.7899	2827.3288	0.0000	1.0000
SYMIG	62	-1238.5210	2829.4239	0.0000	1.0000
SYMI	61	-1241.9566	2830.6115	0.0000	1.0000
JC69G	56	-1269.3569	2856.9943	0.0000	1.0000
F81G	59	-1262.4989	2860.3290	0.0000	1.0000
JC69IG	57	-1268.8676	2861.6992	0.0000	1.0000
JC69I	56	-1272.4082	2863.0969	0.0000	1.0000
F81IG	60	-1261.9625	2864.9398	0.0000	1.0000
F81I	59	-1265.7732	2866.8777	0.0000	1.0000
HKY	59	-1267.7860	2870.9032	0.0000	1.0000
K2P	56	-1277.3873	2873.0552	0.0000	1.0000
GTR	63	-1263.8426	2885.7506	0.0000	1.0000
SYM	60	-1274.4525	2889.9197	0.0000	1.0000
JC69	55	-1303.8759	2920.3486	0.0000	1.0000
F81	58	-1296.8157	2923.2789	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 183-phylip/chr1_15480.phylip = 1 - 294;
 Partition Dummy = 1:183-phylip/chr1_15480.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr1_15480.phylip = 1 - 294;
 Partition Dummy = 1:183-phylip/chr1_15480.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 183-phylip/chr1_15480.phylip = 1 - 294;
 Partition Dummy = 1:183-phylip/chr1_15480.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

