
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:22 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/chr17_2965.phylip" (19 taxa, 316 characters)

Minimum AIC  model: HKYG
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	40	-786.1823	1664.2918	0.4595	0.4595
HKYI	40	-786.8845	1665.6963	0.2277	0.6871
HKYIG	41	-785.7690	1666.1074	0.1854	0.8725
GTRG	44	-782.7169	1668.0464	0.0703	0.9428
GTRIG	45	-782.2028	1669.7388	0.0302	0.9729
GTRI	44	-783.6743	1669.9612	0.0270	0.9999
HKY	39	-796.4846	1682.2736	0.0001	1.0000
GTR	43	-793.5719	1687.0555	0.0000	1.0000
F81G	39	-808.8943	1707.0930	0.0000	1.0000
F81I	39	-809.4139	1708.1321	0.0000	1.0000
F81IG	40	-808.6010	1709.1294	0.0000	1.0000
K2PG	37	-815.5770	1715.2690	0.0000	1.0000
K2PI	37	-816.0236	1716.1622	0.0000	1.0000
K2PIG	38	-815.1337	1716.9678	0.0000	1.0000
SYMG	41	-814.5474	1723.6641	0.0000	1.0000
F81	38	-818.6900	1724.0803	0.0000	1.0000
SYMI	41	-815.0064	1724.5822	0.0000	1.0000
SYMIG	42	-814.0905	1725.4118	0.0000	1.0000
K2P	36	-826.6611	1734.8705	0.0000	1.0000
SYM	40	-825.7193	1743.3659	0.0000	1.0000
JC69G	36	-834.3697	1750.2878	0.0000	1.0000
JC69I	36	-834.8108	1751.1700	0.0000	1.0000
JC69IG	37	-834.0588	1752.2327	0.0000	1.0000
JC69	35	-844.6553	1768.3107	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-786.1823	1652.3646	0.3103	0.3103
GTRG	44	-782.7169	1653.4339	0.1818	0.4921
HKYIG	41	-785.7690	1653.5381	0.1726	0.6646
HKYI	40	-786.8845	1653.7690	0.1537	0.8183
GTRIG	45	-782.2028	1654.4055	0.1118	0.9302
GTRI	44	-783.6743	1655.3486	0.0698	1.0000
HKY	39	-796.4846	1670.9693	0.0000	1.0000
GTR	43	-793.5719	1673.1438	0.0000	1.0000
F81G	39	-808.8943	1695.7886	0.0000	1.0000
F81I	39	-809.4139	1696.8278	0.0000	1.0000
F81IG	40	-808.6010	1697.2021	0.0000	1.0000
K2PG	37	-815.5770	1705.1539	0.0000	1.0000
K2PI	37	-816.0236	1706.0471	0.0000	1.0000
K2PIG	38	-815.1337	1706.2675	0.0000	1.0000
SYMG	41	-814.5474	1711.0948	0.0000	1.0000
SYMI	41	-815.0064	1712.0128	0.0000	1.0000
SYMIG	42	-814.0905	1712.1811	0.0000	1.0000
F81	38	-818.6900	1713.3800	0.0000	1.0000
K2P	36	-826.6611	1725.3222	0.0000	1.0000
SYM	40	-825.7193	1731.4386	0.0000	1.0000
JC69G	36	-834.3697	1740.7394	0.0000	1.0000
JC69I	36	-834.8108	1741.6216	0.0000	1.0000
JC69IG	37	-834.0588	1742.1176	0.0000	1.0000
JC69	35	-844.6553	1759.3107	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-786.1823	1802.5942	0.6323	0.6323
HKYI	40	-786.8845	1803.9987	0.3133	0.9456
HKYIG	41	-785.7690	1807.5235	0.0538	0.9993
HKY	39	-796.4846	1817.4432	0.0004	0.9997
GTRG	44	-782.7169	1818.6865	0.0002	0.9999
GTRI	44	-783.6743	1820.6013	0.0001	1.0000
GTRIG	45	-782.2028	1823.4139	0.0000	1.0000
GTR	43	-793.5719	1834.6407	0.0000	1.0000
F81G	39	-808.8943	1842.2626	0.0000	1.0000
F81I	39	-809.4139	1843.3017	0.0000	1.0000
K2PG	37	-815.5770	1844.1164	0.0000	1.0000
K2PI	37	-816.0236	1845.0096	0.0000	1.0000
F81IG	40	-808.6010	1847.4318	0.0000	1.0000
K2PIG	38	-815.1337	1848.9857	0.0000	1.0000
F81	38	-818.6900	1856.0982	0.0000	1.0000
K2P	36	-826.6611	1860.5289	0.0000	1.0000
SYMG	41	-814.5474	1865.0802	0.0000	1.0000
SYMI	41	-815.0064	1865.9983	0.0000	1.0000
SYMIG	42	-814.0905	1869.9222	0.0000	1.0000
JC69G	36	-834.3697	1875.9461	0.0000	1.0000
JC69I	36	-834.8108	1876.8283	0.0000	1.0000
JC69IG	37	-834.0588	1881.0800	0.0000	1.0000
SYM	40	-825.7193	1881.6683	0.0000	1.0000
JC69	35	-844.6553	1890.7617	0.0000	1.0000



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

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

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

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

