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

Input data from file "261-aggressive-phylip/3367_hg18.phylip" (29 taxa, 225 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	-736.8811	1638.3964	0.6539	0.6539
HKYI	60	-737.8354	1640.3049	0.2518	0.9057
HKYIG	61	-737.0723	1642.5495	0.0820	0.9877
GTRG	64	-733.5771	1647.1542	0.0082	0.9959
GTRI	64	-734.5241	1649.0483	0.0032	0.9991
GTRIG	65	-733.8371	1651.6365	0.0009	0.9999
K2PG	57	-751.8434	1657.2796	0.0001	1.0000
K2PI	57	-753.2632	1660.1192	0.0000	1.0000
K2PIG	58	-752.6028	1662.4346	0.0000	1.0000
SYMG	61	-749.4732	1667.3514	0.0000	1.0000
HKY	59	-753.6459	1668.2009	0.0000	1.0000
SYMI	61	-750.7176	1669.8400	0.0000	1.0000
SYMIG	62	-750.0521	1672.3264	0.0000	1.0000
F81I	59	-756.1241	1673.1572	0.0000	1.0000
F81G	59	-756.1911	1673.2912	0.0000	1.0000
GTR	63	-749.4944	1675.0757	0.0000	1.0000
F81IG	60	-755.5480	1675.7301	0.0000	1.0000
JC69G	56	-768.8310	1687.6620	0.0000	1.0000
JC69I	56	-768.8417	1687.6833	0.0000	1.0000
K2P	56	-768.9592	1687.9185	0.0000	1.0000
JC69IG	57	-768.2441	1690.0810	0.0000	1.0000
SYM	60	-766.2723	1697.1788	0.0000	1.0000
F81	58	-772.5240	1702.2768	0.0000	1.0000
JC69	55	-785.4567	1717.3632	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-736.8811	1593.7623	0.3964	0.3964
GTRG	64	-733.5771	1595.1542	0.1977	0.5941
HKYI	60	-737.8354	1595.6707	0.1527	0.7468
HKYIG	61	-737.0723	1596.1446	0.1205	0.8673
GTRI	64	-734.5241	1597.0483	0.0767	0.9439
GTRIG	65	-733.8371	1597.6742	0.0561	1.0000
K2PG	57	-751.8434	1617.6868	0.0000	1.0000
K2PI	57	-753.2632	1620.5264	0.0000	1.0000
SYMG	61	-749.4732	1620.9465	0.0000	1.0000
K2PIG	58	-752.6028	1621.2057	0.0000	1.0000
SYMI	61	-750.7176	1623.4351	0.0000	1.0000
SYMIG	62	-750.0521	1624.1042	0.0000	1.0000
GTR	63	-749.4944	1624.9888	0.0000	1.0000
HKY	59	-753.6459	1625.2918	0.0000	1.0000
F81I	59	-756.1241	1630.2481	0.0000	1.0000
F81G	59	-756.1911	1630.3821	0.0000	1.0000
F81IG	60	-755.5480	1631.0960	0.0000	1.0000
JC69G	56	-768.8310	1649.6620	0.0000	1.0000
JC69I	56	-768.8417	1649.6833	0.0000	1.0000
K2P	56	-768.9592	1649.9185	0.0000	1.0000
JC69IG	57	-768.2441	1650.4882	0.0000	1.0000
SYM	60	-766.2723	1652.5446	0.0000	1.0000
F81	58	-772.5240	1661.0479	0.0000	1.0000
JC69	55	-785.4567	1680.9135	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-736.8811	1798.7283	0.6933	0.6933
HKYI	60	-737.8354	1800.6368	0.2670	0.9603
HKYIG	61	-737.0723	1804.5267	0.0382	0.9985
K2PG	57	-751.8434	1812.4045	0.0007	0.9993
GTRG	64	-733.5771	1813.7846	0.0004	0.9996
K2PI	57	-753.2632	1815.2441	0.0002	0.9998
GTRI	64	-734.5241	1815.6787	0.0001	1.0000
K2PIG	58	-752.6028	1819.3395	0.0000	1.0000
GTRIG	65	-733.8371	1819.7207	0.0000	1.0000
HKY	59	-753.6459	1826.8417	0.0000	1.0000
SYMG	61	-749.4732	1829.3286	0.0000	1.0000
F81I	59	-756.1241	1831.7981	0.0000	1.0000
SYMI	61	-750.7176	1831.8172	0.0000	1.0000
F81G	59	-756.1911	1831.9320	0.0000	1.0000
SYMIG	62	-750.0521	1835.9024	0.0000	1.0000
F81IG	60	-755.5480	1836.0620	0.0000	1.0000
GTR	63	-749.4944	1840.2031	0.0000	1.0000
JC69G	56	-768.8310	1840.9637	0.0000	1.0000
JC69I	56	-768.8417	1840.9849	0.0000	1.0000
K2P	56	-768.9592	1841.2201	0.0000	1.0000
JC69IG	57	-768.2441	1845.2059	0.0000	1.0000
SYM	60	-766.2723	1857.5107	0.0000	1.0000
F81	58	-772.5240	1859.1817	0.0000	1.0000
JC69	55	-785.4567	1868.7990	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/3367_hg18.phylip = 1 - 225;
 Partition Dummy = 1:261-aggressive-phylip/3367_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/3367_hg18.phylip = 1 - 225;
 Partition Dummy = 1:261-aggressive-phylip/3367_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/3367_hg18.phylip = 1 - 225;
 Partition Dummy = 1:261-aggressive-phylip/3367_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

