
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:35 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/chr1_15480.phylip" (19 taxa, 319 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	-949.8979	1991.5944	0.4622	0.4622
HKYG	40	-950.1163	1992.0311	0.3715	0.8337
HKYIG	41	-950.1160	1994.6652	0.0995	0.9332
GTRI	44	-947.1010	1996.6546	0.0368	0.9700
GTRG	44	-947.5359	1997.5244	0.0238	0.9938
GTRIG	45	-947.5356	2000.2361	0.0061	1.0000
HKY	39	-961.4569	2012.0965	0.0000	1.0000
K2PI	37	-965.2727	2014.5526	0.0000	1.0000
K2PG	37	-965.3108	2014.6288	0.0000	1.0000
K2PIG	38	-965.3105	2017.2068	0.0000	1.0000
GTR	43	-959.0832	2017.9265	0.0000	1.0000
SYMI	41	-962.7038	2019.8407	0.0000	1.0000
SYMG	41	-962.9014	2020.2360	0.0000	1.0000
SYMIG	42	-962.9011	2022.8891	0.0000	1.0000
F81I	39	-971.5205	2032.2238	0.0000	1.0000
F81G	39	-971.7994	2032.7817	0.0000	1.0000
K2P	36	-976.6504	2034.7475	0.0000	1.0000
F81IG	40	-971.7992	2035.3969	0.0000	1.0000
SYM	40	-974.5510	2040.9006	0.0000	1.0000
JC69I	36	-984.8117	2051.0702	0.0000	1.0000
JC69G	36	-984.9535	2051.3538	0.0000	1.0000
F81	38	-982.4563	2051.4982	0.0000	1.0000
JC69IG	37	-984.9532	2053.9135	0.0000	1.0000
JC69	35	-995.8059	2070.5165	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-949.8979	1979.7958	0.3751	0.3751
HKYG	40	-950.1163	1980.2325	0.3015	0.6767
GTRI	44	-947.1010	1982.2020	0.1126	0.7893
HKYIG	41	-950.1160	1982.2320	0.1110	0.9002
GTRG	44	-947.5359	1983.0718	0.0729	0.9732
GTRIG	45	-947.5356	1985.0713	0.0268	1.0000
HKY	39	-961.4569	2000.9137	0.0000	1.0000
GTR	43	-959.0832	2004.1665	0.0000	1.0000
K2PI	37	-965.2727	2004.5455	0.0000	1.0000
K2PG	37	-965.3108	2004.6217	0.0000	1.0000
K2PIG	38	-965.3105	2006.6211	0.0000	1.0000
SYMI	41	-962.7038	2007.4075	0.0000	1.0000
SYMG	41	-962.9014	2007.8028	0.0000	1.0000
SYMIG	42	-962.9011	2009.8021	0.0000	1.0000
F81I	39	-971.5205	2021.0410	0.0000	1.0000
F81G	39	-971.7994	2021.5989	0.0000	1.0000
F81IG	40	-971.7992	2023.5983	0.0000	1.0000
K2P	36	-976.6504	2025.3007	0.0000	1.0000
SYM	40	-974.5510	2029.1020	0.0000	1.0000
F81	38	-982.4563	2040.9125	0.0000	1.0000
JC69I	36	-984.8117	2041.6234	0.0000	1.0000
JC69G	36	-984.9535	2041.9070	0.0000	1.0000
JC69IG	37	-984.9532	2043.9064	0.0000	1.0000
JC69	35	-995.8059	2061.6119	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-949.8979	2130.4035	0.5400	0.5400
HKYG	40	-950.1163	2130.8402	0.4341	0.9741
HKYIG	41	-950.1160	2136.6049	0.0243	0.9985
K2PI	37	-965.2727	2143.8576	0.0006	0.9991
K2PG	37	-965.3108	2143.9337	0.0006	0.9997
HKY	39	-961.4569	2147.7562	0.0001	0.9998
GTRI	44	-947.1010	2147.8704	0.0001	0.9999
GTRG	44	-947.5359	2148.7402	0.0001	1.0000
K2PIG	38	-965.3105	2149.6983	0.0000	1.0000
GTRIG	45	-947.5356	2154.5049	0.0000	1.0000
K2P	36	-976.6504	2160.8476	0.0000	1.0000
SYMI	41	-962.7038	2161.7804	0.0000	1.0000
SYMG	41	-962.9014	2162.1756	0.0000	1.0000
GTR	43	-959.0832	2166.0697	0.0000	1.0000
F81I	39	-971.5205	2167.8835	0.0000	1.0000
SYMIG	42	-962.9011	2167.9401	0.0000	1.0000
F81G	39	-971.7994	2168.4414	0.0000	1.0000
F81IG	40	-971.7992	2174.2060	0.0000	1.0000
JC69I	36	-984.8117	2177.1703	0.0000	1.0000
JC69G	36	-984.9535	2177.4539	0.0000	1.0000
SYM	40	-974.5510	2179.7097	0.0000	1.0000
JC69IG	37	-984.9532	2183.2185	0.0000	1.0000
F81	38	-982.4563	2183.9898	0.0000	1.0000
JC69	35	-995.8059	2193.3936	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/chr1_15480.phylip = 1 - 319;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_15480.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/chr1_15480.phylip = 1 - 319;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_15480.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/chr1_15480.phylip = 1 - 319;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_15480.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

