
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:55 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/chr7_10430.phylip" (19 taxa, 928 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	-2190.7506	4465.1990	0.4575	0.4575
HKYI	40	-2191.3621	4466.4220	0.2482	0.7058
GTRG	44	-2187.8820	4468.2488	0.0996	0.8053
GTRIG	45	-2186.9948	4468.6835	0.0801	0.8855
HKYIG	41	-2191.5203	4468.9277	0.0709	0.9564
GTRI	44	-2188.7077	4469.9001	0.0436	1.0000
HKY	39	-2212.4983	4506.5101	0.0000	1.0000
GTR	43	-2208.8512	4507.9830	0.0000	1.0000
F81G	39	-2220.2345	4521.9825	0.0000	1.0000
F81IG	40	-2219.6399	4522.9777	0.0000	1.0000
F81I	39	-2220.8060	4523.1256	0.0000	1.0000
F81	38	-2240.0215	4559.3772	0.0000	1.0000
K2PG	37	-2252.9384	4583.0364	0.0000	1.0000
K2PIG	38	-2252.3034	4583.9409	0.0000	1.0000
SYMG	41	-2249.5261	4584.9393	0.0000	1.0000
SYMI	41	-2250.2392	4586.3655	0.0000	1.0000
K2PI	37	-2255.8675	4588.8945	0.0000	1.0000
SYMIG	42	-2251.2660	4590.6134	0.0000	1.0000
K2P	36	-2272.5877	4620.1653	0.0000	1.0000
SYM	40	-2268.2845	4620.2668	0.0000	1.0000
JC69G	36	-2275.5928	4626.1755	0.0000	1.0000
JC69IG	37	-2275.0777	4627.3151	0.0000	1.0000
JC69I	36	-2278.5493	4632.0886	0.0000	1.0000
JC69	35	-2294.4552	4661.7355	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2190.7506	4461.5011	0.4057	0.4057
HKYI	40	-2191.3621	4462.7241	0.2201	0.6258
GTRG	44	-2187.8820	4463.7641	0.1309	0.7567
GTRIG	45	-2186.9948	4463.9897	0.1169	0.8736
HKYIG	41	-2191.5203	4465.0405	0.0691	0.9427
GTRI	44	-2188.7077	4465.4154	0.0573	1.0000
HKY	39	-2212.4983	4502.9965	0.0000	1.0000
GTR	43	-2208.8512	4503.7024	0.0000	1.0000
F81G	39	-2220.2345	4518.4690	0.0000	1.0000
F81IG	40	-2219.6399	4519.2798	0.0000	1.0000
F81I	39	-2220.8060	4519.6121	0.0000	1.0000
F81	38	-2240.0215	4556.0431	0.0000	1.0000
K2PG	37	-2252.9384	4579.8768	0.0000	1.0000
K2PIG	38	-2252.3034	4580.6068	0.0000	1.0000
SYMG	41	-2249.5261	4581.0522	0.0000	1.0000
SYMI	41	-2250.2392	4582.4784	0.0000	1.0000
K2PI	37	-2255.8675	4585.7350	0.0000	1.0000
SYMIG	42	-2251.2660	4586.5320	0.0000	1.0000
SYM	40	-2268.2845	4616.5689	0.0000	1.0000
K2P	36	-2272.5877	4617.1754	0.0000	1.0000
JC69G	36	-2275.5928	4623.1856	0.0000	1.0000
JC69IG	37	-2275.0777	4624.1555	0.0000	1.0000
JC69I	36	-2278.5493	4629.0987	0.0000	1.0000
JC69	35	-2294.4552	4658.9104	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2190.7506	4654.8224	0.6419	0.6419
HKYI	40	-2191.3621	4656.0454	0.3483	0.9902
HKYIG	41	-2191.5203	4663.1948	0.0098	1.0000
GTRG	44	-2187.8820	4676.4175	0.0000	1.0000
GTRI	44	-2188.7077	4678.0688	0.0000	1.0000
GTRIG	45	-2186.9948	4681.4761	0.0000	1.0000
HKY	39	-2212.4983	4691.4848	0.0000	1.0000
F81G	39	-2220.2345	4706.9573	0.0000	1.0000
F81I	39	-2220.8060	4708.1003	0.0000	1.0000
GTR	43	-2208.8512	4711.5228	0.0000	1.0000
F81IG	40	-2219.6399	4712.6011	0.0000	1.0000
F81	38	-2240.0215	4739.6983	0.0000	1.0000
K2PG	37	-2252.9384	4758.6990	0.0000	1.0000
K2PIG	38	-2252.3034	4764.2620	0.0000	1.0000
K2PI	37	-2255.8675	4764.5571	0.0000	1.0000
SYMG	41	-2249.5261	4779.2065	0.0000	1.0000
SYMI	41	-2250.2392	4780.6327	0.0000	1.0000
SYMIG	42	-2251.2660	4789.5193	0.0000	1.0000
K2P	36	-2272.5877	4791.1645	0.0000	1.0000
JC69G	36	-2275.5928	4797.1747	0.0000	1.0000
JC69IG	37	-2275.0777	4802.9777	0.0000	1.0000
JC69I	36	-2278.5493	4803.0878	0.0000	1.0000
SYM	40	-2268.2845	4809.8902	0.0000	1.0000
JC69	35	-2294.4552	4828.0665	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/chr7_10430.phylip = 1 - 928;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10430.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/chr7_10430.phylip = 1 - 928;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10430.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/chr7_10430.phylip = 1 - 928;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10430.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

