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

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

Minimum AIC  model: HKYG
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	-783.1227	1720.8099	0.3720	0.3720
HKYG	60	-777.8654	1721.4808	0.2660	0.6380
K2PIG	58	-782.4052	1723.0573	0.1209	0.7589
K2PI	57	-784.4639	1723.4922	0.0973	0.8562
HKYIG	61	-777.2082	1723.9888	0.0759	0.9321
HKYI	60	-779.3518	1724.4536	0.0602	0.9922
SYMG	61	-780.2420	1730.0564	0.0037	0.9959
GTRG	64	-775.1963	1731.7260	0.0016	0.9975
SYMI	61	-781.4547	1732.4817	0.0011	0.9986
SYMIG	62	-779.9086	1733.2603	0.0007	0.9993
GTRI	64	-776.5396	1734.4126	0.0004	0.9997
GTRIG	65	-774.9001	1735.1550	0.0003	1.0000
HKY	59	-796.4147	1754.8046	0.0000	1.0000
K2P	56	-802.6900	1756.3068	0.0000	1.0000
SYM	60	-798.1153	1761.9806	0.0000	1.0000
GTR	63	-792.8039	1762.9709	0.0000	1.0000
F81G	59	-804.5587	1771.0925	0.0000	1.0000
JC69G	56	-810.3998	1771.7264	0.0000	1.0000
F81I	59	-805.1356	1772.2464	0.0000	1.0000
JC69I	56	-810.9717	1772.8702	0.0000	1.0000
F81IG	60	-804.4788	1774.7077	0.0000	1.0000
JC69IG	57	-810.3067	1775.1779	0.0000	1.0000
F81	58	-823.2376	1804.7221	0.0000	1.0000
JC69	55	-828.8333	1805.0000	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-777.8654	1675.7308	0.3742	0.3742
HKYIG	61	-777.2082	1676.4164	0.2656	0.6397
GTRG	64	-775.1963	1678.3927	0.0989	0.7386
HKYI	60	-779.3518	1678.7036	0.0846	0.8232
GTRIG	65	-774.9001	1679.8001	0.0489	0.8721
K2PG	57	-783.1227	1680.2455	0.0391	0.9113
K2PIG	58	-782.4052	1680.8104	0.0295	0.9408
GTRI	64	-776.5396	1681.0792	0.0258	0.9666
SYMG	61	-780.2420	1682.4841	0.0128	0.9794
K2PI	57	-784.4639	1682.9278	0.0102	0.9896
SYMIG	62	-779.9086	1683.8173	0.0066	0.9962
SYMI	61	-781.4547	1684.9094	0.0038	1.0000
HKY	59	-796.4147	1710.8295	0.0000	1.0000
GTR	63	-792.8039	1711.6078	0.0000	1.0000
SYM	60	-798.1153	1716.2306	0.0000	1.0000
K2P	56	-802.6900	1717.3799	0.0000	1.0000
F81G	59	-804.5587	1727.1173	0.0000	1.0000
F81I	59	-805.1356	1728.2713	0.0000	1.0000
F81IG	60	-804.4788	1728.9577	0.0000	1.0000
JC69G	56	-810.3998	1732.7996	0.0000	1.0000
JC69I	56	-810.9717	1733.9434	0.0000	1.0000
JC69IG	57	-810.3067	1734.6135	0.0000	1.0000
F81	58	-823.2376	1762.4752	0.0000	1.0000
JC69	55	-828.8333	1767.6667	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-783.1227	1873.9407	0.6761	0.6761
K2PI	57	-784.4639	1876.6230	0.1768	0.8529
K2PIG	58	-782.4052	1877.9039	0.0932	0.9461
HKYG	60	-777.8654	1879.6205	0.0395	0.9856
HKYI	60	-779.3518	1882.5934	0.0089	0.9945
HKYIG	61	-777.2082	1883.7044	0.0051	0.9996
SYMG	61	-780.2420	1889.7720	0.0002	0.9999
SYMI	61	-781.4547	1892.1973	0.0001	1.0000
SYMIG	62	-779.9086	1894.5034	0.0000	1.0000
GTRG	64	-775.1963	1895.8751	0.0000	1.0000
GTRI	64	-776.5396	1898.5616	0.0000	1.0000
GTRIG	65	-774.9001	1900.6807	0.0000	1.0000
K2P	56	-802.6900	1907.6771	0.0000	1.0000
HKY	59	-796.4147	1911.3211	0.0000	1.0000
SYM	60	-798.1153	1920.1204	0.0000	1.0000
JC69G	56	-810.3998	1923.0967	0.0000	1.0000
JC69I	56	-810.9717	1924.2405	0.0000	1.0000
GTR	63	-792.8039	1925.6921	0.0000	1.0000
F81G	59	-804.5587	1927.6089	0.0000	1.0000
JC69IG	57	-810.3067	1928.3087	0.0000	1.0000
F81I	59	-805.1356	1928.7629	0.0000	1.0000
F81IG	60	-804.4788	1932.8475	0.0000	1.0000
JC69	55	-828.8333	1954.5656	0.0000	1.0000
F81	58	-823.2376	1959.5686	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11592_hg18.phylip = 1 - 221;
 Partition Dummy = 1:261-aggressive-phylip/11592_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 AICc model (K2PG)]

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

