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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-899.3794	1961.4255	0.3752	0.3752
HKYG	60	-906.1055	1962.2111	0.2533	0.6285
GTRI	64	-900.5014	1963.6694	0.1222	0.7507
HKYI	60	-907.0955	1964.1910	0.0941	0.8448
GTRIG	65	-899.2172	1964.3340	0.0876	0.9324
HKYIG	61	-905.8638	1964.8552	0.0675	1.0000
GTR	63	-910.8107	1981.0819	0.0000	1.0000
HKY	59	-919.3667	1985.6313	0.0000	1.0000
K2PG	57	-925.1123	1990.9938	0.0000	1.0000
K2PI	57	-925.7593	1992.2878	0.0000	1.0000
K2PIG	58	-924.8168	1993.4548	0.0000	1.0000
SYMG	61	-921.2405	1995.6086	0.0000	1.0000
SYMI	61	-921.9026	1996.9329	0.0000	1.0000
SYMIG	62	-921.0160	1998.3129	0.0000	1.0000
K2P	56	-940.9109	2019.5637	0.0000	1.0000
F81G	59	-936.8151	2020.5282	0.0000	1.0000
SYM	60	-935.4199	2020.8398	0.0000	1.0000
F81I	59	-937.5386	2021.9751	0.0000	1.0000
F81IG	60	-936.6037	2023.2074	0.0000	1.0000
JC69G	56	-953.2424	2044.2267	0.0000	1.0000
JC69I	56	-953.8531	2045.4481	0.0000	1.0000
F81	58	-951.0096	2045.8403	0.0000	1.0000
JC69IG	57	-953.0032	2046.7756	0.0000	1.0000
JC69	55	-968.8215	2072.3819	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-899.3794	1926.7588	0.5323	0.5323
GTRIG	65	-899.2172	1928.4344	0.2303	0.7625
GTRI	64	-900.5014	1929.0027	0.1733	0.9359
HKYG	60	-906.1055	1932.2111	0.0348	0.9707
HKYIG	61	-905.8638	1933.7277	0.0163	0.9870
HKYI	60	-907.0955	1934.1910	0.0129	1.0000
GTR	63	-910.8107	1947.6213	0.0000	1.0000
HKY	59	-919.3667	1956.7334	0.0000	1.0000
K2PG	57	-925.1123	1964.2245	0.0000	1.0000
SYMG	61	-921.2405	1964.4810	0.0000	1.0000
K2PI	57	-925.7593	1965.5186	0.0000	1.0000
K2PIG	58	-924.8168	1965.6337	0.0000	1.0000
SYMI	61	-921.9026	1965.8053	0.0000	1.0000
SYMIG	62	-921.0160	1966.0319	0.0000	1.0000
SYM	60	-935.4199	1990.8398	0.0000	1.0000
F81G	59	-936.8151	1991.6303	0.0000	1.0000
F81I	59	-937.5386	1993.0772	0.0000	1.0000
F81IG	60	-936.6037	1993.2074	0.0000	1.0000
K2P	56	-940.9109	1993.8218	0.0000	1.0000
F81	58	-951.0096	2018.0191	0.0000	1.0000
JC69G	56	-953.2424	2018.4847	0.0000	1.0000
JC69I	56	-953.8531	2019.7062	0.0000	1.0000
JC69IG	57	-953.0032	2020.0063	0.0000	1.0000
JC69	55	-968.8215	2047.6430	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-906.1055	2155.4298	0.6863	0.6863
HKYI	60	-907.0955	2157.4097	0.2550	0.9413
HKYIG	61	-905.8638	2160.6667	0.0500	0.9914
GTRG	64	-899.3794	2164.8588	0.0062	0.9975
GTRI	64	-900.5014	2167.1027	0.0020	0.9995
GTRIG	65	-899.2172	2170.2547	0.0004	0.9999
HKY	59	-919.3667	2176.2318	0.0000	1.0000
K2PG	57	-925.1123	2176.2823	0.0000	1.0000
K2PI	57	-925.7593	2177.5763	0.0000	1.0000
K2PIG	58	-924.8168	2181.4117	0.0000	1.0000
GTR	63	-910.8107	2182.0010	0.0000	1.0000
SYMG	61	-921.2405	2191.4200	0.0000	1.0000
SYMI	61	-921.9026	2192.7443	0.0000	1.0000
SYMIG	62	-921.0160	2196.6913	0.0000	1.0000
K2P	56	-940.9109	2202.1592	0.0000	1.0000
F81G	59	-936.8151	2211.1287	0.0000	1.0000
F81I	59	-937.5386	2212.5756	0.0000	1.0000
SYM	60	-935.4199	2214.0585	0.0000	1.0000
F81IG	60	-936.6037	2216.4261	0.0000	1.0000
JC69G	56	-953.2424	2226.8222	0.0000	1.0000
JC69I	56	-953.8531	2228.0437	0.0000	1.0000
JC69IG	57	-953.0032	2232.0641	0.0000	1.0000
F81	58	-951.0096	2233.7972	0.0000	1.0000
JC69	55	-968.8215	2252.2601	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/7777_hg18.phylip = 1 - 305;
 Partition Dummy = 1:261-aggressive-phylip/7777_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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7777_hg18.phylip = 1 - 305;
 Partition Dummy = 1:261-aggressive-phylip/7777_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 BIC model (HKYG)]

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

