
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:12 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_3515.phylip" (19 taxa, 733 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-3401.3909	6896.5377	0.5377	0.5377
GTRIG	45	-3400.4591	6896.9444	0.4387	0.9764
HKYG	40	-3409.5151	6903.7702	0.0145	0.9908
HKYIG	41	-3408.8471	6904.6783	0.0092	1.0000
GTRI	44	-3415.7786	6925.3130	0.0000	1.0000
HKYI	40	-3421.6880	6928.1159	0.0000	1.0000
K2PIG	38	-3463.9595	7008.1898	0.0000	1.0000
K2PG	37	-3465.2112	7008.4685	0.0000	1.0000
SYMIG	42	-3460.0001	7009.2349	0.0000	1.0000
SYMG	41	-3461.9822	7010.9484	0.0000	1.0000
SYMI	41	-3469.4314	7025.8468	0.0000	1.0000
K2PI	37	-3474.1595	7026.3651	0.0000	1.0000
F81IG	40	-3480.0649	7044.8696	0.0000	1.0000
F81G	39	-3482.0924	7046.6869	0.0000	1.0000
F81I	39	-3490.1698	7062.8418	0.0000	1.0000
GTR	43	-3510.1362	7111.7645	0.0000	1.0000
JC69IG	37	-3520.2040	7118.4539	0.0000	1.0000
HKY	39	-3520.2252	7122.9526	0.0000	1.0000
JC69G	36	-3523.8242	7123.4760	0.0000	1.0000
JC69I	36	-3530.0563	7135.9401	0.0000	1.0000
SYM	40	-3573.4703	7231.6806	0.0000	1.0000
K2P	36	-3579.7394	7235.3064	0.0000	1.0000
F81	38	-3587.5529	7255.3766	0.0000	1.0000
JC69	35	-3630.9792	7335.5739	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3401.3909	6890.7818	0.5097	0.5097
GTRIG	45	-3400.4591	6890.9182	0.4761	0.9858
HKYG	40	-3409.5151	6899.0303	0.0082	0.9941
HKYIG	41	-3408.8471	6899.6942	0.0059	1.0000
GTRI	44	-3415.7786	6919.5572	0.0000	1.0000
HKYI	40	-3421.6880	6923.3760	0.0000	1.0000
K2PIG	38	-3463.9595	7003.9189	0.0000	1.0000
SYMIG	42	-3460.0001	7004.0001	0.0000	1.0000
K2PG	37	-3465.2112	7004.4224	0.0000	1.0000
SYMG	41	-3461.9822	7005.9643	0.0000	1.0000
SYMI	41	-3469.4314	7020.8628	0.0000	1.0000
K2PI	37	-3474.1595	7022.3190	0.0000	1.0000
F81IG	40	-3480.0649	7040.1297	0.0000	1.0000
F81G	39	-3482.0924	7042.1847	0.0000	1.0000
F81I	39	-3490.1698	7058.3397	0.0000	1.0000
GTR	43	-3510.1362	7106.2725	0.0000	1.0000
JC69IG	37	-3520.2040	7114.4079	0.0000	1.0000
HKY	39	-3520.2252	7118.4504	0.0000	1.0000
JC69G	36	-3523.8242	7119.6484	0.0000	1.0000
JC69I	36	-3530.0563	7132.1125	0.0000	1.0000
SYM	40	-3573.4703	7226.9407	0.0000	1.0000
K2P	36	-3579.7394	7231.4788	0.0000	1.0000
F81	38	-3587.5529	7251.1057	0.0000	1.0000
JC69	35	-3630.9792	7331.9584	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3409.5151	7082.9161	0.9269	0.9269
HKYIG	41	-3408.8471	7088.1772	0.0668	0.9936
GTRG	44	-3401.3909	7093.0563	0.0058	0.9994
GTRIG	45	-3400.4591	7097.7898	0.0005	1.0000
HKYI	40	-3421.6880	7107.2618	0.0000	1.0000
GTRI	44	-3415.7786	7121.8316	0.0000	1.0000
K2PG	37	-3465.2112	7174.5168	0.0000	1.0000
K2PIG	38	-3463.9595	7178.6105	0.0000	1.0000
K2PI	37	-3474.1595	7192.4134	0.0000	1.0000
SYMG	41	-3461.9822	7194.4473	0.0000	1.0000
SYMIG	42	-3460.0001	7197.0802	0.0000	1.0000
SYMI	41	-3469.4314	7209.3457	0.0000	1.0000
F81G	39	-3482.0924	7221.4734	0.0000	1.0000
F81IG	40	-3480.0649	7224.0155	0.0000	1.0000
F81I	39	-3490.1698	7237.6284	0.0000	1.0000
JC69IG	37	-3520.2040	7284.5023	0.0000	1.0000
JC69G	36	-3523.8242	7285.1457	0.0000	1.0000
JC69I	36	-3530.0563	7297.6097	0.0000	1.0000
HKY	39	-3520.2252	7297.7391	0.0000	1.0000
GTR	43	-3510.1362	7303.9497	0.0000	1.0000
K2P	36	-3579.7394	7396.9761	0.0000	1.0000
SYM	40	-3573.4703	7410.8265	0.0000	1.0000
F81	38	-3587.5529	7425.7973	0.0000	1.0000
JC69	35	-3630.9792	7492.8585	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/chrZ_3515.phylip = 1 - 733;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_3515.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/chrZ_3515.phylip = 1 - 733;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_3515.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 (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_3515.phylip = 1 - 733;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_3515.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

