On Computing the Distribution Function for the Sum of Independent and Non-identical Random Indicators Yili Hong

The Poisson binomial distribution is the distribution of the sum of independent and non-identical random indicators. Each indicator follows a Bernoulli distribution with individual success probability. When all success probabilities are equal, the Poisson binomial distribution is a binomial distribution. The Poisson binomial distribution has many applications in different areas such as reliability, survival analysis, survey sampling, econometrics, etc. The computing of the cumulative distribution function (cdf) of the Poisson binomial distribution, however, is not straightforward. Approximation methods such as the Poisson approximation and normal approximations have been used in literature. Recursive formulae also have been used to compute the cdf in some areas. In this paper, we present a simple derivation for an exact formula with a closed-form expression for the cdf of the Poisson binomial distribution. The derivation uses the discrete Fourier transform of the characteristic function of the distribution. We develop an algorithm for efficient implementation of the exact formula. Numerical studies were conducted to study the accuracy of the developed algorithm and the accuracy of approximation methods. We also studied the computational efficiency of different methods. The paper is concluded with a discussion on the use of different methods in practice and some suggestions for practitioners.

[1]  A. Volkova.  A Refinement of the Central Limit Theorem for Sums of Independent Random Indicators , 1996 .

[2]  R. E. Barlow,et al.  Computing k-out-of-n System Reliability , 1984, IEEE Transactions on Reliability.

[3]  K. Athreya,et al.  Measure Theory and Probability Theory , 2006 .

[4]  Yili Hong,et al.  Prediction of remaining life of power transformers based on left truncated and right censored lifetime data , 2009, 0908.2901.

[5]  William S. Griffith,et al.  Optimal Reliability Modeling: Principles and Applications , 2004, Technometrics.

[6]  L. L. Cam,et al.  An approximation theorem for the Poisson binomial distribution. , 1960 .

[7]  M. Fernandez,et al.  Closed-Form Expression for the Poisson-Binomial Probability Density Function , 2010, IEEE Transactions on Aerospace and Electronic Systems.

[8]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[9]  S. M. Samuels.  On the Number of Successes in Independent Trials , 1965 .

[10]  C. Pipper,et al.  [''R"--project for statistical computing]. , 2008, Ugeskrift for laeger.

[11]  Yili Hong,et al.  Field-Failure and Warranty Prediction Based on Auxiliary Use-Rate Information , 2010, Technometrics.

[12]  Jun S. Liu,et al.  Weighted finite population sampling to maximize entropy , 1994 .

[13]  R. Singleton.  An algorithm for computing the mixed radix fast Fourier transform , 1969 .

[14]  Ke Wang,et al.  Multi-Period Corporate Default Prediction with Stochastic Covariates , 2005 .

[15]  Jun S. Liu,et al.  STATISTICAL APPLICATIONS OF THE POISSON-BINOMIAL AND CONDITIONAL BERNOULLI DISTRIBUTIONS , 1997 .

[16]  V. Mikhailov.  On a Refinement of the Central Limit Theorem for Sums of Independent Random Indicators , 1994 .