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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	44	-2050.6801	4195.9712	0.6231	0.6231
GTRIG	45	-2050.3060	4197.5352	0.2851	0.9081
GTRI	44	-2052.5942	4199.7994	0.0919	1.0000
HKYG	40	-2069.8237	4225.0868	0.0000	1.0000
HKYIG	41	-2069.4548	4226.6305	0.0000	1.0000
HKYI	40	-2071.9357	4229.3109	0.0000	1.0000
SYMG	41	-2084.0083	4255.7375	0.0000	1.0000
GTR	43	-2081.8969	4256.1004	0.0000	1.0000
SYMIG	42	-2083.6958	4257.4017	0.0000	1.0000
SYMI	41	-2085.8648	4259.4506	0.0000	1.0000
K2PG	37	-2092.8377	4264.3156	0.0000	1.0000
K2PIG	38	-2092.5374	4265.9740	0.0000	1.0000
K2PI	37	-2095.0822	4268.8047	0.0000	1.0000
HKY	39	-2095.7116	4274.5887	0.0000	1.0000
SYM	40	-2109.2885	4304.0165	0.0000	1.0000
K2P	36	-2116.8276	4310.0440	0.0000	1.0000
F81G	39	-2118.7990	4320.7636	0.0000	1.0000
F81IG	40	-2118.5071	4322.4536	0.0000	1.0000
F81I	39	-2120.8542	4324.8740	0.0000	1.0000
JC69G	36	-2135.6154	4347.6197	0.0000	1.0000
JC69IG	37	-2135.3377	4349.3156	0.0000	1.0000
JC69I	36	-2137.7316	4351.8520	0.0000	1.0000
F81	38	-2142.6779	4366.2549	0.0000	1.0000
JC69	35	-2159.0490	4392.2427	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2050.6801	4189.3602	0.5944	0.5944
GTRIG	45	-2050.3060	4190.6121	0.3179	0.9123
GTRI	44	-2052.5942	4193.1884	0.0877	1.0000
HKYG	40	-2069.8237	4219.6473	0.0000	1.0000
HKYIG	41	-2069.4548	4220.9096	0.0000	1.0000
HKYI	40	-2071.9357	4223.8714	0.0000	1.0000
GTR	43	-2081.8969	4249.7938	0.0000	1.0000
SYMG	41	-2084.0083	4250.0166	0.0000	1.0000
SYMIG	42	-2083.6958	4251.3917	0.0000	1.0000
SYMI	41	-2085.8648	4253.7297	0.0000	1.0000
K2PG	37	-2092.8377	4259.6753	0.0000	1.0000
K2PIG	38	-2092.5374	4261.0748	0.0000	1.0000
K2PI	37	-2095.0822	4264.1644	0.0000	1.0000
HKY	39	-2095.7116	4269.4231	0.0000	1.0000
SYM	40	-2109.2885	4298.5770	0.0000	1.0000
K2P	36	-2116.8276	4305.6552	0.0000	1.0000
F81G	39	-2118.7990	4315.5980	0.0000	1.0000
F81IG	40	-2118.5071	4317.0141	0.0000	1.0000
F81I	39	-2120.8542	4319.7084	0.0000	1.0000
JC69G	36	-2135.6154	4343.2309	0.0000	1.0000
JC69IG	37	-2135.3377	4344.6753	0.0000	1.0000
JC69I	36	-2137.7316	4347.4632	0.0000	1.0000
F81	38	-2142.6779	4361.3558	0.0000	1.0000
JC69	35	-2159.0490	4388.0980	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2050.6801	4385.9389	0.8284	0.8284
GTRI	44	-2052.5942	4389.7671	0.1222	0.9506
GTRIG	45	-2050.3060	4391.6585	0.0475	0.9980
HKYG	40	-2069.8237	4398.3553	0.0017	0.9997
HKYI	40	-2071.9357	4402.5794	0.0002	0.9999
HKYIG	41	-2069.4548	4404.0852	0.0001	1.0000
K2PG	37	-2092.8377	4424.9802	0.0000	1.0000
K2PI	37	-2095.0822	4429.4693	0.0000	1.0000
K2PIG	38	-2092.5374	4430.8474	0.0000	1.0000
SYMG	41	-2084.0083	4433.1922	0.0000	1.0000
SYMI	41	-2085.8648	4436.9053	0.0000	1.0000
SYMIG	42	-2083.6958	4439.0350	0.0000	1.0000
GTR	43	-2081.8969	4441.9048	0.0000	1.0000
HKY	39	-2095.7116	4443.6634	0.0000	1.0000
K2P	36	-2116.8276	4466.4923	0.0000	1.0000
SYM	40	-2109.2885	4477.2850	0.0000	1.0000
F81G	39	-2118.7990	4489.8383	0.0000	1.0000
F81I	39	-2120.8542	4493.9487	0.0000	1.0000
F81IG	40	-2118.5071	4495.7220	0.0000	1.0000
JC69G	36	-2135.6154	4504.0681	0.0000	1.0000
JC69I	36	-2137.7316	4508.3003	0.0000	1.0000
JC69IG	37	-2135.3377	4509.9802	0.0000	1.0000
F81	38	-2142.6779	4531.1283	0.0000	1.0000
JC69	35	-2159.0490	4544.4674	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18665.phylip = 1 - 644;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18665.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18665.phylip = 1 - 644;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18665.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18665.phylip = 1 - 644;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18665.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

