
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:11 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/chr3_457.phylip" (19 taxa, 774 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	-2025.8775	4136.2297	0.3389	0.3389
HKYI	40	-2026.1409	4136.7565	0.2604	0.5993
GTRG	44	-2022.2495	4137.9312	0.1447	0.7441
HKYIG	41	-2025.8756	4138.4562	0.1113	0.8554
GTRI	44	-2022.6599	4138.7519	0.0960	0.9514
GTRIG	45	-2022.2144	4140.1156	0.0486	1.0000
SYMG	41	-2051.0694	4188.8438	0.0000	1.0000
SYMI	41	-2051.3800	4189.4650	0.0000	1.0000
SYMIG	42	-2051.0562	4191.0535	0.0000	1.0000
HKY	39	-2054.4727	4191.1961	0.0000	1.0000
GTR	43	-2050.5760	4192.3356	0.0000	1.0000
K2PG	37	-2058.0748	4193.9702	0.0000	1.0000
K2PI	37	-2058.2692	4194.3591	0.0000	1.0000
K2PIG	38	-2058.0748	4196.1824	0.0000	1.0000
SYM	40	-2079.5459	4243.5665	0.0000	1.0000
F81G	39	-2083.0834	4248.4174	0.0000	1.0000
F81I	39	-2083.3975	4249.0456	0.0000	1.0000
F81IG	40	-2083.0408	4250.5563	0.0000	1.0000
K2P	36	-2087.5929	4250.8005	0.0000	1.0000
F81	38	-2110.8782	4301.7890	0.0000	1.0000
JC69G	36	-2118.6173	4312.8493	0.0000	1.0000
JC69I	36	-2118.7816	4313.1779	0.0000	1.0000
JC69IG	37	-2118.6239	4315.0684	0.0000	1.0000
JC69	35	-2145.2045	4363.8235	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2025.8775	4131.7549	0.2820	0.2820
HKYI	40	-2026.1409	4132.2817	0.2167	0.4987
GTRG	44	-2022.2495	4132.4991	0.1944	0.6931
GTRI	44	-2022.6599	4133.3198	0.1289	0.8220
HKYIG	41	-2025.8756	4133.7513	0.1039	0.9259
GTRIG	45	-2022.2144	4134.4288	0.0741	1.0000
SYMG	41	-2051.0694	4184.1388	0.0000	1.0000
SYMI	41	-2051.3800	4184.7601	0.0000	1.0000
SYMIG	42	-2051.0562	4186.1123	0.0000	1.0000
HKY	39	-2054.4727	4186.9454	0.0000	1.0000
GTR	43	-2050.5760	4187.1520	0.0000	1.0000
K2PG	37	-2058.0748	4190.1496	0.0000	1.0000
K2PI	37	-2058.2692	4190.5384	0.0000	1.0000
K2PIG	38	-2058.0748	4192.1497	0.0000	1.0000
SYM	40	-2079.5459	4239.0917	0.0000	1.0000
F81G	39	-2083.0834	4244.1667	0.0000	1.0000
F81I	39	-2083.3975	4244.7949	0.0000	1.0000
F81IG	40	-2083.0408	4246.0816	0.0000	1.0000
K2P	36	-2087.5929	4247.1858	0.0000	1.0000
F81	38	-2110.8782	4297.7564	0.0000	1.0000
JC69G	36	-2118.6173	4309.2347	0.0000	1.0000
JC69I	36	-2118.7816	4309.5632	0.0000	1.0000
JC69IG	37	-2118.6239	4311.2477	0.0000	1.0000
JC69	35	-2145.2045	4360.4089	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2025.8775	4317.8178	0.5542	0.5542
HKYI	40	-2026.1409	4318.3446	0.4258	0.9800
HKYIG	41	-2025.8756	4324.4657	0.0200	0.9999
GTRG	44	-2022.2495	4337.1682	0.0000	1.0000
GTRI	44	-2022.6599	4337.9890	0.0000	1.0000
GTRIG	45	-2022.2144	4343.7495	0.0000	1.0000
K2PG	37	-2058.0748	4362.2577	0.0000	1.0000
K2PI	37	-2058.2692	4362.6466	0.0000	1.0000
HKY	39	-2054.4727	4368.3567	0.0000	1.0000
K2PIG	38	-2058.0748	4368.9094	0.0000	1.0000
SYMG	41	-2051.0694	4374.8533	0.0000	1.0000
SYMI	41	-2051.3800	4375.4745	0.0000	1.0000
SYMIG	42	-2051.0562	4381.4783	0.0000	1.0000
GTR	43	-2050.5760	4387.1696	0.0000	1.0000
K2P	36	-2087.5929	4414.6424	0.0000	1.0000
SYM	40	-2079.5459	4425.1546	0.0000	1.0000
F81G	39	-2083.0834	4425.5780	0.0000	1.0000
F81I	39	-2083.3975	4426.2062	0.0000	1.0000
F81IG	40	-2083.0408	4432.1445	0.0000	1.0000
F81	38	-2110.8782	4474.5161	0.0000	1.0000
JC69G	36	-2118.6173	4476.6912	0.0000	1.0000
JC69I	36	-2118.7816	4477.0198	0.0000	1.0000
JC69IG	37	-2118.6239	4483.3559	0.0000	1.0000
JC69	35	-2145.2045	4523.2139	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/chr3_457.phylip = 1 - 774;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_457.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/chr3_457.phylip = 1 - 774;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_457.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/chr3_457.phylip = 1 - 774;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_457.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

