
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:07 02/24/2011

Input data from file "/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14414.phylip" (19 taxa, 229 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-781.6856	1660.8181	0.5199	0.5199
HKYG	40	-782.1069	1661.6605	0.3412	0.8611
HKYIG	41	-781.6959	1663.8090	0.1165	0.9776
GTRI	44	-779.8330	1669.1877	0.0079	0.9856
GTRG	44	-780.1299	1669.7815	0.0059	0.9914
K2PI	37	-791.2539	1671.2302	0.0029	0.9943
HKY	39	-788.8376	1672.1832	0.0018	0.9961
GTRIG	45	-779.8405	1672.3039	0.0017	0.9977
K2PG	37	-791.8699	1672.4624	0.0015	0.9993
K2PIG	38	-791.2699	1674.1398	0.0007	0.9999
GTR	43	-786.4407	1679.3354	0.0000	1.0000
SYMI	41	-791.5440	1683.5051	0.0000	1.0000
F81I	39	-794.8385	1684.1850	0.0000	1.0000
SYMG	41	-792.2482	1684.9136	0.0000	1.0000
F81G	39	-795.3468	1685.2016	0.0000	1.0000
K2P	36	-799.6694	1685.2139	0.0000	1.0000
SYMIG	42	-791.5617	1686.5428	0.0000	1.0000
F81IG	40	-794.7909	1687.0285	0.0000	1.0000
JC69I	36	-803.3606	1692.5963	0.0000	1.0000
JC69G	36	-804.0256	1693.9263	0.0000	1.0000
JC69IG	37	-803.3774	1695.4774	0.0000	1.0000
F81	38	-802.6833	1696.9665	0.0000	1.0000
SYM	40	-799.8519	1697.1507	0.0000	1.0000
JC69	35	-811.9826	1707.0223	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-781.6856	1643.3713	0.4402	0.4402
HKYG	40	-782.1069	1644.2137	0.2889	0.7290
HKYIG	41	-781.6959	1645.3918	0.1603	0.8893
GTRI	44	-779.8330	1647.6660	0.0514	0.9407
GTRG	44	-780.1299	1648.2598	0.0382	0.9789
GTRIG	45	-779.8405	1649.6810	0.0188	0.9977
HKY	39	-788.8376	1655.6753	0.0009	0.9986
K2PI	37	-791.2539	1656.5077	0.0006	0.9992
K2PG	37	-791.8699	1657.7399	0.0003	0.9996
K2PIG	38	-791.2699	1658.5398	0.0002	0.9998
GTR	43	-786.4407	1658.8813	0.0002	1.0000
SYMI	41	-791.5440	1665.0880	0.0000	1.0000
SYMG	41	-792.2482	1666.4965	0.0000	1.0000
SYMIG	42	-791.5617	1667.1234	0.0000	1.0000
F81I	39	-794.8385	1667.6770	0.0000	1.0000
F81G	39	-795.3468	1668.6936	0.0000	1.0000
F81IG	40	-794.7909	1669.5817	0.0000	1.0000
K2P	36	-799.6694	1671.3389	0.0000	1.0000
JC69I	36	-803.3606	1678.7213	0.0000	1.0000
SYM	40	-799.8519	1679.7039	0.0000	1.0000
JC69G	36	-804.0256	1680.0513	0.0000	1.0000
JC69IG	37	-803.3774	1680.7548	0.0000	1.0000
F81	38	-802.6833	1681.3665	0.0000	1.0000
JC69	35	-811.9826	1693.9653	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-781.6856	1780.7201	0.4709	0.4709
HKYG	40	-782.1069	1781.5626	0.3090	0.7800
K2PI	37	-791.2539	1783.5554	0.1141	0.8941
K2PG	37	-791.8699	1784.7876	0.0616	0.9557
HKYIG	41	-781.6959	1786.1744	0.0308	0.9865
K2PIG	38	-791.2699	1789.0213	0.0074	0.9939
HKY	39	-788.8376	1789.5904	0.0056	0.9995
K2P	36	-799.6694	1794.9529	0.0004	0.9999
GTRI	44	-779.8330	1798.7497	0.0001	0.9999
GTRG	44	-780.1299	1799.3436	0.0000	1.0000
F81I	39	-794.8385	1801.5922	0.0000	1.0000
JC69I	36	-803.3606	1802.3353	0.0000	1.0000
F81G	39	-795.3468	1802.6088	0.0000	1.0000
JC69G	36	-804.0256	1803.6653	0.0000	1.0000
GTRIG	45	-779.8405	1804.1985	0.0000	1.0000
SYMI	41	-791.5440	1805.8706	0.0000	1.0000
GTR	43	-786.4407	1806.5314	0.0000	1.0000
F81IG	40	-794.7909	1806.9306	0.0000	1.0000
SYMG	41	-792.2482	1807.2791	0.0000	1.0000
JC69IG	37	-803.3774	1807.8026	0.0000	1.0000
SYMIG	42	-791.5617	1811.3397	0.0000	1.0000
F81	38	-802.6833	1811.8479	0.0000	1.0000
JC69	35	-811.9826	1814.1455	0.0000	1.0000
SYM	40	-799.8519	1817.0527	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14414.phylip = 1 - 229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14414.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14414.phylip = 1 - 229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14414.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14414.phylip = 1 - 229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14414.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

