
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:24 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/chr1_14908.phylip" (19 taxa, 709 characters)

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: HKYIG


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
HKYIG	41	-2969.0694	6025.3022	0.6450	0.6450
GTRIG	45	-2965.4493	6027.1429	0.2570	0.9020
HKYG	40	-2972.3763	6029.6627	0.0729	0.9749
GTRG	44	-2969.3200	6032.6038	0.0168	0.9916
HKYI	40	-2974.7142	6034.3386	0.0070	0.9987
GTRI	44	-2971.8657	6037.6953	0.0013	1.0000
K2PIG	38	-2992.0719	6064.5677	0.0000	1.0000
K2PG	37	-2995.5470	6069.2847	0.0000	1.0000
SYMIG	42	-2991.7241	6072.8715	0.0000	1.0000
K2PI	37	-2997.5675	6073.3257	0.0000	1.0000
SYMG	41	-2994.9892	6077.1419	0.0000	1.0000
SYMI	41	-2997.4142	6081.9917	0.0000	1.0000
F81IG	40	-3029.8193	6144.5488	0.0000	1.0000
F81G	39	-3033.3141	6149.2918	0.0000	1.0000
F81I	39	-3034.5776	6151.8189	0.0000	1.0000
JC69IG	37	-3046.5569	6171.3045	0.0000	1.0000
JC69G	36	-3049.6748	6175.3140	0.0000	1.0000
JC69I	36	-3049.8182	6175.6007	0.0000	1.0000
HKY	39	-3137.9333	6358.5303	0.0000	1.0000
GTR	43	-3134.0481	6359.7864	0.0000	1.0000
K2P	36	-3159.3509	6394.6661	0.0000	1.0000
SYM	40	-3157.9776	6400.8654	0.0000	1.0000
F81	38	-3192.0069	6464.4378	0.0000	1.0000
JC69	35	-3208.0734	6489.8911	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-2969.0694	6020.1388	0.5450	0.5450
GTRIG	45	-2965.4493	6020.8986	0.3727	0.9177
HKYG	40	-2972.3763	6024.7525	0.0543	0.9720
GTRG	44	-2969.3200	6026.6400	0.0211	0.9931
HKYI	40	-2974.7142	6029.4284	0.0052	0.9983
GTRI	44	-2971.8657	6031.7314	0.0017	1.0000
K2PIG	38	-2992.0719	6060.1439	0.0000	1.0000
K2PG	37	-2995.5470	6065.0940	0.0000	1.0000
SYMIG	42	-2991.7241	6067.4481	0.0000	1.0000
K2PI	37	-2997.5675	6069.1350	0.0000	1.0000
SYMG	41	-2994.9892	6071.9784	0.0000	1.0000
SYMI	41	-2997.4142	6076.8283	0.0000	1.0000
F81IG	40	-3029.8193	6139.6386	0.0000	1.0000
F81G	39	-3033.3141	6144.6281	0.0000	1.0000
F81I	39	-3034.5776	6147.1552	0.0000	1.0000
JC69IG	37	-3046.5569	6167.1138	0.0000	1.0000
JC69G	36	-3049.6748	6171.3497	0.0000	1.0000
JC69I	36	-3049.8182	6171.6364	0.0000	1.0000
HKY	39	-3137.9333	6353.8666	0.0000	1.0000
GTR	43	-3134.0481	6354.0962	0.0000	1.0000
K2P	36	-3159.3509	6390.7018	0.0000	1.0000
SYM	40	-3157.9776	6395.9552	0.0000	1.0000
F81	38	-3192.0069	6460.0139	0.0000	1.0000
JC69	35	-3208.0734	6486.1467	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	41	-2969.0694	6207.2569	0.4832	0.4832
HKYG	40	-2972.3763	6207.3067	0.4713	0.9544
HKYI	40	-2974.7142	6211.9826	0.0455	0.9999
GTRIG	45	-2965.4493	6226.2721	0.0000	1.0000
GTRG	44	-2969.3200	6227.4496	0.0000	1.0000
GTRI	44	-2971.8657	6232.5411	0.0000	1.0000
K2PIG	38	-2992.0719	6233.5704	0.0000	1.0000
K2PG	37	-2995.5470	6233.9566	0.0000	1.0000
K2PI	37	-2997.5675	6237.9976	0.0000	1.0000
SYMG	41	-2994.9892	6259.0965	0.0000	1.0000
SYMIG	42	-2991.7241	6259.1301	0.0000	1.0000
SYMI	41	-2997.4142	6263.9464	0.0000	1.0000
F81IG	40	-3029.8193	6322.1928	0.0000	1.0000
F81G	39	-3033.3141	6322.6185	0.0000	1.0000
F81I	39	-3034.5776	6325.1456	0.0000	1.0000
JC69G	36	-3049.6748	6335.6485	0.0000	1.0000
JC69I	36	-3049.8182	6335.9352	0.0000	1.0000
JC69IG	37	-3046.5569	6335.9764	0.0000	1.0000
HKY	39	-3137.9333	6531.8570	0.0000	1.0000
GTR	43	-3134.0481	6550.3420	0.0000	1.0000
K2P	36	-3159.3509	6555.0006	0.0000	1.0000
SYM	40	-3157.9776	6578.5095	0.0000	1.0000
F81	38	-3192.0069	6633.4404	0.0000	1.0000
JC69	35	-3208.0734	6645.8817	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_14908.phylip = 1 - 709;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_14908.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_14908.phylip = 1 - 709;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_14908.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_14908.phylip = 1 - 709;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_14908.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

