Erratum to: Uncovering rate variation of lateral gene transfer during bacterial genome evolution

On page 11 of our publication [1], the correction for missing data in a discrete Γ distribution was conducted as

$$Q_{+}={\displaystyle \prod _{i=1}^N{L}_{+}^i}={\displaystyle \prod _{i=1}^N{\displaystyle \sum _{j=1}^M{p}_j\frac{L^i({\mu }_j)}{1-L_{-}^i({\mu }_j)}}.}$$

(1)

This correction, equation (1), is only correct when the category probabilities p _j have no effect on observability. Because the category probabilities have an effect on observability, the appropriate correction is

$$Q_{+}={\displaystyle \prod _{i=1}^N{L}_{+}^i}={\displaystyle \prod _{i=1}^N\frac{{\displaystyle {\sum }_{j=1}^M{p}_j{L}^i({\mu }_j)}}{1-{\displaystyle {\sum }_{j=1}^M{p}_j{L}_{-}^i({\mu }_j)}}.}$$

(2)

Here N is the total number of gene families and M is the number of rate categories. Analyses based on equation (2) as opposed to equation (1) revealed remarkably similar results and all conclusions made based on the old method remain qualitatively unchanged. The additional analyses and the altered tables are provided in Additional File 1 and also available at http://evol.mcmaster.ca/~weilong/research.html. We thank Dr. Matt Spencer for pointing out the difference in averaging methods and thank Drs. David Bryant, Ofir Cohen, Miklós Csűrös, Tal Pupko, and Ed Susko for communication of ideas.