
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:07 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/chr10_5715.phylip" (19 taxa, 341 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-908.2310	1907.3954	0.4062	0.4062
HKYG	40	-908.2649	1907.4632	0.3926	0.7988
HKYIG	41	-907.7726	1909.0636	0.1764	0.9751
GTRI	44	-906.6968	1914.7720	0.0102	0.9853
GTRG	44	-906.7126	1914.8037	0.0100	0.9953
GTRIG	45	-906.2342	1916.5023	0.0043	0.9996
K2PI	37	-919.8144	1922.9093	0.0002	0.9998
K2PG	37	-919.8729	1923.0264	0.0002	0.9999
K2PIG	38	-919.7016	1925.2177	0.0001	1.0000
SYMI	41	-917.5925	1928.7035	0.0000	1.0000
SYMG	41	-917.5993	1928.7169	0.0000	1.0000
SYMIG	42	-917.0882	1930.2972	0.0000	1.0000
HKY	39	-933.0918	1954.5491	0.0000	1.0000
F81I	39	-935.2794	1958.9242	0.0000	1.0000
F81G	39	-935.3016	1958.9685	0.0000	1.0000
F81IG	40	-935.2745	1961.4824	0.0000	1.0000
GTR	43	-931.4274	1961.5955	0.0000	1.0000
K2P	36	-944.8352	1970.4337	0.0000	1.0000
SYM	40	-941.7780	1974.4893	0.0000	1.0000
JC69I	36	-948.5535	1977.8702	0.0000	1.0000
JC69G	36	-948.5585	1977.8802	0.0000	1.0000
JC69IG	37	-948.4914	1980.2634	0.0000	1.0000
F81	38	-959.7396	2005.2937	0.0000	1.0000
JC69	35	-972.1252	2022.5127	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-908.2310	1896.4621	0.3614	0.3614
HKYG	40	-908.2649	1896.5299	0.3493	0.7107
HKYIG	41	-907.7726	1897.5452	0.2103	0.9210
GTRI	44	-906.6968	1901.3937	0.0307	0.9517
GTRG	44	-906.7126	1901.4253	0.0302	0.9819
GTRIG	45	-906.2342	1902.4684	0.0179	0.9998
K2PI	37	-919.8144	1913.6288	0.0001	0.9999
K2PG	37	-919.8729	1913.7459	0.0001	0.9999
K2PIG	38	-919.7016	1915.4032	0.0000	1.0000
SYMI	41	-917.5925	1917.1851	0.0000	1.0000
SYMG	41	-917.5993	1917.1985	0.0000	1.0000
SYMIG	42	-917.0882	1918.1764	0.0000	1.0000
HKY	39	-933.0918	1944.1836	0.0000	1.0000
F81I	39	-935.2794	1948.5587	0.0000	1.0000
F81G	39	-935.3016	1948.6031	0.0000	1.0000
GTR	43	-931.4274	1948.8548	0.0000	1.0000
F81IG	40	-935.2745	1950.5490	0.0000	1.0000
K2P	36	-944.8352	1961.6705	0.0000	1.0000
SYM	40	-941.7780	1963.5559	0.0000	1.0000
JC69I	36	-948.5535	1969.1071	0.0000	1.0000
JC69G	36	-948.5585	1969.1171	0.0000	1.0000
JC69IG	37	-948.4914	1970.9829	0.0000	1.0000
F81	38	-959.7396	1995.4791	0.0000	1.0000
JC69	35	-972.1252	2014.2504	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-908.2310	2049.7374	0.4608	0.4608
HKYG	40	-908.2649	2049.8052	0.4455	0.9063
HKYIG	41	-907.7726	2054.6524	0.0395	0.9458
K2PI	37	-919.8144	2055.4085	0.0270	0.9728
K2PG	37	-919.8729	2055.5256	0.0255	0.9983
K2PIG	38	-919.7016	2061.0147	0.0016	1.0000
GTRI	44	-906.6968	2069.9965	0.0000	1.0000
GTRG	44	-906.7126	2070.0281	0.0000	1.0000
SYMI	41	-917.5925	2074.2923	0.0000	1.0000
SYMG	41	-917.5993	2074.3057	0.0000	1.0000
GTRIG	45	-906.2342	2074.9031	0.0000	1.0000
SYMIG	42	-917.0882	2079.1154	0.0000	1.0000
HKY	39	-933.0918	2093.6271	0.0000	1.0000
F81I	39	-935.2794	2098.0021	0.0000	1.0000
F81G	39	-935.3016	2098.0465	0.0000	1.0000
K2P	36	-944.8352	2099.6183	0.0000	1.0000
F81IG	40	-935.2745	2103.8243	0.0000	1.0000
JC69I	36	-948.5535	2107.0548	0.0000	1.0000
JC69G	36	-948.5585	2107.0648	0.0000	1.0000
JC69IG	37	-948.4914	2112.7625	0.0000	1.0000
GTR	43	-931.4274	2113.6257	0.0000	1.0000
SYM	40	-941.7780	2116.8312	0.0000	1.0000
F81	38	-959.7396	2141.0907	0.0000	1.0000
JC69	35	-972.1252	2148.3662	0.0000	1.0000



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

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

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


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

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


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

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


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

