
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:59 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/chr5_10816.phylip" (19 taxa, 717 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	-2172.6665	4439.2259	0.8145	0.8145
GTRI	44	-2174.2060	4442.3048	0.1747	0.9892
GTRIG	45	-2175.8566	4447.8832	0.0107	0.9999
HKYIG	41	-2186.3702	4459.8427	0.0000	1.0000
HKYG	40	-2187.6632	4460.1785	0.0000	1.0000
HKYI	40	-2188.2519	4461.3560	0.0000	1.0000
SYMG	41	-2198.6555	4484.4132	0.0000	1.0000
SYMIG	42	-2197.9234	4485.2058	0.0000	1.0000
SYMI	41	-2199.8840	4486.8701	0.0000	1.0000
K2PG	37	-2205.6566	4489.4546	0.0000	1.0000
K2PIG	38	-2204.8684	4490.1085	0.0000	1.0000
K2PI	37	-2206.6383	4491.4181	0.0000	1.0000
GTR	43	-2208.7163	4509.0551	0.0000	1.0000
F81G	39	-2218.0834	4518.7754	0.0000	1.0000
F81IG	40	-2217.8317	4520.5155	0.0000	1.0000
F81I	39	-2219.2958	4521.2001	0.0000	1.0000
HKY	39	-2221.1717	4524.9520	0.0000	1.0000
JC69G	36	-2233.7450	4543.4077	0.0000	1.0000
JC69I	36	-2234.4846	4544.8868	0.0000	1.0000
JC69IG	37	-2233.5649	4545.2711	0.0000	1.0000
SYM	40	-2231.8991	4548.6503	0.0000	1.0000
K2P	36	-2238.5086	4552.9349	0.0000	1.0000
F81	38	-2249.6394	4579.6505	0.0000	1.0000
JC69	35	-2264.8777	4603.4559	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2172.6665	4433.3331	0.8132	0.8132
GTRI	44	-2174.2060	4436.4119	0.1744	0.9876
GTRIG	45	-2175.8566	4441.7133	0.0123	1.0000
HKYIG	41	-2186.3702	4454.7404	0.0000	1.0000
HKYG	40	-2187.6632	4455.3264	0.0000	1.0000
HKYI	40	-2188.2519	4456.5039	0.0000	1.0000
SYMG	41	-2198.6555	4479.3110	0.0000	1.0000
SYMIG	42	-2197.9234	4479.8468	0.0000	1.0000
SYMI	41	-2199.8840	4481.7679	0.0000	1.0000
K2PG	37	-2205.6566	4485.3132	0.0000	1.0000
K2PIG	38	-2204.8684	4485.7368	0.0000	1.0000
K2PI	37	-2206.6383	4487.2767	0.0000	1.0000
GTR	43	-2208.7163	4503.4325	0.0000	1.0000
F81G	39	-2218.0834	4514.1669	0.0000	1.0000
F81IG	40	-2217.8317	4515.6634	0.0000	1.0000
F81I	39	-2219.2958	4516.5915	0.0000	1.0000
HKY	39	-2221.1717	4520.3435	0.0000	1.0000
JC69G	36	-2233.7450	4539.4901	0.0000	1.0000
JC69I	36	-2234.4846	4540.9692	0.0000	1.0000
JC69IG	37	-2233.5649	4541.1298	0.0000	1.0000
SYM	40	-2231.8991	4543.7982	0.0000	1.0000
K2P	36	-2238.5086	4549.0173	0.0000	1.0000
F81	38	-2249.6394	4575.2788	0.0000	1.0000
JC69	35	-2264.8777	4599.7555	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2172.6665	4634.6364	0.6743	0.6743
GTRI	44	-2174.2060	4637.7152	0.1446	0.8190
HKYG	40	-2187.6632	4638.3294	0.1064	0.9254
HKYI	40	-2188.2519	4639.5069	0.0591	0.9844
HKYIG	41	-2186.3702	4642.3185	0.0145	0.9989
GTRIG	45	-2175.8566	4647.5917	0.0010	1.0000
K2PG	37	-2205.6566	4654.5910	0.0000	1.0000
K2PI	37	-2206.6383	4656.5545	0.0000	1.0000
K2PIG	38	-2204.8684	4659.5897	0.0000	1.0000
SYMG	41	-2198.6555	4666.8891	0.0000	1.0000
SYMI	41	-2199.8840	4669.3460	0.0000	1.0000
SYMIG	42	-2197.9234	4672.0000	0.0000	1.0000
F81G	39	-2218.0834	4692.5948	0.0000	1.0000
F81I	39	-2219.2958	4695.0195	0.0000	1.0000
F81IG	40	-2217.8317	4698.6665	0.0000	1.0000
HKY	39	-2221.1717	4698.7714	0.0000	1.0000
GTR	43	-2208.7163	4700.1608	0.0000	1.0000
JC69G	36	-2233.7450	4704.1928	0.0000	1.0000
JC69I	36	-2234.4846	4705.6719	0.0000	1.0000
JC69IG	37	-2233.5649	4710.4076	0.0000	1.0000
K2P	36	-2238.5086	4713.7200	0.0000	1.0000
SYM	40	-2231.8991	4726.8012	0.0000	1.0000
F81	38	-2249.6394	4749.1317	0.0000	1.0000
JC69	35	-2264.8777	4759.8831	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/chr5_10816.phylip = 1 - 717;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10816.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/chr5_10816.phylip = 1 - 717;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10816.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/chr5_10816.phylip = 1 - 717;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10816.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

