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.     

Testing overdispersion in the zero-inflated Poisson model

The zero-inflated negative binomial (ZINB) model is used to account for commonly occurring overdispersion detected in data that are initially analyzed under the zero-inflated Poisson (ZIP) model. Tests for overdispersion (Wald test, likelihood ratio test [LRT], and score test) based on ZINB model for use in ZIP regression models have been developed. Due to similarity to the ZINB model, we consider the zero-inflated generalized Poisson (ZIGP) model as an alternate model for overdispersed zero-inflated count data. The score test has an advantage over the LRT and the Wald test in that the score test only requires that the parameter of interest be estimated under the null hypothesis. This paper proposes score tests for overdispersion based on the ZIGP model and illustrates that the derived score statistics are exactly the same as the score statistics under the ZINB model. A simulation study indicates the proposed score statistics are preferred to other tests for higher empirical power. In practice, based on the approximate mean–variance relationship in the data, the ZINB or ZIGP model can be considered, and a formal score test based on asymptotic standard normal distribution can be employed for assessing overdispersion in the ZIP model. We provide an example to illustrate the procedures for data analysis.     

A note on LeCam’s bound for the distance between the Poisson binomial and the Poisson distribution

n possibly different success probabilities p ₁, p ₂, ..., p _n is frequently approximated by a Poisson distribution with parameter λ = p ₁ + p ₂ + ... + p _n . LeCam's bound p ² ₁ + p ² ₂ + ... + p _n ² for the total variation distance between both distributions is particularly useful provided the success probabilities are small. The paper presents an improved version of LeCam's bound if a generalized d-dimensional Poisson binomial distribution is to be approximated by a compound Poisson distribution. Received: May 10, 2000; revised version: January 15, 2001     

Survival functions for the frailty models based on the discrete compound Poisson process

Frailty models are often used to model heterogeneity in survival analysis. The distribution of the frailty is generally assumed to be continuous. In some circumstances, it is appropriate to consider discrete frailty distributions. Having zero frailty can be interpreted as being immune, and population heterogeneity may be analysed using discrete frailty models. In this paper, survival functions are derived for the frailty models based on the discrete compound Poisson process. Maximum likelihood estimation procedures for the parameters are studied. We examine the fit of the models to earthquake and the traffic accidents’ data sets from Turkey.     

Semi-parametric extended Poisson process models for count data

A general framework for the analysis of count data (with covariates) is proposed using formulations for the transition rates of a state-dependent birth process. The form for the transition rates incorporates covariates proportionally, with the residual distribution determined from a smooth non-parametric state-dependent form. Computation of the resulting probabilities is discussed, leading to model estimation using a penalized likelihood function. Two data sets are used as illustrative examples, one representing underdispersed Poisson-like data and the other overdispersed binomial-like data.     

Time series count data regression

The count data model studied in the paper extends the Poisson model by al-lowing for overdispersion and serial correlation. Alternative approaches to esti-mate nuisance parameters, required for the correction of the Poisson maximum likelihood covariance matrix estimator and for a quasi-likelihood estimator, are studied. The estimators are evaluated by finite sample Monte Carlo experi-mentation. It is found that the Poisson maximum likelihood estimator with corrected covariance matrix estimators provide reliable inferences for longer time series. Overdispersion test statistics are wellbehaved, while conventional portmanteau statistics for white noise have too large sizes. Two empirical illustrations are included.     

Interval Estimations for the Mean of the Generalized Poisson and Generalized Negative Binomial Distributions

In this article, we extend the Wald, score, skewness-corrected score, likelihood ratio, and mid-P intervals for the means of the generalized Poisson and generalized negative binomial distributions. These distributions are the members of the discrete version of the natural exponential family (NEF) with cubic variance function (CVF). Also, the coverage probabilities, the distal and mesial noncoverage probabilities, and the lengths of the proposed confidence intervals are estimated by means of a Monte Carlo simulation study. Finally, some practical examples are provided to show the applicability of the proposed intervals in applied studies.     

On some inflated generalized discrete distributions

Generalized discrete distributions such as the double Poisson and the double binomial family of Lagrange distributions are considered when the probabilities are inflated by a constant λ (0 < λ < 1). In each of the above cases, the effect of inflation on the variance is discussed. Also, the Bayesian estimate of inflation as well as those of the parameters are attempted. A maximum likelihood method is also suggested.     

Analysis of over-dispersed count data with extra zeros using the Poisson log-skew-normal distribution

ABSTRACT

Mixed Poisson distributions are widely used in various applications of count data mainly when extra variation is present. This paper introduces an extension in terms of a mixed strategy to jointly deal with extra-Poisson variation and zero-inflated counts. In particular, we propose the Poisson log-skew-normal distribution which utilizes the log-skew-normal as a mixing prior and present its main properties. This is directly done through additional hierarchy level to the lognormal prior and includes the Poisson lognormal distribution as its special case. Two numerical methods are developed for the evaluation of associated likelihoods based on the Gauss–Hermite quadrature and the Lambert's W function. By conducting simulation studies, we show that the proposed distribution performs better than several commonly used distributions that allow for over-dispersion or zero inflation. The usefulness of the proposed distribution in empirical work is highlighted by the analysis of a real data set taken from health economics contexts.     

Characterizations of discrete compound Poisson distributions

ABSTRACT

The aim of this paper is to give some new characterizations of discrete compound Poisson distributions. Firstly, we give a characterization by the Lévy–Khintchine formula of infinitely divisible distributions under some conditions. The second characterization need to present by row sum of random triangular arrays converges in distribution. And we give an application in probabilistic number theory, the strongly additive function converging to a discrete compound Poisson in distribution. The next characterization, is an extension of Watanabe’s theorem of characterization of homogeneous Poisson process. The last characterization will be illustrated by waiting time distributions, especially the matrix-exponential representation.     

Incomplete moments of modified power series districutions with applications

The recurrence relations between the incomplete moments and the factorial incomplete moments of the modified power series distributions (MPSD) are derived. These relations are employed to obtain the experessions for the incomplete moments and the incomplete factorial moments of some particular members of the MPSD class such as the generalized negative binomial, the generalized Poisson, the generalized logrithmic series, the lost game distribution and the distribution of the number of customers served in a busy period. An application of the incomplete moments of the generalized Poisson distribution is provided in the economic selection of a manufactured product. A numerical example is provided using the Poisson distribution and the Generalized Poisson distribution. The example illustrates the difference in results using the two models     

The analysis of incontinence episodes and other count data in patients with overactive bladder by Poisson and negative binomial regression

Clinical studies in overactive bladder have traditionally used analysis of covariance or nonparametric methods to analyse the number of incontinence episodes and other count data. It is known that if the underlying distributional assumptions of a particular parametric method do not hold, an alternative parametric method may be more efficient than a nonparametric one, which makes no assumptions regarding the underlying distribution of the data. Therefore, there are advantages in using methods based on the Poisson distribution or extensions of that method, which incorporate specific features that provide a modelling framework for count data. One challenge with count data is overdispersion, but methods are available that can account for this through the introduction of random effect terms in the modelling, and it is this modelling framework that leads to the negative binomial distribution. These models can also provide clinicians with a clearer and more appropriate interpretation of treatment effects in terms of rate ratios. In this paper, the previously used parametric and non‐parametric approaches are contrasted with those based on Poisson regression and various extensions in trials evaluating solifenacin and mirabegron in patients with overactive bladder. In these applications, negative binomial models are seen to fit the data well. Copyright © 2014 John Wiley & Sons, Ltd.     

Generalized Self-Consistency: Multinomial logit model and Poisson likelihood

A generalized self-consistency approach to maximum likelihood estimation (MLE) and model building was developed in Tsodikov [2003. Semiparametric models: a generalized self-consistency approach. J. Roy. Statist. Soc. Ser. B Statist. Methodology 65(3), 759–774] and applied to a survival analysis problem. We extend the framework to obtain second-order results such as information matrix and properties of the variance. Multinomial model motivates the paper and is used throughout as an example. Computational challenges with the multinomial likelihood motivated Baker [1994. The Multinomial–Poisson transformation. The Statist. 43, 495–504] to develop the Multinomial–Poisson (MP) transformation for a large variety of regression models with multinomial likelihood kernel. Multinomial regression is transformed into a Poisson regression at the cost of augmenting model parameters and restricting the problem to discrete covariates. Imposing normalization restrictions by means of Lagrange multipliers [Lang, J., 1996. On the comparison of multinomial and Poisson log-linear models. J. Roy. Statist. Soc. Ser. B Statist. Methodology 58, 253–266] justifies the approach. Using the self-consistency framework we develop an alternative solution to multinomial model fitting that does not require augmenting parameters while allowing for a Poisson likelihood and arbitrary covariate structures. Normalization restrictions are imposed by averaging over artificial “missing data” (fake mixture). Lack of probabilistic interpretation at the “complete-data” level makes the use of the generalized self-consistency machinery essential.     

Zero-truncated compound Poisson integer-valued GARCH models for time series

Starting from the compound Poisson INGARCH models, we introduce in this paper a new family of integer-valued models suitable to describe count data without zeros that we name zero-truncated CP-INGARCH processes. For such class of models, a probabilistic study concerning moments existence, stationarity and ergodicity is developed. The conditional quasi-maximum likelihood method is introduced to consistently estimate the parameters of a wide zero-truncated compound Poisson subclass of models. The conditional maximum likelihood method is also used to estimate the parameters of ZTCP-INGARCH processes associated with well-specified conditional laws. A simulation study that compares some of those estimators and illustrates their finite distance behaviour as well as a real-data application conclude the paper.     

Bayesian analysis of Cannabis offences using generalized Poisson geometric process model with flexible dispersion

ABSTRACT

This paper derives models to analyse Cannabis offences count series from New South Wales, Australia. The data display substantial overdispersion as well as underdispersion for a subset, trend movement and population heterogeneity. To describe the trend dynamic in the data, the Poisson geometric process model is first adopted and is extended to the generalized Poisson geometric process model to capture both over- and underdispersion. By further incorporating mixture effect, the model accommodates population heterogeneity and enables classification of homogeneous units. The model is implemented using Markov chain Monte Carlo algorithms via the user-friendly WinBUGS software and its performance is evaluated through a simulation study.     

Pairwise likelihood estimation for multivariate mixed Poisson models generated by Gamma intensities

Estimating the parameters of multivariate mixed Poisson models is an important problem in image processing applications, especially for active imaging or astronomy. The classical maximum likelihood approach cannot be used for these models since the corresponding masses cannot be expressed in a simple closed form. This paper studies a maximum pairwise likelihood approach to estimate the parameters of multivariate mixed Poisson models when the mixing distribution is a multivariate Gamma distribution. The consistency and asymptotic normality of this estimator are derived. Simulations conducted on synthetic data illustrate these results and show that the proposed estimator outperforms classical estimators based on the method of moments. An application to change detection in low-flux images is also investigated.     

Efficient regression modeling for correlated and overdispersed count data

Abstract

The objective of this paper is to propose an efficient estimation procedure in a marginal mean regression model for longitudinal count data and to develop a hypothesis test for detecting the presence of overdispersion. We extend the matrix expansion idea of quadratic inference functions to the negative binomial regression framework that entails accommodating both the within-subject correlation and overdispersion issue. Theoretical and numerical results show that the proposed procedure yields a more efficient estimator asymptotically than the one ignoring either the within-subject correlation or overdispersion. When the overdispersion is absent in data, the proposed method might hinder the estimation efficiency in practice, yet the Poisson regression based regression model is fitted to the data sufficiently well. Therefore, we construct the hypothesis test that recommends an appropriate model for the analysis of the correlated count data. Extensive simulation studies indicate that the proposed test can identify the effective model consistently. The proposed procedure is also applied to a transportation safety study and recommends the proposed negative binomial regression model.     

Score test for homogeneity of dispersion in generalized Poisson mixed models with excess zeros

In many applications, the clustered count data often contain excess zeros and the zero-inflated generalized Poisson mixed (ZIGPM) regression model may be suitable. However, dispersion in ZIGPM is often treated as fixed unknown parameter, and this assumption may be not appropriate in some situations. In this article, a score test for homogeneity of dispersion parameter in ZIGPM regression model is developed and corresponding test statistic is obtained. Sampling distribution and power of the score test statistic are investigated through Monte Carlo simulation. Finally, results from a biological example illustrate the usefulness of the diagnostic statistic.