
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:15 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_4763.phylip" (19 taxa, 274 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	-910.9517	1907.8187	0.7955	0.7955
K2PI	37	-913.0821	1912.0794	0.0945	0.8900
K2PIG	38	-912.4700	1913.5528	0.0452	0.9352
HKYG	40	-910.3399	1914.7571	0.0248	0.9600
HKYIG	41	-909.2914	1915.4276	0.0177	0.9777
SYMG	41	-909.3730	1915.5909	0.0163	0.9940
HKYI	40	-912.5178	1919.1128	0.0028	0.9968
SYMI	41	-911.3568	1919.5583	0.0022	0.9990
SYMIG	42	-911.3207	1922.2779	0.0006	0.9996
GTRG	44	-909.0808	1923.4542	0.0003	0.9999
GTRI	44	-911.1029	1927.4983	0.0000	1.0000
GTRIG	45	-910.9922	1930.1423	0.0000	1.0000
K2P	36	-938.9847	1961.2099	0.0000	1.0000
HKY	39	-937.9245	1967.1824	0.0000	1.0000
SYM	40	-939.1213	1972.3198	0.0000	1.0000
GTR	43	-938.3849	1979.2219	0.0000	1.0000
JC69G	36	-963.8285	2010.8974	0.0000	1.0000
JC69IG	37	-962.9909	2011.8970	0.0000	1.0000
JC69I	36	-965.4390	2014.1186	0.0000	1.0000
F81G	39	-963.0082	2017.3498	0.0000	1.0000
F81IG	40	-962.1163	2018.3098	0.0000	1.0000
F81I	39	-964.6718	2020.6769	0.0000	1.0000
JC69	35	-989.2223	2059.0328	0.0000	1.0000
F81	38	-988.2160	2065.0447	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-910.9517	1895.9034	0.6620	0.6620
K2PI	37	-913.0821	1900.1641	0.0786	0.7406
HKYIG	41	-909.2914	1900.5828	0.0638	0.8044
HKYG	40	-910.3399	1900.6798	0.0608	0.8652
SYMG	41	-909.3730	1900.7461	0.0588	0.9239
K2PIG	38	-912.4700	1900.9401	0.0534	0.9773
SYMI	41	-911.3568	1904.7135	0.0081	0.9854
HKYI	40	-912.5178	1905.0356	0.0069	0.9923
GTRG	44	-909.0808	1906.1616	0.0039	0.9962
SYMIG	42	-911.3207	1906.6415	0.0031	0.9993
GTRI	44	-911.1029	1910.2057	0.0005	0.9998
GTRIG	45	-910.9922	1911.9844	0.0002	1.0000
K2P	36	-938.9847	1949.9694	0.0000	1.0000
HKY	39	-937.9245	1953.8490	0.0000	1.0000
SYM	40	-939.1213	1958.2425	0.0000	1.0000
GTR	43	-938.3849	1962.7697	0.0000	1.0000
JC69G	36	-963.8285	1999.6569	0.0000	1.0000
JC69IG	37	-962.9909	1999.9817	0.0000	1.0000
JC69I	36	-965.4390	2002.8781	0.0000	1.0000
F81G	39	-963.0082	2004.0164	0.0000	1.0000
F81IG	40	-962.1163	2004.2325	0.0000	1.0000
F81I	39	-964.6718	2007.3435	0.0000	1.0000
JC69	35	-989.2223	2048.4445	0.0000	1.0000
F81	38	-988.2160	2052.4319	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-910.9517	2029.5891	0.8829	0.8829
K2PI	37	-913.0821	2033.8499	0.1049	0.9878
K2PIG	38	-912.4700	2038.2389	0.0117	0.9995
HKYG	40	-910.3399	2045.2050	0.0004	0.9998
HKYIG	41	-909.2914	2048.7210	0.0001	0.9999
SYMG	41	-909.3730	2048.8843	0.0001	1.0000
HKYI	40	-912.5178	2049.5607	0.0000	1.0000
SYMI	41	-911.3568	2052.8518	0.0000	1.0000
SYMIG	42	-911.3207	2058.3929	0.0000	1.0000
GTRG	44	-909.0808	2065.1392	0.0000	1.0000
GTRI	44	-911.1029	2069.1834	0.0000	1.0000
GTRIG	45	-910.9922	2074.5752	0.0000	1.0000
K2P	36	-938.9847	2080.0420	0.0000	1.0000
HKY	39	-937.9245	2094.7610	0.0000	1.0000
SYM	40	-939.1213	2102.7676	0.0000	1.0000
GTR	43	-938.3849	2118.1342	0.0000	1.0000
JC69G	36	-963.8285	2129.7295	0.0000	1.0000
JC69I	36	-965.4390	2132.9507	0.0000	1.0000
JC69IG	37	-962.9909	2133.6675	0.0000	1.0000
F81G	39	-963.0082	2144.9284	0.0000	1.0000
F81I	39	-964.6718	2148.2555	0.0000	1.0000
F81IG	40	-962.1163	2148.7576	0.0000	1.0000
JC69	35	-989.2223	2174.9040	0.0000	1.0000
F81	38	-988.2160	2189.7308	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/chrZ_4763.phylip = 1 - 274;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_4763.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/chrZ_4763.phylip = 1 - 274;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_4763.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/chrZ_4763.phylip = 1 - 274;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_4763.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

