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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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	-964.2997	2017.4403	0.4513	0.4513
HKYI	40	-964.4123	2017.6656	0.4032	0.8546
HKYIG	41	-964.2634	2019.8349	0.1363	0.9909
GTRG	44	-964.0352	2026.8606	0.0041	0.9949
GTRI	44	-964.1354	2027.0610	0.0037	0.9986
GTRIG	45	-964.0041	2029.3198	0.0012	0.9998
HKY	39	-973.5027	2033.3925	0.0002	1.0000
F81G	39	-975.8968	2038.1808	0.0000	1.0000
F81I	39	-975.9870	2038.3611	0.0000	1.0000
K2PG	37	-978.9757	2039.4700	0.0000	1.0000
K2PI	37	-979.1255	2039.7698	0.0000	1.0000
F81IG	40	-975.8810	2040.6029	0.0000	1.0000
K2PIG	38	-978.9309	2041.8082	0.0000	1.0000
GTR	43	-973.3653	2043.0132	0.0000	1.0000
SYMG	41	-978.1686	2047.6453	0.0000	1.0000
SYMI	41	-978.3552	2048.0185	0.0000	1.0000
SYMIG	42	-978.1037	2049.9961	0.0000	1.0000
F81	38	-985.0094	2053.9653	0.0000	1.0000
K2P	36	-987.7739	2054.6517	0.0000	1.0000
JC69G	36	-989.2276	2057.5593	0.0000	1.0000
JC69I	36	-989.3426	2057.7892	0.0000	1.0000
JC69IG	37	-989.2031	2059.9249	0.0000	1.0000
SYM	40	-987.0264	2062.8937	0.0000	1.0000
JC69	35	-998.1613	2073.0246	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-964.2997	2008.5994	0.4292	0.4292
HKYI	40	-964.4123	2008.8246	0.3835	0.8127
HKYIG	41	-964.2634	2010.5268	0.1637	0.9765
GTRG	44	-964.0352	2016.0704	0.0102	0.9867
GTRI	44	-964.1354	2016.2708	0.0093	0.9960
GTRIG	45	-964.0041	2018.0083	0.0039	0.9998
HKY	39	-973.5027	2025.0054	0.0001	1.0000
F81G	39	-975.8968	2029.7937	0.0000	1.0000
F81I	39	-975.9870	2029.9740	0.0000	1.0000
F81IG	40	-975.8810	2031.7619	0.0000	1.0000
K2PG	37	-978.9757	2031.9513	0.0000	1.0000
K2PI	37	-979.1255	2032.2511	0.0000	1.0000
GTR	43	-973.3653	2032.7306	0.0000	1.0000
K2PIG	38	-978.9309	2033.8619	0.0000	1.0000
SYMG	41	-978.1686	2038.3372	0.0000	1.0000
SYMI	41	-978.3552	2038.7104	0.0000	1.0000
SYMIG	42	-978.1037	2040.2074	0.0000	1.0000
F81	38	-985.0094	2046.0189	0.0000	1.0000
K2P	36	-987.7739	2047.5477	0.0000	1.0000
JC69G	36	-989.2276	2050.4553	0.0000	1.0000
JC69I	36	-989.3426	2050.6852	0.0000	1.0000
JC69IG	37	-989.2031	2052.4061	0.0000	1.0000
SYM	40	-987.0264	2054.0528	0.0000	1.0000
JC69	35	-998.1613	2066.3225	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-964.2997	2169.4403	0.5118	0.5118
HKYI	40	-964.4123	2169.6655	0.4573	0.9691
HKYIG	41	-964.2634	2175.3887	0.0261	0.9953
K2PG	37	-978.9757	2180.7292	0.0018	0.9971
K2PI	37	-979.1255	2181.0289	0.0016	0.9987
HKY	39	-973.5027	2181.8253	0.0010	0.9997
F81G	39	-975.8968	2186.6136	0.0001	0.9998
K2PIG	38	-978.9309	2186.6607	0.0001	0.9999
F81I	39	-975.9870	2186.7939	0.0001	1.0000
K2P	36	-987.7739	2192.3046	0.0000	1.0000
F81IG	40	-975.8810	2192.6029	0.0000	1.0000
GTRG	44	-964.0352	2192.9954	0.0000	1.0000
GTRI	44	-964.1354	2193.1958	0.0000	1.0000
JC69G	36	-989.2276	2195.2121	0.0000	1.0000
JC69I	36	-989.3426	2195.4420	0.0000	1.0000
F81	38	-985.0094	2198.8178	0.0000	1.0000
GTRIG	45	-964.0041	2198.9544	0.0000	1.0000
JC69IG	37	-989.2031	2201.1840	0.0000	1.0000
SYMG	41	-978.1686	2203.1992	0.0000	1.0000
SYMI	41	-978.3552	2203.5723	0.0000	1.0000
GTR	43	-973.3653	2205.6346	0.0000	1.0000
JC69	35	-998.1613	2207.0583	0.0000	1.0000
SYMIG	42	-978.1037	2209.0904	0.0000	1.0000
SYM	40	-987.0264	2214.8937	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_10443.phylip = 1 - 412;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10443.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_10443.phylip = 1 - 412;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10443.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 (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10443.phylip = 1 - 412;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10443.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;


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

