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

Input data from file "261-aggressive-phylip/11448_hg18.phylip" (29 taxa, 238 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	-794.4623	1739.6579	0.5470	0.5470
HKYG	60	-790.4954	1742.3467	0.1426	0.6896
K2PI	57	-795.8396	1742.4125	0.1380	0.8276
K2PIG	58	-794.3767	1742.9880	0.1035	0.9311
HKYI	60	-791.8883	1745.1325	0.0354	0.9665
HKYIG	61	-790.4003	1745.7780	0.0256	0.9922
GTRG	64	-786.9147	1749.9219	0.0032	0.9954
SYMG	61	-793.1090	1751.1953	0.0017	0.9971
K2P	56	-802.7605	1752.7917	0.0008	0.9979
GTRI	64	-788.6509	1753.3943	0.0006	0.9985
GTRIG	65	-786.7872	1753.4581	0.0006	0.9990
SYMI	61	-794.6417	1754.2607	0.0004	0.9994
SYMIG	62	-793.0182	1754.6764	0.0003	0.9997
HKY	59	-799.0044	1755.7841	0.0002	0.9999
JC69G	56	-805.0864	1757.4436	0.0001	0.9999
JC69I	56	-806.3419	1759.9544	0.0000	1.0000
F81G	59	-801.2898	1760.3549	0.0000	1.0000
JC69IG	57	-805.0180	1760.7693	0.0000	1.0000
F81I	59	-802.6225	1763.0203	0.0000	1.0000
F81IG	60	-801.2262	1763.8084	0.0000	1.0000
SYM	60	-801.3204	1763.9967	0.0000	1.0000
GTR	63	-795.8515	1764.0478	0.0000	1.0000
JC69	55	-813.0996	1770.0454	0.0000	1.0000
F81	58	-809.2639	1772.7624	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-790.4954	1700.9907	0.2963	0.2963
GTRG	64	-786.9147	1701.8294	0.1948	0.4912
HKYIG	61	-790.4003	1702.8007	0.1199	0.6111
K2PG	57	-794.4623	1702.9246	0.1127	0.7238
GTRIG	65	-786.7872	1703.5744	0.0814	0.8052
HKYI	60	-791.8883	1703.7765	0.0736	0.8788
K2PIG	58	-794.3767	1704.7534	0.0452	0.9240
GTRI	64	-788.6509	1705.3018	0.0343	0.9583
K2PI	57	-795.8396	1705.6792	0.0284	0.9867
SYMG	61	-793.1090	1708.2180	0.0080	0.9947
SYMIG	62	-793.0182	1710.0364	0.0032	0.9979
SYMI	61	-794.6417	1711.2834	0.0017	0.9996
HKY	59	-799.0044	1716.0088	0.0002	0.9998
K2P	56	-802.7605	1717.5210	0.0001	0.9999
GTR	63	-795.8515	1717.7030	0.0001	1.0000
F81G	59	-801.2898	1720.5797	0.0000	1.0000
JC69G	56	-805.0864	1722.1729	0.0000	1.0000
F81IG	60	-801.2262	1722.4524	0.0000	1.0000
SYM	60	-801.3204	1722.6408	0.0000	1.0000
F81I	59	-802.6225	1723.2450	0.0000	1.0000
JC69IG	57	-805.0180	1724.0360	0.0000	1.0000
JC69I	56	-806.3419	1724.6837	0.0000	1.0000
F81	58	-809.2639	1734.5277	0.0000	1.0000
JC69	55	-813.0996	1736.1992	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-794.4623	1900.8440	0.7427	0.7427
K2PI	57	-795.8396	1903.5986	0.1874	0.9301
K2PIG	58	-794.3767	1906.1451	0.0524	0.9826
HKYG	60	-790.4954	1909.3270	0.0107	0.9932
K2P	56	-802.7605	1911.9682	0.0029	0.9961
HKYI	60	-791.8883	1912.1128	0.0027	0.9987
HKYIG	61	-790.4003	1914.6092	0.0008	0.9995
JC69G	56	-805.0864	1916.6200	0.0003	0.9998
JC69I	56	-806.3419	1919.1309	0.0001	0.9999
SYMG	61	-793.1090	1920.0265	0.0001	0.9999
HKY	59	-799.0044	1920.8727	0.0000	1.0000
JC69IG	57	-805.0180	1921.9554	0.0000	1.0000
SYMI	61	-794.6417	1923.0919	0.0000	1.0000
GTRG	64	-786.9147	1924.0547	0.0000	1.0000
SYMIG	62	-793.0182	1925.3172	0.0000	1.0000
F81G	59	-801.2898	1925.4436	0.0000	1.0000
JC69	55	-813.0996	1927.1741	0.0000	1.0000
GTRI	64	-788.6509	1927.5271	0.0000	1.0000
F81I	59	-802.6225	1928.1089	0.0000	1.0000
GTRIG	65	-786.7872	1929.2720	0.0000	1.0000
F81IG	60	-801.2262	1930.7887	0.0000	1.0000
SYM	60	-801.3204	1930.9770	0.0000	1.0000
F81	58	-809.2639	1935.9194	0.0000	1.0000
GTR	63	-795.8515	1936.4561	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/11448_hg18.phylip = 1 - 238;
 Partition Dummy = 1:261-aggressive-phylip/11448_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/11448_hg18.phylip = 1 - 238;
 Partition Dummy = 1:261-aggressive-phylip/11448_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/11448_hg18.phylip = 1 - 238;
 Partition Dummy = 1:261-aggressive-phylip/11448_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

