Statistical control charts for shifted generalized poisson distribution

Summary The use of shifted (or zero-truncated) generalized Poisson distribution to describe the occurrence of events in production processes is considered. The methods of moments and maximum likelihood are proposed for estimating the parameters of shifted generalized Poisson distribution. Control charts for the total number of events and for the average number of events are developed. Finally, a numerical example is used to illustrate the construction of control charts.     

A surplus process involving a compound Poisson counting process and applications

Abstract

In this paper, we introduce a surplus process involving a compound Poisson counting process, which is a generalization of the classical ruin model where the claim-counting process is a homogeneous Poisson process. The incentive is to model batch arrival of claims using a counting process that is based on a compound distribution. This reduces the difficulty of modeling claim amounts and is consistent with industrial data. Recursive formula, some properties and relevant main ruin theory results are provided. Further, we consider applications involving zero-truncated negative binomial and zero-truncated binomial batch arrivals when the claim amounts follow exponential or Erlang distribution.     

A class of observation-driven random coefficient INAR(1) processes based on negative binomial thinning

Integer-valued time series models and their applications have attracted a lot of attention over the last years. In this paper, we introduce a class of observation-driven random coefficient integer-valued autoregressive processes based on negative binomial thinning, where the autoregressive parameter depends on the observed values of the previous moment. Basic probability and statistics properties of the process are established. The unknown parameters are estimated by the conditional least squares and empirical likelihood methods. Specially, we consider three aspects of the empirical likelihood method: maximum empirical likelihood estimate, confidence region and EL test. The performance of the two estimation methods is compared through simulation studies. Finally, an application to a real data example is provided.     

Practical Maximum Pseudolikelihood for Spatial Point Patterns(with Discussion)

This paper describes a technique for computing approximate maximum pseudolikelihood estimates of the parameters of a spatial point process. The method is an extension of Berman & Turner's (1992) device for maximizing the likelihoods of inhomogeneous spatial Poisson processes. For a very wide class of spatial point process models the likelihood is intractable, while the pseudolikelihood is known explicitly, except for the computation of an integral over the sampling region. Approximation of this integral by a finite sum in a special way yields an approximate pseudolikelihood which is formally equivalent to the (weighted) likelihood of a loglinear model with Poisson responses. This can be maximized using standard statistical software for generalized linear or additive models, provided the conditional intensity of the process takes an 'exponential family' form. Using this approach a wide variety of spatial point process models of Gibbs type can be fitted rapidly, incorporating spatial trends, interaction between points, dependence on spatial covariates, and mark information.     

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.     

Conditional Akaike Information Criteria for a Class of Poisson Mixture Models with Random Effects

Focusing on the model selection problems in the family of Poisson mixture models (including the Poisson mixture regression model with random effects and zero‐inflated Poisson regression model with random effects), the current paper derives two conditional Akaike information criteria. The criteria are the unbiased estimators of the conditional Akaike information based on the conditional log‐likelihood and the conditional Akaike information based on the joint log‐likelihood, respectively. The derivation is free from the specific parametric assumptions about the conditional mean of the true data‐generating model and applies to different types of estimation methods. Additionally, the derivation is not based on the asymptotic argument. Simulations show that the proposed criteria have promising estimation accuracy. In addition, it is found that the criterion based on the conditional log‐likelihood demonstrates good model selection performance under different scenarios. Two sets of real data are used to illustrate the proposed method.     

Extended Poisson INAR(1) processes with equidispersion,underdispersion and overdispersion

Real count data time series often show the phenomenon of the underdispersion and overdispersion. In this paper, we develop two extensions of the first-order integer-valued autoregressive process with Poisson innovations, based on binomial thinning, for modeling integer-valued time series with equidispersion, underdispersion, and overdispersion. The main properties of the models are derived. The methods of conditional maximum likelihood, Yule–Walker, and conditional least squares are used for estimating the parameters, and their asymptotic properties are established. We also use a test based on our processes for checking if the count time series considered is overdispersed or underdispersed. The proposed models are fitted to time series of the weekly number of syphilis cases and monthly counts of family violence illustrating its capabilities in challenging the overdispersed and underdispersed count data.     

Bivariate regression models based on compound Poisson distribution

Based on a compound Poisson distribution, new bivariate regression models are introduced and studied. The parameters of the bivariate regression models are estimated by using the maximum likelihood method. Two applications on real datasets are presented to illustrate the models. The results show that these models are compatible to other bivariate Poisson models.     

Score Tests for Testing Independence in the Zero-Truncated Bivariate Poisson Models

In this study, score test statistics for testing independence in the zero-truncated bivariate Poisson distributions are proposed. The Monte Carlo study shows that the score tests proposed in this article keep the significance level close to the nominal one, but the LR and Wald tests over-reject the null hypothesis when it is true. The score tests for testing independence in the zero-truncated bivariate Poisson regression models are also derived in this study.     

Zero-inflated compound Poisson distributions in integer-valued GARCH models

In this paper we introduce a wide class of integer-valued stochastic processes that allows to take into consideration, simultaneously, relevant characteristics observed in count data namely zero inflation, overdispersion and conditional heteroscedasticity. This class includes, in particular, the compound Poisson, the zero-inflated Poisson and the zero-inflated negative binomial INGARCH models, recently proposed in literature. The main probabilistic analysis of this class of processes is here developed. Precisely, first- and second-order stationarity conditions are derived, the autocorrelation function is deduced and the strict stationarity is established in a large subclass. We also analyse in a particular model the existence of higher-order moments and deduce the explicit form for the first four cumulants, as well as its skewness and kurtosis.     

Insurance Premium Prediction via Gradient Tree-Boosted Tweedie Compound Poisson Models

The Tweedie GLM is a widely used method for predicting insurance premiums. However, the structure of the logarithmic mean is restricted to a linear form in the Tweedie GLM, which can be too rigid for many applications. As a better alternative, we propose a gradient tree-boosting algorithm and apply it to Tweedie compound Poisson models for pure premiums. We use a profile likelihood approach to estimate the index and dispersion parameters. Our method is capable of fitting a flexible nonlinear Tweedie model and capturing complex interactions among predictors. A simulation study confirms the excellent prediction performance of our method. As an application, we apply our method to an auto-insurance claim data and show that the new method is superior to the existing methods in the sense that it generates more accurate premium predictions, thus helping solve the adverse selection issue. We have implemented our method in a user-friendly R package that also includes a nice visualization tool for interpreting the fitted model.     

On compounded geometric distributions and their applications

Here, we introduce two-parameter compounded geometric distributions with monotone failure rates. These distributions are derived by compounding geometric distribution and zero-truncated Poisson distribution. Some statistical and reliability properties of the distributions are investigated. Parameters of the proposed distributions are estimated by the maximum likelihood method as well as through the minimum distance method of estimation. Performance of the estimates by both the methods of estimation is compared based on Monte Carlo simulations. An illustration with Air Crash casualties demonstrates that the distributions can be considered as a suitable model under several real situations.     

Monitoring parameter shift with Poisson integer-valued GARCH models

This study examines the statistical process control chart used to detect a parameter shift with Poisson integer-valued GARCH (INGARCH) models and zero-inflated Poisson INGARCH models. INGARCH models have a conditional mean structure similar to GARCH models and are well known to be appropriate to analyzing count data that feature overdispersion. Special attention is paid in this study to conditional and general likelihood ratio-based (CLR and GLR) CUSUM charts and the score function-based CUSUM (SFCUSUM) chart. The performance of each of the proposed methods is evaluated through a simulation study, by calculating their average run length. Our findings show that the proposed methods perform adequately, and that the CLR chart outperforms the GLR chart when there is an increased shift of parameters. Moreover, the use of the SFCUSUM chart in particular is found to lead to a lower false alarm rate than the use of the CLR chart.     

Signed compound poisson integer-valued GARCH processes

Abstract

We propose signed compound Poisson integer-valued GARCH processes for the modeling of the difference of count time series data. We investigate the theoretical properties of these processes and we state their ergodicity and stationarity under mild conditions. We discuss the conditional maximum likelihood estimator when the series appearing in the difference are INGARCH with geometric distribution and explore its finite sample properties in a simulation study. Two real data examples illustrate this methodology.     

Estimating the parameters of a BINMA Poisson model for a non-stationary bivariate time series

This article proposes a novel non-stationary BINMA time series model by extending two INMA processes where their innovation series follow the bivariate Poisson under time-varying moment assumptions. This article also demonstrates, through simulation studies, the use and superiority of the generalized quasi-likelihood (GQL) approach to estimate the regression effects, which is computationally less complicated as compared to conditional maximum likelihood estimation (CMLE) and the feasible generalized least squares (FGLS). The serial and bivariate dependence correlations are estimated by a robust method of moments.     

Regression analysis of zero-inflated time-series counts: application to air pollution related emergency room visit data

Time-series count data with excessive zeros frequently occur in environmental, medical and biological studies. These data have been traditionally handled by conditional and marginal modeling approaches separately in the literature. The conditional modeling approaches are computationally much simpler, whereas marginal modeling approaches can link the overall mean with covariates directly. In this paper, we propose new models that can have conditional and marginal modeling interpretations for zero-inflated time-series counts using compound Poisson distributed random effects. We also develop a computationally efficient estimation method for our models using a quasi-likelihood approach. The proposed method is illustrated with an application to air pollution-related emergency room visits. We also evaluate the performance of our method through simulation studies.     

Likelihood-based and Bayesian methods for Tweedie compound Poisson linear mixed models

The Tweedie compound Poisson distribution is a subclass of the exponential dispersion family with a power variance function, in which the value of the power index lies in the interval (1,2). It is well known that the Tweedie compound Poisson density function is not analytically tractable, and numerical procedures that allow the density to be accurately and fast evaluated did not appear until fairly recently. Unsurprisingly, there has been little statistical literature devoted to full maximum likelihood inference for Tweedie compound Poisson mixed models. To date, the focus has been on estimation methods in the quasi-likelihood framework. Further, Tweedie compound Poisson mixed models involve an unknown variance function, which has a significant impact on hypothesis tests and predictive uncertainty measures. The estimation of the unknown variance function is thus of independent interest in many applications. However, quasi-likelihood-based methods are not well suited to this task. This paper presents several likelihood-based inferential methods for the Tweedie compound Poisson mixed model that enable estimation of the variance function from the data. These algorithms include the likelihood approximation method, in which both the integral over the random effects and the compound Poisson density function are evaluated numerically; and the latent variable approach, in which maximum likelihood estimation is carried out via the Monte Carlo EM algorithm, without the need for approximating the density function. In addition, we derive the corresponding Markov Chain Monte Carlo algorithm for a Bayesian formulation of the mixed model. We demonstrate the use of the various methods through a numerical example, and conduct an array of simulation studies to evaluate the statistical properties of the proposed estimators.     

Estimation Under Inequality Constraints: Semiparametric Estimation of Conditional Duration Models

This article proposes a semiparametric estimator of the parameter in a conditional duration model when there are inequality constraints on some parameters and the error distribution may be unknown. We propose to estimate the parameter by a constrained version of an unrestricted semiparametrically efficient estimator. The main requirement for applying this method is that the initial unrestricted estimator converges in distribution. Apart from this, additional regularity conditions on the data generating process or the likelihood function, are not required. Hence the method is applicable to a broad range of models where the parameter space is constrained by inequality constraints, such as the conditional duration models. In a simulation study involving conditional duration models, the overall performance of the constrained estimator was better than its competitors, in terms of mean squared error. A data example is used to illustrate the method.     

A weighted Poisson distribution and its application to cure rate models

The family of weighted Poisson distributions offers great flexibility in modeling discrete data due to its potential to capture over/under-dispersion by an appropriate selection of the weight function. In this paper, we introduce a flexible weighted Poisson distribution and further study its properties by using it in the context of cure rate modeling under a competing cause scenario. A special case of the new distribution is the COM-Poisson distribution which in turn encompasses the Bernoulli, Poisson, and geometric distributions; hence, many of the well-studied cure rate models may be seen as special cases of the proposed model. We focus on the estimation, through the maximum likelihood method, of the cured proportion and the properties of the failure time of the susceptibles/non cured individuals; a profile likelihood approach is also adopted for estimating the parameters of the weighted Poisson distribution. A Monte Carlo simulation study demonstrates the accuracy of the proposed inferential method. Finally, as an illustration, we fit the proposed model to a cutaneous melanoma data set.     

零浮动Poisson项目计数法在敏感数据抽样调查中的应用

文章将Poisson-Poisson项目计数法进行推广,提出零浮动Poisson项目计数法,其中,非敏感辅助变量来自于一个参数已知的零浮动Poisson分布。并给出了该模型下敏感参数极大似然估计的EM算法以及构造其置信区间的bootstrap方法。此外,还对该模型保护受访者隐私的能力加以讨论,发现该模型的隐私保护要优于Poisson-Poisson项目计数法。最后,从随机模拟的结果表明在该模型下利用本文所介绍的分析方法可以得到敏感参数的较为准确的估计。