
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:22 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_7827.phylip" (19 taxa, 739 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	-3083.2130	6260.1321	0.6903	0.6903
GTRIG	45	-3082.8804	6261.7349	0.3097	1.0000
GTRI	44	-3094.8315	6283.3690	0.0000	1.0000
HKYG	40	-3100.0683	6284.8357	0.0000	1.0000
HKYIG	41	-3099.7528	6286.4469	0.0000	1.0000
HKYI	40	-3109.7944	6304.2880	0.0000	1.0000
SYMG	41	-3109.8029	6306.5470	0.0000	1.0000
SYMIG	42	-3109.4730	6308.1356	0.0000	1.0000
K2PG	37	-3116.7988	6311.6090	0.0000	1.0000
K2PIG	38	-3116.4630	6313.1603	0.0000	1.0000
SYMI	41	-3119.6925	6326.3261	0.0000	1.0000
K2PI	37	-3126.5225	6331.0564	0.0000	1.0000
GTR	43	-3150.7527	6392.9500	0.0000	1.0000
HKY	39	-3166.3070	6415.0775	0.0000	1.0000
SYM	40	-3176.8509	6438.4010	0.0000	1.0000
K2P	36	-3185.1044	6446.0038	0.0000	1.0000
F81IG	40	-3195.8447	6476.3885	0.0000	1.0000
F81G	39	-3198.4637	6479.3908	0.0000	1.0000
F81I	39	-3204.5741	6491.6117	0.0000	1.0000
JC69IG	37	-3208.1944	6494.4002	0.0000	1.0000
JC69G	36	-3210.3149	6496.4246	0.0000	1.0000
JC69I	36	-3217.0036	6509.8020	0.0000	1.0000
F81	38	-3259.1722	6598.5786	0.0000	1.0000
JC69	35	-3272.0755	6617.7357	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3083.2130	6254.4261	0.6609	0.6609
GTRIG	45	-3082.8804	6255.7608	0.3391	1.0000
GTRI	44	-3094.8315	6277.6630	0.0000	1.0000
HKYG	40	-3100.0683	6280.1366	0.0000	1.0000
HKYIG	41	-3099.7528	6281.5057	0.0000	1.0000
HKYI	40	-3109.7944	6299.5888	0.0000	1.0000
SYMG	41	-3109.8029	6301.6058	0.0000	1.0000
SYMIG	42	-3109.4730	6302.9460	0.0000	1.0000
K2PG	37	-3116.7988	6307.5976	0.0000	1.0000
K2PIG	38	-3116.4630	6308.9260	0.0000	1.0000
SYMI	41	-3119.6925	6321.3849	0.0000	1.0000
K2PI	37	-3126.5225	6327.0450	0.0000	1.0000
GTR	43	-3150.7527	6387.5054	0.0000	1.0000
HKY	39	-3166.3070	6410.6139	0.0000	1.0000
SYM	40	-3176.8509	6433.7019	0.0000	1.0000
K2P	36	-3185.1044	6442.2089	0.0000	1.0000
F81IG	40	-3195.8447	6471.6894	0.0000	1.0000
F81G	39	-3198.4637	6474.9273	0.0000	1.0000
F81I	39	-3204.5741	6487.1482	0.0000	1.0000
JC69IG	37	-3208.1944	6490.3888	0.0000	1.0000
JC69G	36	-3210.3149	6492.6297	0.0000	1.0000
JC69I	36	-3217.0036	6506.0071	0.0000	1.0000
F81	38	-3259.1722	6594.3443	0.0000	1.0000
JC69	35	-3272.0755	6614.1511	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3083.2130	6457.0592	0.9270	0.9270
GTRIG	45	-3082.8804	6462.9992	0.0476	0.9745
HKYG	40	-3100.0683	6464.3485	0.0242	0.9987
HKYIG	41	-3099.7528	6470.3229	0.0012	1.0000
K2PG	37	-3116.7988	6477.9936	0.0000	1.0000
GTRI	44	-3094.8315	6480.2961	0.0000	1.0000
HKYI	40	-3109.7944	6483.8007	0.0000	1.0000
K2PIG	38	-3116.4630	6483.9273	0.0000	1.0000
SYMG	41	-3109.8029	6490.4230	0.0000	1.0000
SYMIG	42	-3109.4730	6496.3685	0.0000	1.0000
K2PI	37	-3126.5225	6497.4410	0.0000	1.0000
SYMI	41	-3119.6925	6510.2021	0.0000	1.0000
GTR	43	-3150.7527	6585.5332	0.0000	1.0000
HKY	39	-3166.3070	6590.2206	0.0000	1.0000
K2P	36	-3185.1044	6607.9996	0.0000	1.0000
SYM	40	-3176.8509	6617.9138	0.0000	1.0000
F81G	39	-3198.4637	6654.5339	0.0000	1.0000
F81IG	40	-3195.8447	6655.9013	0.0000	1.0000
JC69G	36	-3210.3149	6658.4205	0.0000	1.0000
JC69IG	37	-3208.1944	6660.7848	0.0000	1.0000
F81I	39	-3204.5741	6666.7548	0.0000	1.0000
JC69I	36	-3217.0036	6671.7979	0.0000	1.0000
F81	38	-3259.1722	6769.3456	0.0000	1.0000
JC69	35	-3272.0755	6775.3365	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_7827.phylip = 1 - 739;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7827.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_7827.phylip = 1 - 739;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7827.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/chrZ_7827.phylip = 1 - 739;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7827.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

