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

Input data from file "261-aggressive-phylip/5622_hg18.phylip" (29 taxa, 242 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	-851.9775	1853.8898	0.4787	0.4787
K2PIG	58	-851.2772	1855.9532	0.1706	0.6493
SYMG	61	-846.0865	1856.1953	0.1512	0.8004
SYMIG	62	-845.0700	1857.7825	0.0684	0.8688
K2PI	57	-854.1088	1858.1524	0.0568	0.9256
HKYIG	61	-847.7001	1859.4224	0.0301	0.9557
SYMI	61	-848.3744	1860.7710	0.0153	0.9711
GTRG	64	-842.9701	1860.9457	0.0141	0.9851
HKYG	60	-850.7984	1862.0388	0.0081	0.9933
GTRIG	65	-842.0406	1862.8313	0.0055	0.9987
HKYI	60	-853.2324	1866.9069	0.0007	0.9994
GTRI	64	-846.2038	1867.4133	0.0006	1.0000
JC69IG	57	-891.7252	1933.3853	0.0000	1.0000
F81IG	60	-888.8732	1938.1883	0.0000	1.0000
JC69G	56	-897.6308	1941.7697	0.0000	1.0000
K2P	56	-898.0695	1942.6470	0.0000	1.0000
F81G	59	-893.2293	1943.3598	0.0000	1.0000
SYM	60	-892.4705	1945.3830	0.0000	1.0000
JC69I	56	-899.6192	1945.7465	0.0000	1.0000
HKY	59	-895.2186	1947.3383	0.0000	1.0000
F81I	59	-895.6024	1948.1060	0.0000	1.0000
GTR	63	-889.3049	1949.9131	0.0000	1.0000
JC69	55	-937.8169	2018.7521	0.0000	1.0000
F81	58	-932.9496	2019.2981	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-842.9701	1813.9401	0.2322	0.2322
GTRIG	65	-842.0406	1814.0813	0.2163	0.4485
SYMIG	62	-845.0700	1814.1400	0.2101	0.6586
SYMG	61	-846.0865	1814.1731	0.2066	0.8652
HKYIG	61	-847.7001	1817.4002	0.0412	0.9064
K2PG	57	-851.9775	1817.9550	0.0312	0.9376
K2PIG	58	-851.2772	1818.5543	0.0231	0.9607
SYMI	61	-848.3744	1818.7488	0.0210	0.9817
GTRI	64	-846.2038	1820.4076	0.0091	0.9908
HKYG	60	-850.7984	1821.5968	0.0050	0.9959
K2PI	57	-854.1088	1822.2176	0.0037	0.9996
HKYI	60	-853.2324	1826.4649	0.0004	1.0000
JC69IG	57	-891.7252	1897.4505	0.0000	1.0000
F81IG	60	-888.8732	1897.7463	0.0000	1.0000
F81G	59	-893.2293	1904.4587	0.0000	1.0000
GTR	63	-889.3049	1904.6097	0.0000	1.0000
SYM	60	-892.4705	1904.9411	0.0000	1.0000
JC69G	56	-897.6308	1907.2616	0.0000	1.0000
K2P	56	-898.0695	1908.1389	0.0000	1.0000
HKY	59	-895.2186	1908.4372	0.0000	1.0000
F81I	59	-895.6024	1909.2049	0.0000	1.0000
JC69I	56	-899.6192	1911.2384	0.0000	1.0000
F81	58	-932.9496	1981.8992	0.0000	1.0000
JC69	55	-937.8169	1985.6338	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-851.9775	2016.8245	0.7947	0.7947
K2PIG	58	-851.2772	2020.9127	0.1029	0.8976
K2PI	57	-854.1088	2021.0871	0.0943	0.9920
SYMG	61	-846.0865	2026.9983	0.0049	0.9969
HKYIG	61	-847.7001	2030.2254	0.0010	0.9978
SYMIG	62	-845.0700	2030.4541	0.0009	0.9987
HKYG	60	-850.7984	2030.9331	0.0007	0.9994
SYMI	61	-848.3744	2031.5740	0.0005	0.9999
HKYI	60	-853.2324	2035.8012	0.0001	1.0000
GTRG	64	-842.9701	2037.2321	0.0000	1.0000
GTRIG	65	-842.0406	2040.8622	0.0000	1.0000
GTRI	64	-846.2038	2043.6997	0.0000	1.0000
JC69IG	57	-891.7252	2096.3199	0.0000	1.0000
JC69G	56	-897.6308	2102.6421	0.0000	1.0000
K2P	56	-898.0695	2103.5194	0.0000	1.0000
JC69I	56	-899.6192	2106.6189	0.0000	1.0000
F81IG	60	-888.8732	2107.0826	0.0000	1.0000
F81G	59	-893.2293	2110.3060	0.0000	1.0000
SYM	60	-892.4705	2114.2773	0.0000	1.0000
HKY	59	-895.2186	2114.2846	0.0000	1.0000
F81I	59	-895.6024	2115.0522	0.0000	1.0000
GTR	63	-889.3049	2124.4128	0.0000	1.0000
JC69	55	-937.8169	2177.5254	0.0000	1.0000
F81	58	-932.9496	2184.2576	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/5622_hg18.phylip = 1 - 242;
 Partition Dummy = 1:261-aggressive-phylip/5622_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/5622_hg18.phylip = 1 - 242;
 Partition Dummy = 1:261-aggressive-phylip/5622_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/5622_hg18.phylip = 1 - 242;
 Partition Dummy = 1:261-aggressive-phylip/5622_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

