
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:41 08/01/2011

Input data from file "261-aggressive-phylip/7434_hg18.phylip" (29 taxa, 408 characters)

Minimum AIC  model: SYMG
Minimum AICc model: K2PG
Minimum BIC  model: K2P


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
K2PG	57	-1284.6705	2702.2325	0.3104	0.3104
K2PI	57	-1284.9209	2702.7331	0.2416	0.5520
SYMG	61	-1280.3320	2704.5252	0.0986	0.6506
K2PIG	58	-1284.6951	2705.0006	0.0778	0.7284
SYMI	61	-1280.6366	2705.1344	0.0727	0.8011
K2P	56	-1287.6579	2705.5038	0.0605	0.8616
HKYG	60	-1282.8039	2706.7029	0.0332	0.8948
HKYI	60	-1283.0451	2707.1852	0.0261	0.9209
SYMIG	62	-1280.3470	2707.3375	0.0242	0.9451
SYM	60	-1283.2614	2707.6179	0.0210	0.9661
GTRG	64	-1278.4920	2709.2406	0.0093	0.9754
HKYIG	61	-1282.8266	2709.5145	0.0081	0.9836
GTRI	64	-1278.7866	2709.8297	0.0070	0.9905
HKY	59	-1285.9454	2710.2356	0.0057	0.9962
GTRIG	65	-1278.5008	2712.0893	0.0022	0.9984
GTR	63	-1281.6834	2712.8087	0.0016	1.0000
JC69G	56	-1297.8078	2725.8037	0.0000	1.0000
JC69I	56	-1298.1064	2726.4009	0.0000	1.0000
JC69IG	57	-1297.8268	2728.5450	0.0000	1.0000
JC69	55	-1300.8164	2729.1328	0.0000	1.0000
F81G	59	-1296.1545	2730.6538	0.0000	1.0000
F81I	59	-1296.4376	2731.2200	0.0000	1.0000
F81IG	60	-1296.1766	2733.4484	0.0000	1.0000
F81	58	-1299.2496	2734.1094	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	61	-1280.3320	2682.6640	0.1968	0.1968
SYMI	61	-1280.6366	2683.2732	0.1451	0.3419
K2PG	57	-1284.6705	2683.3411	0.1403	0.4821
K2PI	57	-1284.9209	2683.8417	0.1092	0.5913
SYMIG	62	-1280.3470	2684.6940	0.0713	0.6627
GTRG	64	-1278.4920	2684.9840	0.0617	0.7243
K2PIG	58	-1284.6951	2685.3903	0.0503	0.7747
GTRI	64	-1278.7866	2685.5732	0.0459	0.8206
HKYG	60	-1282.8039	2685.6078	0.0452	0.8658
HKYI	60	-1283.0451	2686.0901	0.0355	0.9013
SYM	60	-1283.2614	2686.5228	0.0286	0.9298
GTRIG	65	-1278.5008	2687.0016	0.0225	0.9523
K2P	56	-1287.6579	2687.3158	0.0192	0.9716
HKYIG	61	-1282.8266	2687.6532	0.0162	0.9878
GTR	63	-1281.6834	2689.3668	0.0069	0.9947
HKY	59	-1285.9454	2689.8908	0.0053	1.0000
JC69G	56	-1297.8078	2707.6157	0.0000	1.0000
JC69I	56	-1298.1064	2708.2128	0.0000	1.0000
JC69IG	57	-1297.8268	2709.6535	0.0000	1.0000
F81G	59	-1296.1545	2710.3090	0.0000	1.0000
F81I	59	-1296.4376	2710.8752	0.0000	1.0000
JC69	55	-1300.8164	2711.6328	0.0000	1.0000
F81IG	60	-1296.1766	2712.3533	0.0000	1.0000
F81	58	-1299.2496	2714.4991	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-1287.6579	2911.9467	0.3575	0.3575
K2PG	57	-1284.6705	2911.9833	0.3510	0.7085
K2PI	57	-1284.9209	2912.4839	0.2733	0.9818
K2PIG	58	-1284.6951	2918.0438	0.0170	0.9988
HKYG	60	-1282.8039	2926.2838	0.0003	0.9990
HKY	59	-1285.9454	2926.5556	0.0002	0.9993
HKYI	60	-1283.0451	2926.7662	0.0002	0.9995
SYM	60	-1283.2614	2927.1988	0.0002	0.9997
SYMG	61	-1280.3320	2927.3513	0.0002	0.9998
SYMI	61	-1280.6366	2927.9605	0.0001	0.9999
JC69G	56	-1297.8078	2932.2466	0.0000	1.0000
JC69	55	-1300.8164	2932.2525	0.0000	1.0000
HKYIG	61	-1282.8266	2932.3405	0.0000	1.0000
JC69I	56	-1298.1064	2932.8438	0.0000	1.0000
SYMIG	62	-1280.3470	2933.3926	0.0000	1.0000
JC69IG	57	-1297.8268	2938.2958	0.0000	1.0000
GTRG	64	-1278.4920	2941.7051	0.0000	1.0000
GTR	63	-1281.6834	2942.0767	0.0000	1.0000
GTRI	64	-1278.7866	2942.2943	0.0000	1.0000
F81G	59	-1296.1545	2946.9737	0.0000	1.0000
F81	58	-1299.2496	2947.1526	0.0000	1.0000
F81I	59	-1296.4376	2947.5399	0.0000	1.0000
GTRIG	65	-1278.5008	2947.7340	0.0000	1.0000
F81IG	60	-1296.1766	2953.0293	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7434_hg18.phylip = 1 - 408;
 Partition Dummy = 1:261-aggressive-phylip/7434_hg18.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7434_hg18.phylip = 1 - 408;
 Partition Dummy = 1:261-aggressive-phylip/7434_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7434_hg18.phylip = 1 - 408;
 Partition Dummy = 1:261-aggressive-phylip/7434_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

