
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:01 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_1350.phylip" (19 taxa, 693 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	-2630.5475	5355.2060	0.6781	0.6781
GTRIG	45	-2630.1507	5356.7001	0.3212	0.9993
GTRI	44	-2637.4361	5368.9832	0.0007	1.0000
HKYG	40	-2648.1994	5381.4295	0.0000	1.0000
HKYIG	41	-2647.8291	5382.9486	0.0000	1.0000
HKYI	40	-2654.6721	5394.3749	0.0000	1.0000
SYMG	41	-2655.5193	5398.3290	0.0000	1.0000
SYMIG	42	-2655.2202	5399.9972	0.0000	1.0000
K2PG	37	-2665.9310	5410.1552	0.0000	1.0000
K2PIG	38	-2665.6072	5411.7465	0.0000	1.0000
SYMI	41	-2662.3006	5411.8915	0.0000	1.0000
K2PI	37	-2672.9609	5424.2150	0.0000	1.0000
GTR	43	-2675.7899	5443.4102	0.0000	1.0000
HKY	39	-2690.8913	5464.5605	0.0000	1.0000
SYM	40	-2699.7491	5484.5288	0.0000	1.0000
K2P	36	-2709.2498	5494.5607	0.0000	1.0000
F81G	39	-2712.1202	5507.0184	0.0000	1.0000
F81IG	40	-2714.5090	5514.0487	0.0000	1.0000
F81I	39	-2717.8922	5518.5624	0.0000	1.0000
JC69G	36	-2724.3014	5524.6637	0.0000	1.0000
JC69IG	37	-2724.0542	5526.4016	0.0000	1.0000
JC69I	36	-2730.6909	5537.4427	0.0000	1.0000
F81	38	-2749.9609	5580.4538	0.0000	1.0000
JC69	35	-2763.1548	5600.1451	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2630.5475	5349.0949	0.6460	0.6460
GTRIG	45	-2630.1507	5350.3013	0.3534	0.9993
GTRI	44	-2637.4361	5362.8721	0.0007	1.0000
HKYG	40	-2648.1994	5376.3988	0.0000	1.0000
HKYIG	41	-2647.8291	5377.6582	0.0000	1.0000
HKYI	40	-2654.6721	5389.3442	0.0000	1.0000
SYMG	41	-2655.5193	5393.0387	0.0000	1.0000
SYMIG	42	-2655.2202	5394.4403	0.0000	1.0000
K2PG	37	-2665.9310	5405.8620	0.0000	1.0000
SYMI	41	-2662.3006	5406.6012	0.0000	1.0000
K2PIG	38	-2665.6072	5407.2144	0.0000	1.0000
K2PI	37	-2672.9609	5419.9218	0.0000	1.0000
GTR	43	-2675.7899	5437.5797	0.0000	1.0000
HKY	39	-2690.8913	5459.7826	0.0000	1.0000
SYM	40	-2699.7491	5479.4982	0.0000	1.0000
K2P	36	-2709.2498	5490.4997	0.0000	1.0000
F81G	39	-2712.1202	5502.2404	0.0000	1.0000
F81IG	40	-2714.5090	5509.0180	0.0000	1.0000
F81I	39	-2717.8922	5513.7844	0.0000	1.0000
JC69G	36	-2724.3014	5520.6028	0.0000	1.0000
JC69IG	37	-2724.0542	5522.1084	0.0000	1.0000
JC69I	36	-2730.6909	5533.3817	0.0000	1.0000
F81	38	-2749.9609	5575.9217	0.0000	1.0000
JC69	35	-2763.1548	5596.3095	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2630.5475	5548.9002	0.9359	0.9359
GTRIG	45	-2630.1507	5554.6477	0.0529	0.9888
HKYG	40	-2648.1994	5558.0400	0.0097	0.9985
GTRI	44	-2637.4361	5562.6774	0.0010	0.9994
HKYIG	41	-2647.8291	5563.8405	0.0005	1.0000
HKYI	40	-2654.6721	5570.9854	0.0000	1.0000
K2PG	37	-2665.9310	5573.8801	0.0000	1.0000
SYMG	41	-2655.5193	5579.2209	0.0000	1.0000
K2PIG	38	-2665.6072	5579.7735	0.0000	1.0000
SYMIG	42	-2655.2202	5585.1636	0.0000	1.0000
K2PI	37	-2672.9609	5587.9399	0.0000	1.0000
SYMI	41	-2662.3006	5592.7834	0.0000	1.0000
GTR	43	-2675.7899	5632.8440	0.0000	1.0000
HKY	39	-2690.8913	5636.8827	0.0000	1.0000
K2P	36	-2709.2498	5653.9768	0.0000	1.0000
SYM	40	-2699.7491	5661.1394	0.0000	1.0000
F81G	39	-2712.1202	5679.3406	0.0000	1.0000
JC69G	36	-2724.3014	5684.0798	0.0000	1.0000
JC69IG	37	-2724.0542	5690.1265	0.0000	1.0000
F81IG	40	-2714.5090	5690.6592	0.0000	1.0000
F81I	39	-2717.8922	5690.8846	0.0000	1.0000
JC69I	36	-2730.6909	5696.8588	0.0000	1.0000
F81	38	-2749.9609	5748.4809	0.0000	1.0000
JC69	35	-2763.1548	5755.2456	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/chr6_1350.phylip = 1 - 693;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_1350.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/chr6_1350.phylip = 1 - 693;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_1350.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_1350.phylip = 1 - 693;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_1350.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

