
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:38 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_6615.phylip" (19 taxa, 347 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	-950.7766	1992.2721	0.3356	0.3356
HKYI	40	-950.9492	1992.6174	0.2824	0.6180
K2PG	37	-955.1167	1993.3336	0.1974	0.8154
HKYIG	41	-950.7713	1994.8344	0.0932	0.9087
K2PI	37	-956.6484	1996.3971	0.0427	0.9513
K2PIG	38	-956.4148	1998.4530	0.0153	0.9666
GTRG	44	-949.0742	1999.2610	0.0102	0.9768
GTRI	44	-949.2628	1999.6381	0.0084	0.9852
SYMG	41	-953.5003	2000.2923	0.0061	0.9913
GTRIG	45	-949.0667	2001.8875	0.0027	0.9940
HKY	39	-956.8982	2001.9594	0.0026	0.9967
K2P	36	-961.3006	2003.1947	0.0014	0.9981
SYMI	41	-955.0558	2003.4035	0.0013	0.9994
SYMIG	42	-954.8069	2005.4953	0.0005	0.9998
GTR	43	-955.2838	2009.0559	0.0001	0.9999
SYM	40	-959.6901	2010.0992	0.0000	1.0000
F81G	39	-962.5376	2013.2380	0.0000	1.0000
JC69G	36	-966.4803	2013.5542	0.0000	1.0000
F81I	39	-962.7206	2013.6041	0.0000	1.0000
F81IG	40	-962.5319	2015.7827	0.0000	1.0000
JC69I	36	-968.0232	2016.6400	0.0000	1.0000
JC69IG	37	-967.7982	2018.6967	0.0000	1.0000
F81	38	-968.4001	2022.4236	0.0000	1.0000
JC69	35	-972.3747	2022.8524	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-950.7766	1981.5531	0.3550	0.3550
HKYI	40	-950.9492	1981.8984	0.2987	0.6537
HKYIG	41	-950.7713	1983.5426	0.1313	0.7850
K2PG	37	-955.1167	1984.2333	0.0930	0.8780
GTRG	44	-949.0742	1986.1484	0.0357	0.9137
GTRI	44	-949.2628	1986.5256	0.0295	0.9432
K2PI	37	-956.6484	1987.2968	0.0201	0.9633
GTRIG	45	-949.0667	1988.1334	0.0132	0.9765
K2PIG	38	-956.4148	1988.8296	0.0093	0.9859
SYMG	41	-953.5003	1989.0005	0.0086	0.9944
HKY	39	-956.8982	1991.7965	0.0021	0.9966
SYMI	41	-955.0558	1992.1117	0.0018	0.9984
SYMIG	42	-954.8069	1993.6137	0.0009	0.9992
K2P	36	-961.3006	1994.6011	0.0005	0.9997
GTR	43	-955.2838	1996.5675	0.0002	0.9999
SYM	40	-959.6901	1999.3803	0.0000	1.0000
F81G	39	-962.5376	2003.0751	0.0000	1.0000
F81I	39	-962.7206	2003.4413	0.0000	1.0000
JC69G	36	-966.4803	2004.9606	0.0000	1.0000
F81IG	40	-962.5319	2005.0637	0.0000	1.0000
JC69I	36	-968.0232	2008.0464	0.0000	1.0000
JC69IG	37	-967.7982	2009.5964	0.0000	1.0000
F81	38	-968.4001	2012.8002	0.0000	1.0000
JC69	35	-972.3747	2014.7495	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-955.1167	2126.6583	0.7737	0.7737
K2PI	37	-956.6484	2129.7218	0.1672	0.9409
K2P	36	-961.3006	2133.1768	0.0297	0.9706
K2PIG	38	-956.4148	2135.1040	0.0113	0.9820
HKYG	40	-950.7766	2135.5261	0.0092	0.9911
HKYI	40	-950.9492	2135.8714	0.0077	0.9989
HKYIG	41	-950.7713	2141.3649	0.0005	0.9994
HKY	39	-956.8982	2141.9202	0.0004	0.9997
JC69G	36	-966.4803	2143.5363	0.0002	0.9999
JC69I	36	-968.0232	2146.6221	0.0000	0.9999
SYMG	41	-953.5003	2146.8228	0.0000	1.0000
JC69	35	-972.3747	2149.4759	0.0000	1.0000
SYMI	41	-955.0558	2149.9340	0.0000	1.0000
JC69IG	37	-967.7982	2152.0214	0.0000	1.0000
F81G	39	-962.5376	2153.1988	0.0000	1.0000
SYM	40	-959.6901	2153.3533	0.0000	1.0000
F81I	39	-962.7206	2153.5649	0.0000	1.0000
SYMIG	42	-954.8069	2155.2854	0.0000	1.0000
GTRG	44	-949.0742	2155.5187	0.0000	1.0000
GTRI	44	-949.2628	2155.8959	0.0000	1.0000
F81IG	40	-962.5319	2159.0367	0.0000	1.0000
F81	38	-968.4001	2159.0746	0.0000	1.0000
GTRIG	45	-949.0667	2161.3530	0.0000	1.0000
GTR	43	-955.2838	2162.0885	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/chrZ_6615.phylip = 1 - 347;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6615.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/chrZ_6615.phylip = 1 - 347;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6615.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/chrZ_6615.phylip = 1 - 347;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6615.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

