
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:05 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/chr8_3358.phylip" (19 taxa, 336 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	-995.0249	2081.1685	0.6470	0.6470
HKYIG	41	-994.9656	2083.6456	0.1875	0.8345
HKYI	40	-996.9094	2084.9374	0.0983	0.9328
GTRG	44	-992.4279	2086.4640	0.0458	0.9786
GTRIG	45	-992.3720	2089.0198	0.0128	0.9913
GTRI	44	-994.3070	2090.2223	0.0070	0.9983
HKY	39	-1002.3415	2093.2235	0.0016	0.9999
GTR	43	-999.9202	2098.7993	0.0001	1.0000
K2PG	37	-1014.9234	2113.2830	0.0000	1.0000
K2PIG	38	-1014.8012	2115.5822	0.0000	1.0000
K2PI	37	-1016.7355	2116.9071	0.0000	1.0000
SYMG	41	-1013.4509	2120.6160	0.0000	1.0000
SYMIG	42	-1013.3276	2122.9828	0.0000	1.0000
SYMI	41	-1015.3493	2124.4128	0.0000	1.0000
K2P	36	-1023.0260	2126.9617	0.0000	1.0000
F81G	39	-1022.2864	2133.1133	0.0000	1.0000
SYM	40	-1021.5902	2134.2990	0.0000	1.0000
F81IG	40	-1022.2217	2135.5621	0.0000	1.0000
F81I	39	-1024.1003	2136.7412	0.0000	1.0000
F81	38	-1029.4347	2144.8492	0.0000	1.0000
JC69G	36	-1039.5541	2160.0180	0.0000	1.0000
JC69IG	37	-1039.4576	2162.3514	0.0000	1.0000
JC69I	36	-1041.4266	2163.7629	0.0000	1.0000
JC69	35	-1047.1558	2172.7117	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-995.0249	2070.0499	0.5199	0.5199
HKYIG	41	-994.9656	2071.9313	0.2029	0.7228
GTRG	44	-992.4279	2072.8558	0.1278	0.8506
HKYI	40	-996.9094	2073.8188	0.0790	0.9296
GTRIG	45	-992.3720	2074.7440	0.0497	0.9793
GTRI	44	-994.3070	2076.6141	0.0195	0.9989
HKY	39	-1002.3415	2082.6830	0.0009	0.9998
GTR	43	-999.9202	2085.8404	0.0002	1.0000
K2PG	37	-1014.9234	2103.8468	0.0000	1.0000
K2PIG	38	-1014.8012	2105.6024	0.0000	1.0000
K2PI	37	-1016.7355	2107.4709	0.0000	1.0000
SYMG	41	-1013.4509	2108.9018	0.0000	1.0000
SYMIG	42	-1013.3276	2110.6552	0.0000	1.0000
SYMI	41	-1015.3493	2112.6985	0.0000	1.0000
K2P	36	-1023.0260	2118.0520	0.0000	1.0000
F81G	39	-1022.2864	2122.5728	0.0000	1.0000
SYM	40	-1021.5902	2123.1803	0.0000	1.0000
F81IG	40	-1022.2217	2124.4435	0.0000	1.0000
F81I	39	-1024.1003	2126.2006	0.0000	1.0000
F81	38	-1029.4347	2134.8694	0.0000	1.0000
JC69G	36	-1039.5541	2151.1083	0.0000	1.0000
JC69IG	37	-1039.4576	2152.9152	0.0000	1.0000
JC69I	36	-1041.4266	2154.8532	0.0000	1.0000
JC69	35	-1047.1558	2164.3117	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-995.0249	2222.7343	0.8182	0.8182
HKYI	40	-996.9094	2226.5032	0.1243	0.9425
HKYIG	41	-994.9656	2228.4328	0.0474	0.9899
HKY	39	-1002.3415	2231.5503	0.0100	0.9999
GTRG	44	-992.4279	2240.8087	0.0001	1.0000
GTRI	44	-994.3070	2244.5670	0.0000	1.0000
K2PG	37	-1014.9234	2245.0799	0.0000	1.0000
GTRIG	45	-992.3720	2246.5140	0.0000	1.0000
K2PI	37	-1016.7355	2248.7040	0.0000	1.0000
GTR	43	-999.9202	2249.9762	0.0000	1.0000
K2PIG	38	-1014.8012	2250.6526	0.0000	1.0000
K2P	36	-1023.0260	2255.4680	0.0000	1.0000
SYMG	41	-1013.4509	2265.4033	0.0000	1.0000
SYMI	41	-1015.3493	2269.2001	0.0000	1.0000
SYMIG	42	-1013.3276	2270.9738	0.0000	1.0000
F81G	39	-1022.2864	2271.4401	0.0000	1.0000
F81I	39	-1024.1003	2275.0680	0.0000	1.0000
SYM	40	-1021.5902	2275.8648	0.0000	1.0000
F81IG	40	-1022.2217	2277.1279	0.0000	1.0000
F81	38	-1029.4347	2279.9196	0.0000	1.0000
JC69G	36	-1039.5541	2288.5243	0.0000	1.0000
JC69I	36	-1041.4266	2292.2692	0.0000	1.0000
JC69IG	37	-1039.4576	2294.1483	0.0000	1.0000
JC69	35	-1047.1558	2297.9106	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/chr8_3358.phylip = 1 - 336;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_3358.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/chr8_3358.phylip = 1 - 336;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_3358.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/chr8_3358.phylip = 1 - 336;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_3358.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

