A distribution function from population genetics statistics using Stirling numbers of the first kind: asymptotics, inversion and numerical evaluation

DOI10.1090/MCOM/3711MaRDI QIDQ5029477zbMATH OpenOpenAlexFDO

Authors Swaine L. Chen, Nico M. Temme

Publication date 14 February 2022

Published in Mathematics of Computation (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2111.11221

zbMATH Keywords

cumulative distribution functionasymptotic analysisnumerical algorithmsStirling numbers of the first kindevolutionary inference from sequence alignmentspopulation genetics statistics

Mathematics Subject Classification ID

Asymptotic properties of parametric estimators (62F12) Asymptotic approximations, asymptotic expansions (steepest descent, etc.) (41A60) Bell and Stirling numbers (11B73)

Abstract: Stirling numbers of the first kind are common in number theory and combinatorics; through Ewen's sampling formula, these numbers enter into the calculation of several population genetics statistics, such as Fu's Fs. In previous papers we have considered an asymptotic estimator for a finite sum of Stirling numbers, which enables rapid and accurate calculation of Fu's Fs. These sums can also be viewed as a cumulative distribution function; this formulation leads directly to an inversion problem, where, given a value for Fu's Fs, the goal is to solve for one of the input parameters. We solve this inversion using Newton iteration for small parameters. For large parameters we need to extend the earlier obtained asymptotic results to handle the inversion problem asymptotically. Numerical experiments are given to show the efficiency of both solving the inversion problem and the expanded estimator for the statistical quantities.

Recommendations

• Asymptotic expansions for the Stirling numbers of the first kind
• Approximation of inverse moments of discrete distributions
• Review of the stirling numbers, their generalizations and Statistical Applications
• scientific article; zbMATH DE number 837073
• scientific article; zbMATH DE number 1829175

Cites work

• Asymptotic Estimates of Stirling Numbers
• Asymptotic inversion of the incomplete beta function
• Asymptotic methods for integrals
• Efficient algorithms for the inversion of the cumulative central beta distribution
• Numerical Methods for Special Functions

Cited in

(2)

• The genetic sums' representation for the moments of a system of Stieltjes functions and its application
• Mode and Edgeworth expansion for the Ewens distribution and the Stirling numbers

Describes a project that uses

Uses Software

• DLMF

This page was built for publication: A distribution function from population genetics statistics using Stirling numbers of the first kind: asymptotics, inversion and numerical evaluation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5029477)