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

Input data from file "261-aggressive-phylip/10312_hg18.phylip" (29 taxa, 375 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	-882.1304	1907.5730	0.3571	0.3571
HKYI	60	-882.3964	1908.1049	0.2737	0.6308
GTRG	64	-877.4182	1909.6751	0.1248	0.7556
HKYIG	61	-881.8290	1909.8241	0.1159	0.8715
GTRI	64	-877.8380	1910.5147	0.0820	0.9535
GTRIG	65	-876.9464	1911.6598	0.0463	0.9998
HKY	59	-891.0907	1922.6577	0.0002	1.0000
GTR	63	-886.9278	1925.7849	0.0000	1.0000
K2PG	57	-910.9156	1956.6893	0.0000	1.0000
K2PI	57	-911.1146	1957.0873	0.0000	1.0000
K2PIG	58	-910.5497	1958.7576	0.0000	1.0000
SYMG	61	-907.4436	1961.0534	0.0000	1.0000
SYMI	61	-907.6581	1961.4823	0.0000	1.0000
F81G	59	-911.2228	1962.9218	0.0000	1.0000
SYMIG	62	-906.9809	1963.0003	0.0000	1.0000
F81I	59	-911.3451	1963.1664	0.0000	1.0000
F81IG	60	-910.7887	1964.8895	0.0000	1.0000
K2P	56	-920.4804	1973.0363	0.0000	1.0000
F81	58	-919.6064	1976.8709	0.0000	1.0000
SYM	60	-917.1954	1977.7030	0.0000	1.0000
JC69G	56	-936.6164	2005.3082	0.0000	1.0000
JC69I	56	-936.8482	2005.7718	0.0000	1.0000
JC69IG	57	-936.2896	2007.4372	0.0000	1.0000
JC69	55	-945.6827	2020.6757	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-877.4182	1882.8364	0.2978	0.2978
GTRI	64	-877.8380	1883.6760	0.1957	0.4935
GTRIG	65	-876.9464	1883.8928	0.1756	0.6692
HKYG	60	-882.1304	1884.2609	0.1461	0.8153
HKYI	60	-882.3964	1884.7928	0.1120	0.9272
HKYIG	61	-881.8290	1885.6579	0.0727	0.9999
GTR	63	-886.9278	1899.8556	0.0001	0.9999
HKY	59	-891.0907	1900.1815	0.0001	1.0000
K2PG	57	-910.9156	1935.8312	0.0000	1.0000
K2PI	57	-911.1146	1936.2292	0.0000	1.0000
SYMG	61	-907.4436	1936.8872	0.0000	1.0000
K2PIG	58	-910.5497	1937.0993	0.0000	1.0000
SYMI	61	-907.6581	1937.3161	0.0000	1.0000
SYMIG	62	-906.9809	1937.9618	0.0000	1.0000
F81G	59	-911.2228	1940.4456	0.0000	1.0000
F81I	59	-911.3451	1940.6902	0.0000	1.0000
F81IG	60	-910.7887	1941.5774	0.0000	1.0000
K2P	56	-920.4804	1952.9609	0.0000	1.0000
SYM	60	-917.1954	1954.3909	0.0000	1.0000
F81	58	-919.6064	1955.2127	0.0000	1.0000
JC69G	56	-936.6164	1985.2327	0.0000	1.0000
JC69I	56	-936.8482	1985.6963	0.0000	1.0000
JC69IG	57	-936.2896	1986.5791	0.0000	1.0000
JC69	55	-945.6827	2001.3654	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-882.1304	2119.8764	0.5434	0.5434
HKYI	60	-882.3964	2120.4083	0.4165	0.9600
HKYIG	61	-881.8290	2125.2004	0.0379	0.9979
HKY	59	-891.0907	2131.8701	0.0014	0.9993
GTRG	64	-877.4182	2134.1597	0.0004	0.9997
GTRI	64	-877.8380	2134.9992	0.0003	1.0000
GTRIG	65	-876.9464	2139.1430	0.0000	1.0000
GTR	63	-886.9278	2147.2520	0.0000	1.0000
K2PG	57	-910.9156	2159.6660	0.0000	1.0000
K2PI	57	-911.1146	2160.0640	0.0000	1.0000
K2PIG	58	-910.5497	2164.8610	0.0000	1.0000
F81G	59	-911.2228	2172.1343	0.0000	1.0000
F81I	59	-911.3451	2172.3789	0.0000	1.0000
K2P	56	-920.4804	2172.8687	0.0000	1.0000
SYMG	61	-907.4436	2176.4297	0.0000	1.0000
SYMI	61	-907.6581	2176.8586	0.0000	1.0000
F81IG	60	-910.7887	2177.1929	0.0000	1.0000
SYMIG	62	-906.9809	2181.4312	0.0000	1.0000
F81	58	-919.6064	2182.9744	0.0000	1.0000
SYM	60	-917.1954	2190.0064	0.0000	1.0000
JC69G	56	-936.6164	2205.1406	0.0000	1.0000
JC69I	56	-936.8482	2205.6042	0.0000	1.0000
JC69IG	57	-936.2896	2210.4139	0.0000	1.0000
JC69	55	-945.6827	2217.3463	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/10312_hg18.phylip = 1 - 375;
 Partition Dummy = 1:261-aggressive-phylip/10312_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/10312_hg18.phylip = 1 - 375;
 Partition Dummy = 1:261-aggressive-phylip/10312_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/10312_hg18.phylip = 1 - 375;
 Partition Dummy = 1:261-aggressive-phylip/10312_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

