
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:40 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/chr5_14682.phylip" (19 taxa, 289 characters)

Minimum AIC  model: K2PG
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	37	-789.1044	1663.4119	0.3660	0.3660
K2PIG	38	-788.2546	1664.3651	0.2272	0.5932
K2PI	37	-789.7267	1664.6566	0.1964	0.7897
HKYG	40	-786.7110	1666.6479	0.0726	0.8623
HKYIG	41	-785.7702	1667.4838	0.0478	0.9100
HKYI	40	-787.3656	1667.9571	0.0377	0.9478
SYMG	41	-786.5909	1669.1251	0.0210	0.9688
SYMIG	42	-785.6839	1670.0507	0.0132	0.9820
SYMI	41	-787.4408	1670.8250	0.0090	0.9910
GTRG	44	-784.0424	1672.3143	0.0043	0.9953
GTRIG	45	-782.9989	1673.0349	0.0030	0.9983
GTRI	44	-784.9427	1674.1149	0.0017	1.0000
K2P	36	-803.3803	1689.3319	0.0000	1.0000
HKY	39	-801.1860	1692.9022	0.0000	1.0000
SYM	40	-800.3789	1693.9836	0.0000	1.0000
JC69G	36	-806.9005	1696.3724	0.0000	1.0000
JC69IG	37	-806.1045	1697.4122	0.0000	1.0000
JC69I	36	-807.5455	1697.6623	0.0000	1.0000
GTR	43	-798.1646	1697.7741	0.0000	1.0000
F81G	39	-805.0727	1700.6756	0.0000	1.0000
F81IG	40	-804.2405	1701.7068	0.0000	1.0000
F81I	39	-805.7213	1701.9727	0.0000	1.0000
JC69	35	-820.8990	1721.7584	0.0000	1.0000
F81	38	-819.1926	1726.2411	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-789.1044	1652.2087	0.2164	0.2164
K2PIG	38	-788.2546	1652.5091	0.1862	0.4026
HKYG	40	-786.7110	1653.4221	0.1180	0.5206
K2PI	37	-789.7267	1653.4534	0.1161	0.6367
HKYIG	41	-785.7702	1653.5405	0.1112	0.7479
HKYI	40	-787.3656	1654.7313	0.0613	0.8092
SYMG	41	-786.5909	1655.1818	0.0489	0.8581
SYMIG	42	-785.6839	1655.3678	0.0446	0.9027
GTRIG	45	-782.9989	1655.9979	0.0325	0.9353
GTRG	44	-784.0424	1656.0848	0.0312	0.9664
SYMI	41	-787.4408	1656.8817	0.0209	0.9873
GTRI	44	-784.9427	1657.8854	0.0127	1.0000
K2P	36	-803.3803	1678.7605	0.0000	1.0000
HKY	39	-801.1860	1680.3720	0.0000	1.0000
SYM	40	-800.3789	1680.7578	0.0000	1.0000
GTR	43	-798.1646	1682.3292	0.0000	1.0000
JC69G	36	-806.9005	1685.8010	0.0000	1.0000
JC69IG	37	-806.1045	1686.2090	0.0000	1.0000
JC69I	36	-807.5455	1687.0909	0.0000	1.0000
F81G	39	-805.0727	1688.1455	0.0000	1.0000
F81IG	40	-804.2405	1688.4810	0.0000	1.0000
F81I	39	-805.7213	1689.4426	0.0000	1.0000
JC69	35	-820.8990	1711.7979	0.0000	1.0000
F81	38	-819.1926	1714.3851	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-789.1044	1787.8665	0.5959	0.5959
K2PI	37	-789.7267	1789.1112	0.3198	0.9156
K2PIG	38	-788.2546	1791.8334	0.0820	0.9976
HKYG	40	-786.7110	1800.0791	0.0013	0.9990
HKYI	40	-787.3656	1801.3884	0.0007	0.9997
HKYIG	41	-785.7702	1803.8640	0.0002	0.9999
SYMG	41	-786.5909	1805.5053	0.0001	0.9999
SYMI	41	-787.4408	1807.2052	0.0000	1.0000
SYMIG	42	-785.6839	1809.3577	0.0000	1.0000
K2P	36	-803.3803	1810.7519	0.0000	1.0000
GTRG	44	-784.0424	1817.4075	0.0000	1.0000
JC69G	36	-806.9005	1817.7923	0.0000	1.0000
JC69I	36	-807.5455	1819.0823	0.0000	1.0000
GTRI	44	-784.9427	1819.2082	0.0000	1.0000
GTRIG	45	-782.9989	1820.9871	0.0000	1.0000
JC69IG	37	-806.1045	1821.8668	0.0000	1.0000
HKY	39	-801.1860	1823.3627	0.0000	1.0000
SYM	40	-800.3789	1827.4149	0.0000	1.0000
F81G	39	-805.0727	1831.1361	0.0000	1.0000
F81I	39	-805.7213	1832.4332	0.0000	1.0000
F81IG	40	-804.2405	1835.1381	0.0000	1.0000
GTR	43	-798.1646	1839.9856	0.0000	1.0000
JC69	35	-820.8990	1840.1229	0.0000	1.0000
F81	38	-819.1926	1853.7093	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14682.phylip = 1 - 289;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14682.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 AICc model (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14682.phylip = 1 - 289;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14682.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14682.phylip = 1 - 289;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14682.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

