
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:58 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/chr5_10794.phylip" (19 taxa, 709 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	-2528.0750	5141.0602	0.5247	0.5247
HKYIG	41	-2527.8904	5142.9442	0.2046	0.7293
GTRG	44	-2525.0049	5143.9737	0.1223	0.8516
HKYI	40	-2529.9978	5144.9058	0.0767	0.9283
GTRIG	45	-2524.7737	5145.7917	0.0493	0.9775
GTRI	44	-2526.6988	5147.3615	0.0225	1.0000
K2PG	37	-2558.6132	5195.4172	0.0000	1.0000
K2PIG	38	-2558.4541	5197.3320	0.0000	1.0000
SYMG	41	-2555.7027	5198.5688	0.0000	1.0000
K2PI	37	-2560.7641	5199.7190	0.0000	1.0000
SYMIG	42	-2555.5857	5200.5948	0.0000	1.0000
SYMI	41	-2557.6666	5202.4967	0.0000	1.0000
HKY	39	-2593.6232	5269.9101	0.0000	1.0000
F81G	39	-2594.8186	5272.3010	0.0000	1.0000
F81IG	40	-2594.6667	5274.2436	0.0000	1.0000
GTR	43	-2591.3835	5274.4573	0.0000	1.0000
F81I	39	-2596.0699	5274.8034	0.0000	1.0000
JC69G	36	-2618.6749	5313.3140	0.0000	1.0000
JC69IG	37	-2618.5438	5315.2784	0.0000	1.0000
JC69I	36	-2620.0141	5315.9925	0.0000	1.0000
K2P	36	-2622.7503	5321.4648	0.0000	1.0000
SYM	40	-2620.6896	5326.2894	0.0000	1.0000
F81	38	-2653.5397	5387.5033	0.0000	1.0000
JC69	35	-2676.8302	5427.4048	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2528.0750	5136.1500	0.4467	0.4467
HKYIG	41	-2527.8904	5137.7808	0.1976	0.6443
GTRG	44	-2525.0049	5138.0099	0.1763	0.8206
GTRIG	45	-2524.7737	5139.5473	0.0817	0.9023
HKYI	40	-2529.9978	5139.9956	0.0653	0.9676
GTRI	44	-2526.6988	5141.3977	0.0324	1.0000
K2PG	37	-2558.6132	5191.2265	0.0000	1.0000
K2PIG	38	-2558.4541	5192.9081	0.0000	1.0000
SYMG	41	-2555.7027	5193.4054	0.0000	1.0000
SYMIG	42	-2555.5857	5195.1714	0.0000	1.0000
K2PI	37	-2560.7641	5195.5283	0.0000	1.0000
SYMI	41	-2557.6666	5197.3333	0.0000	1.0000
HKY	39	-2593.6232	5265.2464	0.0000	1.0000
F81G	39	-2594.8186	5267.6373	0.0000	1.0000
GTR	43	-2591.3835	5268.7670	0.0000	1.0000
F81IG	40	-2594.6667	5269.3334	0.0000	1.0000
F81I	39	-2596.0699	5270.1397	0.0000	1.0000
JC69G	36	-2618.6749	5309.3497	0.0000	1.0000
JC69IG	37	-2618.5438	5311.0877	0.0000	1.0000
JC69I	36	-2620.0141	5312.0282	0.0000	1.0000
K2P	36	-2622.7503	5317.5005	0.0000	1.0000
SYM	40	-2620.6896	5321.3793	0.0000	1.0000
F81	38	-2653.5397	5383.0794	0.0000	1.0000
JC69	35	-2676.8302	5423.6603	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2528.0750	5318.7042	0.8393	0.8393
HKYI	40	-2529.9978	5322.5498	0.1227	0.9620
HKYIG	41	-2527.8904	5324.8989	0.0379	1.0000
GTRG	44	-2525.0049	5338.8195	0.0000	1.0000
GTRI	44	-2526.6988	5342.2073	0.0000	1.0000
GTRIG	45	-2524.7737	5344.9208	0.0000	1.0000
K2PG	37	-2558.6132	5360.0891	0.0000	1.0000
K2PI	37	-2560.7641	5364.3909	0.0000	1.0000
K2PIG	38	-2558.4541	5366.3346	0.0000	1.0000
SYMG	41	-2555.7027	5380.5235	0.0000	1.0000
SYMI	41	-2557.6666	5384.4513	0.0000	1.0000
SYMIG	42	-2555.5857	5386.8533	0.0000	1.0000
HKY	39	-2593.6232	5443.2367	0.0000	1.0000
F81G	39	-2594.8186	5445.6276	0.0000	1.0000
F81I	39	-2596.0699	5448.1301	0.0000	1.0000
F81IG	40	-2594.6667	5451.8876	0.0000	1.0000
GTR	43	-2591.3835	5465.0128	0.0000	1.0000
JC69G	36	-2618.6749	5473.6485	0.0000	1.0000
JC69I	36	-2620.0141	5476.3270	0.0000	1.0000
JC69IG	37	-2618.5438	5479.9503	0.0000	1.0000
K2P	36	-2622.7503	5481.7993	0.0000	1.0000
SYM	40	-2620.6896	5503.9335	0.0000	1.0000
F81	38	-2653.5397	5556.5059	0.0000	1.0000
JC69	35	-2676.8302	5583.3953	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/chr5_10794.phylip = 1 - 709;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10794.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/chr5_10794.phylip = 1 - 709;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10794.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/chr5_10794.phylip = 1 - 709;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10794.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

