
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:52 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_26677.phylip" (19 taxa, 771 characters)

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	41	-2315.9603	4718.6448	0.5362	0.5362
GTRIG	45	-2312.7297	4721.1698	0.1517	0.6879
HKYI	40	-2318.3522	4721.1977	0.1496	0.8375
HKYG	40	-2318.8952	4722.2835	0.0869	0.9244
GTRG	44	-2315.2158	4723.8862	0.0390	0.9634
GTRI	44	-2315.2807	4724.0159	0.0366	1.0000
HKY	39	-2344.5648	4771.3977	0.0000	1.0000
GTR	43	-2340.9752	4773.1554	0.0000	1.0000
K2PIG	38	-2371.3816	4822.8124	0.0000	1.0000
K2PI	37	-2373.7720	4825.3802	0.0000	1.0000
K2PG	37	-2373.8125	4825.4612	0.0000	1.0000
F81IG	40	-2371.2832	4827.0595	0.0000	1.0000
F81I	39	-2373.7683	4829.8047	0.0000	1.0000
F81G	39	-2374.0159	4830.2999	0.0000	1.0000
SYMIG	42	-2371.1711	4831.3037	0.0000	1.0000
SYMI	41	-2373.5071	4833.7386	0.0000	1.0000
SYMG	41	-2373.5973	4833.9189	0.0000	1.0000
K2P	36	-2400.7296	4877.0887	0.0000	1.0000
F81	38	-2399.6448	4879.3388	0.0000	1.0000
SYM	40	-2400.4002	4885.2936	0.0000	1.0000
JC69IG	37	-2418.8568	4915.5498	0.0000	1.0000
JC69G	36	-2420.9206	4917.4705	0.0000	1.0000
JC69I	36	-2421.2750	4918.1795	0.0000	1.0000
JC69	35	-2447.0175	4967.4635	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-2315.9603	4713.9205	0.4856	0.4856
GTRIG	45	-2312.7297	4715.4595	0.2250	0.7106
HKYI	40	-2318.3522	4716.7045	0.1207	0.8313
HKYG	40	-2318.8952	4717.7904	0.0701	0.9014
GTRG	44	-2315.2158	4718.4316	0.0509	0.9523
GTRI	44	-2315.2807	4718.5614	0.0477	1.0000
HKY	39	-2344.5648	4767.1296	0.0000	1.0000
GTR	43	-2340.9752	4767.9505	0.0000	1.0000
K2PIG	38	-2371.3816	4818.7633	0.0000	1.0000
K2PI	37	-2373.7720	4821.5440	0.0000	1.0000
K2PG	37	-2373.8125	4821.6249	0.0000	1.0000
F81IG	40	-2371.2832	4822.5663	0.0000	1.0000
F81I	39	-2373.7683	4825.5366	0.0000	1.0000
F81G	39	-2374.0159	4826.0318	0.0000	1.0000
SYMIG	42	-2371.1711	4826.3422	0.0000	1.0000
SYMI	41	-2373.5071	4829.0143	0.0000	1.0000
SYMG	41	-2373.5973	4829.1946	0.0000	1.0000
K2P	36	-2400.7296	4873.4593	0.0000	1.0000
F81	38	-2399.6448	4875.2896	0.0000	1.0000
SYM	40	-2400.4002	4880.8004	0.0000	1.0000
JC69IG	37	-2418.8568	4911.7135	0.0000	1.0000
JC69G	36	-2420.9206	4913.8411	0.0000	1.0000
JC69I	36	-2421.2750	4914.5501	0.0000	1.0000
JC69	35	-2447.0175	4964.0349	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-2318.3522	4902.6120	0.5063	0.5063
HKYG	40	-2318.8952	4903.6979	0.2942	0.8005
HKYIG	41	-2315.9603	4904.4757	0.1994	1.0000
GTRG	44	-2315.2158	4922.9299	0.0000	1.0000
GTRI	44	-2315.2807	4923.0597	0.0000	1.0000
GTRIG	45	-2312.7297	4924.6055	0.0000	1.0000
HKY	39	-2344.5648	4948.3894	0.0000	1.0000
GTR	43	-2340.9752	4967.8011	0.0000	1.0000
K2PI	37	-2373.7720	4993.5084	0.0000	1.0000
K2PG	37	-2373.8125	4993.5894	0.0000	1.0000
K2PIG	38	-2371.3816	4995.3754	0.0000	1.0000
F81I	39	-2373.7683	5006.7964	0.0000	1.0000
F81G	39	-2374.0159	5007.2916	0.0000	1.0000
F81IG	40	-2371.2832	5008.4738	0.0000	1.0000
SYMI	41	-2373.5071	5019.5695	0.0000	1.0000
SYMG	41	-2373.5973	5019.7498	0.0000	1.0000
SYMIG	42	-2371.1711	5021.5451	0.0000	1.0000
K2P	36	-2400.7296	5040.7760	0.0000	1.0000
F81	38	-2399.6448	5051.9018	0.0000	1.0000
SYM	40	-2400.4002	5066.7080	0.0000	1.0000
JC69G	36	-2420.9206	5081.1579	0.0000	1.0000
JC69I	36	-2421.2750	5081.8668	0.0000	1.0000
JC69IG	37	-2418.8568	5083.6780	0.0000	1.0000
JC69	35	-2447.0175	5126.7040	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26677.phylip = 1 - 771;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26677.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26677.phylip = 1 - 771;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26677.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) 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/chr2_26677.phylip = 1 - 771;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26677.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

