
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:26 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/chrZ_7924.phylip" (19 taxa, 694 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	-2013.0491	4111.1212	0.4295	0.4295
HKYIG	41	-2012.2660	4111.8142	0.3037	0.7332
GTRG	44	-2009.6257	4113.3532	0.1407	0.8738
GTRIG	45	-2008.9183	4114.2254	0.0910	0.9648
HKYI	40	-2015.7919	4116.6067	0.0277	0.9925
GTRI	44	-2012.5518	4119.2054	0.0075	1.0000
K2PIG	38	-2027.0600	4134.6451	0.0000	1.0000
K2PG	37	-2028.2901	4134.8668	0.0000	1.0000
K2PI	37	-2030.7070	4139.7007	0.0000	1.0000
SYMG	41	-2026.7338	4140.7498	0.0000	1.0000
SYMIG	42	-2025.9342	4141.4167	0.0000	1.0000
SYMI	41	-2029.6389	4146.5601	0.0000	1.0000
HKY	39	-2037.0171	4156.8048	0.0000	1.0000
GTR	43	-2032.9813	4157.7841	0.0000	1.0000
K2P	36	-2051.8439	4179.7426	0.0000	1.0000
SYM	40	-2050.1300	4185.2830	0.0000	1.0000
F81G	39	-2056.3646	4195.4999	0.0000	1.0000
F81IG	40	-2055.7287	4196.4805	0.0000	1.0000
F81I	39	-2058.6769	4200.1244	0.0000	1.0000
JC69G	36	-2070.4211	4216.8970	0.0000	1.0000
JC69IG	37	-2069.7715	4217.8296	0.0000	1.0000
JC69I	36	-2072.6948	4221.4443	0.0000	1.0000
F81	38	-2079.2489	4239.0230	0.0000	1.0000
JC69	35	-2093.4834	4260.7966	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2013.0491	4106.0982	0.3471	0.3471
HKYIG	41	-2012.2660	4106.5320	0.2795	0.6266
GTRG	44	-2009.6257	4107.2515	0.1950	0.8216
GTRIG	45	-2008.9183	4107.8365	0.1456	0.9672
HKYI	40	-2015.7919	4111.5837	0.0224	0.9895
GTRI	44	-2012.5518	4113.1037	0.0105	1.0000
K2PIG	38	-2027.0600	4130.1199	0.0000	1.0000
K2PG	37	-2028.2901	4130.5803	0.0000	1.0000
K2PI	37	-2030.7070	4135.4141	0.0000	1.0000
SYMG	41	-2026.7338	4135.4676	0.0000	1.0000
SYMIG	42	-2025.9342	4135.8683	0.0000	1.0000
SYMI	41	-2029.6389	4141.2779	0.0000	1.0000
GTR	43	-2032.9813	4151.9625	0.0000	1.0000
HKY	39	-2037.0171	4152.0341	0.0000	1.0000
K2P	36	-2051.8439	4175.6878	0.0000	1.0000
SYM	40	-2050.1300	4180.2600	0.0000	1.0000
F81G	39	-2056.3646	4190.7293	0.0000	1.0000
F81IG	40	-2055.7287	4191.4575	0.0000	1.0000
F81I	39	-2058.6769	4195.3537	0.0000	1.0000
JC69G	36	-2070.4211	4212.8422	0.0000	1.0000
JC69IG	37	-2069.7715	4213.5431	0.0000	1.0000
JC69I	36	-2072.6948	4217.3895	0.0000	1.0000
F81	38	-2079.2489	4234.4978	0.0000	1.0000
JC69	35	-2093.4834	4256.9669	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2013.0491	4287.7971	0.8674	0.8674
HKYIG	41	-2012.2660	4292.7734	0.0720	0.9394
HKYI	40	-2015.7919	4293.2826	0.0559	0.9953
K2PG	37	-2028.2901	4298.6517	0.0038	0.9991
K2PIG	38	-2027.0600	4302.7339	0.0005	0.9996
K2PI	37	-2030.7070	4303.4855	0.0003	0.9999
GTRG	44	-2009.6257	4307.1202	0.0001	1.0000
GTRIG	45	-2008.9183	4312.2478	0.0000	1.0000
GTRI	44	-2012.5518	4312.9724	0.0000	1.0000
SYMG	41	-2026.7338	4321.7090	0.0000	1.0000
SYMIG	42	-2025.9342	4326.6521	0.0000	1.0000
SYMI	41	-2029.6389	4327.5192	0.0000	1.0000
HKY	39	-2037.0171	4329.1905	0.0000	1.0000
K2P	36	-2051.8439	4339.2168	0.0000	1.0000
GTR	43	-2032.9813	4347.2888	0.0000	1.0000
SYM	40	-2050.1300	4361.9589	0.0000	1.0000
F81G	39	-2056.3646	4367.8857	0.0000	1.0000
F81I	39	-2058.6769	4372.5101	0.0000	1.0000
F81IG	40	-2055.7287	4373.1564	0.0000	1.0000
JC69G	36	-2070.4211	4376.3712	0.0000	1.0000
JC69I	36	-2072.6948	4380.9185	0.0000	1.0000
JC69IG	37	-2069.7715	4381.6145	0.0000	1.0000
F81	38	-2079.2489	4407.1117	0.0000	1.0000
JC69	35	-2093.4834	4415.9534	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/chrZ_7924.phylip = 1 - 694;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7924.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/chrZ_7924.phylip = 1 - 694;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7924.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/chrZ_7924.phylip = 1 - 694;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7924.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

