
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:53 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_10694.phylip" (19 taxa, 847 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	-2397.7227	4879.5149	0.5324	0.5324
HKYIG	41	-2397.3517	4880.9816	0.2557	0.7881
HKYI	40	-2398.8551	4881.7796	0.1716	0.9597
GTRG	44	-2396.5007	4885.9391	0.0214	0.9812
GTRIG	45	-2396.1256	4887.4198	0.0102	0.9914
GTRI	44	-2397.4123	4887.7622	0.0086	1.0000
HKY	39	-2436.0927	4954.0516	0.0000	1.0000
GTR	43	-2434.9034	4960.5192	0.0000	1.0000
K2PG	37	-2456.1106	4989.6971	0.0000	1.0000
K2PIG	38	-2455.5666	4990.8015	0.0000	1.0000
K2PI	37	-2457.4237	4992.3233	0.0000	1.0000
SYMG	41	-2453.5903	4993.4588	0.0000	1.0000
SYMIG	42	-2452.9554	4994.4033	0.0000	1.0000
F81G	39	-2456.3630	4994.5923	0.0000	1.0000
F81IG	40	-2455.9832	4996.0359	0.0000	1.0000
SYMI	41	-2455.0574	4996.3931	0.0000	1.0000
F81I	39	-2457.5556	4996.9774	0.0000	1.0000
F81	38	-2491.2660	5062.2003	0.0000	1.0000
K2P	36	-2493.8608	5063.0105	0.0000	1.0000
SYM	40	-2491.3926	5066.8548	0.0000	1.0000
JC69G	36	-2505.3988	5086.0864	0.0000	1.0000
JC69IG	37	-2504.9273	5087.3306	0.0000	1.0000
JC69I	36	-2506.7413	5088.7715	0.0000	1.0000
JC69	35	-2540.6948	5154.4970	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2397.7227	4875.4454	0.5061	0.5061
HKYIG	41	-2397.3517	4876.7033	0.2698	0.7759
HKYI	40	-2398.8551	4877.7101	0.1631	0.9391
GTRG	44	-2396.5007	4881.0014	0.0315	0.9705
GTRIG	45	-2396.1256	4882.2512	0.0168	0.9874
GTRI	44	-2397.4123	4882.8246	0.0126	1.0000
HKY	39	-2436.0927	4950.1854	0.0000	1.0000
GTR	43	-2434.9034	4955.8069	0.0000	1.0000
K2PG	37	-2456.1106	4986.2212	0.0000	1.0000
K2PIG	38	-2455.5666	4987.1331	0.0000	1.0000
K2PI	37	-2457.4237	4988.8474	0.0000	1.0000
SYMG	41	-2453.5903	4989.1806	0.0000	1.0000
SYMIG	42	-2452.9554	4989.9107	0.0000	1.0000
F81G	39	-2456.3630	4990.7261	0.0000	1.0000
F81IG	40	-2455.9832	4991.9664	0.0000	1.0000
SYMI	41	-2455.0574	4992.1149	0.0000	1.0000
F81I	39	-2457.5556	4993.1113	0.0000	1.0000
F81	38	-2491.2660	5058.5320	0.0000	1.0000
K2P	36	-2493.8608	5059.7216	0.0000	1.0000
SYM	40	-2491.3926	5062.7853	0.0000	1.0000
JC69G	36	-2505.3988	5082.7975	0.0000	1.0000
JC69IG	37	-2504.9273	5083.8547	0.0000	1.0000
JC69I	36	-2506.7413	5085.4826	0.0000	1.0000
JC69	35	-2540.6948	5151.3897	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2397.7227	5065.1134	0.7288	0.7288
HKYI	40	-2398.8551	5067.3781	0.2349	0.9637
HKYIG	41	-2397.3517	5071.1130	0.0363	1.0000
GTRG	44	-2396.5007	5089.6363	0.0000	1.0000
GTRI	44	-2397.4123	5091.4594	0.0000	1.0000
GTRIG	45	-2396.1256	5095.6278	0.0000	1.0000
HKY	39	-2436.0927	5135.1117	0.0000	1.0000
GTR	43	-2434.9034	5159.7000	0.0000	1.0000
K2PG	37	-2456.1106	5161.6641	0.0000	1.0000
K2PI	37	-2457.4237	5164.2903	0.0000	1.0000
K2PIG	38	-2455.5666	5167.3178	0.0000	1.0000
F81G	39	-2456.3630	5175.6524	0.0000	1.0000
F81I	39	-2457.5556	5178.0376	0.0000	1.0000
F81IG	40	-2455.9832	5181.6344	0.0000	1.0000
SYMG	41	-2453.5903	5183.5903	0.0000	1.0000
SYMI	41	-2455.0574	5186.5246	0.0000	1.0000
SYMIG	42	-2452.9554	5189.0622	0.0000	1.0000
K2P	36	-2493.8608	5230.4229	0.0000	1.0000
F81	38	-2491.2660	5238.7166	0.0000	1.0000
SYM	40	-2491.3926	5252.4533	0.0000	1.0000
JC69G	36	-2505.3988	5253.4988	0.0000	1.0000
JC69I	36	-2506.7413	5256.1838	0.0000	1.0000
JC69IG	37	-2504.9273	5259.2976	0.0000	1.0000
JC69	35	-2540.6948	5317.3492	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_10694.phylip = 1 - 847;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10694.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_10694.phylip = 1 - 847;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10694.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_10694.phylip = 1 - 847;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10694.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

