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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: SYMG


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	-2941.8812	5979.6508	0.8908	0.8908
GTRIG	45	-2943.7881	5985.8396	0.0404	0.9312
SYMIG	42	-2947.4295	5986.0257	0.0368	0.9680
SYMG	41	-2948.7424	5986.3045	0.0320	1.0000
HKYG	40	-2956.6634	5999.8090	0.0000	1.0000
HKYIG	41	-2956.6650	6002.1498	0.0000	1.0000
K2PG	37	-2965.7420	6011.0086	0.0000	1.0000
K2PIG	38	-2965.7419	6013.3184	0.0000	1.0000
GTRI	44	-2961.0952	6018.0789	0.0000	1.0000
SYMI	41	-2964.7367	6018.2932	0.0000	1.0000
HKYI	40	-2972.7064	6031.8950	0.0000	1.0000
K2PI	37	-2980.6997	6040.9240	0.0000	1.0000
F81IG	40	-3044.6587	6175.7996	0.0000	1.0000
F81G	39	-3051.0587	6186.2712	0.0000	1.0000
JC69G	36	-3060.2634	6197.7502	0.0000	1.0000
JC69IG	37	-3060.2671	6200.0588	0.0000	1.0000
F81I	39	-3060.9239	6206.0016	0.0000	1.0000
JC69I	36	-3074.8130	6226.8494	0.0000	1.0000
SYM	40	-3080.8345	6248.1513	0.0000	1.0000
GTR	43	-3077.4109	6248.3447	0.0000	1.0000
HKY	39	-3091.8967	6267.9472	0.0000	1.0000
K2P	36	-3101.3796	6279.9826	0.0000	1.0000
F81	38	-3173.5009	6428.8365	0.0000	1.0000
JC69	35	-3181.9223	6438.7760	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2941.8812	5971.7624	0.9054	0.9054
GTRIG	45	-2943.7881	5977.5762	0.0495	0.9549
SYMIG	42	-2947.4295	5978.8590	0.0261	0.9809
SYMG	41	-2948.7424	5979.4847	0.0191	1.0000
HKYG	40	-2956.6634	5993.3268	0.0000	1.0000
HKYIG	41	-2956.6650	5995.3300	0.0000	1.0000
K2PG	37	-2965.7420	6005.4841	0.0000	1.0000
K2PIG	38	-2965.7419	6007.4837	0.0000	1.0000
GTRI	44	-2961.0952	6010.1905	0.0000	1.0000
SYMI	41	-2964.7367	6011.4734	0.0000	1.0000
HKYI	40	-2972.7064	6025.4128	0.0000	1.0000
K2PI	37	-2980.6997	6035.3994	0.0000	1.0000
F81IG	40	-3044.6587	6169.3174	0.0000	1.0000
F81G	39	-3051.0587	6180.1173	0.0000	1.0000
JC69G	36	-3060.2634	6192.5267	0.0000	1.0000
JC69IG	37	-3060.2671	6194.5343	0.0000	1.0000
F81I	39	-3060.9239	6199.8477	0.0000	1.0000
JC69I	36	-3074.8130	6221.6259	0.0000	1.0000
GTR	43	-3077.4109	6240.8218	0.0000	1.0000
SYM	40	-3080.8345	6241.6691	0.0000	1.0000
HKY	39	-3091.8967	6261.7934	0.0000	1.0000
K2P	36	-3101.3796	6274.7591	0.0000	1.0000
F81	38	-3173.5009	6423.0018	0.0000	1.0000
JC69	35	-3181.9223	6433.8445	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	41	-2948.7424	6155.9671	0.7963	0.7963
SYMIG	42	-2947.4295	6159.6459	0.1265	0.9229
GTRG	44	-2941.8812	6161.1581	0.0594	0.9823
K2PG	37	-2965.7420	6164.7487	0.0099	0.9922
HKYG	40	-2956.6634	6165.5048	0.0068	0.9989
K2PIG	38	-2965.7419	6171.0528	0.0004	0.9993
GTRIG	45	-2943.7881	6171.2764	0.0004	0.9997
HKYIG	41	-2956.6650	6171.8124	0.0003	1.0000
SYMI	41	-2964.7367	6187.9558	0.0000	1.0000
K2PI	37	-2980.6997	6194.6640	0.0000	1.0000
HKYI	40	-2972.7064	6197.5908	0.0000	1.0000
GTRI	44	-2961.0952	6199.5862	0.0000	1.0000
F81IG	40	-3044.6587	6341.4954	0.0000	1.0000
JC69G	36	-3060.2634	6347.4869	0.0000	1.0000
F81G	39	-3051.0587	6347.9908	0.0000	1.0000
JC69IG	37	-3060.2671	6353.7989	0.0000	1.0000
F81I	39	-3060.9239	6367.7212	0.0000	1.0000
JC69I	36	-3074.8130	6376.5861	0.0000	1.0000
SYM	40	-3080.8345	6413.8470	0.0000	1.0000
GTR	43	-3077.4109	6425.9131	0.0000	1.0000
HKY	39	-3091.8967	6429.6669	0.0000	1.0000
K2P	36	-3101.3796	6429.7193	0.0000	1.0000
JC69	35	-3181.9223	6584.5002	0.0000	1.0000
F81	38	-3173.5009	6586.5709	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/chr7_1396.phylip = 1 - 547;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_1396.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/chr7_1396.phylip = 1 - 547;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_1396.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 (SYMG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_1396.phylip = 1 - 547;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_1396.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=Fixed(Equal);
END;


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

