Abstract
The random variable ‘substitution number’ N t , i.e.the number of mutations that have accumulated in a sequence under neutral evolution during a time t, is a cornerstone of the field of molecular evolution. We show here that a complete solution of the moments of this random variable can be obtained explicitly by recurrence, using a simple matrix shift method. This result leads to simple expressions for both the short and long time limits of the moments that can be computed directly from the substitution matrix used to model the neutral evolution. The method developed here is also used to compute the moments of the complementary variable, T n , the time it takes to accumulate n mutations. The method we develop here necessitates only elementary operations on the substitution matrix and does not involve spectral decomposition.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
