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

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

Minimum AIC  model: GTRG
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	-826.9687	1813.2922	0.5430	0.5430
HKYI	60	-827.8470	1815.0488	0.2256	0.7687
HKYIG	61	-826.9715	1816.8294	0.0926	0.8613
GTRG	64	-821.9431	1817.6004	0.0630	0.9243
GTRI	64	-823.0048	1819.7240	0.0218	0.9461
K2PG	57	-835.3785	1819.7412	0.0216	0.9677
GTRIG	65	-821.9448	1821.2930	0.0099	0.9776
K2PI	57	-836.2751	1821.5343	0.0088	0.9864
HKY	59	-832.9693	1821.7996	0.0077	0.9941
K2PIG	58	-835.3681	1823.1404	0.0039	0.9981
GTR	63	-828.0137	1826.0930	0.0009	0.9990
SYMG	61	-832.0946	1827.0757	0.0006	0.9996
SYMI	61	-832.9955	1828.8775	0.0002	0.9998
K2P	56	-842.1974	1829.9948	0.0001	0.9999
SYMIG	62	-832.0835	1830.6236	0.0001	1.0000
SYM	60	-838.8163	1836.9874	0.0000	1.0000
F81G	59	-846.3519	1848.5649	0.0000	1.0000
F81I	59	-847.2619	1850.3847	0.0000	1.0000
F81IG	60	-846.3564	1852.0676	0.0000	1.0000
JC69G	56	-853.5799	1852.7597	0.0000	1.0000
JC69I	56	-854.5496	1854.6992	0.0000	1.0000
JC69IG	57	-853.5719	1856.1279	0.0000	1.0000
F81	58	-852.3597	1857.1237	0.0000	1.0000
JC69	55	-860.0094	1862.2701	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-821.9431	1771.8861	0.4228	0.4228
GTRIG	65	-821.9448	1773.8897	0.1553	0.5781
HKYG	60	-826.9687	1773.9373	0.1516	0.7297
GTRI	64	-823.0048	1774.0097	0.1462	0.8759
HKYI	60	-827.8470	1775.6939	0.0630	0.9389
HKYIG	61	-826.9715	1775.9429	0.0556	0.9945
GTR	63	-828.0137	1782.0274	0.0027	0.9972
HKY	59	-832.9693	1783.9386	0.0010	0.9982
K2PG	57	-835.3785	1784.7570	0.0007	0.9989
SYMG	61	-832.0946	1786.1892	0.0003	0.9992
K2PI	57	-836.2751	1786.5502	0.0003	0.9995
K2PIG	58	-835.3681	1786.7361	0.0003	0.9997
SYMI	61	-832.9955	1787.9910	0.0001	0.9999
SYMIG	62	-832.0835	1788.1671	0.0001	1.0000
K2P	56	-842.1974	1796.3948	0.0000	1.0000
SYM	60	-838.8163	1797.6326	0.0000	1.0000
F81G	59	-846.3519	1810.7039	0.0000	1.0000
F81I	59	-847.2619	1812.5237	0.0000	1.0000
F81IG	60	-846.3564	1812.7128	0.0000	1.0000
JC69G	56	-853.5799	1819.1597	0.0000	1.0000
F81	58	-852.3597	1820.7194	0.0000	1.0000
JC69I	56	-854.5496	1821.0992	0.0000	1.0000
JC69IG	57	-853.5719	1821.1438	0.0000	1.0000
JC69	55	-860.0094	1830.0188	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-826.9687	1984.5006	0.3559	0.3559
K2PG	57	-835.3785	1984.7922	0.3077	0.6636
HKYI	60	-827.8470	1986.2572	0.1479	0.8115
K2PI	57	-836.2751	1986.5853	0.1255	0.9370
HKYIG	61	-826.9715	1990.0156	0.0226	0.9596
K2PIG	58	-835.3681	1990.2806	0.0198	0.9794
HKY	59	-832.9693	1990.9925	0.0139	0.9932
K2P	56	-842.1974	1992.9205	0.0053	0.9985
GTRG	64	-821.9431	1996.4870	0.0009	0.9994
GTRI	64	-823.0048	1998.6105	0.0003	0.9997
SYMG	61	-832.0946	2000.2619	0.0001	0.9998
GTRIG	65	-821.9448	2001.9999	0.0001	0.9999
SYMI	61	-832.9955	2002.0637	0.0001	1.0000
GTR	63	-828.0137	2003.1189	0.0000	1.0000
SYMIG	62	-832.0835	2005.7492	0.0000	1.0000
SYM	60	-838.8163	2008.1959	0.0000	1.0000
JC69G	56	-853.5799	2015.6855	0.0000	1.0000
JC69I	56	-854.5496	2017.6249	0.0000	1.0000
F81G	59	-846.3519	2017.7578	0.0000	1.0000
F81I	59	-847.2619	2019.5777	0.0000	1.0000
JC69IG	57	-853.5719	2021.1789	0.0000	1.0000
JC69	55	-860.0094	2023.0352	0.0000	1.0000
F81IG	60	-846.3564	2023.2761	0.0000	1.0000
F81	58	-852.3597	2024.2639	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/12095_hg18.phylip = 1 - 247;
 Partition Dummy = 1:261-aggressive-phylip/12095_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 (HKYG)]

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

