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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-804.8837	1767.4817	0.4355	0.4355
GTRI	64	-805.0705	1767.8553	0.3613	0.7968
GTRIG	65	-804.8747	1770.5020	0.0962	0.8930
GTR	63	-808.1398	1770.9771	0.0759	0.9689
HKYG	60	-814.3015	1774.3777	0.0139	0.9827
HKYI	60	-814.3642	1774.5030	0.0130	0.9957
HKYIG	61	-814.2990	1777.3260	0.0032	0.9989
HKY	59	-818.3639	1779.5699	0.0010	1.0000
K2PG	57	-826.1829	1789.4041	0.0000	1.0000
SYMG	61	-820.3520	1789.4319	0.0000	1.0000
K2PI	57	-826.2337	1789.5057	0.0000	1.0000
SYMI	61	-820.4826	1789.6931	0.0000	1.0000
F81G	59	-823.6646	1790.1712	0.0000	1.0000
F81I	59	-823.7146	1790.2713	0.0000	1.0000
K2PIG	58	-826.1775	1792.2850	0.0000	1.0000
SYMIG	62	-820.3382	1792.3785	0.0000	1.0000
F81IG	60	-823.6597	1793.0940	0.0000	1.0000
SYM	60	-823.9821	1793.7388	0.0000	1.0000
F81	58	-827.8300	1795.5900	0.0000	1.0000
K2P	56	-830.8090	1795.7847	0.0000	1.0000
JC69G	56	-835.0457	1804.2580	0.0000	1.0000
JC69I	56	-835.0982	1804.3630	0.0000	1.0000
JC69IG	57	-835.0377	1807.1138	0.0000	1.0000
JC69	55	-839.5015	1810.3180	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-804.8837	1737.7675	0.4318	0.4318
GTRI	64	-805.0705	1738.1410	0.3582	0.7900
GTRIG	65	-804.8747	1739.7494	0.1603	0.9503
GTR	63	-808.1398	1742.2796	0.0452	0.9955
HKYG	60	-814.3015	1748.6030	0.0019	0.9974
HKYI	60	-814.3642	1748.7283	0.0018	0.9992
HKYIG	61	-814.2990	1750.5981	0.0007	0.9999
HKY	59	-818.3639	1754.7278	0.0001	1.0000
SYMG	61	-820.3520	1762.7040	0.0000	1.0000
SYMI	61	-820.4826	1762.9652	0.0000	1.0000
SYMIG	62	-820.3382	1764.6764	0.0000	1.0000
F81G	59	-823.6646	1765.3291	0.0000	1.0000
F81I	59	-823.7146	1765.4292	0.0000	1.0000
K2PG	57	-826.1829	1766.3658	0.0000	1.0000
K2PI	57	-826.2337	1766.4674	0.0000	1.0000
F81IG	60	-823.6597	1767.3193	0.0000	1.0000
SYM	60	-823.9821	1767.9642	0.0000	1.0000
K2PIG	58	-826.1775	1768.3550	0.0000	1.0000
F81	58	-827.8300	1771.6599	0.0000	1.0000
K2P	56	-830.8090	1773.6180	0.0000	1.0000
JC69G	56	-835.0457	1782.0913	0.0000	1.0000
JC69I	56	-835.0982	1782.1963	0.0000	1.0000
JC69IG	57	-835.0377	1784.0755	0.0000	1.0000
JC69	55	-839.5015	1789.0031	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-814.3015	1979.2157	0.3727	0.3727
HKYI	60	-814.3642	1979.3410	0.3501	0.7229
HKY	59	-818.3639	1981.4969	0.1191	0.8420
GTRG	64	-804.8837	1983.7543	0.0385	0.8805
GTRI	64	-805.0705	1984.1278	0.0320	0.9125
GTR	63	-808.1398	1984.4229	0.0276	0.9401
HKYIG	61	-814.2990	1985.0543	0.0201	0.9602
K2PG	57	-826.1829	1985.4478	0.0165	0.9767
K2PI	57	-826.2337	1985.5494	0.0157	0.9924
K2P	56	-830.8090	1988.8565	0.0030	0.9954
GTRIG	65	-804.8747	1989.5797	0.0021	0.9975
K2PIG	58	-826.1775	1991.2806	0.0009	0.9984
F81G	59	-823.6646	1992.0982	0.0006	0.9990
F81I	59	-823.7146	1992.1984	0.0006	0.9996
F81	58	-827.8300	1994.5855	0.0002	0.9998
SYMG	61	-820.3520	1997.1602	0.0000	0.9998
JC69G	56	-835.0457	1997.3298	0.0000	0.9998
SYMI	61	-820.4826	1997.4214	0.0000	0.9999
JC69I	56	-835.0982	1997.4348	0.0000	0.9999
F81IG	60	-823.6597	1997.9320	0.0000	1.0000
SYM	60	-823.9821	1998.5768	0.0000	1.0000
JC69	55	-839.5015	2000.3980	0.0000	1.0000
SYMIG	62	-820.3382	2002.9761	0.0000	1.0000
JC69IG	57	-835.0377	2003.1575	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8865_hg18.phylip = 1 - 345;
 Partition Dummy = 1:261-aggressive-phylip/8865_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8865_hg18.phylip = 1 - 345;
 Partition Dummy = 1:261-aggressive-phylip/8865_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,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/8865_hg18.phylip = 1 - 345;
 Partition Dummy = 1:261-aggressive-phylip/8865_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

