
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:39 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_10494.phylip" (19 taxa, 798 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	-2593.1484	5279.5557	0.6493	0.6493
GTRIG	45	-2592.7009	5280.9071	0.3303	0.9796
GTRI	44	-2596.6099	5286.4787	0.0204	1.0000
HKYG	40	-2608.2157	5300.7642	0.0000	1.0000
HKYIG	41	-2607.8400	5302.2356	0.0000	1.0000
HKYI	40	-2611.6406	5307.6142	0.0000	1.0000
SYMG	41	-2626.1585	5338.8725	0.0000	1.0000
SYMIG	42	-2625.8475	5340.4790	0.0000	1.0000
SYMI	41	-2629.2201	5344.9957	0.0000	1.0000
K2PG	37	-2635.0245	5347.7490	0.0000	1.0000
GTR	43	-2628.9061	5348.8307	0.0000	1.0000
K2PIG	38	-2634.7153	5349.3357	0.0000	1.0000
K2PI	37	-2638.1189	5353.9378	0.0000	1.0000
HKY	39	-2641.9972	5366.1104	0.0000	1.0000
F81G	39	-2655.3444	5392.8050	0.0000	1.0000
F81IG	40	-2655.0835	5394.4999	0.0000	1.0000
F81I	39	-2658.1562	5398.4285	0.0000	1.0000
SYM	40	-2658.8496	5402.0322	0.0000	1.0000
K2P	36	-2667.9763	5411.4532	0.0000	1.0000
JC69G	36	-2680.3064	5436.1135	0.0000	1.0000
JC69IG	37	-2680.0533	5437.8066	0.0000	1.0000
JC69I	36	-2683.1768	5441.8542	0.0000	1.0000
F81	38	-2687.2748	5454.4547	0.0000	1.0000
JC69	35	-2711.8079	5496.9228	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2593.1484	5274.2967	0.6223	0.6223
GTRIG	45	-2592.7009	5275.4017	0.3581	0.9805
GTRI	44	-2596.6099	5281.2197	0.0195	1.0000
HKYG	40	-2608.2157	5296.4313	0.0000	1.0000
HKYIG	41	-2607.8400	5297.6801	0.0000	1.0000
HKYI	40	-2611.6406	5303.2813	0.0000	1.0000
SYMG	41	-2626.1585	5334.3169	0.0000	1.0000
SYMIG	42	-2625.8475	5335.6949	0.0000	1.0000
SYMI	41	-2629.2201	5340.4401	0.0000	1.0000
GTR	43	-2628.9061	5343.8121	0.0000	1.0000
K2PG	37	-2635.0245	5344.0490	0.0000	1.0000
K2PIG	38	-2634.7153	5345.4305	0.0000	1.0000
K2PI	37	-2638.1189	5350.2378	0.0000	1.0000
HKY	39	-2641.9972	5361.9943	0.0000	1.0000
F81G	39	-2655.3444	5388.6889	0.0000	1.0000
F81IG	40	-2655.0835	5390.1670	0.0000	1.0000
F81I	39	-2658.1562	5394.3124	0.0000	1.0000
SYM	40	-2658.8496	5397.6993	0.0000	1.0000
K2P	36	-2667.9763	5407.9526	0.0000	1.0000
JC69G	36	-2680.3064	5432.6129	0.0000	1.0000
JC69IG	37	-2680.0533	5434.1066	0.0000	1.0000
JC69I	36	-2683.1768	5438.3535	0.0000	1.0000
F81	38	-2687.2748	5450.5496	0.0000	1.0000
JC69	35	-2711.8079	5493.6157	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2593.1484	5480.3095	0.7787	0.7787
HKYG	40	-2608.2157	5483.7157	0.1418	0.9205
GTRIG	45	-2592.7009	5486.0966	0.0431	0.9636
GTRI	44	-2596.6099	5487.2325	0.0244	0.9881
HKYIG	41	-2607.8400	5489.6465	0.0073	0.9954
HKYI	40	-2611.6406	5490.5656	0.0046	1.0000
K2PG	37	-2635.0245	5517.2870	0.0000	1.0000
K2PIG	38	-2634.7153	5523.3507	0.0000	1.0000
K2PI	37	-2638.1189	5523.4758	0.0000	1.0000
SYMG	41	-2626.1585	5526.2834	0.0000	1.0000
SYMIG	42	-2625.8475	5532.3435	0.0000	1.0000
SYMI	41	-2629.2201	5532.4066	0.0000	1.0000
HKY	39	-2641.9972	5544.5966	0.0000	1.0000
GTR	43	-2628.9061	5545.1428	0.0000	1.0000
F81G	39	-2655.3444	5571.2911	0.0000	1.0000
K2P	36	-2667.9763	5576.5085	0.0000	1.0000
F81I	39	-2658.1562	5576.9146	0.0000	1.0000
F81IG	40	-2655.0835	5577.4513	0.0000	1.0000
SYM	40	-2658.8496	5584.9836	0.0000	1.0000
JC69G	36	-2680.3064	5601.1688	0.0000	1.0000
JC69I	36	-2683.1768	5606.9094	0.0000	1.0000
JC69IG	37	-2680.0533	5607.3447	0.0000	1.0000
F81	38	-2687.2748	5628.4697	0.0000	1.0000
JC69	35	-2711.8079	5657.4895	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_10494.phylip = 1 - 798;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10494.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_10494.phylip = 1 - 798;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10494.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/chr7_10494.phylip = 1 - 798;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10494.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

