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

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

Minimum AIC  model: GTRG
Minimum AICc model: K2PG
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
K2PG	57	-796.4400	1747.6948	0.6431	0.6431
K2PIG	58	-796.1024	1750.7142	0.1421	0.7853
K2PI	57	-798.1868	1751.1883	0.1121	0.8974
HKYG	60	-792.9594	1751.9566	0.0764	0.9737
HKYIG	61	-792.9640	1755.8014	0.0112	0.9849
HKYI	60	-795.0393	1756.1162	0.0095	0.9944
SYMG	61	-794.6036	1759.0807	0.0022	0.9966
GTRG	64	-788.8293	1759.3360	0.0019	0.9985
SYMIG	62	-794.3044	1762.3668	0.0004	0.9989
SYMI	61	-796.3756	1762.6246	0.0004	0.9993
GTRIG	65	-788.9478	1763.6099	0.0002	0.9995
GTRI	64	-791.0955	1763.8685	0.0002	0.9997
JC69G	56	-806.4592	1764.0840	0.0002	0.9999
JC69IG	57	-806.1708	1767.1565	0.0000	1.0000
JC69I	56	-808.1756	1767.5168	0.0000	1.0000
F81G	59	-803.4686	1769.1871	0.0000	1.0000
F81IG	60	-803.4539	1772.9455	0.0000	1.0000
F81I	59	-805.4405	1773.1309	0.0000	1.0000
K2P	56	-818.2233	1787.6122	0.0000	1.0000
HKY	59	-813.7688	1789.7876	0.0000	1.0000
GTR	63	-808.8174	1795.3272	0.0000	1.0000
SYM	60	-815.3568	1796.7512	0.0000	1.0000
JC69	55	-827.7355	1803.0319	0.0000	1.0000
F81	58	-824.1747	1806.8587	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-788.8293	1705.6585	0.2659	0.2659
HKYG	60	-792.9594	1705.9189	0.2334	0.4993
K2PG	57	-796.4400	1706.8800	0.1444	0.6436
GTRIG	65	-788.9478	1707.8956	0.0869	0.7305
HKYIG	61	-792.9640	1707.9280	0.0855	0.8160
K2PIG	58	-796.1024	1708.2049	0.0744	0.8904
HKYI	60	-795.0393	1710.0785	0.0292	0.9196
GTRI	64	-791.0955	1710.1911	0.0276	0.9471
K2PI	57	-798.1868	1710.3735	0.0252	0.9723
SYMG	61	-794.6036	1711.2072	0.0166	0.9889
SYMIG	62	-794.3044	1712.6088	0.0082	0.9971
SYMI	61	-796.3756	1714.7512	0.0028	0.9999
JC69G	56	-806.4592	1724.9184	0.0000	1.0000
F81G	59	-803.4686	1724.9371	0.0000	1.0000
JC69IG	57	-806.1708	1726.3417	0.0000	1.0000
F81IG	60	-803.4539	1726.9078	0.0000	1.0000
JC69I	56	-808.1756	1728.3511	0.0000	1.0000
F81I	59	-805.4405	1728.8809	0.0000	1.0000
GTR	63	-808.8174	1743.6349	0.0000	1.0000
HKY	59	-813.7688	1745.5376	0.0000	1.0000
K2P	56	-818.2233	1748.4466	0.0000	1.0000
SYM	60	-815.3568	1750.7135	0.0000	1.0000
F81	58	-824.1747	1764.3494	0.0000	1.0000
JC69	55	-827.7355	1765.4710	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-796.4400	1900.3168	0.7802	0.7802
K2PI	57	-798.1868	1903.8103	0.1360	0.9162
K2PIG	58	-796.1024	1905.0353	0.0737	0.9899
HKYG	60	-792.9594	1909.5365	0.0078	0.9977
HKYI	60	-795.0393	1913.6962	0.0010	0.9987
HKYIG	61	-792.9640	1914.9393	0.0005	0.9992
JC69G	56	-806.4592	1914.9615	0.0005	0.9997
SYMG	61	-794.6036	1918.2185	0.0001	0.9998
JC69I	56	-808.1756	1918.3943	0.0001	0.9999
JC69IG	57	-806.1708	1919.7784	0.0000	1.0000
SYMI	61	-796.3756	1921.7625	0.0000	1.0000
GTRG	64	-788.8293	1922.8507	0.0000	1.0000
SYMIG	62	-794.3044	1923.0137	0.0000	1.0000
F81G	59	-803.4686	1925.1611	0.0000	1.0000
GTRI	64	-791.0955	1927.3832	0.0000	1.0000
GTRIG	65	-788.9478	1928.4814	0.0000	1.0000
F81I	59	-805.4405	1929.1049	0.0000	1.0000
F81IG	60	-803.4539	1930.5254	0.0000	1.0000
K2P	56	-818.2233	1938.4897	0.0000	1.0000
HKY	59	-813.7688	1945.7616	0.0000	1.0000
JC69	55	-827.7355	1952.1205	0.0000	1.0000
SYM	60	-815.3568	1954.3312	0.0000	1.0000
GTR	63	-808.8174	1957.4334	0.0000	1.0000
F81	58	-824.1747	1961.1798	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/7793_hg18.phylip = 1 - 220;
 Partition Dummy = 1:261-aggressive-phylip/7793_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 (K2PG)]

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


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

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

