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

Input data from file "261-aggressive-phylip/8466_hg18.phylip" (29 taxa, 283 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	-821.9366	1787.2599	0.5424	0.5424
HKYG	60	-817.8945	1788.7620	0.2560	0.7984
K2PIG	58	-821.9367	1790.4270	0.1113	0.9098
HKYIG	61	-817.8945	1792.0153	0.0503	0.9601
K2PI	57	-825.4970	1794.3806	0.0154	0.9755
SYMG	61	-819.6070	1795.4403	0.0091	0.9846
HKYI	60	-821.5601	1796.0931	0.0066	0.9911
K2P	56	-828.7317	1797.7112	0.0029	0.9940
GTRG	64	-816.0826	1798.3303	0.0021	0.9962
SYMIG	62	-819.6071	1798.7233	0.0018	0.9979
HKY	59	-824.8473	1799.4436	0.0012	0.9992
GTRIG	65	-816.0825	1801.7042	0.0004	0.9996
SYMI	61	-823.0974	1802.4211	0.0003	0.9998
SYM	60	-825.8509	1804.6747	0.0001	0.9999
GTRI	64	-819.6910	1805.5472	0.0001	1.0000
GTR	63	-822.6309	1808.0838	0.0000	1.0000
JC69G	56	-841.6534	1823.5546	0.0000	1.0000
F81G	59	-837.7131	1825.1751	0.0000	1.0000
JC69IG	57	-841.6535	1826.6936	0.0000	1.0000
F81IG	60	-837.7130	1828.3990	0.0000	1.0000
JC69I	56	-845.1875	1830.6228	0.0000	1.0000
F81I	59	-841.3300	1832.4088	0.0000	1.0000
JC69	55	-848.3104	1833.7574	0.0000	1.0000
F81	58	-844.6753	1835.9042	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-817.8945	1755.7890	0.4674	0.4674
HKYIG	61	-817.8945	1757.7891	0.1719	0.6393
K2PG	57	-821.9366	1757.8732	0.1648	0.8041
K2PIG	58	-821.9367	1759.8734	0.0606	0.8648
GTRG	64	-816.0826	1760.1652	0.0524	0.9172
SYMG	61	-819.6070	1761.2141	0.0310	0.9482
GTRIG	65	-816.0825	1762.1650	0.0193	0.9675
HKYI	60	-821.5601	1763.1201	0.0120	0.9795
SYMIG	62	-819.6071	1763.2142	0.0114	0.9909
K2PI	57	-825.4970	1764.9940	0.0047	0.9956
GTRI	64	-819.6910	1767.3820	0.0014	0.9970
HKY	59	-824.8473	1767.6947	0.0012	0.9982
SYMI	61	-823.0974	1768.1948	0.0009	0.9991
K2P	56	-828.7317	1769.4634	0.0005	0.9996
GTR	63	-822.6309	1771.2619	0.0002	0.9998
SYM	60	-825.8509	1771.7017	0.0002	1.0000
F81G	59	-837.7131	1793.4262	0.0000	1.0000
JC69G	56	-841.6534	1795.3068	0.0000	1.0000
F81IG	60	-837.7130	1795.4260	0.0000	1.0000
JC69IG	57	-841.6535	1797.3070	0.0000	1.0000
F81I	59	-841.3300	1800.6599	0.0000	1.0000
JC69I	56	-845.1875	1802.3751	0.0000	1.0000
F81	58	-844.6753	1805.3506	0.0000	1.0000
JC69	55	-848.3104	1806.6208	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-821.9366	1965.6637	0.8928	0.8928
K2PIG	58	-821.9367	1971.3093	0.0531	0.9459
K2PI	57	-825.4970	1972.7844	0.0254	0.9713
K2P	56	-828.7317	1973.6085	0.0168	0.9881
HKYG	60	-817.8945	1974.5158	0.0107	0.9988
HKYIG	61	-817.8945	1980.1613	0.0006	0.9994
HKYI	60	-821.5601	1981.8469	0.0003	0.9997
HKY	59	-824.8473	1982.7760	0.0002	0.9999
SYMG	61	-819.6070	1983.5864	0.0001	1.0000
SYMIG	62	-819.6071	1989.2319	0.0000	1.0000
SYM	60	-825.8509	1990.4285	0.0000	1.0000
SYMI	61	-823.0974	1990.5671	0.0000	1.0000
GTRG	64	-816.0826	1993.4738	0.0000	1.0000
GTRIG	65	-816.0825	1999.1190	0.0000	1.0000
JC69G	56	-841.6534	1999.4519	0.0000	1.0000
GTRI	64	-819.6910	2000.6906	0.0000	1.0000
GTR	63	-822.6309	2000.9250	0.0000	1.0000
JC69IG	57	-841.6535	2005.0975	0.0000	1.0000
JC69I	56	-845.1875	2006.5201	0.0000	1.0000
JC69	55	-848.3104	2007.1204	0.0000	1.0000
F81G	59	-837.7131	2008.5076	0.0000	1.0000
F81IG	60	-837.7130	2014.1528	0.0000	1.0000
F81I	59	-841.3300	2015.7413	0.0000	1.0000
F81	58	-844.6753	2016.7865	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/8466_hg18.phylip = 1 - 283;
 Partition Dummy = 1:261-aggressive-phylip/8466_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/8466_hg18.phylip = 1 - 283;
 Partition Dummy = 1:261-aggressive-phylip/8466_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/8466_hg18.phylip = 1 - 283;
 Partition Dummy = 1:261-aggressive-phylip/8466_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

