
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:28 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_8059.phylip" (19 taxa, 409 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2PI


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
GTRI	44	-850.8866	1800.6523	0.3113	0.3113
HKYI	40	-855.9107	1800.7343	0.2988	0.6101
GTRG	44	-851.9041	1802.6873	0.1125	0.7226
GTRIG	45	-850.8879	1803.1808	0.0879	0.8105
HKYIG	41	-855.9110	1803.2062	0.0868	0.8974
HKYG	40	-857.1492	1803.2114	0.0866	0.9839
K2PI	37	-863.5805	1808.7405	0.0055	0.9894
SYMI	41	-858.9703	1809.3248	0.0041	0.9935
K2PG	37	-864.6386	1810.8568	0.0019	0.9954
SYMG	41	-859.7666	1810.9175	0.0018	0.9972
K2PIG	38	-863.5811	1811.1729	0.0016	0.9988
SYMIG	42	-858.9663	1811.8014	0.0012	1.0000
GTR	43	-864.4523	1825.2717	0.0000	1.0000
HKY	39	-869.9028	1826.2609	0.0000	1.0000
F81I	39	-872.6043	1831.6638	0.0000	1.0000
SYM	40	-871.7559	1832.4248	0.0000	1.0000
K2P	36	-877.0415	1833.2444	0.0000	1.0000
F81G	39	-873.6514	1833.7581	0.0000	1.0000
F81IG	40	-872.6045	1834.1221	0.0000	1.0000
JC69I	36	-879.4938	1838.1490	0.0000	1.0000
JC69G	36	-880.4473	1840.0558	0.0000	1.0000
JC69IG	37	-879.4940	1840.5676	0.0000	1.0000
F81	38	-887.5908	1859.1925	0.0000	1.0000
JC69	35	-892.6580	1862.0721	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	44	-850.8866	1789.7731	0.4271	0.4271
GTRIG	45	-850.8879	1791.7758	0.1569	0.5841
GTRG	44	-851.9041	1791.8082	0.1544	0.7385
HKYI	40	-855.9107	1791.8213	0.1534	0.8918
HKYIG	41	-855.9110	1793.8220	0.0564	0.9483
HKYG	40	-857.1492	1794.2983	0.0445	0.9927
SYMI	41	-858.9703	1799.9406	0.0026	0.9954
K2PI	37	-863.5805	1801.1609	0.0014	0.9968
SYMG	41	-859.7666	1801.5333	0.0012	0.9980
SYMIG	42	-858.9663	1801.9326	0.0010	0.9990
K2PIG	38	-863.5811	1803.1621	0.0005	0.9995
K2PG	37	-864.6386	1803.2773	0.0005	1.0000
GTR	43	-864.4523	1814.9045	0.0000	1.0000
HKY	39	-869.9028	1817.8056	0.0000	1.0000
F81I	39	-872.6043	1823.2085	0.0000	1.0000
SYM	40	-871.7559	1823.5117	0.0000	1.0000
F81IG	40	-872.6045	1825.2090	0.0000	1.0000
F81G	39	-873.6514	1825.3028	0.0000	1.0000
K2P	36	-877.0415	1826.0831	0.0000	1.0000
JC69I	36	-879.4938	1830.9877	0.0000	1.0000
JC69G	36	-880.4473	1832.8945	0.0000	1.0000
JC69IG	37	-879.4940	1832.9881	0.0000	1.0000
F81	38	-887.5908	1851.1816	0.0000	1.0000
JC69	35	-892.6580	1855.3160	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-863.5805	1949.6684	0.5732	0.5732
K2PG	37	-864.6386	1951.7848	0.1989	0.7721
HKYI	40	-855.9107	1952.3699	0.1485	0.9206
HKYG	40	-857.1492	1954.8469	0.0430	0.9636
K2PIG	38	-863.5811	1955.6833	0.0283	0.9919
HKYIG	41	-855.9110	1958.3843	0.0073	0.9993
SYMI	41	-858.9703	1964.5029	0.0003	0.9996
SYMG	41	-859.7666	1966.0956	0.0002	0.9998
GTRI	44	-850.8866	1966.3766	0.0001	0.9999
GTRG	44	-851.9041	1968.4116	0.0000	1.0000
SYMIG	42	-858.9663	1970.5086	0.0000	1.0000
K2P	36	-877.0415	1970.5768	0.0000	1.0000
GTRIG	45	-850.8879	1972.3930	0.0000	1.0000
HKY	39	-869.9028	1974.3405	0.0000	1.0000
JC69I	36	-879.4938	1975.4814	0.0000	1.0000
JC69G	36	-880.4473	1977.3883	0.0000	1.0000
F81I	39	-872.6043	1979.7434	0.0000	1.0000
JC69IG	37	-879.4940	1981.4955	0.0000	1.0000
F81G	39	-873.6514	1981.8377	0.0000	1.0000
SYM	40	-871.7559	1984.0603	0.0000	1.0000
F81IG	40	-872.6045	1985.7576	0.0000	1.0000
GTR	43	-864.4523	1987.4943	0.0000	1.0000
JC69	35	-892.6580	1995.7961	0.0000	1.0000
F81	38	-887.5908	2003.7028	0.0000	1.0000



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

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

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


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

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


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

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


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

