On Stein's method, smoothing estimates in total variation distance and mixture distributions

In many settings it is useful to have bounds on the total variation distance between some random variable Z and its shifted version Z+1. For example, such quantities are often needed when applying Stein's method for probability approximation. This note considers one way in which such bounds can be derived, in cases where Z is either the equilibrium distribution of some birth-death process or the mixture of such a distribution. Applications of these bounds are given to translated Poisson and compound Poisson approximations for Poisson mixtures and the Pólya distribution.     

Extended Poisson process modelling of dilution series data

Summary.  Data comprising colony counts, or a binary variable representing fertile (or sterile) samples, as a dilution series of the containing medium are analysed by using extended Poisson process modelling. These models form a class of flexible probability distributions that are widely applicable to count and grouped binary data. Standard distributions such as Poisson and binomial, and those representing overdispersion and underdispersion relative to these distributions can be expressed within this class. For all the models in the class, likelihoods can be obtained. These models have not been widely used because of the perceived difficulty of performing the calculations and the lack of associated software. Exact calculation of the probabilities that are involved can be time consuming although accurate approximations that use considerably less computational time are available. Although dilution series data are the focus here, the models are applicable to any count or binary data. A benefit of the approach is the ability to draw likelihood-based inferences from the data.     

Multi-Center Clinical Trials with Random Enrollment: Theoretical Approximations

ABSTRACT

We consider the problem of analyzing multi-center clinical trials when the number of patients at each center and on each treatment arm is random and follows the Poisson distribution. Theoretical approximations are made for the first two moments of the mean square errors (MSE's) for three different estimators of treatment effect difference that are commonly used in multi-center clinical trials. To construct these approximations, approximations are needed for the harmonic mean and negative moments of the Poisson distribution. This is achieved through the use of recurrence relations. The accuracy of the approximations for the moments of the MSE's were then validated through comparing the theoretical values to those obtained from a simulation study under two different enrollment environments.     

The probability function of a geometric Poisson distribution

The geometric Poisson (also called Pólya–Aeppli) distribution is a particular case of the compound Poisson distribution. In this study, the explicit probability function of the geometric Poisson distribution is derived and a straightforward proof for this function is given. By means of a proposed algorithm, some numerical examples and an application on traffic accidents are also given to illustrate the usage of the probability function and proposed algorithm.     

On Infinitely Divisible Exponential Dispersion Model Related to Poisson-Exponential Distribution

We construct a univariate exponential dispersion model comprised of discrete infinitely divisible distributions. This model emerges in the theory of branching processes. We obtain a representation for the Lévy measure of relevant distributions and characterize their laws as Poisson mixtures and/or compound Poisson distributions. The regularity of the unit variance function of this model is employed for the derivation of approximations by the Poisson-exponential model. We emphasize the role of the latter class. We construct local approximations relating them to properties of special functions and branching diffusions.     

Approximations for Multiple Scan Statistics on the Circle

Multiple scan statistic is usually used by epidemiologist to test the uniformity or clustering of data. In this article, we extend the work of Lin (1999) to give a general expression for the moments of multiple scan statistic on a circle, and use these moments to approximate its distribution using Markov chain and compound Poisson approximations proposed by Huffer and Lin (1997a) and Lin (1993). Numerical results are presented to evaluate the performance of these approximations.     

Some Utilizations of Lambert W Function in Distribution Theory

We introduce a new class of positive infinitely divisible probability laws calling them 𝔏_γ distributions. Their cumulant-generating functions (cgf) are expressed in terms of the principal branch of the Lambert W function. The probability density functions (pdfs) of 𝔏_γ laws are bounded resembling pdf of a Lévy stable distribution. The exponential dispersion model constructed starting from an 𝔏_γ distribution admits the inverse Gaussian approximation. The natural exponential family constructed starting from an 𝔏_γ distribution constitutes the reciprocal of the natural exponential family generated by a spectrally negative stable law with α = 1. We derive new results on 𝔏_γ laws and the related exponential dispersion models, including their convolution and scaling closure properties. We generate another exponential dispersion model starting from an exponentially compounded 𝔏_γ law. This distribution emerges in the Poisson mixture representation of a generalized Poisson law. We extend the Poisson approximation for the scaled Neyman type A exponential dispersion model. We derive saddlepoint-type approximations for some of these exponential dispersion models. The role of the Lambert W function is emphasized.     

SMOOTH TESTS FOR THE BIVARIATE POISSON DISTRIBUTION

A theorem of Rayner & Best (1989) is generalised to permit the construction of smooth tests of goodness of fit without requiring a set of orthonormal functions on the hypothesised distribution. This result is used to construct smooth tests for the bivariate Poisson distribution. The test due to Crockett (1979) is similar to a smooth test that assesses the variance structure under the bivariate Poisson model; the test due to Loukas & Kemp (1986) is related to a smooth test that seeks to detect a particular linear relationship between the variances and covariance under the bivariate Poisson model. Using focused smooth tests may be more informative than using previously suggested tests. The distribution of the Loukas & Kemp (1986) statistic is not well approximated by the x²distribution for larger correlations, and a revised statistic is suggested.     

More Damned Lies and Statistics: How Numbers Confuse Public Issues and Chance: A Guide to Gambling,Love, the Stock Market, & Just About Everything Else

In this article we review two historical approximations to the Poisson and binomial cumulative distribution functions (CDFs); that is, the Wilson–Hilferty and Camp–Paulson approximations. Both of these approximations reduce to standard normal formulas that produce very accurate estimates of the Poisson and binomial CDFs, and are thus quite simple to implement. Additionally, in an upper-division undergraduate or master’s level probability and inference course, the derivation of these approximations presents a nice opportunity to introduce and study the distributional relationships between the gamma and Poisson CDFs, and the binomial, beta, and F CDFs. This article presents the basic theorems and lemmas needed to derive each approximation, along with some relevant examples that compare and contrast the precision of these approximations with their large-sample, limiting normal counterparts.     

LOG-LINEAR MODELS FOR MEAN AND DISPERSION IN MIXED POISSON REGRESSION MODELS

This paper is concerned with the analysis of repeated measures count data overdispersed relative to a Poisson distribution, with the overdispersion possibly heterogeneous. To accommodate the overdispersion, the Poisson random variable is compounded with a gamma random variable, and both the mean of the Poisson and the variance of the gamma are modelled using log linear models. Maximum likelihood estimates (MLE) are then obtained. The paper also gives extended quasi-likelihood estimates for a more general class of compounding distributions which are shown to be approximations to the MLEs obtained for the gamma case. The theory is illustrated by modelling the determination of asbestos fibre intensity on membrane filters mounted on microscope slides.     

Abel series distributions with applications to fluctuations of sample functions of stochastic processes

A two-parameter class of discrete distributions, Abel series distributions, generated by expanding a suitable pa,rametric function into a series of Abel polynomials is discussed. An Abel series distribution occurs in fluctuations of sample functions of stochastic processes and has applications in insurance risk, queueing, dam and storage processes. The probability generating function and the factorial moments of the Abel series distributions are obtained in closed forms. It is pointed out that the name of the generalized Poisson distribution of Consul and Jain is justified by the form of its generating function. Finally it is shown that this generalized Poisson distribution is the only member of the Abel series distributions which is closed under convolution.     

Some Theoretical Comparisons of Negative Binomial and Zero-Inflated Poisson Distributions

In this article, we compare the zero-inflated Poisson (ZIP) and negative binomial (NB) distributions based on three most important criteria: the probability of zero, the mean value, and the variance. Our results show that with same mean value and variance, the ZIP distribution always has a larger probability of zeros; with same mean value and probability of zeros, the NB distribution always has a larger variance; and with same variance and probability of zeros, the ZIP distribution always has a larger mean value. We also study the properties of Vuong test in model selection in three cases by simulations.     

Poisson approximation and dna sequence matching

The Poisson distribution is commonly used to model the number of occurrences of independent rare events. However, many instances arise where dependence exists, for example, in counting the length of long head runs in coin tossing, or matches between two DNA sequences. The Chen-Stein method of Poisson approximation yields bounds on the error incurred when approximating the number of occurrences of possibly dependent events by a Poisson random variable of the same mean. In addition to the problems related to the motivating examples from molecular biology involving runs and matches, the method may be applied to questions as varied as calculating probabilities involving extremes of sequences of random variables and approximating the probability of general birthday coincidences.     

Exact short poisson confidence intervals

The authors propose a new method for constructing a confidence interval for the expectation θ of a Poisson random variable. The interval they obtain cannot be shortened without the infimum over θ of the coverage probability falling below 1 ‐ α. In addition, the endpoints of the interval are strictly increasing functions of the observed variable. An easy‐to‐program algorithm is provided for computing this interval.     

Zero-inflated models and estimation in zero-inflated Poisson distribution

In this paper, we briefly overview different zero-inflated probability distributions. We compare the performance of the estimates of Poisson, Generalized Poisson, ZIP, ZIGP and ZINB models through Mean square error (MSE), bias and Standard error (SE) when the samples are generated from ZIP distribution. We propose a new estimator referred to as probability estimator (PE) of inflation parameter of ZIP distribution based on moment estimator (ME) of the mean parameter and compare its performance with ME and maximum likelihood estimator (MLE) through a simulation study. We use the PE along with ME and MLE to fit ZIP distribution to various zero-inflated datasets and observe that the results do not differ significantly. We recommend using PE in place of MLE since it is easy to calculate and the simulation study in this paper demonstrates that the PE performs as good as MLE irrespective of the sample size.     

On Conditional and Partial Correlation

A graphical technique, similar in spirit to probability plotting, can be used to judge whether a Poisson model is appropriate for an observed frequency distribution. This “Poissonness plot” can equally be applied to truncated Poisson situations. It provides a type of robustness for detecting isolated discrepancies in otherwise well-behaved frequency distributions.     

Hypotheses testing for a multidimensional parameter of inhomogeneous Poisson processes

We observe a realization of an inhomogeneous Poisson process whose intensity function depends on an unknown multidimensional parameter. We consider the asymptotic behaviour of the Rao score test for a simple null hypothesis against the multilateral alternative. By using the Edgeworth type expansion (under the null hypothesis) for a vector of stochastic integrals with respect to the Poisson process, we refine the (classic) threshold of the test (obtained by the central limit theorem), which improves the first type probability of error. The expansion allows us to describe the power of the test under the local alternative, i.e. a sequence of alternatives, which converge to the null hypothesis with a certain rate. The rates can be different for components of the parameter.     

Use of an empirical probability generating function for testing a Poisson model

We present a test of the fit to a Poisson model based on the empirical probability generating function (epgf). We derive the limiting distribution of the test under the Poisson hypothesis and show that a rescaling of it is approximately independent of the mean parameter in the Poisson distribution. When inspected under a simulation study over a range of alternative distributions, we find that this test shows reasonable behaviour compared to other goodness-of-fit tests like the Poisson index of dispersion and smooth test applied to the Poisson model. These results illustrate that epgf-based methods for anlyzing count data are promising.     

A new technique for calculating fourier integrals based on the poisson summation formula

This paper derives a general method for the calculation of Fourier integrals based on the Poisson summation formula. The approach is an extension of that due to Medhurst (1960) and can be thought of as an extension of the Mobius inversion formula. Useful applications are to the calculation of probability density funtions from their Fourier transforms.