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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	40	-3461.3567	7007.1220	0.5973	0.5973
HKYIG	41	-3461.3566	7009.3485	0.1962	0.7936
GTRG	44	-3458.2288	7009.8090	0.1559	0.9494
GTRIG	45	-3458.2289	7012.0599	0.0506	1.0000
HKYI	40	-3475.3893	7035.1871	0.0000	1.0000
GTRI	44	-3472.9349	7039.2213	0.0000	1.0000
SYMG	41	-3512.0687	7110.7727	0.0000	1.0000
GTR	43	-3510.6861	7112.4788	0.0000	1.0000
SYMIG	42	-3512.0688	7113.0056	0.0000	1.0000
HKY	39	-3516.3306	7114.8491	0.0000	1.0000
SYMI	41	-3529.9073	7146.4498	0.0000	1.0000
K2PG	37	-3539.4464	7156.6572	0.0000	1.0000
K2PIG	38	-3539.4463	7158.8658	0.0000	1.0000
K2PI	37	-3553.2436	7184.2515	0.0000	1.0000
SYM	40	-3565.7775	7215.9635	0.0000	1.0000
F81G	39	-3573.0381	7228.2640	0.0000	1.0000
F81IG	40	-3573.0380	7230.4846	0.0000	1.0000
F81I	39	-3585.4381	7253.0641	0.0000	1.0000
K2P	36	-3594.1943	7263.9501	0.0000	1.0000
F81	38	-3630.8562	7341.6856	0.0000	1.0000
JC69G	36	-3691.1146	7457.7907	0.0000	1.0000
JC69IG	37	-3691.1147	7459.9938	0.0000	1.0000
JC69I	36	-3703.1965	7481.9544	0.0000	1.0000
JC69	35	-3739.4475	7552.2595	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-3461.3567	7002.7134	0.5155	0.5155
GTRG	44	-3458.2288	7004.4577	0.2155	0.7311
HKYIG	41	-3461.3566	7004.7132	0.1897	0.9207
GTRIG	45	-3458.2289	7006.4578	0.0793	1.0000
HKYI	40	-3475.3893	7030.7785	0.0000	1.0000
GTRI	44	-3472.9349	7033.8699	0.0000	1.0000
SYMG	41	-3512.0687	7106.1375	0.0000	1.0000
GTR	43	-3510.6861	7107.3722	0.0000	1.0000
SYMIG	42	-3512.0688	7108.1377	0.0000	1.0000
HKY	39	-3516.3306	7110.6611	0.0000	1.0000
SYMI	41	-3529.9073	7141.8146	0.0000	1.0000
K2PG	37	-3539.4464	7152.8928	0.0000	1.0000
K2PIG	38	-3539.4463	7154.8926	0.0000	1.0000
K2PI	37	-3553.2436	7180.4871	0.0000	1.0000
SYM	40	-3565.7775	7211.5549	0.0000	1.0000
F81G	39	-3573.0381	7224.0761	0.0000	1.0000
F81IG	40	-3573.0380	7226.0760	0.0000	1.0000
F81I	39	-3585.4381	7248.8761	0.0000	1.0000
K2P	36	-3594.1943	7260.3886	0.0000	1.0000
F81	38	-3630.8562	7337.7124	0.0000	1.0000
JC69G	36	-3691.1146	7454.2292	0.0000	1.0000
JC69IG	37	-3691.1147	7456.2294	0.0000	1.0000
JC69I	36	-3703.1965	7478.3929	0.0000	1.0000
JC69	35	-3739.4475	7548.8950	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3461.3567	7189.3407	0.9655	0.9655
HKYIG	41	-3461.3566	7196.0063	0.0345	1.0000
GTRG	44	-3458.2288	7209.7477	0.0000	1.0000
GTRIG	45	-3458.2289	7216.4135	0.0000	1.0000
HKYI	40	-3475.3893	7217.4058	0.0000	1.0000
GTRI	44	-3472.9349	7239.1600	0.0000	1.0000
HKY	39	-3516.3306	7292.6228	0.0000	1.0000
SYMG	41	-3512.0687	7297.4305	0.0000	1.0000
SYMIG	42	-3512.0688	7304.0964	0.0000	1.0000
GTR	43	-3510.6861	7307.9966	0.0000	1.0000
K2PG	37	-3539.4464	7325.5231	0.0000	1.0000
K2PIG	38	-3539.4463	7332.1885	0.0000	1.0000
SYMI	41	-3529.9073	7333.1076	0.0000	1.0000
K2PI	37	-3553.2436	7353.1174	0.0000	1.0000
SYM	40	-3565.7775	7398.1822	0.0000	1.0000
F81G	39	-3573.0381	7406.0378	0.0000	1.0000
F81IG	40	-3573.0380	7412.7034	0.0000	1.0000
K2P	36	-3594.1943	7428.3532	0.0000	1.0000
F81I	39	-3585.4381	7430.8378	0.0000	1.0000
F81	38	-3630.8562	7515.0084	0.0000	1.0000
JC69G	36	-3691.1146	7622.1938	0.0000	1.0000
JC69IG	37	-3691.1147	7628.8597	0.0000	1.0000
JC69I	36	-3703.1965	7646.3576	0.0000	1.0000
JC69	35	-3739.4475	7712.1940	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7262.phylip = 1 - 785;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7262.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7262.phylip = 1 - 785;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7262.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7262.phylip = 1 - 785;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7262.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

