
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:23 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_11214.phylip" (19 taxa, 267 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	-791.6021	1669.4838	0.4185	0.4185
K2PI	37	-791.6926	1669.6647	0.3823	0.8008
K2PIG	38	-791.5676	1672.1352	0.1111	0.9119
HKYG	40	-789.8485	1674.2103	0.0394	0.9513
SYMG	41	-789.2386	1675.7839	0.0179	0.9692
SYMI	41	-789.3313	1675.9692	0.0163	0.9856
HKYI	40	-791.7961	1678.1054	0.0056	0.9912
SYMIG	42	-789.2470	1678.6190	0.0043	0.9955
HKYIG	41	-791.4942	1680.2950	0.0019	0.9974
GTRG	44	-787.4011	1680.6400	0.0016	0.9990
GTRIG	45	-787.3987	1683.5305	0.0004	0.9994
K2P	36	-800.0169	1683.6165	0.0004	0.9997
GTRI	44	-789.2545	1684.3468	0.0002	1.0000
HKY	39	-798.3651	1688.4746	0.0000	1.0000
SYM	40	-798.4910	1691.4954	0.0000	1.0000
GTR	43	-796.8086	1696.5858	0.0000	1.0000
JC69I	36	-814.7423	1713.0673	0.0000	1.0000
JC69G	36	-814.8416	1713.2659	0.0000	1.0000
JC69IG	37	-814.8576	1715.9946	0.0000	1.0000
F81I	39	-812.7877	1717.3199	0.0000	1.0000
F81G	39	-812.8646	1717.4738	0.0000	1.0000
F81IG	40	-812.7701	1720.0535	0.0000	1.0000
JC69	35	-822.5517	1726.0125	0.0000	1.0000
F81	38	-820.5822	1730.1643	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-791.6021	1657.2043	0.3145	0.3145
K2PI	37	-791.6926	1657.3852	0.2873	0.6017
K2PIG	38	-791.5676	1659.1352	0.1197	0.7215
HKYG	40	-789.8485	1659.6971	0.0904	0.8119
SYMG	41	-789.2386	1660.4773	0.0612	0.8731
SYMI	41	-789.3313	1660.6626	0.0558	0.9289
SYMIG	42	-789.2470	1662.4940	0.0223	0.9512
GTRG	44	-787.4011	1662.8022	0.0191	0.9704
HKYI	40	-791.7961	1663.5921	0.0129	0.9833
GTRIG	45	-787.3987	1664.7974	0.0071	0.9903
HKYIG	41	-791.4942	1664.9884	0.0064	0.9967
GTRI	44	-789.2545	1666.5089	0.0030	0.9997
K2P	36	-800.0169	1672.0339	0.0002	0.9999
HKY	39	-798.3651	1674.7301	0.0000	1.0000
SYM	40	-798.4910	1676.9821	0.0000	1.0000
GTR	43	-796.8086	1679.6172	0.0000	1.0000
JC69I	36	-814.7423	1701.4847	0.0000	1.0000
JC69G	36	-814.8416	1701.6833	0.0000	1.0000
F81I	39	-812.7877	1703.5754	0.0000	1.0000
JC69IG	37	-814.8576	1703.7152	0.0000	1.0000
F81G	39	-812.8646	1703.7293	0.0000	1.0000
F81IG	40	-812.7701	1705.5402	0.0000	1.0000
JC69	35	-822.5517	1715.1034	0.0000	1.0000
F81	38	-820.5822	1717.1643	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-791.6021	1789.9325	0.5045	0.5045
K2PI	37	-791.6926	1790.1134	0.4608	0.9653
K2PIG	38	-791.5676	1795.4507	0.0320	0.9973
K2P	36	-800.0169	1801.1749	0.0018	0.9991
HKYG	40	-789.8485	1803.1870	0.0007	0.9997
HKYI	40	-791.7961	1807.0820	0.0001	0.9998
SYMG	41	-789.2386	1807.5545	0.0001	0.9999
SYMI	41	-789.3313	1807.7398	0.0001	1.0000
HKYIG	41	-791.4942	1812.0656	0.0000	1.0000
SYMIG	42	-789.2470	1813.1585	0.0000	1.0000
HKY	39	-798.3651	1814.6328	0.0000	1.0000
SYM	40	-798.4910	1820.4720	0.0000	1.0000
GTRG	44	-787.4011	1820.6411	0.0000	1.0000
GTRI	44	-789.2545	1824.3479	0.0000	1.0000
GTRIG	45	-787.3987	1826.2236	0.0000	1.0000
JC69I	36	-814.7423	1830.6256	0.0000	1.0000
JC69G	36	-814.8416	1830.8243	0.0000	1.0000
GTR	43	-796.8086	1833.8689	0.0000	1.0000
JC69IG	37	-814.8576	1836.4434	0.0000	1.0000
JC69	35	-822.5517	1840.6571	0.0000	1.0000
F81I	39	-812.7877	1843.4781	0.0000	1.0000
F81G	39	-812.8646	1843.6320	0.0000	1.0000
F81IG	40	-812.7701	1849.0301	0.0000	1.0000
F81	38	-820.5822	1853.4798	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_11214.phylip = 1 - 267;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11214.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_11214.phylip = 1 - 267;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11214.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_11214.phylip = 1 - 267;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11214.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

