
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:54 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/chr7_6598.phylip" (19 taxa, 311 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	40	-922.0762	1936.3006	0.4348	0.4348
HKYI	40	-922.1917	1936.5315	0.3874	0.8223
HKYIG	41	-922.4203	1939.6437	0.0817	0.9040
GTRI	44	-919.1420	1941.1712	0.0381	0.9421
GTRG	44	-919.3590	1941.6053	0.0307	0.9728
K2PG	37	-929.8454	1943.9911	0.0093	0.9820
K2PI	37	-929.9611	1944.2225	0.0083	0.9903
GTRIG	45	-919.3750	1944.3727	0.0077	0.9980
K2PIG	38	-930.5346	1947.9662	0.0013	0.9993
SYMG	41	-928.0895	1950.9821	0.0003	0.9996
SYMI	41	-928.2319	1951.2667	0.0002	0.9998
F81G	39	-932.1846	1953.8822	0.0001	0.9999
F81I	39	-932.3738	1954.2605	0.0001	0.9999
SYMIG	42	-928.7260	1954.9296	0.0000	1.0000
F81IG	40	-932.0036	1956.1554	0.0000	1.0000
JC69G	36	-939.3146	1960.3517	0.0000	1.0000
JC69I	36	-939.4744	1960.6714	0.0000	1.0000
JC69IG	37	-940.1282	1964.5569	0.0000	1.0000
HKY	39	-942.3881	1974.2890	0.0000	1.0000
GTR	43	-939.7173	1979.6069	0.0000	1.0000
K2P	36	-951.0379	1983.7984	0.0000	1.0000
SYM	40	-949.5670	1991.2821	0.0000	1.0000
F81	38	-952.2929	1991.4829	0.0000	1.0000
JC69	35	-959.9082	1998.9801	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-922.0762	1924.1525	0.3456	0.3456
HKYI	40	-922.1917	1924.3833	0.3079	0.6535
GTRI	44	-919.1420	1926.2839	0.1190	0.7726
GTRG	44	-919.3590	1926.7181	0.0958	0.8684
HKYIG	41	-922.4203	1926.8407	0.0901	0.9585
GTRIG	45	-919.3750	1928.7500	0.0347	0.9932
K2PG	37	-929.8454	1933.6907	0.0029	0.9961
K2PI	37	-929.9611	1933.9221	0.0026	0.9987
K2PIG	38	-930.5346	1937.0692	0.0005	0.9993
SYMG	41	-928.0895	1938.1791	0.0003	0.9996
SYMI	41	-928.2319	1938.4637	0.0003	0.9999
SYMIG	42	-928.7260	1941.4520	0.0001	0.9999
F81G	39	-932.1846	1942.3693	0.0000	1.0000
F81I	39	-932.3738	1942.7476	0.0000	1.0000
F81IG	40	-932.0036	1944.0072	0.0000	1.0000
JC69G	36	-939.3146	1950.6291	0.0000	1.0000
JC69I	36	-939.4744	1950.9488	0.0000	1.0000
JC69IG	37	-940.1282	1954.2565	0.0000	1.0000
HKY	39	-942.3881	1962.7761	0.0000	1.0000
GTR	43	-939.7173	1965.4347	0.0000	1.0000
K2P	36	-951.0379	1974.0758	0.0000	1.0000
SYM	40	-949.5670	1979.1340	0.0000	1.0000
F81	38	-952.2929	1980.5858	0.0000	1.0000
JC69	35	-959.9082	1989.8164	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-929.8454	2072.0631	0.3629	0.3629
K2PI	37	-929.9611	2072.2945	0.3232	0.6861
HKYG	40	-922.0762	2073.7442	0.1566	0.8426
HKYI	40	-922.1917	2073.9750	0.1395	0.9821
K2PIG	38	-930.5346	2079.1813	0.0103	0.9925
HKYIG	41	-922.4203	2080.1722	0.0063	0.9988
JC69G	36	-939.3146	2085.2617	0.0005	0.9993
JC69I	36	-939.4744	2085.5813	0.0004	0.9997
F81G	39	-932.1846	2088.2212	0.0001	0.9998
F81I	39	-932.3738	2088.5995	0.0001	0.9999
GTRI	44	-919.1420	2090.8348	0.0000	0.9999
GTRG	44	-919.3590	2091.2689	0.0000	0.9999
SYMG	41	-928.0895	2091.5106	0.0000	1.0000
SYMI	41	-928.2319	2091.7952	0.0000	1.0000
JC69IG	37	-940.1282	2092.6288	0.0000	1.0000
F81IG	40	-932.0036	2093.5989	0.0000	1.0000
GTRIG	45	-919.3750	2097.0407	0.0000	1.0000
SYMIG	42	-928.7260	2098.5233	0.0000	1.0000
HKY	39	-942.3881	2108.6280	0.0000	1.0000
K2P	36	-951.0379	2108.7083	0.0000	1.0000
JC69	35	-959.9082	2120.7092	0.0000	1.0000
F81	38	-952.2929	2122.6980	0.0000	1.0000
GTR	43	-939.7173	2126.2458	0.0000	1.0000
SYM	40	-949.5670	2128.7257	0.0000	1.0000



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

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

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


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6598.phylip = 1 - 311;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6598.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) 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/chr7_6598.phylip = 1 - 311;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6598.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

