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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	60	-757.1327	1674.9320	0.6158	0.6158
HKYI	60	-758.0379	1676.7425	0.2491	0.8649
HKYIG	61	-757.0087	1678.2743	0.1158	0.9807
GTRG	64	-754.0721	1683.4168	0.0089	0.9895
GTRI	64	-755.1262	1685.5251	0.0031	0.9926
K2PG	57	-768.0550	1686.2412	0.0022	0.9948
HKY	59	-764.8851	1686.8862	0.0016	0.9963
GTRIG	65	-753.9419	1686.9123	0.0015	0.9979
K2PI	57	-768.7005	1687.5321	0.0011	0.9990
K2PIG	58	-767.9453	1689.4950	0.0004	0.9994
F81G	59	-766.6045	1690.3250	0.0003	0.9997
F81I	59	-767.3901	1691.8962	0.0001	0.9999
F81IG	60	-766.4948	1693.6563	0.0001	0.9999
SYMG	61	-765.0490	1694.3550	0.0000	0.9999
GTR	63	-761.8509	1695.2610	0.0000	1.0000
SYMI	61	-765.7021	1695.6613	0.0000	1.0000
SYMIG	62	-764.9574	1697.8025	0.0000	1.0000
K2P	56	-775.6941	1698.0838	0.0000	1.0000
F81	58	-774.0001	1701.6047	0.0000	1.0000
JC69G	56	-777.5170	1701.7296	0.0000	1.0000
JC69I	56	-778.2044	1703.1045	0.0000	1.0000
JC69IG	57	-777.4207	1704.9725	0.0000	1.0000
SYM	60	-773.0418	1706.7503	0.0000	1.0000
JC69	55	-784.9122	1713.1218	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-757.1327	1634.2653	0.3977	0.3977
HKYIG	61	-757.0087	1636.0174	0.1656	0.5633
HKYI	60	-758.0379	1636.0759	0.1608	0.7241
GTRG	64	-754.0721	1636.1441	0.1554	0.8795
GTRIG	65	-753.9419	1637.8837	0.0651	0.9447
GTRI	64	-755.1262	1638.2524	0.0542	0.9988
HKY	59	-764.8851	1647.7702	0.0005	0.9993
GTR	63	-761.8509	1649.7017	0.0002	0.9995
K2PG	57	-768.0550	1650.1101	0.0001	0.9996
F81G	59	-766.6045	1651.2090	0.0001	0.9997
K2PI	57	-768.7005	1651.4009	0.0001	0.9998
K2PIG	58	-767.9453	1651.8906	0.0001	0.9998
SYMG	61	-765.0490	1652.0980	0.0001	0.9999
F81I	59	-767.3901	1652.7802	0.0000	0.9999
F81IG	60	-766.4948	1652.9897	0.0000	1.0000
SYMI	61	-765.7021	1653.4043	0.0000	1.0000
SYMIG	62	-764.9574	1653.9149	0.0000	1.0000
K2P	56	-775.6941	1663.3882	0.0000	1.0000
F81	58	-774.0001	1664.0003	0.0000	1.0000
SYM	60	-773.0418	1666.0836	0.0000	1.0000
JC69G	56	-777.5170	1667.0340	0.0000	1.0000
JC69I	56	-778.2044	1668.4089	0.0000	1.0000
JC69IG	57	-777.4207	1668.8414	0.0000	1.0000
JC69	55	-784.9122	1679.8245	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-757.1327	1843.3531	0.6270	0.6270
HKYI	60	-758.0379	1845.1637	0.2536	0.8806
HKYIG	61	-757.0087	1848.5900	0.0457	0.9263
K2PG	57	-768.0550	1848.7435	0.0423	0.9687
K2PI	57	-768.7005	1850.0343	0.0222	0.9909
HKY	59	-764.8851	1853.3732	0.0042	0.9951
K2PIG	58	-767.9453	1854.0088	0.0030	0.9981
F81G	59	-766.6045	1856.8120	0.0007	0.9988
F81I	59	-767.3901	1858.3832	0.0003	0.9992
K2P	56	-775.6941	1858.5368	0.0003	0.9995
GTRG	64	-754.0721	1859.1711	0.0002	0.9997
GTRI	64	-755.1262	1861.2794	0.0001	0.9998
F81IG	60	-766.4948	1862.0775	0.0001	0.9999
JC69G	56	-777.5170	1862.1826	0.0001	0.9999
JC69I	56	-778.2044	1863.5575	0.0000	0.9999
GTRIG	65	-753.9419	1864.3955	0.0000	1.0000
SYMG	61	-765.0490	1864.6706	0.0000	1.0000
SYMI	61	-765.7021	1865.9769	0.0000	1.0000
F81	58	-774.0001	1866.1185	0.0000	1.0000
JC69IG	57	-777.4207	1867.4748	0.0000	1.0000
GTR	63	-761.8509	1869.2439	0.0000	1.0000
SYMIG	62	-764.9574	1869.9723	0.0000	1.0000
JC69	55	-784.9122	1871.4883	0.0000	1.0000
SYM	60	-773.0418	1875.1715	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/3390_hg18.phylip = 1 - 241;
 Partition Dummy = 1:261-aggressive-phylip/3390_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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/3390_hg18.phylip = 1 - 241;
 Partition Dummy = 1:261-aggressive-phylip/3390_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 BIC model (HKYG)]

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


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

