
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:18 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/chrZ_11457.phylip" (19 taxa, 728 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	44	-2157.8188	4409.4355	0.5417	0.5417
GTRI	44	-2158.4690	4410.7360	0.2827	0.8244
GTRIG	45	-2157.8098	4411.6900	0.1755	0.9999
HKYG	40	-2171.1963	4427.1669	0.0001	0.9999
HKYI	40	-2171.9510	4428.6764	0.0000	1.0000
HKYIG	41	-2171.1920	4429.4043	0.0000	1.0000
GTR	43	-2173.5089	4438.5500	0.0000	1.0000
HKY	39	-2185.4615	4453.4578	0.0000	1.0000
SYMG	41	-2223.2752	4533.5709	0.0000	1.0000
SYMIG	42	-2223.2359	4535.7449	0.0000	1.0000
SYMI	41	-2224.4503	4535.9209	0.0000	1.0000
K2PG	37	-2229.3148	4536.7049	0.0000	1.0000
K2PIG	38	-2229.2468	4538.7954	0.0000	1.0000
K2PI	37	-2230.7371	4539.5495	0.0000	1.0000
F81G	39	-2231.4252	4545.3854	0.0000	1.0000
F81IG	40	-2231.4324	4547.6391	0.0000	1.0000
F81I	39	-2236.6172	4555.7693	0.0000	1.0000
SYM	40	-2239.8560	4564.4864	0.0000	1.0000
K2P	36	-2245.6715	4567.1983	0.0000	1.0000
F81	38	-2244.3651	4569.0321	0.0000	1.0000
JC69G	36	-2280.5058	4636.8669	0.0000	1.0000
JC69I	36	-2281.2547	4638.3647	0.0000	1.0000
JC69IG	37	-2280.5079	4639.0912	0.0000	1.0000
JC69	35	-2294.2124	4662.0663	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2157.8188	4403.6376	0.5282	0.5282
GTRI	44	-2158.4690	4404.9380	0.2757	0.8039
GTRIG	45	-2157.8098	4405.6196	0.1961	0.9999
HKYG	40	-2171.1963	4422.3926	0.0000	1.0000
HKYI	40	-2171.9510	4423.9020	0.0000	1.0000
HKYIG	41	-2171.1920	4424.3839	0.0000	1.0000
GTR	43	-2173.5089	4433.0178	0.0000	1.0000
HKY	39	-2185.4615	4448.9230	0.0000	1.0000
SYMG	41	-2223.2752	4528.5505	0.0000	1.0000
SYMIG	42	-2223.2359	4530.4719	0.0000	1.0000
SYMI	41	-2224.4503	4530.9005	0.0000	1.0000
K2PG	37	-2229.3148	4532.6296	0.0000	1.0000
K2PIG	38	-2229.2468	4534.4936	0.0000	1.0000
K2PI	37	-2230.7371	4535.4742	0.0000	1.0000
F81G	39	-2231.4252	4540.8505	0.0000	1.0000
F81IG	40	-2231.4324	4542.8647	0.0000	1.0000
F81I	39	-2236.6172	4551.2344	0.0000	1.0000
SYM	40	-2239.8560	4559.7120	0.0000	1.0000
K2P	36	-2245.6715	4563.3430	0.0000	1.0000
F81	38	-2244.3651	4564.7302	0.0000	1.0000
JC69G	36	-2280.5058	4633.0116	0.0000	1.0000
JC69I	36	-2281.2547	4634.5094	0.0000	1.0000
JC69IG	37	-2280.5079	4635.0159	0.0000	1.0000
JC69	35	-2294.2124	4658.4247	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2157.8188	4605.6108	0.3575	0.3575
HKYG	40	-2171.1963	4606.0046	0.2936	0.6511
GTRI	44	-2158.4690	4606.9112	0.1866	0.8377
HKYI	40	-2171.9510	4607.5141	0.1380	0.9757
GTRIG	45	-2157.8098	4612.1831	0.0134	0.9891
HKYIG	41	-2171.1920	4612.5863	0.0109	1.0000
HKY	39	-2185.4615	4627.9447	0.0000	1.0000
GTR	43	-2173.5089	4630.4008	0.0000	1.0000
K2PG	37	-2229.3148	4702.4707	0.0000	1.0000
K2PI	37	-2230.7371	4705.3153	0.0000	1.0000
K2PIG	38	-2229.2468	4708.9250	0.0000	1.0000
SYMG	41	-2223.2752	4716.7528	0.0000	1.0000
SYMI	41	-2224.4503	4719.1029	0.0000	1.0000
F81G	39	-2231.4252	4719.8722	0.0000	1.0000
SYMIG	42	-2223.2359	4723.2645	0.0000	1.0000
F81IG	40	-2231.4324	4726.4768	0.0000	1.0000
K2P	36	-2245.6715	4728.5938	0.0000	1.0000
F81I	39	-2236.6172	4730.2562	0.0000	1.0000
F81	38	-2244.3651	4739.1616	0.0000	1.0000
SYM	40	-2239.8560	4743.3241	0.0000	1.0000
JC69G	36	-2280.5058	4798.2624	0.0000	1.0000
JC69I	36	-2281.2547	4799.7602	0.0000	1.0000
JC69IG	37	-2280.5079	4804.8570	0.0000	1.0000
JC69	35	-2294.2124	4819.0853	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_11457.phylip = 1 - 728;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_11457.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_11457.phylip = 1 - 728;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_11457.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_11457.phylip = 1 - 728;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_11457.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

