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

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

Minimum AIC  model: HKYI
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	57	-775.0951	1692.8136	0.2598	0.2598
HKYI	60	-770.4778	1693.0609	0.2295	0.4893
K2PG	57	-775.2595	1693.1423	0.2204	0.7097
HKYG	60	-770.6867	1693.4787	0.1863	0.8960
K2PIG	58	-775.2592	1696.2749	0.0460	0.9420
HKYIG	61	-770.7587	1696.8391	0.0347	0.9767
GTRI	64	-767.1802	1699.5033	0.0092	0.9859
GTRG	64	-767.4101	1699.9631	0.0073	0.9931
SYMI	61	-773.2680	1701.8577	0.0028	0.9960
SYMG	61	-773.4664	1702.2544	0.0023	0.9983
GTRIG	65	-767.5515	1703.5783	0.0012	0.9995
SYMIG	62	-773.4662	1705.4988	0.0005	0.9999
HKY	59	-780.9276	1710.7722	0.0000	1.0000
K2P	56	-785.9717	1711.4606	0.0000	1.0000
GTR	63	-776.3847	1714.6094	0.0000	1.0000
SYM	60	-783.1181	1718.3414	0.0000	1.0000
JC69I	56	-790.2222	1719.9615	0.0000	1.0000
JC69G	56	-790.4265	1720.3703	0.0000	1.0000
F81I	59	-785.8661	1720.6492	0.0000	1.0000
F81G	59	-786.0999	1721.1169	0.0000	1.0000
JC69IG	57	-790.4263	1723.4760	0.0000	1.0000
F81IG	60	-786.1881	1724.4814	0.0000	1.0000
F81	58	-795.8854	1737.5272	0.0000	1.0000
JC69	55	-800.6262	1737.6901	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-770.4778	1660.9557	0.2793	0.2793
HKYG	60	-770.6867	1661.3734	0.2266	0.5059
GTRI	64	-767.1802	1662.3604	0.1383	0.6442
GTRG	64	-767.4101	1662.8202	0.1099	0.7541
HKYIG	61	-770.7587	1663.5175	0.0776	0.8317
K2PI	57	-775.0951	1664.1902	0.0554	0.8871
K2PG	57	-775.2595	1664.5189	0.0470	0.9341
GTRIG	65	-767.5515	1665.1030	0.0351	0.9692
K2PIG	58	-775.2592	1666.5184	0.0173	0.9866
SYMI	61	-773.2680	1668.5361	0.0063	0.9929
SYMG	61	-773.4664	1668.9328	0.0052	0.9980
SYMIG	62	-773.4662	1670.9324	0.0019	0.9999
GTR	63	-776.3847	1678.7694	0.0000	1.0000
HKY	59	-780.9276	1679.8552	0.0000	1.0000
K2P	56	-785.9717	1683.9433	0.0000	1.0000
SYM	60	-783.1181	1686.2361	0.0000	1.0000
F81I	59	-785.8661	1689.7322	0.0000	1.0000
F81G	59	-786.0999	1690.1998	0.0000	1.0000
F81IG	60	-786.1881	1692.3762	0.0000	1.0000
JC69I	56	-790.2222	1692.4443	0.0000	1.0000
JC69G	56	-790.4265	1692.8530	0.0000	1.0000
JC69IG	57	-790.4263	1694.8526	0.0000	1.0000
F81	58	-795.8854	1707.7707	0.0000	1.0000
JC69	55	-800.6262	1711.2523	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-775.0951	1873.1765	0.5162	0.5162
K2PG	57	-775.2595	1873.5053	0.4380	0.9542
K2PIG	58	-775.2592	1879.1711	0.0258	0.9800
HKYI	60	-770.4778	1880.9413	0.0106	0.9906
HKYG	60	-770.6867	1881.3590	0.0086	0.9993
HKYIG	61	-770.7587	1887.1695	0.0005	0.9997
K2P	56	-785.9717	1889.2632	0.0002	0.9999
SYMI	61	-773.2680	1892.1881	0.0000	1.0000
SYMG	61	-773.4664	1892.5849	0.0000	1.0000
HKY	59	-780.9276	1896.1744	0.0000	1.0000
GTRI	64	-767.1802	1897.0117	0.0000	1.0000
GTRG	64	-767.4101	1897.4715	0.0000	1.0000
JC69I	56	-790.2222	1897.7642	0.0000	1.0000
JC69G	56	-790.4265	1898.1729	0.0000	1.0000
SYMIG	62	-773.4662	1898.2508	0.0000	1.0000
GTRIG	65	-767.5515	1903.4207	0.0000	1.0000
JC69IG	57	-790.4263	1903.8389	0.0000	1.0000
F81I	59	-785.8661	1906.0513	0.0000	1.0000
SYM	60	-783.1181	1906.2217	0.0000	1.0000
F81G	59	-786.0999	1906.5190	0.0000	1.0000
GTR	63	-776.3847	1909.7542	0.0000	1.0000
F81IG	60	-786.1881	1912.3618	0.0000	1.0000
JC69	55	-800.6262	1912.9058	0.0000	1.0000
F81	58	-795.8854	1920.4235	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10090_hg18.phylip = 1 - 289;
 Partition Dummy = 1:261-aggressive-phylip/10090_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (K2PI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10090_hg18.phylip = 1 - 289;
 Partition Dummy = 1:261-aggressive-phylip/10090_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10090_hg18.phylip = 1 - 289;
 Partition Dummy = 1:261-aggressive-phylip/10090_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

