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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	40	-814.1233	1721.3666	0.6614	0.6614
HKYIG	41	-813.9510	1723.7333	0.2025	0.8639
HKYI	40	-815.9265	1724.9730	0.1090	0.9729
GTRIG	45	-811.3480	1729.5939	0.0108	0.9837
K2PG	37	-822.5559	1730.2264	0.0079	0.9916
GTRI	44	-814.0710	1732.2395	0.0029	0.9944
GTRG	44	-814.2083	1732.5142	0.0025	0.9969
K2PI	37	-824.1976	1733.5098	0.0015	0.9985
K2PIG	38	-822.9023	1733.5665	0.0015	1.0000
SYMI	41	-823.4308	1742.6928	0.0000	1.0000
SYMIG	42	-822.1616	1742.8878	0.0000	1.0000
SYMG	41	-823.9783	1743.7878	0.0000	1.0000
F81IG	40	-826.2051	1745.5301	0.0000	1.0000
F81I	39	-827.7223	1745.8749	0.0000	1.0000
F81G	39	-828.4120	1747.2542	0.0000	1.0000
JC69IG	37	-834.5717	1754.2580	0.0000	1.0000
JC69I	36	-835.9506	1754.3893	0.0000	1.0000
JC69G	36	-836.4166	1755.3213	0.0000	1.0000
HKY	39	-850.2982	1791.0266	0.0000	1.0000
GTR	43	-847.2598	1795.8394	0.0000	1.0000
K2P	36	-859.6931	1801.8744	0.0000	1.0000
SYM	40	-858.5545	1810.2290	0.0000	1.0000
F81	38	-861.9010	1811.5640	0.0000	1.0000
JC69	35	-870.4790	1820.8403	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-814.1233	1708.2466	0.5706	0.5706
HKYIG	41	-813.9510	1709.9020	0.2494	0.8200
HKYI	40	-815.9265	1711.8530	0.0940	0.9140
GTRIG	45	-811.3480	1712.6960	0.0617	0.9757
GTRI	44	-814.0710	1716.1419	0.0110	0.9867
GTRG	44	-814.2083	1716.4166	0.0096	0.9963
K2PG	37	-822.5559	1719.1118	0.0025	0.9988
K2PIG	38	-822.9023	1721.8046	0.0006	0.9995
K2PI	37	-824.1976	1722.3952	0.0005	0.9999
SYMIG	42	-822.1616	1728.3233	0.0000	1.0000
SYMI	41	-823.4308	1728.8615	0.0000	1.0000
SYMG	41	-823.9783	1729.9565	0.0000	1.0000
F81IG	40	-826.2051	1732.4101	0.0000	1.0000
F81I	39	-827.7223	1733.4446	0.0000	1.0000
F81G	39	-828.4120	1734.8239	0.0000	1.0000
JC69IG	37	-834.5717	1743.1434	0.0000	1.0000
JC69I	36	-835.9506	1743.9011	0.0000	1.0000
JC69G	36	-836.4166	1744.8331	0.0000	1.0000
HKY	39	-850.2982	1778.5963	0.0000	1.0000
GTR	43	-847.2598	1780.5196	0.0000	1.0000
K2P	36	-859.6931	1791.3862	0.0000	1.0000
SYM	40	-858.5545	1797.1090	0.0000	1.0000
F81	38	-861.9010	1799.8021	0.0000	1.0000
JC69	35	-870.4790	1810.9580	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-822.5559	1855.0248	0.4206	0.4206
HKYG	40	-814.1233	1855.1795	0.3893	0.8098
K2PI	37	-824.1976	1858.3082	0.0814	0.8913
HKYI	40	-815.9265	1858.7859	0.0641	0.9554
HKYIG	41	-813.9510	1860.5083	0.0271	0.9825
K2PIG	38	-822.9023	1861.3909	0.0174	0.9999
JC69I	36	-835.9506	1876.1408	0.0000	1.0000
F81I	39	-827.7223	1876.7042	0.0000	1.0000
JC69G	36	-836.4166	1877.0728	0.0000	1.0000
GTRI	44	-814.0710	1877.7681	0.0000	1.0000
GTRIG	45	-811.3480	1877.9955	0.0000	1.0000
GTRG	44	-814.2083	1878.0429	0.0000	1.0000
F81G	39	-828.4120	1878.0835	0.0000	1.0000
JC69IG	37	-834.5717	1879.0564	0.0000	1.0000
F81IG	40	-826.2051	1879.3431	0.0000	1.0000
SYMI	41	-823.4308	1879.4678	0.0000	1.0000
SYMG	41	-823.9783	1880.5628	0.0000	1.0000
SYMIG	42	-822.1616	1882.6029	0.0000	1.0000
HKY	39	-850.2982	1921.8559	0.0000	1.0000
K2P	36	-859.6931	1923.6259	0.0000	1.0000
GTR	43	-847.2598	1938.4725	0.0000	1.0000
F81	38	-861.9010	1939.3883	0.0000	1.0000
JC69	35	-870.4790	1939.5243	0.0000	1.0000
SYM	40	-858.5545	1944.0420	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/chr24_1122.phylip = 1 - 291;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr24_1122.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/chr24_1122.phylip = 1 - 291;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr24_1122.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 (K2PG)]

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

