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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-2037.4619	4169.0349	0.5532	0.5532
GTRIG	45	-2036.8378	4170.0743	0.3290	0.8822
GTRI	44	-2039.1860	4172.4831	0.0987	0.9809
SYMG	41	-2044.8118	4176.9139	0.0108	0.9916
SYMIG	42	-2044.2621	4178.0812	0.0060	0.9977
SYMI	41	-2046.4110	4180.1123	0.0022	0.9998
HKYG	40	-2050.8908	4186.8122	0.0001	0.9999
HKYIG	41	-2050.3890	4188.0683	0.0000	0.9999
HKYI	40	-2052.2469	4189.5245	0.0000	1.0000
K2PG	37	-2055.6429	4189.5790	0.0000	1.0000
K2PIG	38	-2055.1702	4190.8726	0.0000	1.0000
K2PI	37	-2056.9529	4192.1990	0.0000	1.0000
GTR	43	-2072.9664	4237.7633	0.0000	1.0000
SYM	40	-2078.9025	4242.8357	0.0000	1.0000
HKY	39	-2084.7601	4252.2981	0.0000	1.0000
K2P	36	-2088.8712	4253.8033	0.0000	1.0000
JC69G	36	-2107.8765	4291.8139	0.0000	1.0000
F81G	39	-2104.6590	4292.0959	0.0000	1.0000
JC69IG	37	-2107.3194	4292.9319	0.0000	1.0000
F81IG	40	-2104.0807	4293.1921	0.0000	1.0000
JC69I	36	-2109.1240	4294.3090	0.0000	1.0000
F81I	39	-2105.9448	4294.6676	0.0000	1.0000
JC69	35	-2140.0928	4354.0212	0.0000	1.0000
F81	38	-2137.0175	4354.5670	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2037.4619	4162.9237	0.5294	0.5294
GTRIG	45	-2036.8378	4163.6756	0.3635	0.8929
GTRI	44	-2039.1860	4166.3719	0.0944	0.9873
SYMG	41	-2044.8118	4171.6235	0.0068	0.9942
SYMIG	42	-2044.2621	4172.5242	0.0044	0.9985
SYMI	41	-2046.4110	4174.8219	0.0014	0.9999
HKYG	40	-2050.8908	4181.7815	0.0000	0.9999
HKYIG	41	-2050.3890	4182.7780	0.0000	1.0000
HKYI	40	-2052.2469	4184.4938	0.0000	1.0000
K2PG	37	-2055.6429	4185.2859	0.0000	1.0000
K2PIG	38	-2055.1702	4186.3404	0.0000	1.0000
K2PI	37	-2056.9529	4187.9059	0.0000	1.0000
GTR	43	-2072.9664	4231.9328	0.0000	1.0000
SYM	40	-2078.9025	4237.8050	0.0000	1.0000
HKY	39	-2084.7601	4247.5202	0.0000	1.0000
K2P	36	-2088.8712	4249.7424	0.0000	1.0000
F81G	39	-2104.6590	4287.3180	0.0000	1.0000
JC69G	36	-2107.8765	4287.7529	0.0000	1.0000
F81IG	40	-2104.0807	4288.1614	0.0000	1.0000
JC69IG	37	-2107.3194	4288.6388	0.0000	1.0000
F81I	39	-2105.9448	4289.8896	0.0000	1.0000
JC69I	36	-2109.1240	4290.2480	0.0000	1.0000
F81	38	-2137.0175	4350.0349	0.0000	1.0000
JC69	35	-2140.0928	4350.1856	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-2055.6429	4353.3040	0.6739	0.6739
K2PI	37	-2056.9529	4355.9240	0.1818	0.8557
SYMG	41	-2044.8118	4357.8058	0.0710	0.9267
K2PIG	38	-2055.1702	4358.8996	0.0411	0.9678
SYMI	41	-2046.4110	4361.0042	0.0143	0.9821
GTRG	44	-2037.4619	4362.7291	0.0061	0.9882
SYMIG	42	-2044.2621	4363.2475	0.0047	0.9928
HKYG	40	-2050.8908	4363.4227	0.0043	0.9971
HKYI	40	-2052.2469	4366.1350	0.0011	0.9982
GTRI	44	-2039.1860	4366.1773	0.0011	0.9993
GTRIG	45	-2036.8378	4368.0219	0.0004	0.9997
HKYIG	41	-2050.3890	4368.9602	0.0003	1.0000
K2P	36	-2088.8712	4413.2194	0.0000	1.0000
SYM	40	-2078.9025	4419.4462	0.0000	1.0000
HKY	39	-2084.7601	4424.6203	0.0000	1.0000
GTR	43	-2072.9664	4427.1971	0.0000	1.0000
JC69G	36	-2107.8765	4451.2300	0.0000	1.0000
JC69I	36	-2109.1240	4453.7251	0.0000	1.0000
JC69IG	37	-2107.3194	4456.6569	0.0000	1.0000
F81G	39	-2104.6590	4464.4181	0.0000	1.0000
F81I	39	-2105.9448	4466.9898	0.0000	1.0000
F81IG	40	-2104.0807	4469.8026	0.0000	1.0000
JC69	35	-2140.0928	4509.1217	0.0000	1.0000
F81	38	-2137.0175	4522.5941	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/chr7_11574.phylip = 1 - 693;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_11574.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/chr7_11574.phylip = 1 - 693;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_11574.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 (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_11574.phylip = 1 - 693;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_11574.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

