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

Minimum AIC  model: GTRIG
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	-2949.3628	5992.3267	0.5110	0.5110
GTRIG	45	-2948.2755	5992.4150	0.4890	1.0000
GTRI	44	-2962.1655	6017.9321	0.0000	1.0000
HKYG	40	-2966.7016	6018.0164	0.0000	1.0000
HKYIG	41	-2965.6552	6018.1612	0.0000	1.0000
SYMG	41	-2974.8956	6036.6418	0.0000	1.0000
SYMIG	42	-2973.8950	6036.8845	0.0000	1.0000
HKYI	40	-2976.1692	6036.9517	0.0000	1.0000
K2PG	37	-2982.6286	6043.1956	0.0000	1.0000
K2PIG	38	-2981.6060	6043.3691	0.0000	1.0000
SYMI	41	-2984.1474	6055.1455	0.0000	1.0000
K2PI	37	-2991.1049	6060.1482	0.0000	1.0000
F81G	39	-3037.3533	6157.0886	0.0000	1.0000
F81IG	40	-3036.4636	6157.5404	0.0000	1.0000
GTR	43	-3033.5766	6158.4978	0.0000	1.0000
F81I	39	-3045.0559	6172.4938	0.0000	1.0000
JC69G	36	-3049.6950	6175.1158	0.0000	1.0000
JC69IG	37	-3048.9255	6175.7894	0.0000	1.0000
HKY	39	-3046.9385	6176.2590	0.0000	1.0000
JC69I	36	-3056.6871	6189.1000	0.0000	1.0000
SYM	40	-3055.8062	6196.2256	0.0000	1.0000
K2P	36	-3062.8158	6201.3575	0.0000	1.0000
F81	38	-3112.5207	6305.1984	0.0000	1.0000
JC69	35	-3124.2267	6321.9730	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-2948.2755	5986.5510	0.5218	0.5218
GTRG	44	-2949.3628	5986.7256	0.4782	1.0000
GTRI	44	-2962.1655	6012.3309	0.0000	1.0000
HKYIG	41	-2965.6552	6013.3105	0.0000	1.0000
HKYG	40	-2966.7016	6013.4032	0.0000	1.0000
SYMIG	42	-2973.8950	6031.7900	0.0000	1.0000
SYMG	41	-2974.8956	6031.7911	0.0000	1.0000
HKYI	40	-2976.1692	6032.3385	0.0000	1.0000
K2PIG	38	-2981.6060	6039.2120	0.0000	1.0000
K2PG	37	-2982.6286	6039.2572	0.0000	1.0000
SYMI	41	-2984.1474	6050.2948	0.0000	1.0000
K2PI	37	-2991.1049	6056.2099	0.0000	1.0000
F81G	39	-3037.3533	6152.7066	0.0000	1.0000
F81IG	40	-3036.4636	6152.9272	0.0000	1.0000
GTR	43	-3033.5766	6153.1531	0.0000	1.0000
F81I	39	-3045.0559	6168.1118	0.0000	1.0000
JC69G	36	-3049.6950	6171.3899	0.0000	1.0000
JC69IG	37	-3048.9255	6171.8510	0.0000	1.0000
HKY	39	-3046.9385	6171.8770	0.0000	1.0000
JC69I	36	-3056.6871	6185.3741	0.0000	1.0000
SYM	40	-3055.8062	6191.6124	0.0000	1.0000
K2P	36	-3062.8158	6197.6316	0.0000	1.0000
F81	38	-3112.5207	6301.0413	0.0000	1.0000
JC69	35	-3124.2267	6318.4535	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2949.3628	6190.1260	0.8876	0.8876
GTRIG	45	-2948.2755	6194.5741	0.0960	0.9836
HKYG	40	-2966.7016	6198.3126	0.0148	0.9984
HKYIG	41	-2965.6552	6202.8426	0.0015	1.0000
K2PG	37	-2982.6286	6210.2985	0.0000	1.0000
K2PIG	38	-2981.6060	6214.8760	0.0000	1.0000
GTRI	44	-2962.1655	6215.7313	0.0000	1.0000
HKYI	40	-2976.1692	6217.2480	0.0000	1.0000
SYMG	41	-2974.8956	6221.3233	0.0000	1.0000
SYMIG	42	-2973.8950	6225.9449	0.0000	1.0000
K2PI	37	-2991.1049	6227.2511	0.0000	1.0000
SYMI	41	-2984.1474	6239.8270	0.0000	1.0000
F81G	39	-3037.3533	6332.9933	0.0000	1.0000
JC69G	36	-3049.6950	6337.8084	0.0000	1.0000
F81IG	40	-3036.4636	6337.8366	0.0000	1.0000
JC69IG	37	-3048.9255	6342.8922	0.0000	1.0000
F81I	39	-3045.0559	6348.3985	0.0000	1.0000
JC69I	36	-3056.6871	6351.7926	0.0000	1.0000
GTR	43	-3033.5766	6351.9308	0.0000	1.0000
HKY	39	-3046.9385	6352.1637	0.0000	1.0000
K2P	36	-3062.8158	6364.0501	0.0000	1.0000
SYM	40	-3055.8062	6376.5219	0.0000	1.0000
F81	38	-3112.5207	6476.7053	0.0000	1.0000
JC69	35	-3124.2267	6480.2492	0.0000	1.0000



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

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

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

