
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:09 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/chr19_372.phylip" (19 taxa, 280 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	-810.1402	1714.0042	0.4933	0.4933
HKYG	40	-810.6446	1715.0131	0.2979	0.7912
GTRI	44	-806.3137	1717.4784	0.0868	0.8780
HKYIG	41	-810.6443	1717.7592	0.0755	0.9535
GTRG	44	-807.2828	1719.4167	0.0329	0.9864
GTRIG	45	-807.2830	1722.2582	0.0080	0.9944
K2PI	37	-819.3727	1724.3652	0.0028	0.9972
K2PG	37	-820.0309	1725.6816	0.0014	0.9986
SYMI	41	-815.3399	1727.1504	0.0007	0.9993
K2PIG	38	-820.0307	1728.3602	0.0004	0.9997
SYMG	41	-816.3156	1729.1018	0.0003	0.9999
SYMIG	42	-816.3161	1731.8727	0.0001	1.0000
HKY	39	-827.8244	1746.6489	0.0000	1.0000
F81I	39	-830.1536	1751.3073	0.0000	1.0000
F81G	39	-830.2532	1751.5065	0.0000	1.0000
GTR	43	-825.8571	1753.7481	0.0000	1.0000
F81IG	40	-830.2529	1754.2296	0.0000	1.0000
K2P	36	-836.4040	1755.7709	0.0000	1.0000
JC69I	36	-837.3962	1757.7553	0.0000	1.0000
JC69G	36	-837.6055	1758.1739	0.0000	1.0000
JC69IG	37	-837.6052	1760.8302	0.0000	1.0000
SYM	40	-834.0896	1761.9031	0.0000	1.0000
F81	38	-844.6801	1777.6590	0.0000	1.0000
JC69	35	-851.8030	1783.9339	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-810.1402	1700.2803	0.3216	0.3216
GTRI	44	-806.3137	1700.6273	0.2704	0.5920
HKYG	40	-810.6446	1701.2892	0.1942	0.7862
GTRG	44	-807.2828	1702.5656	0.1026	0.8887
HKYIG	41	-810.6443	1703.2886	0.0715	0.9602
GTRIG	45	-807.2830	1704.5659	0.0377	0.9979
SYMI	41	-815.3399	1712.6798	0.0007	0.9986
K2PI	37	-819.3727	1712.7454	0.0006	0.9992
K2PG	37	-820.0309	1714.0617	0.0003	0.9995
SYMG	41	-816.3156	1714.6312	0.0002	0.9998
K2PIG	38	-820.0307	1716.0614	0.0001	0.9999
SYMIG	42	-816.3161	1716.6322	0.0001	1.0000
HKY	39	-827.8244	1733.6489	0.0000	1.0000
GTR	43	-825.8571	1737.7142	0.0000	1.0000
F81I	39	-830.1536	1738.3073	0.0000	1.0000
F81G	39	-830.2532	1738.5065	0.0000	1.0000
F81IG	40	-830.2529	1740.5057	0.0000	1.0000
K2P	36	-836.4040	1744.8079	0.0000	1.0000
JC69I	36	-837.3962	1746.7923	0.0000	1.0000
JC69G	36	-837.6055	1747.2110	0.0000	1.0000
SYM	40	-834.0896	1748.1793	0.0000	1.0000
JC69IG	37	-837.6052	1749.2104	0.0000	1.0000
F81	38	-844.6801	1765.3602	0.0000	1.0000
JC69	35	-851.8030	1773.6060	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-810.1402	1845.6719	0.4254	0.4254
HKYG	40	-810.6446	1846.6808	0.2569	0.6822
K2PI	37	-819.3727	1847.2326	0.1949	0.8771
K2PG	37	-820.0309	1848.5490	0.1009	0.9781
HKYIG	41	-810.6443	1852.3150	0.0154	0.9934
K2PIG	38	-820.0307	1854.1834	0.0060	0.9995
GTRI	44	-806.3137	1860.5581	0.0002	0.9997
SYMI	41	-815.3399	1861.7062	0.0001	0.9998
GTRG	44	-807.2828	1862.4964	0.0001	0.9999
SYMG	41	-816.3156	1863.6576	0.0001	1.0000
GTRIG	45	-807.2830	1868.1315	0.0000	1.0000
SYMIG	42	-816.3161	1869.2933	0.0000	1.0000
HKY	39	-827.8244	1875.4057	0.0000	1.0000
K2P	36	-836.4040	1875.6603	0.0000	1.0000
JC69I	36	-837.3962	1877.6447	0.0000	1.0000
JC69G	36	-837.6055	1878.0634	0.0000	1.0000
F81I	39	-830.1536	1880.0641	0.0000	1.0000
F81G	39	-830.2532	1880.2633	0.0000	1.0000
JC69IG	37	-837.6052	1883.6976	0.0000	1.0000
F81IG	40	-830.2529	1885.8973	0.0000	1.0000
SYM	40	-834.0896	1893.5708	0.0000	1.0000
GTR	43	-825.8571	1894.0102	0.0000	1.0000
JC69	35	-851.8030	1900.8236	0.0000	1.0000
F81	38	-844.6801	1903.4822	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/chr19_372.phylip = 1 - 280;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr19_372.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/chr19_372.phylip = 1 - 280;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr19_372.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/chr19_372.phylip = 1 - 280;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr19_372.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

