
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:16 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_6769.phylip" (19 taxa, 827 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	-3795.0090	7674.1911	0.6600	0.6600
HKYIG	41	-3794.7078	7675.8028	0.2948	0.9548
GTRG	44	-3793.6122	7680.2884	0.0313	0.9861
GTRIG	45	-3793.3163	7681.9335	0.0138	0.9999
HKYI	40	-3803.5730	7691.3190	0.0001	1.0000
GTRI	44	-3801.7108	7696.4856	0.0000	1.0000
SYMG	41	-3816.3585	7719.1042	0.0000	1.0000
SYMIG	42	-3816.2100	7721.0271	0.0000	1.0000
K2PG	37	-3824.3401	7726.2443	0.0000	1.0000
K2PIG	38	-3824.1825	7728.1265	0.0000	1.0000
SYMI	41	-3825.0113	7736.4098	0.0000	1.0000
K2PI	37	-3832.6176	7742.7993	0.0000	1.0000
F81G	39	-3864.9721	7811.9085	0.0000	1.0000
F81IG	40	-3864.8173	7813.8077	0.0000	1.0000
F81I	39	-3872.6611	7827.2866	0.0000	1.0000
JC69G	36	-3889.0315	7853.4352	0.0000	1.0000
JC69IG	37	-3888.9453	7855.4546	0.0000	1.0000
GTR	43	-3888.2448	7867.3224	0.0000	1.0000
HKY	39	-3892.9101	7867.7847	0.0000	1.0000
JC69I	36	-3896.2182	7867.8086	0.0000	1.0000
SYM	40	-3912.1408	7908.4546	0.0000	1.0000
K2P	36	-3920.4251	7916.2224	0.0000	1.0000
F81	38	-3956.9275	7993.6164	0.0000	1.0000
JC69	35	-3980.4970	8034.1798	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-3795.0090	7670.0181	0.6219	0.6219
HKYIG	41	-3794.7078	7671.4156	0.3092	0.9311
GTRG	44	-3793.6122	7675.2245	0.0460	0.9771
GTRIG	45	-3793.3163	7676.6326	0.0228	0.9999
HKYI	40	-3803.5730	7687.1460	0.0001	1.0000
GTRI	44	-3801.7108	7691.4217	0.0000	1.0000
SYMG	41	-3816.3585	7714.7170	0.0000	1.0000
SYMIG	42	-3816.2100	7716.4199	0.0000	1.0000
K2PG	37	-3824.3401	7722.6803	0.0000	1.0000
K2PIG	38	-3824.1825	7724.3650	0.0000	1.0000
SYMI	41	-3825.0113	7732.0226	0.0000	1.0000
K2PI	37	-3832.6176	7739.2353	0.0000	1.0000
F81G	39	-3864.9721	7807.9441	0.0000	1.0000
F81IG	40	-3864.8173	7809.6346	0.0000	1.0000
F81I	39	-3872.6611	7823.3222	0.0000	1.0000
JC69G	36	-3889.0315	7850.0630	0.0000	1.0000
JC69IG	37	-3888.9453	7851.8906	0.0000	1.0000
GTR	43	-3888.2448	7862.4897	0.0000	1.0000
HKY	39	-3892.9101	7863.8203	0.0000	1.0000
JC69I	36	-3896.2182	7864.4364	0.0000	1.0000
SYM	40	-3912.1408	7904.2816	0.0000	1.0000
K2P	36	-3920.4251	7912.8502	0.0000	1.0000
F81	38	-3956.9275	7989.8550	0.0000	1.0000
JC69	35	-3980.4970	8030.9939	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3795.0090	7858.7303	0.9549	0.9549
HKYIG	41	-3794.7078	7864.8456	0.0449	0.9998
HKYI	40	-3803.5730	7875.8582	0.0002	1.0000
GTRG	44	-3793.6122	7882.8079	0.0000	1.0000
GTRIG	45	-3793.3163	7888.9338	0.0000	1.0000
K2PG	37	-3824.3401	7897.2390	0.0000	1.0000
GTRI	44	-3801.7108	7899.0051	0.0000	1.0000
K2PIG	38	-3824.1825	7903.6416	0.0000	1.0000
SYMG	41	-3816.3585	7908.1470	0.0000	1.0000
K2PI	37	-3832.6176	7913.7941	0.0000	1.0000
SYMIG	42	-3816.2100	7914.5677	0.0000	1.0000
SYMI	41	-3825.0113	7925.4526	0.0000	1.0000
F81G	39	-3864.9721	7991.9385	0.0000	1.0000
F81IG	40	-3864.8173	7998.3468	0.0000	1.0000
F81I	39	-3872.6611	8007.3166	0.0000	1.0000
JC69G	36	-3889.0315	8019.9040	0.0000	1.0000
JC69IG	37	-3888.9453	8026.4494	0.0000	1.0000
JC69I	36	-3896.2182	8034.2774	0.0000	1.0000
HKY	39	-3892.9101	8047.8147	0.0000	1.0000
GTR	43	-3888.2448	8065.3553	0.0000	1.0000
K2P	36	-3920.4251	8082.6912	0.0000	1.0000
SYM	40	-3912.1408	8092.9938	0.0000	1.0000
F81	38	-3956.9275	8169.1316	0.0000	1.0000
JC69	35	-3980.4970	8196.1171	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_6769.phylip = 1 - 827;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6769.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_6769.phylip = 1 - 827;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6769.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_6769.phylip = 1 - 827;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6769.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

