
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:39 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/chr24_967.phylip" (19 taxa, 266 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	-785.2194	1665.0166	0.5034	0.5034
HKYG	40	-785.8060	1666.1898	0.2800	0.7833
HKYIG	41	-785.2354	1667.8459	0.1223	0.9057
GTRI	44	-782.1271	1670.1727	0.0382	0.9439
GTRG	44	-782.5864	1671.0914	0.0241	0.9680
HKY	39	-789.8574	1671.5202	0.0195	0.9875
GTRIG	45	-782.1387	1673.0957	0.0089	0.9964
GTR	43	-785.9390	1674.9231	0.0036	0.9999
K2PI	37	-798.7797	1683.8928	0.0000	1.0000
K2PG	37	-799.2418	1684.8168	0.0000	1.0000
K2PIG	38	-798.7924	1686.6420	0.0000	1.0000
K2P	36	-802.9998	1689.6327	0.0000	1.0000
F81I	39	-799.2500	1690.3052	0.0000	1.0000
SYMI	41	-796.8642	1691.1033	0.0000	1.0000
F81G	39	-799.7595	1691.3243	0.0000	1.0000
SYMG	41	-797.3445	1692.0641	0.0000	1.0000
F81IG	40	-799.2640	1693.1059	0.0000	1.0000
SYMIG	42	-796.8772	1693.9518	0.0000	1.0000
SYM	40	-800.3164	1695.2105	0.0000	1.0000
F81	38	-803.1701	1695.3975	0.0000	1.0000
JC69I	36	-810.2830	1704.1991	0.0000	1.0000
JC69G	36	-810.7417	1705.1167	0.0000	1.0000
JC69IG	37	-810.2956	1706.9245	0.0000	1.0000
JC69	35	-814.1663	1709.2891	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-785.2194	1650.4388	0.3605	0.3605
HKYG	40	-785.8060	1651.6120	0.2005	0.5610
GTRI	44	-782.1271	1652.2541	0.1455	0.7065
HKYIG	41	-785.2354	1652.4709	0.1305	0.8370
GTRG	44	-782.5864	1653.1728	0.0919	0.9289
GTRIG	45	-782.1387	1654.2775	0.0529	0.9818
HKY	39	-789.8574	1657.7149	0.0095	0.9912
GTR	43	-785.9390	1657.8781	0.0087	1.0000
K2PI	37	-798.7797	1671.5594	0.0000	1.0000
K2PG	37	-799.2418	1672.4835	0.0000	1.0000
K2PIG	38	-798.7924	1673.5848	0.0000	1.0000
SYMI	41	-796.8642	1675.7283	0.0000	1.0000
F81I	39	-799.2500	1676.4999	0.0000	1.0000
SYMG	41	-797.3445	1676.6891	0.0000	1.0000
F81G	39	-799.7595	1677.5190	0.0000	1.0000
SYMIG	42	-796.8772	1677.7545	0.0000	1.0000
K2P	36	-802.9998	1677.9996	0.0000	1.0000
F81IG	40	-799.2640	1678.5281	0.0000	1.0000
SYM	40	-800.3164	1680.6328	0.0000	1.0000
F81	38	-803.1701	1682.3402	0.0000	1.0000
JC69I	36	-810.2830	1692.5659	0.0000	1.0000
JC69G	36	-810.7417	1693.4835	0.0000	1.0000
JC69IG	37	-810.2956	1694.5911	0.0000	1.0000
JC69	35	-814.1663	1698.3325	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-785.2194	1793.7787	0.5600	0.5600
HKYG	40	-785.8060	1794.9519	0.3115	0.8714
HKY	39	-789.8574	1797.4712	0.0884	0.9598
HKYIG	41	-785.2354	1799.3942	0.0338	0.9936
K2PI	37	-798.7797	1804.1488	0.0031	0.9967
K2PG	37	-799.2418	1805.0729	0.0020	0.9987
K2P	36	-802.9998	1807.0054	0.0008	0.9994
K2PIG	38	-798.7924	1809.7576	0.0002	0.9996
GTRI	44	-782.1271	1809.9279	0.0002	0.9998
GTRG	44	-782.5864	1810.8467	0.0001	0.9999
GTR	43	-785.9390	1811.9684	0.0001	1.0000
GTRIG	45	-782.1387	1815.5348	0.0000	1.0000
F81I	39	-799.2500	1816.2563	0.0000	1.0000
F81G	39	-799.7595	1817.2753	0.0000	1.0000
F81	38	-803.1701	1818.5131	0.0000	1.0000
JC69I	36	-810.2830	1821.5718	0.0000	1.0000
F81IG	40	-799.2640	1821.8679	0.0000	1.0000
JC69G	36	-810.7417	1822.4893	0.0000	1.0000
SYMI	41	-796.8642	1822.6517	0.0000	1.0000
SYMG	41	-797.3445	1823.6124	0.0000	1.0000
JC69	35	-814.1663	1823.7549	0.0000	1.0000
SYM	40	-800.3164	1823.9726	0.0000	1.0000
JC69IG	37	-810.2956	1827.1805	0.0000	1.0000
SYMIG	42	-796.8772	1828.2613	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/chr24_967.phylip = 1 - 266;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr24_967.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/chr24_967.phylip = 1 - 266;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr24_967.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/chr24_967.phylip = 1 - 266;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr24_967.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

