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

Minimum AIC  model: K2PI
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	37	-891.7469	1868.8325	0.4121	0.4121
K2PG	37	-891.7985	1868.9357	0.3914	0.8035
K2PIG	38	-891.7981	1871.5962	0.1035	0.9070
HKYG	40	-890.2812	1873.9501	0.0319	0.9389
HKYI	40	-890.3437	1874.0751	0.0300	0.9689
SYMG	41	-889.7719	1875.6586	0.0136	0.9825
HKYIG	41	-890.3236	1876.7619	0.0078	0.9903
SYMI	41	-890.5701	1877.2549	0.0061	0.9964
SYMIG	42	-890.6104	1880.0850	0.0015	0.9979
GTRG	44	-887.8537	1880.1388	0.0014	0.9993
GTRI	44	-888.9126	1882.2567	0.0005	0.9998
GTRIG	45	-888.4877	1884.2253	0.0002	1.0000
JC69I	36	-911.4042	1905.5072	0.0000	1.0000
JC69G	36	-911.4645	1905.6277	0.0000	1.0000
JC69IG	37	-911.4641	1908.2669	0.0000	1.0000
F81I	39	-910.0232	1910.7294	0.0000	1.0000
F81G	39	-910.0417	1910.7664	0.0000	1.0000
F81IG	40	-910.0411	1913.4700	0.0000	1.0000
K2P	36	-917.2920	1917.2827	0.0000	1.0000
SYM	40	-912.3252	1918.0382	0.0000	1.0000
HKY	39	-915.5610	1921.8049	0.0000	1.0000
GTR	43	-910.4035	1922.4433	0.0000	1.0000
JC69	35	-936.7744	1953.6289	0.0000	1.0000
F81	38	-935.3876	1958.7751	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	37	-891.7469	1857.4938	0.3252	0.3252
K2PG	37	-891.7985	1857.5970	0.3089	0.6341
K2PIG	38	-891.7981	1859.5962	0.1137	0.7478
HKYG	40	-890.2812	1860.5624	0.0701	0.8179
HKYI	40	-890.3437	1860.6873	0.0659	0.8838
SYMG	41	-889.7719	1861.5438	0.0429	0.9267
HKYIG	41	-890.3236	1862.6472	0.0247	0.9514
SYMI	41	-890.5701	1863.1402	0.0193	0.9707
GTRG	44	-887.8537	1863.7073	0.0146	0.9853
SYMIG	42	-890.6104	1865.2208	0.0068	0.9921
GTRI	44	-888.9126	1865.8251	0.0050	0.9972
GTRIG	45	-888.4877	1866.9753	0.0028	1.0000
JC69I	36	-911.4042	1894.8084	0.0000	1.0000
JC69G	36	-911.4645	1894.9289	0.0000	1.0000
JC69IG	37	-911.4641	1896.9282	0.0000	1.0000
F81I	39	-910.0232	1898.0465	0.0000	1.0000
F81G	39	-910.0417	1898.0834	0.0000	1.0000
F81IG	40	-910.0411	1900.0823	0.0000	1.0000
SYM	40	-912.3252	1904.6505	0.0000	1.0000
K2P	36	-917.2920	1906.5839	0.0000	1.0000
GTR	43	-910.4035	1906.8070	0.0000	1.0000
HKY	39	-915.5610	1909.1219	0.0000	1.0000
JC69	35	-936.7744	1943.5489	0.0000	1.0000
F81	38	-935.3876	1946.7751	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-891.7469	1992.7655	0.4981	0.4981
K2PG	37	-891.7985	1992.8687	0.4730	0.9711
K2PIG	38	-891.7981	1998.5239	0.0280	0.9990
HKYG	40	-890.2812	2006.8020	0.0004	0.9995
HKYI	40	-890.3437	2006.9270	0.0004	0.9999
SYMG	41	-889.7719	2011.4395	0.0000	1.0000
HKYIG	41	-890.3236	2012.5428	0.0000	1.0000
SYMI	41	-890.5701	2013.0358	0.0000	1.0000
SYMIG	42	-890.6104	2018.7724	0.0000	1.0000
GTRG	44	-887.8537	2024.5709	0.0000	1.0000
JC69I	36	-911.4042	2026.4241	0.0000	1.0000
JC69G	36	-911.4645	2026.5446	0.0000	1.0000
GTRI	44	-888.9126	2026.6888	0.0000	1.0000
GTRIG	45	-888.4877	2031.4949	0.0000	1.0000
JC69IG	37	-911.4641	2032.1999	0.0000	1.0000
K2P	36	-917.2920	2038.1996	0.0000	1.0000
F81I	39	-910.0232	2040.6302	0.0000	1.0000
F81G	39	-910.0417	2040.6671	0.0000	1.0000
F81IG	40	-910.0411	2046.3220	0.0000	1.0000
SYM	40	-912.3252	2050.8902	0.0000	1.0000
HKY	39	-915.5610	2051.7056	0.0000	1.0000
GTR	43	-910.4035	2064.0146	0.0000	1.0000
JC69	35	-936.7744	2071.5086	0.0000	1.0000
F81	38	-935.3876	2085.7028	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_14552.phylip = 1 - 286;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_14552.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_14552.phylip = 1 - 286;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_14552.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_14552.phylip = 1 - 286;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_14552.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

