
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:12 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/chr2_24807.phylip" (19 taxa, 300 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	-999.7641	2092.1923	0.5548	0.5548
HKYI	40	-1000.6978	2094.0597	0.2181	0.7729
HKYIG	41	-999.6963	2094.7415	0.1551	0.9280
GTRG	44	-997.4648	2098.4590	0.0242	0.9521
K2PG	37	-1007.1435	2099.0199	0.0183	0.9704
GTRI	44	-998.2707	2100.0707	0.0108	0.9812
K2PI	37	-1008.1183	2100.9694	0.0069	0.9881
GTRIG	45	-997.4509	2101.2010	0.0061	0.9942
K2PIG	38	-1007.0278	2101.4119	0.0055	0.9998
SYMG	41	-1006.6887	2108.7262	0.0001	0.9999
SYMI	41	-1007.5953	2110.5395	0.0001	1.0000
SYMIG	42	-1006.6026	2111.2597	0.0000	1.0000
F81G	39	-1014.7234	2119.4467	0.0000	1.0000
F81I	39	-1015.5769	2121.1539	0.0000	1.0000
F81IG	40	-1014.6658	2121.9957	0.0000	1.0000
HKY	39	-1016.5058	2123.0117	0.0000	1.0000
JC69G	36	-1021.6427	2125.4147	0.0000	1.0000
JC69I	36	-1022.5831	2127.2955	0.0000	1.0000
GTR	43	-1013.4598	2127.7009	0.0000	1.0000
JC69IG	37	-1021.5456	2127.8240	0.0000	1.0000
K2P	36	-1024.7574	2131.6441	0.0000	1.0000
SYM	40	-1023.7830	2140.2302	0.0000	1.0000
F81	38	-1031.2937	2149.9436	0.0000	1.0000
JC69	35	-1038.8211	2157.1877	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-999.7641	2079.5282	0.4668	0.4668
HKYIG	41	-999.6963	2081.3927	0.1838	0.6506
HKYI	40	-1000.6978	2081.3956	0.1835	0.8341
GTRG	44	-997.4648	2082.9296	0.0852	0.9194
GTRI	44	-998.2707	2084.5413	0.0381	0.9574
GTRIG	45	-997.4509	2084.9018	0.0318	0.9892
K2PG	37	-1007.1435	2088.2870	0.0059	0.9951
K2PIG	38	-1007.0278	2090.0555	0.0024	0.9975
K2PI	37	-1008.1183	2090.2366	0.0022	0.9997
SYMG	41	-1006.6887	2095.3773	0.0002	0.9999
SYMI	41	-1007.5953	2097.1907	0.0001	0.9999
SYMIG	42	-1006.6026	2097.2052	0.0001	1.0000
F81G	39	-1014.7234	2107.4467	0.0000	1.0000
F81I	39	-1015.5769	2109.1539	0.0000	1.0000
F81IG	40	-1014.6658	2109.3316	0.0000	1.0000
HKY	39	-1016.5058	2111.0117	0.0000	1.0000
GTR	43	-1013.4598	2112.9197	0.0000	1.0000
JC69G	36	-1021.6427	2115.2854	0.0000	1.0000
JC69IG	37	-1021.5456	2117.0912	0.0000	1.0000
JC69I	36	-1022.5831	2117.1662	0.0000	1.0000
K2P	36	-1024.7574	2121.5148	0.0000	1.0000
SYM	40	-1023.7830	2127.5661	0.0000	1.0000
F81	38	-1031.2937	2138.5873	0.0000	1.0000
JC69	35	-1038.8211	2147.6422	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1007.1435	2225.3270	0.5288	0.5288
K2PI	37	-1008.1183	2227.2765	0.1995	0.7284
HKYG	40	-999.7641	2227.6795	0.1631	0.8915
HKYI	40	-1000.6978	2229.5469	0.0641	0.9556
K2PIG	38	-1007.0278	2230.7993	0.0343	0.9899
HKYIG	41	-999.6963	2233.2477	0.0101	1.0000
GTRG	44	-997.4648	2245.8960	0.0000	1.0000
SYMG	41	-1006.6887	2247.2324	0.0000	1.0000
GTRI	44	-998.2707	2247.5077	0.0000	1.0000
JC69G	36	-1021.6427	2248.6216	0.0000	1.0000
SYMI	41	-1007.5953	2249.0457	0.0000	1.0000
JC69I	36	-1022.5831	2250.5023	0.0000	1.0000
GTRIG	45	-997.4509	2251.5720	0.0000	1.0000
F81G	39	-1014.7234	2251.8943	0.0000	1.0000
SYMIG	42	-1006.6026	2252.7641	0.0000	1.0000
F81I	39	-1015.5769	2253.6014	0.0000	1.0000
JC69IG	37	-1021.5456	2254.1311	0.0000	1.0000
K2P	36	-1024.7574	2254.8510	0.0000	1.0000
HKY	39	-1016.5058	2255.4592	0.0000	1.0000
F81IG	40	-1014.6658	2257.4829	0.0000	1.0000
GTR	43	-1013.4598	2272.1823	0.0000	1.0000
SYM	40	-1023.7830	2275.7174	0.0000	1.0000
JC69	35	-1038.8211	2277.2746	0.0000	1.0000
F81	38	-1031.2937	2279.3311	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/chr2_24807.phylip = 1 - 300;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24807.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/chr2_24807.phylip = 1 - 300;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24807.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/chr2_24807.phylip = 1 - 300;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24807.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

