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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-2114.3590	4323.6776	0.5532	0.5532
GTRIG	45	-2113.4256	4324.1400	0.4390	0.9922
HKYG	40	-2124.0637	4333.8518	0.0034	0.9956
GTRI	44	-2119.8110	4334.5816	0.0024	0.9979
HKYIG	41	-2123.4510	4334.9230	0.0020	0.9999
HKYI	40	-2128.0717	4341.8676	0.0001	1.0000
F81G	39	-2164.1475	4411.7305	0.0000	1.0000
F81IG	40	-2163.6587	4413.0416	0.0000	1.0000
SYMG	41	-2162.5752	4413.1713	0.0000	1.0000
SYMIG	42	-2161.7819	4413.8895	0.0000	1.0000
K2PG	37	-2169.5226	4417.9272	0.0000	1.0000
K2PIG	38	-2168.7968	4418.7485	0.0000	1.0000
F81I	39	-2167.8735	4419.1825	0.0000	1.0000
SYMI	41	-2166.8473	4421.7155	0.0000	1.0000
GTR	43	-2165.1173	4422.8733	0.0000	1.0000
K2PI	37	-2173.3471	4425.5761	0.0000	1.0000
HKY	39	-2178.7090	4440.8536	0.0000	1.0000
JC69G	36	-2206.4194	4489.4558	0.0000	1.0000
JC69IG	37	-2205.8778	4490.6375	0.0000	1.0000
JC69I	36	-2209.9606	4496.5382	0.0000	1.0000
F81	38	-2213.3759	4507.9065	0.0000	1.0000
SYM	40	-2212.0364	4509.7971	0.0000	1.0000
K2P	36	-2221.2157	4519.0484	0.0000	1.0000
JC69	35	-2255.8666	4586.0931	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2114.3590	4316.7181	0.5140	0.5140
GTRIG	45	-2113.4256	4316.8513	0.4809	0.9949
GTRI	44	-2119.8110	4327.6221	0.0022	0.9971
HKYG	40	-2124.0637	4328.1275	0.0017	0.9988
HKYIG	41	-2123.4510	4328.9020	0.0012	1.0000
HKYI	40	-2128.0717	4336.1433	0.0000	1.0000
F81G	39	-2164.1475	4406.2950	0.0000	1.0000
SYMG	41	-2162.5752	4407.1503	0.0000	1.0000
F81IG	40	-2163.6587	4407.3174	0.0000	1.0000
SYMIG	42	-2161.7819	4407.5637	0.0000	1.0000
K2PG	37	-2169.5226	4413.0453	0.0000	1.0000
K2PIG	38	-2168.7968	4413.5937	0.0000	1.0000
F81I	39	-2167.8735	4413.7469	0.0000	1.0000
SYMI	41	-2166.8473	4415.6945	0.0000	1.0000
GTR	43	-2165.1173	4416.2347	0.0000	1.0000
K2PI	37	-2173.3471	4420.6942	0.0000	1.0000
HKY	39	-2178.7090	4435.4181	0.0000	1.0000
JC69G	36	-2206.4194	4484.8388	0.0000	1.0000
JC69IG	37	-2205.8778	4485.7556	0.0000	1.0000
JC69I	36	-2209.9606	4491.9212	0.0000	1.0000
F81	38	-2213.3759	4502.7517	0.0000	1.0000
SYM	40	-2212.0364	4504.0729	0.0000	1.0000
K2P	36	-2221.2157	4514.4315	0.0000	1.0000
JC69	35	-2255.8666	4581.7332	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2124.0637	4504.9273	0.8766	0.8766
HKYIG	41	-2123.4510	4510.1218	0.0653	0.9419
GTRG	44	-2114.3590	4511.1978	0.0381	0.9800
HKYI	40	-2128.0717	4512.9431	0.0159	0.9959
GTRIG	45	-2113.4256	4515.7510	0.0039	0.9998
GTRI	44	-2119.8110	4522.1018	0.0002	1.0000
K2PG	37	-2169.5226	4576.5851	0.0000	1.0000
F81G	39	-2164.1475	4578.6748	0.0000	1.0000
K2PIG	38	-2168.7968	4581.5535	0.0000	1.0000
F81IG	40	-2163.6587	4584.1172	0.0000	1.0000
K2PI	37	-2173.3471	4584.2340	0.0000	1.0000
F81I	39	-2167.8735	4586.1267	0.0000	1.0000
SYMG	41	-2162.5752	4588.3701	0.0000	1.0000
SYMIG	42	-2161.7819	4593.2035	0.0000	1.0000
SYMI	41	-2166.8473	4596.9143	0.0000	1.0000
GTR	43	-2165.1173	4606.2945	0.0000	1.0000
HKY	39	-2178.7090	4607.7979	0.0000	1.0000
JC69G	36	-2206.4194	4643.9586	0.0000	1.0000
JC69IG	37	-2205.8778	4649.2954	0.0000	1.0000
JC69I	36	-2209.9606	4651.0410	0.0000	1.0000
F81	38	-2213.3759	4670.7115	0.0000	1.0000
K2P	36	-2221.2157	4673.5513	0.0000	1.0000
SYM	40	-2212.0364	4680.8727	0.0000	1.0000
JC69	35	-2255.8666	4736.4330	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_11571.phylip = 1 - 614;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11571.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_11571.phylip = 1 - 614;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11571.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 (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11571.phylip = 1 - 614;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11571.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

