
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:17 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/chr3_24787.phylip" (19 taxa, 536 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	-2225.3984	4537.4230	0.7504	0.7504
GTRG	44	-2222.5542	4541.1736	0.1150	0.8654
HKYIG	41	-2226.3562	4541.6841	0.0891	0.9545
GTRIG	45	-2222.5541	4543.5572	0.0349	0.9895
HKYI	40	-2229.7242	4546.0747	0.0099	0.9994
GTRI	44	-2227.7815	4551.6281	0.0006	1.0000
K2PG	37	-2249.1790	4578.0045	0.0000	1.0000
K2PIG	38	-2249.1843	4580.3324	0.0000	1.0000
SYMG	41	-2247.2737	4583.5190	0.0000	1.0000
SYMIG	42	-2247.2735	4585.8735	0.0000	1.0000
K2PI	37	-2253.4741	4586.5948	0.0000	1.0000
HKY	39	-2251.3960	4587.0823	0.0000	1.0000
GTR	43	-2248.1471	4589.9854	0.0000	1.0000
SYMI	41	-2252.0175	4593.0067	0.0000	1.0000
K2P	36	-2276.4871	4630.3129	0.0000	1.0000
SYM	40	-2274.3326	4635.2915	0.0000	1.0000
F81G	39	-2282.9524	4650.1951	0.0000	1.0000
F81IG	40	-2282.9546	4652.5355	0.0000	1.0000
F81I	39	-2286.4035	4657.0973	0.0000	1.0000
JC69G	36	-2302.7683	4682.8752	0.0000	1.0000
JC69IG	37	-2302.7681	4685.1828	0.0000	1.0000
JC69I	36	-2306.3392	4690.0172	0.0000	1.0000
F81	38	-2306.5376	4695.0391	0.0000	1.0000
JC69	35	-2326.8210	4728.6821	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2225.3984	4530.7968	0.6302	0.6302
GTRG	44	-2222.5542	4533.1084	0.1984	0.8286
HKYIG	41	-2226.3562	4534.7124	0.0890	0.9176
GTRIG	45	-2222.5541	4535.1083	0.0730	0.9906
HKYI	40	-2229.7242	4539.4485	0.0083	0.9989
GTRI	44	-2227.7815	4543.5629	0.0011	1.0000
K2PG	37	-2249.1790	4572.3579	0.0000	1.0000
K2PIG	38	-2249.1843	4574.3687	0.0000	1.0000
SYMG	41	-2247.2737	4576.5473	0.0000	1.0000
SYMIG	42	-2247.2735	4578.5469	0.0000	1.0000
HKY	39	-2251.3960	4580.7920	0.0000	1.0000
K2PI	37	-2253.4741	4580.9482	0.0000	1.0000
GTR	43	-2248.1471	4582.2943	0.0000	1.0000
SYMI	41	-2252.0175	4586.0350	0.0000	1.0000
K2P	36	-2276.4871	4624.9742	0.0000	1.0000
SYM	40	-2274.3326	4628.6653	0.0000	1.0000
F81G	39	-2282.9524	4643.9048	0.0000	1.0000
F81IG	40	-2282.9546	4645.9092	0.0000	1.0000
F81I	39	-2286.4035	4650.8070	0.0000	1.0000
JC69G	36	-2302.7683	4677.5365	0.0000	1.0000
JC69IG	37	-2302.7681	4679.5363	0.0000	1.0000
JC69I	36	-2306.3392	4684.6785	0.0000	1.0000
F81	38	-2306.5376	4689.0753	0.0000	1.0000
JC69	35	-2326.8210	4723.6421	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2225.3984	4702.1621	0.9710	0.9710
HKYIG	41	-2226.3562	4710.3619	0.0161	0.9871
HKYI	40	-2229.7242	4710.8138	0.0128	0.9999
GTRG	44	-2222.5542	4721.6103	0.0001	1.0000
GTRIG	45	-2222.5541	4727.8943	0.0000	1.0000
K2PG	37	-2249.1790	4730.8709	0.0000	1.0000
GTRI	44	-2227.7815	4732.0648	0.0000	1.0000
K2PIG	38	-2249.1843	4737.1658	0.0000	1.0000
K2PI	37	-2253.4741	4739.4611	0.0000	1.0000
HKY	39	-2251.3960	4747.8733	0.0000	1.0000
SYMG	41	-2247.2737	4752.1968	0.0000	1.0000
SYMIG	42	-2247.2735	4758.4806	0.0000	1.0000
SYMI	41	-2252.0175	4761.6845	0.0000	1.0000
GTR	43	-2248.1471	4766.5121	0.0000	1.0000
K2P	36	-2276.4871	4779.2030	0.0000	1.0000
SYM	40	-2274.3326	4800.0306	0.0000	1.0000
F81G	39	-2282.9524	4810.9860	0.0000	1.0000
F81IG	40	-2282.9546	4817.2746	0.0000	1.0000
F81I	39	-2286.4035	4817.8882	0.0000	1.0000
JC69G	36	-2302.7683	4831.7653	0.0000	1.0000
JC69IG	37	-2302.7681	4838.0492	0.0000	1.0000
JC69I	36	-2306.3392	4838.9073	0.0000	1.0000
F81	38	-2306.5376	4851.8724	0.0000	1.0000
JC69	35	-2326.8210	4873.5868	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/chr3_24787.phylip = 1 - 536;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24787.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/chr3_24787.phylip = 1 - 536;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24787.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/chr3_24787.phylip = 1 - 536;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24787.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

