
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:29 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/chr6_7266.phylip" (19 taxa, 467 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	-2150.2799	4388.2594	0.6261	0.6261
GTRG	44	-2146.7301	4390.8441	0.1719	0.7980
HKYIG	41	-2150.4733	4391.0502	0.1551	0.9531
GTRIG	45	-2146.8558	4393.5454	0.0445	0.9977
HKYI	40	-2156.4309	4400.5613	0.0013	0.9990
GTRI	44	-2151.8797	4401.1432	0.0010	1.0000
SYMG	41	-2195.5569	4481.2173	0.0000	1.0000
F81G	39	-2198.4042	4482.1152	0.0000	1.0000
SYMIG	42	-2195.4792	4483.4773	0.0000	1.0000
F81IG	40	-2198.1790	4484.0576	0.0000	1.0000
F81I	39	-2200.9325	4487.1718	0.0000	1.0000
SYMI	41	-2201.6313	4493.3662	0.0000	1.0000
K2PG	37	-2208.5930	4497.7409	0.0000	1.0000
K2PIG	38	-2208.6209	4500.1671	0.0000	1.0000
K2PI	37	-2212.6750	4505.9047	0.0000	1.0000
GTR	43	-2239.8760	4574.6975	0.0000	1.0000
JC69G	36	-2251.0201	4580.2356	0.0000	1.0000
JC69IG	37	-2251.0995	4582.7538	0.0000	1.0000
HKY	39	-2249.5570	4584.4208	0.0000	1.0000
JC69I	36	-2254.0448	4586.2848	0.0000	1.0000
F81	38	-2288.6268	4660.1788	0.0000	1.0000
SYM	40	-2290.9613	4669.6222	0.0000	1.0000
K2P	36	-2301.2443	4680.6839	0.0000	1.0000
JC69	35	-2337.5309	4750.9086	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2150.2799	4380.5598	0.4644	0.4644
GTRG	44	-2146.7301	4381.4602	0.2961	0.7604
HKYIG	41	-2150.4733	4382.9467	0.1408	0.9012
GTRIG	45	-2146.8558	4383.7116	0.0960	0.9973
GTRI	44	-2151.8797	4391.7594	0.0017	0.9990
HKYI	40	-2156.4309	4392.8618	0.0010	1.0000
SYMG	41	-2195.5569	4473.1138	0.0000	1.0000
F81G	39	-2198.4042	4474.8084	0.0000	1.0000
SYMIG	42	-2195.4792	4474.9584	0.0000	1.0000
F81IG	40	-2198.1790	4476.3581	0.0000	1.0000
F81I	39	-2200.9325	4479.8650	0.0000	1.0000
SYMI	41	-2201.6313	4485.2627	0.0000	1.0000
K2PG	37	-2208.5930	4491.1861	0.0000	1.0000
K2PIG	38	-2208.6209	4493.2419	0.0000	1.0000
K2PI	37	-2212.6750	4499.3499	0.0000	1.0000
GTR	43	-2239.8760	4565.7519	0.0000	1.0000
JC69G	36	-2251.0201	4574.0403	0.0000	1.0000
JC69IG	37	-2251.0995	4576.1990	0.0000	1.0000
HKY	39	-2249.5570	4577.1140	0.0000	1.0000
JC69I	36	-2254.0448	4580.0895	0.0000	1.0000
F81	38	-2288.6268	4653.2536	0.0000	1.0000
SYM	40	-2290.9613	4661.9227	0.0000	1.0000
K2P	36	-2301.2443	4674.4886	0.0000	1.0000
JC69	35	-2337.5309	4745.0618	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2150.2799	4546.4130	0.9611	0.9611
HKYIG	41	-2150.4733	4552.9462	0.0367	0.9978
HKYI	40	-2156.4309	4558.7150	0.0020	0.9998
GTRG	44	-2146.7301	4563.8987	0.0002	1.0000
GTRIG	45	-2146.8558	4570.2965	0.0000	1.0000
GTRI	44	-2151.8797	4574.1978	0.0000	1.0000
F81G	39	-2198.4042	4636.5153	0.0000	1.0000
F81I	39	-2200.9325	4641.5718	0.0000	1.0000
F81IG	40	-2198.1790	4642.2113	0.0000	1.0000
SYMG	41	-2195.5569	4643.1133	0.0000	1.0000
K2PG	37	-2208.5930	4644.6003	0.0000	1.0000
SYMIG	42	-2195.4792	4649.1043	0.0000	1.0000
K2PIG	38	-2208.6209	4650.8024	0.0000	1.0000
K2PI	37	-2212.6750	4652.7641	0.0000	1.0000
SYMI	41	-2201.6313	4655.2622	0.0000	1.0000
JC69G	36	-2251.0201	4723.3081	0.0000	1.0000
JC69I	36	-2254.0448	4729.3574	0.0000	1.0000
JC69IG	37	-2251.0995	4729.6132	0.0000	1.0000
HKY	39	-2249.5570	4738.8208	0.0000	1.0000
GTR	43	-2239.8760	4744.0441	0.0000	1.0000
F81	38	-2288.6268	4810.8141	0.0000	1.0000
K2P	36	-2301.2443	4823.7564	0.0000	1.0000
SYM	40	-2290.9613	4827.7759	0.0000	1.0000
JC69	35	-2337.5309	4890.1833	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/chr6_7266.phylip = 1 - 467;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7266.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/chr6_7266.phylip = 1 - 467;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7266.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/chr6_7266.phylip = 1 - 467;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7266.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

