
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:19 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/chr15_4329.phylip" (19 taxa, 240 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	-751.9052	1591.7313	0.4728	0.4728
K2PI	37	-752.3618	1592.6443	0.2995	0.7724
K2PIG	38	-751.8045	1594.3552	0.1273	0.8997
HKYG	40	-749.9735	1596.4293	0.0451	0.9448
HKYI	40	-750.4321	1597.3467	0.0285	0.9734
HKYIG	41	-749.8746	1599.1431	0.0116	0.9850
SYMG	41	-750.4078	1600.2095	0.0068	0.9918
SYMI	41	-750.8410	1601.0759	0.0044	0.9962
SYMIG	42	-750.3244	1602.9838	0.0017	0.9979
GTRG	44	-748.5359	1605.3795	0.0005	0.9984
JC69G	36	-760.1768	1605.4768	0.0005	0.9989
GTRI	44	-749.0041	1606.3159	0.0003	0.9992
JC69I	36	-760.6331	1606.3893	0.0003	0.9996
JC69IG	37	-760.0876	1608.0961	0.0001	0.9997
GTRIG	45	-748.4322	1608.2046	0.0001	0.9998
F81G	39	-757.5581	1608.7163	0.0001	0.9999
F81I	39	-758.0245	1609.6490	0.0001	1.0000
F81IG	40	-757.4613	1611.4051	0.0000	1.0000
K2P	36	-769.1859	1623.4949	0.0000	1.0000
HKY	39	-767.0325	1627.6650	0.0000	1.0000
SYM	40	-767.0676	1630.6176	0.0000	1.0000
GTR	43	-765.0174	1635.3410	0.0000	1.0000
JC69	35	-777.1519	1636.6568	0.0000	1.0000
F81	38	-774.6878	1640.1219	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-751.9052	1577.8105	0.3373	0.3373
K2PI	37	-752.3618	1578.7235	0.2137	0.5510
K2PIG	38	-751.8045	1579.6090	0.1372	0.6882
HKYG	40	-749.9735	1579.9469	0.1159	0.8041
HKYI	40	-750.4321	1580.8643	0.0733	0.8774
HKYIG	41	-749.8746	1581.7492	0.0471	0.9244
SYMG	41	-750.4078	1582.8155	0.0276	0.9521
SYMI	41	-750.8410	1583.6819	0.0179	0.9700
SYMIG	42	-750.3244	1584.6488	0.0110	0.9810
GTRG	44	-748.5359	1585.0718	0.0089	0.9900
GTRI	44	-749.0041	1586.0082	0.0056	0.9955
GTRIG	45	-748.4322	1586.8644	0.0036	0.9992
JC69G	36	-760.1768	1592.3536	0.0002	0.9994
F81G	39	-757.5581	1593.1163	0.0002	0.9996
JC69I	36	-760.6331	1593.2661	0.0001	0.9997
F81I	39	-758.0245	1594.0490	0.0001	0.9998
JC69IG	37	-760.0876	1594.1753	0.0001	0.9999
F81IG	40	-757.4613	1594.9227	0.0001	1.0000
K2P	36	-769.1859	1610.3718	0.0000	1.0000
HKY	39	-767.0325	1612.0650	0.0000	1.0000
SYM	40	-767.0676	1614.1352	0.0000	1.0000
GTR	43	-765.0174	1616.0349	0.0000	1.0000
JC69	35	-777.1519	1624.3039	0.0000	1.0000
F81	38	-774.6878	1625.3757	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-751.9052	1706.5941	0.5831	0.5831
K2PI	37	-752.3618	1707.5071	0.3694	0.9525
K2PIG	38	-751.8045	1711.8733	0.0416	0.9941
JC69G	36	-760.1768	1717.6566	0.0023	0.9964
JC69I	36	-760.6331	1718.5691	0.0015	0.9979
HKYG	40	-749.9735	1719.1725	0.0011	0.9990
HKYI	40	-750.4321	1720.0899	0.0007	0.9997
JC69IG	37	-760.0876	1722.9589	0.0002	0.9998
HKYIG	41	-749.8746	1724.4554	0.0001	0.9999
SYMG	41	-750.4078	1725.5217	0.0000	1.0000
SYMI	41	-750.8410	1726.3881	0.0000	1.0000
F81G	39	-757.5581	1728.8612	0.0000	1.0000
F81I	39	-758.0245	1729.7939	0.0000	1.0000
SYMIG	42	-750.3244	1730.8356	0.0000	1.0000
F81IG	40	-757.4613	1734.1482	0.0000	1.0000
K2P	36	-769.1859	1735.6748	0.0000	1.0000
GTRG	44	-748.5359	1738.2199	0.0000	1.0000
GTRI	44	-749.0041	1739.1563	0.0000	1.0000
GTRIG	45	-748.4322	1743.4931	0.0000	1.0000
JC69	35	-777.1519	1746.1262	0.0000	1.0000
HKY	39	-767.0325	1747.8099	0.0000	1.0000
SYM	40	-767.0676	1753.3607	0.0000	1.0000
F81	38	-774.6878	1757.6399	0.0000	1.0000
GTR	43	-765.0174	1765.7023	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/chr15_4329.phylip = 1 - 240;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_4329.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/chr15_4329.phylip = 1 - 240;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_4329.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/chr15_4329.phylip = 1 - 240;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_4329.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

