Implementation of a Stirling number estimator enables direct calculation of population genetics tests for large sequence datasets

摘要

Motivation Stirling numbers enter into the calculation of several population genetics statistics, including Fu's F-s. However, as alignments become large (>= 50 sequences), the Stirling numbers required rapidly exceed the standard floating point range. Another recursive method for calculating Fu's F-s suffers from floating point underflow issues.