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

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

Minimum AIC  model: SYMIG
Minimum AICc model: SYMIG
Minimum BIC  model: K2PG


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
SYMIG	62	-986.3927	2132.4566	0.4347	0.4347
SYMG	61	-988.2874	2132.9567	0.3386	0.7733
SYMI	61	-989.6302	2135.6422	0.0884	0.8617
K2PIG	58	-994.7301	2136.1507	0.0686	0.9303
K2PG	57	-996.7980	2137.1138	0.0424	0.9726
K2PI	57	-997.8607	2139.2393	0.0146	0.9873
GTRIG	65	-985.9995	2141.7213	0.0042	0.9915
GTRG	64	-987.8080	2141.9569	0.0038	0.9952
HKYIG	61	-993.3411	2143.0640	0.0022	0.9974
HKYG	60	-995.4131	2143.9483	0.0014	0.9988
GTRI	64	-989.4079	2145.1569	0.0008	0.9996
HKYI	60	-996.5480	2146.2182	0.0004	1.0000
JC69IG	57	-1009.2751	2162.0681	0.0000	1.0000
JC69G	56	-1011.2576	2162.8886	0.0000	1.0000
JC69I	56	-1012.2763	2164.9260	0.0000	1.0000
F81G	59	-1009.2175	2168.3270	0.0000	1.0000
F81IG	60	-1009.0110	2171.1441	0.0000	1.0000
F81I	59	-1012.1056	2174.1030	0.0000	1.0000
SYM	60	-1025.8264	2204.7750	0.0000	1.0000
K2P	56	-1035.5191	2211.4116	0.0000	1.0000
GTR	63	-1024.8463	2212.6835	0.0000	1.0000
HKY	59	-1033.5890	2217.0699	0.0000	1.0000
JC69	55	-1049.6830	2236.6227	0.0000	1.0000
F81	58	-1047.5802	2241.8510	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMIG	62	-986.3927	2096.7854	0.6055	0.6055
SYMG	61	-988.2874	2098.5748	0.2475	0.8529
SYMI	61	-989.6302	2101.2604	0.0646	0.9175
GTRIG	65	-985.9995	2101.9991	0.0447	0.9622
GTRG	64	-987.8080	2103.6159	0.0199	0.9821
K2PIG	58	-994.7301	2105.4601	0.0079	0.9900
GTRI	64	-989.4079	2106.8159	0.0040	0.9940
K2PG	57	-996.7980	2107.5960	0.0027	0.9968
HKYIG	61	-993.3411	2108.6822	0.0016	0.9983
K2PI	57	-997.8607	2109.7214	0.0009	0.9993
HKYG	60	-995.4131	2110.8261	0.0005	0.9998
HKYI	60	-996.5480	2113.0960	0.0002	1.0000
JC69IG	57	-1009.2751	2132.5502	0.0000	1.0000
JC69G	56	-1011.2576	2134.5153	0.0000	1.0000
F81G	59	-1009.2175	2136.4351	0.0000	1.0000
JC69I	56	-1012.2763	2136.5527	0.0000	1.0000
F81IG	60	-1009.0110	2138.0219	0.0000	1.0000
F81I	59	-1012.1056	2142.2111	0.0000	1.0000
SYM	60	-1025.8264	2171.6529	0.0000	1.0000
GTR	63	-1024.8463	2175.6927	0.0000	1.0000
K2P	56	-1035.5191	2183.0382	0.0000	1.0000
HKY	59	-1033.5890	2185.1780	0.0000	1.0000
JC69	55	-1049.6830	2209.3660	0.0000	1.0000
F81	58	-1047.5802	2211.1604	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-996.7980	2315.1847	0.5204	0.5204
K2PIG	58	-994.7301	2316.6907	0.2451	0.7655
K2PI	57	-997.8607	2317.3101	0.1798	0.9453
SYMG	61	-988.2874	2320.7312	0.0325	0.9778
SYMIG	62	-986.3927	2322.5836	0.0129	0.9907
SYMI	61	-989.6302	2323.4167	0.0085	0.9992
HKYG	60	-995.4131	2329.3406	0.0004	0.9996
HKYIG	61	-993.3411	2330.8385	0.0002	0.9998
HKYI	60	-996.5480	2331.6104	0.0001	1.0000
GTRG	64	-987.8080	2336.6980	0.0000	1.0000
JC69G	56	-1011.2576	2338.4621	0.0000	1.0000
GTRIG	65	-985.9995	2338.7230	0.0000	1.0000
GTRI	64	-989.4079	2339.8979	0.0000	1.0000
JC69IG	57	-1009.2751	2340.1389	0.0000	1.0000
JC69I	56	-1012.2763	2340.4995	0.0000	1.0000
F81G	59	-1009.2175	2351.3076	0.0000	1.0000
F81IG	60	-1009.0110	2356.5363	0.0000	1.0000
F81I	59	-1012.1056	2357.0836	0.0000	1.0000
K2P	56	-1035.5191	2386.9850	0.0000	1.0000
SYM	60	-1025.8264	2390.1673	0.0000	1.0000
HKY	59	-1033.5890	2400.0505	0.0000	1.0000
GTR	63	-1024.8463	2405.1328	0.0000	1.0000
JC69	55	-1049.6830	2409.6709	0.0000	1.0000
F81	58	-1047.5802	2422.3910	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11432_hg18.phylip = 1 - 282;
 Partition Dummy = 1:261-aggressive-phylip/11432_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11432_hg18.phylip = 1 - 282;
 Partition Dummy = 1:261-aggressive-phylip/11432_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11432_hg18.phylip = 1 - 282;
 Partition Dummy = 1:261-aggressive-phylip/11432_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;


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

