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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-737.2907	1632.9060	0.4257	0.4257
HKYI	60	-737.7213	1633.7671	0.2767	0.7024
K2PG	57	-743.8541	1635.7907	0.1006	0.8030
K2PI	57	-744.0445	1636.1715	0.0832	0.8862
HKYIG	61	-737.2636	1636.3377	0.0765	0.9627
K2PIG	58	-743.6428	1638.7468	0.0229	0.9857
GTRG	64	-734.3392	1641.1703	0.0068	0.9925
GTRI	64	-734.5634	1641.6188	0.0055	0.9980
GTRIG	65	-734.1170	1644.3630	0.0014	0.9993
SYMI	61	-742.7545	1647.3194	0.0003	0.9997
SYMG	61	-742.9076	1647.6257	0.0003	0.9999
SYMIG	62	-742.4844	1650.3021	0.0001	1.0000
HKY	59	-751.6052	1658.0853	0.0000	1.0000
K2P	56	-758.8264	1662.3912	0.0000	1.0000
GTR	63	-747.8276	1664.5488	0.0000	1.0000
F81G	59	-755.2685	1665.4120	0.0000	1.0000
F81I	59	-755.6979	1666.2708	0.0000	1.0000
F81IG	60	-755.3885	1669.1016	0.0000	1.0000
JC69G	56	-763.8651	1672.4686	0.0000	1.0000
JC69I	56	-764.0481	1672.8346	0.0000	1.0000
SYM	60	-757.4348	1673.1941	0.0000	1.0000
JC69IG	57	-763.7070	1675.4965	0.0000	1.0000
F81	58	-769.5567	1690.5745	0.0000	1.0000
JC69	55	-778.4572	1698.3429	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-737.2907	1594.5814	0.3463	0.3463
HKYI	60	-737.7213	1595.4425	0.2251	0.5714
HKYIG	61	-737.2636	1596.5271	0.1309	0.7023
GTRG	64	-734.3392	1596.6783	0.1214	0.8236
GTRI	64	-734.5634	1597.1268	0.0970	0.9206
GTRIG	65	-734.1170	1598.2340	0.0558	0.9764
K2PG	57	-743.8541	1601.7082	0.0098	0.9862
K2PI	57	-744.0445	1602.0890	0.0081	0.9943
K2PIG	58	-743.6428	1603.2857	0.0045	0.9987
SYMI	61	-742.7545	1607.5089	0.0005	0.9993
SYMG	61	-742.9076	1607.8152	0.0005	0.9997
SYMIG	62	-742.4844	1608.9687	0.0003	1.0000
HKY	59	-751.6052	1621.2103	0.0000	1.0000
GTR	63	-747.8276	1621.6552	0.0000	1.0000
F81G	59	-755.2685	1628.5370	0.0000	1.0000
F81I	59	-755.6979	1629.3958	0.0000	1.0000
K2P	56	-758.8264	1629.6527	0.0000	1.0000
F81IG	60	-755.3885	1630.7770	0.0000	1.0000
SYM	60	-757.4348	1634.8695	0.0000	1.0000
JC69G	56	-763.8651	1639.7301	0.0000	1.0000
JC69I	56	-764.0481	1640.0961	0.0000	1.0000
JC69IG	57	-763.7070	1641.4140	0.0000	1.0000
F81	58	-769.5567	1655.1133	0.0000	1.0000
JC69	55	-778.4572	1666.9144	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-743.8541	1802.8856	0.4528	0.4528
K2PI	57	-744.0445	1803.2664	0.3743	0.8271
HKYG	60	-737.2907	1806.3471	0.0802	0.9073
HKYI	60	-737.7213	1807.2082	0.0522	0.9595
K2PIG	58	-743.6428	1807.9926	0.0352	0.9947
HKYIG	61	-737.2636	1811.8223	0.0052	0.9999
GTRG	64	-734.3392	1822.5618	0.0000	0.9999
SYMI	61	-742.7545	1822.8041	0.0000	1.0000
GTRI	64	-734.5634	1823.0103	0.0000	1.0000
SYMG	61	-742.9076	1823.1104	0.0000	1.0000
K2P	56	-758.8264	1827.3007	0.0000	1.0000
GTRIG	65	-734.1170	1827.6469	0.0000	1.0000
SYMIG	62	-742.4844	1827.7933	0.0000	1.0000
HKY	59	-751.6052	1829.4466	0.0000	1.0000
F81G	59	-755.2685	1836.7733	0.0000	1.0000
JC69G	56	-763.8651	1837.3782	0.0000	1.0000
F81I	59	-755.6979	1837.6321	0.0000	1.0000
JC69I	56	-764.0481	1837.7441	0.0000	1.0000
F81IG	60	-755.3885	1842.5427	0.0000	1.0000
JC69IG	57	-763.7070	1842.5915	0.0000	1.0000
GTR	63	-747.8276	1844.0093	0.0000	1.0000
SYM	60	-757.4348	1846.6353	0.0000	1.0000
F81	58	-769.5567	1859.8202	0.0000	1.0000
JC69	55	-778.4572	1861.0330	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/5391_hg18.phylip = 1 - 252;
 Partition Dummy = 1:261-aggressive-phylip/5391_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/5391_hg18.phylip = 1 - 252;
 Partition Dummy = 1:261-aggressive-phylip/5391_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 (K2PG)]

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

