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

Input data from file "261-aggressive-phylip/11158_hg18.phylip" (29 taxa, 308 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	-910.0150	1969.6657	0.4212	0.4212
HKYI	60	-910.3304	1970.2965	0.3072	0.7284
GTRG	64	-905.2028	1972.6443	0.0950	0.8234
HKYIG	61	-909.9972	1972.7424	0.0904	0.9138
GTRI	64	-905.5591	1973.3569	0.0665	0.9804
GTRIG	65	-905.1972	1975.8490	0.0191	0.9995
HKY	59	-918.6310	1983.8104	0.0004	0.9998
GTR	63	-913.3905	1985.8301	0.0001	1.0000
F81G	59	-922.3943	1991.3370	0.0000	1.0000
F81I	59	-922.6853	1991.9190	0.0000	1.0000
K2PG	57	-926.7046	1993.8572	0.0000	1.0000
F81IG	60	-922.3619	1994.3593	0.0000	1.0000
K2PI	57	-926.9978	1994.4435	0.0000	1.0000
K2PIG	58	-926.6709	1996.8278	0.0000	1.0000
SYMG	61	-923.0409	1998.8298	0.0000	1.0000
SYMI	61	-923.3043	1999.3566	0.0000	1.0000
SYMIG	62	-923.0315	2001.9488	0.0000	1.0000
F81	58	-931.4257	2006.3374	0.0000	1.0000
K2P	56	-936.1525	2009.7392	0.0000	1.0000
SYM	60	-931.7504	2013.1365	0.0000	1.0000
JC69G	56	-938.3278	2014.0899	0.0000	1.0000
JC69I	56	-938.5975	2014.6293	0.0000	1.0000
JC69IG	57	-938.2944	2017.0368	0.0000	1.0000
JC69	55	-947.5863	2029.6171	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-905.2028	1938.4056	0.3328	0.3328
GTRI	64	-905.5591	1939.1182	0.2330	0.5658
HKYG	60	-910.0150	1940.0301	0.1477	0.7135
GTRIG	65	-905.1972	1940.3945	0.1231	0.8366
HKYI	60	-910.3304	1940.6609	0.1078	0.9444
HKYIG	61	-909.9972	1941.9944	0.0553	0.9997
GTR	63	-913.3905	1952.7809	0.0003	0.9999
HKY	59	-918.6310	1955.2620	0.0001	1.0000
F81G	59	-922.3943	1962.7887	0.0000	1.0000
F81I	59	-922.6853	1963.3706	0.0000	1.0000
F81IG	60	-922.3619	1964.7237	0.0000	1.0000
K2PG	57	-926.7046	1967.4092	0.0000	1.0000
K2PI	57	-926.9978	1967.9955	0.0000	1.0000
SYMG	61	-923.0409	1968.0819	0.0000	1.0000
SYMI	61	-923.3043	1968.6087	0.0000	1.0000
K2PIG	58	-926.6709	1969.3419	0.0000	1.0000
SYMIG	62	-923.0315	1970.0631	0.0000	1.0000
F81	58	-931.4257	1978.8515	0.0000	1.0000
SYM	60	-931.7504	1983.5009	0.0000	1.0000
K2P	56	-936.1525	1984.3049	0.0000	1.0000
JC69G	56	-938.3278	1988.6556	0.0000	1.0000
JC69I	56	-938.5975	1989.1950	0.0000	1.0000
JC69IG	57	-938.2944	1990.5888	0.0000	1.0000
JC69	55	-947.5863	2005.1726	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-910.0150	2163.8361	0.5575	0.5575
HKYI	60	-910.3304	2164.4668	0.4067	0.9642
HKYIG	61	-909.9972	2169.5305	0.0323	0.9966
HKY	59	-918.6310	2175.3379	0.0018	0.9983
GTRG	64	-905.2028	2177.1320	0.0007	0.9991
GTRI	64	-905.5591	2177.8446	0.0005	0.9996
K2PG	57	-926.7046	2180.0249	0.0002	0.9997
K2PI	57	-926.9978	2180.6112	0.0001	0.9999
GTRIG	65	-905.1972	2182.8509	0.0000	0.9999
F81G	59	-922.3943	2182.8645	0.0000	1.0000
F81I	59	-922.6853	2183.4465	0.0000	1.0000
K2PIG	58	-926.6709	2185.6876	0.0000	1.0000
GTR	63	-913.3905	2187.7772	0.0000	1.0000
F81IG	60	-922.3619	2188.5297	0.0000	1.0000
K2P	56	-936.1525	2193.1905	0.0000	1.0000
F81	58	-931.4257	2195.1972	0.0000	1.0000
SYMG	61	-923.0409	2195.6179	0.0000	1.0000
SYMI	61	-923.3043	2196.1448	0.0000	1.0000
JC69G	56	-938.3278	2197.5412	0.0000	1.0000
JC69I	56	-938.5975	2198.0806	0.0000	1.0000
SYMIG	62	-923.0315	2201.3292	0.0000	1.0000
JC69IG	57	-938.2944	2203.2045	0.0000	1.0000
SYM	60	-931.7504	2207.3068	0.0000	1.0000
JC69	55	-947.5863	2210.3281	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/11158_hg18.phylip = 1 - 308;
 Partition Dummy = 1:261-aggressive-phylip/11158_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/11158_hg18.phylip = 1 - 308;
 Partition Dummy = 1:261-aggressive-phylip/11158_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/11158_hg18.phylip = 1 - 308;
 Partition Dummy = 1:261-aggressive-phylip/11158_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

