Modelling claim number using a new mixture model: negative binomial gamma distribution

ABSTRACT

In actuarial applications, mixed Poisson distributions are widely used for modelling claim counts as observed data on the number of claims often exhibit a variance noticeably exceeding the mean. In this study, a new claim number distribution is obtained by mixing negative binomial parameter p which is reparameterized as p?=?exp( ?λ) with Gamma distribution. Basic properties of this new distribution are given. Maximum likelihood estimators of the parameters are calculated using the Newton–Raphson and genetic algorithm (GA). We compared the performance of these methods in terms of efficiency by simulation. A numerical example is provided.     

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.     

Inflated modified power series distributions with applications

In certain applications involving discrete data, it is sometimes found that X = 0 is observed with a frequency significantly higher than predicted by the assumed model. Zero inflated Poisson, binomial and negative binomial models have been employed in some clinical trials and in some regression analysis problems.

In this paper, we study the zero inflated modified power series distributions (IMPSD) which include among others the generalized Poisson and the generalized negative binomial distributions and hence the Poisson, binomial and negative binomial distributions. The structural properties along with the distribution of the sum of independent IMPSD variables are studied. The maximum likelihood estimation of the parameters of the model is examined and the variance-covariance matrix of the estimators is obtained. Finally, examples are presented for the generalized Poisson distribution to illustrate the results.     

A note on an integer valued time series model with Poisson–negative binomial marginal distribution

Abstract

This paper proposes a new model for autoregressive time series of counts in terms of a convolution of Poisson and negative binomial random variables, known as Poisson–negative binomial (PNB) distribution. The corresponding first-order integer valued time series models are developed and their properties are discussed. The geometric PNB and the geometric semi PNB distributions are also introduced and studied.     

Tests for detecting overdispersion in poisson models

Collings and Margolin(1985) developed a locally most powerful unbiased test for detecting negative binomial departures from a Poisson model, when the variance is a quadratic function of the mean. Kim and Park(1992) developed a locally most powerful unbiased test, when the variance is a linear function of the mean. It is found that a different mean-variance structure of a negative binomial derives a different locally optimal test statistic.

In this paper Collings and Margolin's and Kim and Park's results are unified and extended by developing a test for overdispersion in Poisson model against Katz family of distributions, Our setup has two extensions: First, Katz family of distributions is employed as an extension of the negative binomial distribution. Second, the mean-variance structure of the mixed Poisson model is given by σ² = μ+cμ^r for arbitrary but fixed r. We derive a local score test for testing H₀ : c = 0. Superiority of a new test is proved by the asymtotic relative efficiency as well as the simulation study.     

Approximate tolerance intervals for the discrete Poisson–Lindley distribution

The Poisson–Lindley distribution is a compound discrete distribution that can be used as an alternative to other discrete distributions, like the negative binomial. This paper develops approximate one-sided and equal-tailed two-sided tolerance intervals for the Poisson–Lindley distribution. Practical applications of the Poisson–Lindley distribution frequently involve large samples, thus we utilize large-sample Wald confidence intervals in the construction of our tolerance intervals. A coverage study is presented to demonstrate the efficacy of the proposed tolerance intervals. The tolerance intervals are also demonstrated using two real data sets. The R code developed for our discussion is briefly highlighted and included in the tolerance package.     

On a class of bivariate mixed Sarmanov distributions

Multivariate distributions are more and more used to model the dependence encountered in many fields. However, classical multivariate distributions can be restrictive by their nature, while Sarmanov's multivariate distribution, by joining different marginals in a flexible and tractable dependence structure, often provides a valuable alternative. In this paper, we introduce some bivariate mixed Sarmanov distributions with the purpose to extend the class of bivariate Sarmanov distributions and to obtain new dependency structures. Special attention is paid to the bivariate mixed Sarmanov distribution with Poisson marginals and, in particular, to the resulting bivariate Sarmanov distributions with negative binomial and with Poisson‐inverse Gaussian marginals; these particular types of mixed distributions have possible applications in, for example modelling bivariate count data. The extension to higher dimensions is also discussed. Moreover, concerning the dependency structure, we also present some correlation formulas.     

Property of bivariate poisson distribution and its application to stochastic processes

Assuming that the random vectors X ₁ and X ₂ have independent bivariate Poisson distributions, the conditional distribution of X ₁ given X ₁?+?X ₂?=?n is obtained. The conditional distribution turns out to be a finite mixture of distributions involving univariate binomial distributions and the mixing proportions are based on a bivariate Poisson (BVP) distribution. The result is used to establish two properties of a bivariate Poisson stochastic process which are the bivariate extensions of the properties for a Poisson process given by Karlin, S. and Taylor, H. M. (1975). A First Course in Stochastic Processes, Academic Press, New York.     

The Least Upper Bound on the Poisson-Negative Binomial Relative Error

In this article, basic mathematical computations are used to determine the least upper bound on the relative error between the negative binomial cumulative distribution function with parameters n and p and the Poisson cumulative distribution function with mean λ =nq = n(1 ? p). Following this bound, it is indicated that the negative binomial cumulative distribution function can be properly approximated by the Poisson cumulative distribution function whenever q is sufficiently small. Five numerical examples are presented to illustrate the obtained result.     

ESTIMATING A RATIO OF MEANS FROM BIVARIATE COUNTS WITH APPLICATIONS IN STEREOLOGY

In survey sampling and in stereology, it is often desirable to estimate the ratio of means θ= E(Y)/E(X) from bivariate count data (X, Y) with unknown joint distribution. We review methods that are available for this problem, with particular reference to stereological applications. We also develop new methods based on explicit statistical models for the data, and associated model diagnostics. The methods are tested on a stereological dataset. For point‐count data, binomial regression and bivariate binomial models are generally adequate. Intercept‐count data are often overdispersed relative to Poisson regression models, but adequately fitted by negative binomial regression.     

Joint (k,l)-hurdle random effects models for mixed longitudinal power series and normal outcomes

A random effects model for analyzing mixed longitudinal normal and count outcomes with and without the possibility of non ignorable missing outcomes is presented. The count response is inflated in two points (k and l) and the (k, l)-Hurdle power series is used as its distribution. The new distribution contains, as special submodels, several important distributions which are discussed, such as (k, l)-Hurdle Poisson and (k, l)-Hurdle negative binomial and (k, l)-Hurdle binomial distributions among others. Random effects are used to take into account the correlation between longitudinal outcomes and inflation parameters. A full likelihood-based approach is used to yield maximum likelihood estimates of the model parameters. A simulation study is performed in which for count outcome (k, l)-Hurdle Poisson, (k, l)-Hurdle negative binomial and (k, l)-Hurdle binomial distributions are considered. To illustrate the application of such modelling the longitudinal data of body mass index and the number of joint damage are analyzed.     

A Bayesian random effects model for analyzing mixed negative binomial and continuous longitudinal responses

ABSTRACT

A general Bayesian random effects model for analyzing longitudinal mixed correlated continuous and negative binomial responses with and without missing data is presented. This Bayesian model, given some random effects, uses a normal distribution for the continuous response and a negative binomial distribution for the count response. A Markov Chain Monte Carlo sampling algorithm is described for estimating the posterior distribution of the parameters. This Bayesian model is illustrated by a simulation study. For sensitivity analysis to investigate the change of parameter estimates with respect to the perturbation from missing at random to not missing at random assumption, the use of posterior curvature is proposed. The model is applied to a medical data, obtained from an observational study on women, where the correlated responses are the negative binomial response of joint damage and continuous response of body mass index. The simultaneous effects of some covariates on both responses are also investigated.     

Expectation identity for the binomial distribution and its application in the calculations of high-order binomial moments

ABSTRACT

For the exponential families normal, gamma, beta, Poisson, and negative binomial, there exists an expectation identity for each of the family. For the binomial family, we discover an expectation identity, which is useful in analytical calculations of its high-order moments.     

Bayesian estimation for non zero inflated modified power series distribution under linex and generalized entropy loss functions

Abstract

In this paper, we derive Bayesian estimators of the parameters of modified power series distributions inflated at any of a support point under linex and general entropy loss function. We assume that the prior information can be summarized by a uniform, Beta, two-sided power, Gamma or generalized Pareto distributions. The obtained results are demonstrated on the generalized Poisson and the generalized negative binomial distribution inflated at a given point.     

The Generalized Waring Process and Its Application

The generalized Waring distribution is a discrete distribution with a wide spectrum of applications in areas such as accident statistics, income analysis, environmental statistics, etc. It has been used as a model that better describes such practical situations as opposed to the Poisson distribution or the negative binomial distribution. Associated to both the Poisson and negative binomial distributions are the well-known Poisson and Pólya processes. In this article, the generalized Waring process is defined. Two models have been shown to lead to the generalized Waring process. One is related to a Cox process, while the other is a compound Poisson process. The defined generalized Waring process is shown to be a stationary, but non homogenous Markov process. Several properties are studied and the intensity, individual intensity, and Chapman–Kolmogorov differential equations of it are obtained. Moreover, the Poisson and Pólya processes are shown to arise as special cases of the generalized Waring process. Using this fact, some known results and some properties of them are obtained.     

A Conway Maxwell Poisson type generalization of the negative hypergeometric distribution

Abstract

Negative hypergeometric distribution arises as a waiting time distribution when we sample without replacement from a finite population. It has applications in many areas such as inspection sampling and estimation of wildlife populations. However, as is well known, the negative hypergeometric distribution is over-dispersed in the sense that its variance is greater than the mean. To make it more flexible and versatile, we propose a modified version of negative hypergeometric distribution called COM-Negative Hypergeometric distribution (COM-NH) by introducing a shape parameter as in the COM-Poisson and COMP-Binomial distributions. It is shown that under some limiting conditions, COM-NH approaches to a distribution that we call the COM-Negative binomial (COMP-NB), which in turn, approaches to the COM Poisson distribution. For the proposed model, we investigate the dispersion characteristics and shape of the probability mass function for different combinations of parameters. We also develop statistical inference for this model including parameter estimation and hypothesis tests. In particular, we investigate some properties such as bias, MSE, and coverage probabilities of the maximum likelihood estimators for its parameters by Monte Carlo simulation and likelihood ratio test to assess shape parameter of the underlying model. We present illustrative data to provide discussion.     

Discriminating between the Poisson and negative binomial distributions:an application to goal scoring in association football

Negative binomial and Poisson distributions are fitted to data on scores in Association Football for the seasons 1946–47 to 1983–84. There are strong grounds for preferring the negative binomial up to 1970; thereafter the Poisson seems adequate. Simplification is achieved by fitting the negative binomial with a common parameter. The analyses are set in the context of previous applications and interpretations in the area. Different models giving rise to the negative binomial or Poisson are investigated and some support found for models not previously advanced in this context. Notwith-standing the success of such exercises some scepticism is expressed about the interpretations placed on previous analyses.     

Negative binomial control limits for count data with extra‐Poisson variation

Traditional techniques for calculating control limits for processes with discrete responses are based on the Poisson distribution. However, for many processes, the assumption of a Poisson distribution is violated. In such cases, use of traditional Poisson control limits may result in an inflated risk of Type I error. The negative binomial distribution is a natural extension of the Poisson distribution and allows for over‐dispersion relative to the Poisson distribution. A simple approach to calculating exact and approximate control limits for count data based on the negative binomial distribution is described. The approach is illustrated by application to water bacteria count data taken from a water purification system. Copyright © 2003 John Wiley & Sons, Ltd.     

On the Bayes estimators of the parameters of size-biased generalized power series distributions

ABSTRACT

In this paper, we derive the Bayes estimators of functions of parameters of the size-biased generalized power series distribution under squared error loss function and weighted square error loss function. The results of size-biased GPSD are then used to obtain particular cases of the size-biased negative binomial, size-biased logarithmic series, and size-biased Poisson distributions. These estimators are better than the classical minimum variance unbiased estimators in the sense that they increase the range of the estimation. Finally, an example is provided to illustrate the results and a goodness of fit test is done using the maximum likelihood and Bayes estimators.     

New polya and inverse polya distributions of order k

New Polya and inverse Polya distributions of order k are derived by means of generalized urn models and by compounding the binomial and negative binomial distributions of order k of Philippou (1986, 1983) with the beta distribution. It i s noted that the present Polpa distribution of order k includes as special cases a new hypergeometric distribution of order k, a negative one,an inverse one, and a discrete uniform of the same order. The probability generating functions, means and variances of the new distributions are obtained, and five asymptotic results are established relating them to the abovedmentioned binomial and negative binomial distributions of order k, and to the Poisson distribution of the same order of Philippou (1983).Moment estimates are also given and applications are indicated.