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

Input data from file "261-aggressive-phylip/4245_hg18.phylip" (29 taxa, 326 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	-915.6304	1978.8835	0.4326	0.4326
HKYIG	61	-914.7088	1980.0691	0.2391	0.6717
HKYI	60	-916.3052	1980.2330	0.2203	0.8920
GTRG	64	-911.4456	1982.7686	0.0620	0.9540
GTRIG	65	-910.8135	1984.6270	0.0245	0.9785
GTRI	64	-912.5084	1984.8941	0.0214	0.9999
F81G	59	-927.1392	1998.8949	0.0000	1.0000
F81IG	60	-925.9309	1999.4845	0.0000	1.0000
F81I	59	-927.8581	2000.3328	0.0000	1.0000
K2PG	57	-931.6024	2001.8763	0.0000	1.0000
K2PIG	58	-930.3870	2002.4070	0.0000	1.0000
K2PI	57	-932.2933	2003.2582	0.0000	1.0000
SYMG	61	-928.3660	2007.3836	0.0000	1.0000
SYMIG	62	-927.2426	2008.1886	0.0000	1.0000
SYMI	61	-929.1155	2008.8825	0.0000	1.0000
HKY	59	-934.6832	2013.9829	0.0000	1.0000
JC69G	56	-939.8553	2015.4429	0.0000	1.0000
JC69IG	57	-938.5848	2015.8413	0.0000	1.0000
GTR	63	-929.5949	2015.9684	0.0000	1.0000
JC69I	56	-940.5735	2016.8794	0.0000	1.0000
F81	58	-945.9017	2033.4363	0.0000	1.0000
K2P	56	-950.8779	2037.4881	0.0000	1.0000
SYM	60	-947.1205	2041.8637	0.0000	1.0000
JC69	55	-958.8651	2050.5450	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-911.4456	1950.8912	0.2462	0.2462
HKYG	60	-915.6304	1951.2608	0.2047	0.4509
HKYIG	61	-914.7088	1951.4176	0.1893	0.6402
GTRIG	65	-910.8135	1951.6270	0.1705	0.8107
HKYI	60	-916.3052	1952.6104	0.1042	0.9149
GTRI	64	-912.5084	1953.0167	0.0851	1.0000
F81IG	60	-925.9309	1971.8619	0.0000	1.0000
F81G	59	-927.1392	1972.2784	0.0000	1.0000
F81I	59	-927.8581	1973.7162	0.0000	1.0000
K2PIG	58	-930.3870	1976.7740	0.0000	1.0000
K2PG	57	-931.6024	1977.2047	0.0000	1.0000
SYMIG	62	-927.2426	1978.4852	0.0000	1.0000
K2PI	57	-932.2933	1978.5865	0.0000	1.0000
SYMG	61	-928.3660	1978.7320	0.0000	1.0000
SYMI	61	-929.1155	1980.2310	0.0000	1.0000
GTR	63	-929.5949	1985.1897	0.0000	1.0000
HKY	59	-934.6832	1987.3663	0.0000	1.0000
JC69IG	57	-938.5848	1991.1697	0.0000	1.0000
JC69G	56	-939.8553	1991.7106	0.0000	1.0000
JC69I	56	-940.5735	1993.1470	0.0000	1.0000
F81	58	-945.9017	2007.8034	0.0000	1.0000
K2P	56	-950.8779	2013.7557	0.0000	1.0000
SYM	60	-947.1205	2014.2411	0.0000	1.0000
JC69	55	-958.8651	2027.7301	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-915.6304	2178.4747	0.6058	0.6058
HKYI	60	-916.3052	2179.8242	0.3085	0.9143
HKYIG	61	-914.7088	2182.4183	0.0843	0.9986
K2PG	57	-931.6024	2193.0579	0.0004	0.9990
GTRG	64	-911.4456	2193.2527	0.0004	0.9994
K2PI	57	-932.2933	2194.4397	0.0002	0.9996
GTRI	64	-912.5084	2195.3782	0.0001	0.9997
F81G	59	-927.1392	2195.7053	0.0001	0.9998
K2PIG	58	-930.3870	2196.4140	0.0001	0.9999
F81I	59	-927.8581	2197.1432	0.0001	0.9999
GTRIG	65	-910.8135	2197.7753	0.0000	1.0000
F81IG	60	-925.9309	2199.0757	0.0000	1.0000
JC69G	56	-939.8553	2203.7768	0.0000	1.0000
JC69I	56	-940.5735	2205.2133	0.0000	1.0000
JC69IG	57	-938.5848	2207.0228	0.0000	1.0000
SYMG	61	-928.3660	2209.7328	0.0000	1.0000
HKY	59	-934.6832	2210.7933	0.0000	1.0000
SYMI	61	-929.1155	2211.2317	0.0000	1.0000
SYMIG	62	-927.2426	2213.2729	0.0000	1.0000
GTR	63	-929.5949	2223.7643	0.0000	1.0000
K2P	56	-950.8779	2225.8220	0.0000	1.0000
F81	58	-945.9017	2227.4434	0.0000	1.0000
JC69	55	-958.8651	2236.0095	0.0000	1.0000
SYM	60	-947.1205	2241.4549	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 261-aggressive-phylip/4245_hg18.phylip = 1 - 326;
 Partition Dummy = 1:261-aggressive-phylip/4245_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=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 261-aggressive-phylip/4245_hg18.phylip = 1 - 326;
 Partition Dummy = 1:261-aggressive-phylip/4245_hg18.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 261-aggressive-phylip/4245_hg18.phylip = 1 - 326;
 Partition Dummy = 1:261-aggressive-phylip/4245_hg18.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

