Experience Rating with Poisson Mixtures

Garfield Brown,Steve Brooks,Winston Buckley

DOI: https://doi.org/10.48550/arXiv.1012.4702

2010-12-21

Abstract:We present a mixture Poisson model for claims counts in which the number of components in the mixture are estimated by reversible jump MCMC methods.

Applications

What problem does this paper attempt to address?

About the Justification of Experience Rating: Bonus Malus System and a new Poisson Mixture Model

Magda Schiegl

DOI: https://doi.org/10.48550/arXiv.1009.4142

2010-09-22

Abstract:The claim experience of the past is a very important information to calculate the fair price of an insurance contract. In a lot of European countries for instance the prices for motor car insurance depend on the number of claims the driver has reported to the insurance company during the last years. Classically these prices are calculated on the basis of a mixed Poisson model with a gamma mixing distribution. The mixing distribution models the car drivers' qualities across the insured portfolio. This is just one example for experience rating. In the classical context the price is equal to the expectation of the Bayesian posterior distribution. In some lines of business (especially third party liability and lines with exposure to extreme weather events) we that the real world data cannot be described well enough by the classical Poisson - gamma model. Therefore we investigate the influence of the mixing distribution on the posterior distribution conditional on the experienced number of claims. This enables the application of other - more risk adequate premium principles than the expectation principle. We introduce the inverse - gamma distribution as a new mixing distribution to model claim numbers and compare it to the classical gamma distribution. In both cases a closed analytic representation of the mixed distribution can be found: In the classic case the well known negative binomial distribution, in our new one a representation using the Bessel functions. Additionally we present numerical results about the tail behaviour of the mixed Poisson - inverse - gamma distribution. Finally we introduce the concept of resolution. It enables us to decide if the classification of risk groups via the number of experienced claims is a risk adequate procedure.

Risk Management,Statistical Finance
Bayesian estimation and prediction for certain mixtures

Aziz LMoudden,Éric Marchand

DOI: https://doi.org/10.48550/arXiv.2005.02879

2020-09-08

Abstract:For two vast families of mixture distributions and a given prior, we provide unified representations of posterior and predictive distributions. Model applications presented include bivariate mixtures of Gamma distributions labelled as Kibble-type, non-central Chi-square and F distributions, the distribution of $R^2$ in multiple regression, variance mixture of normal distributions, and mixtures of location-scale exponential distributions including the multivariate Lomax distribution. An emphasis is also placed on analytical representations and the relationships with a host of existing distributions and several hypergeomtric functions of one or two variables.

Statistics Theory
Bivariate phase-type distributions for experience rating in disability insurance

Christian Furrer,Jacob Juhl Sørensen,Jorge Yslas

2024-05-30

Abstract:In this paper, we consider the problem of experience rating within the classic Markov chain life insurance framework. We begin by investigating various multivariate mixed Poisson models with mixing distributions encompassing independent Gamma, hierarchical Gamma, and multivariate phase-type. In particular, we demonstrate how maximum likelihood estimation for these proposed models can be performed using expectation-maximization algorithms, which might be of independent interest. Subsequently, we establish a link between mixed Poisson distributions and the problem of pricing group disability insurance contracts that exhibit heterogeneity. We focus on shrinkage estimation of disability and recovery rates, taking into account sampling effects such as right-censoring. Finally, we showcase the practicality of these shrinkage estimators through a numerical study based on simulated yet realistic insurance data. Our findings highlight that by allowing for dependency between latent group effects, estimates of recovery and disability rates mutually improve, leading to enhanced predictive performance.

Statistics Theory
Erlang mixture modeling for Poisson process intensities

Hyotae Kim,Athanasios Kottas

DOI: https://doi.org/10.48550/arXiv.2110.12513

2021-10-25

Abstract:We develop a prior probability model for temporal Poisson process intensities through structured mixtures of Erlang densities with common scale parameter, mixing on the integer shape parameters. The mixture weights are constructed through increments of a cumulative intensity function which is modeled nonparametrically with a gamma process prior. Such model specification provides a novel extension of Erlang mixtures for density estimation to the intensity estimation setting. The prior model structure supports general shapes for the point process intensity function, and it also enables effective handling of the Poisson process likelihood normalizing term resulting in efficient posterior simulation. The Erlang mixture modeling approach is further elaborated to develop an inference method for spatial Poisson processes. The methodology is examined relative to existing Bayesian nonparametric modeling approaches, including empirical comparison with Gaussian process prior based models, and is illustrated with synthetic and real data examples.

Methodology
Mixture Models With a Prior on the Number of Components

Jeffrey W. Miller,Matthew T. Harrison

DOI: https://doi.org/10.1080/01621459.2016.1255636

IF: 4.369

2017-11-13

Journal of the American Statistical Association

Abstract:A natural Bayesian approach for mixture models with an unknown number of components is to take the usual finite mixture model with symmetric Dirichlet weights, and put a prior on the number of components-that is, to use a mixture of finite mixtures (MFM). The most commonly-used method of inference for MFMs is reversible jump Markov chain Monte Carlo, but it can be nontrivial to design good reversible jump moves, especially in high-dimensional spaces. Meanwhile, there are samplers for Dirichlet process mixture (DPM) models that are relatively simple and are easily adapted to new applications. It turns out that, in fact, many of the essential properties of DPMs are also exhibited by MFMs-an exchangeable partition distribution, restaurant process, random measure representation, and stick-breaking representation-and crucially, the MFM analogues are simple enough that they can be used much like the corresponding DPM properties. Consequently, many of the powerful methods developed for inference in DPMs can be directly applied to MFMs as well; this simplifies the implementation of MFMs and can substantially improve mixing. We illustrate with real and simulated data, including high-dimensional gene expression data used to discriminate cancer subtypes.

statistics & probability
mixture Poisson exponential–inverse Gaussian regression model

Emilio Gómez-Déniz,Enrique Calderín-Ojeda

DOI: https://doi.org/10.51936/rqpj6167

2016-07-01

Advances in Methodology and Statistics

Abstract:In this paper a mixed Poisson regression model for count data is introduced. This model is derived by mixing the Poisson distribution with the one–parameter continuous exponential–inverse Gaussian distribution. The obtained probability mass function is over-dispersed and unimodal with modal value located at zero. Estimation is performed by maximum likelihood. As an application, the demand for health services among people 65 and over is examined using this regression model since empirical evidence has suggested that the over–dispersion and a large portion of non–users are common features of medical care utilization data.

English Else
Bivariate Gamma Mixture of Experts Models for Joint Insurance Claims Modeling

Sen Hu,T Brendan Murphy,Adrian O'Hagan

DOI: https://doi.org/10.48550/arXiv.1904.04699

2019-04-09

Abstract:In general insurance, risks from different categories are often modeled independently and their sum is regarded as the total risk the insurer takes on in exchange for a premium. The dependence from multiple risks is generally neglected even when correlation could exist, for example a single car accident may result in claims from multiple risk categories. It is desirable to take the covariance of different categories into consideration in modeling in order to better predict future claims and hence allow greater accuracy in ratemaking. In this work multivariate severity models are investigated using mixture of experts models with bivariate gamma distributions, where the dependence structure is modeled directly using a GLM framework, and covariates can be placed in both gating and expert networks. Furthermore, parsimonious parameterisations are considered, which leads to a family of bivariate gamma mixture of experts models. It can be viewed as a model-based clustering approach that clusters policyholders into sub-groups with different dependencies, and the parameters of the mixture models are dependent on the covariates. Clustering is shown to be important in separating the data into sub-groupings where strong dependence is often present, even if the overall data set exhibits only weak dependence. In doing so, the correlation within different components features prominently in the model. It is shown that, by applying to both simulated data and a real-world Irish GI insurer data set, claim predictions can be improved.

Applications
A Mixture Model of Proportional Reversed Hazard Rate

Xiaohu Li,Zhouping Li

DOI: https://doi.org/10.1080/03610920802050935

2008-01-01

Abstract:This article investigates properties of mixture model of proportional reversed hazard rate. Firstly, the mixing random variable and the overall population variable are proved to be positively likelihood dependent. Secondly, lower bounds for the distribution function as well as the conditional distribution are established in the case that the mixing variable belongs to certain nonparametric classes. Finally, some stochastic orders on the mixing (baseline) variables are proved to be translated to the corresponding overall population variables.
Maximum-penalized-likelihood estimation for independent and Markov-dependent mixture models.

M. Puterman,B. Leroux

DOI: https://doi.org/10.2307/2532308

IF: 1.701

1992-06-01

Biometrics

Abstract:This paper concerns the use and implementation of maximum-penalized-likelihood procedures for choosing the number of mixing components and estimating the parameters in independent and Markov-dependent mixture models. Computation of the estimates is achieved via algorithms for the automatic generation of starting values for the EM algorithm. Computation of the information matrix is also discussed. Poisson mixture models are applied to a sequence of counts of movements by a fetal lamb in utero obtained by ultrasound. The resulting estimates are seen to provide plausible mechanisms for the physiological process.

Mathematics,Medicine
Bayesian Ratemaking with Common Effects Modeled by Mixture of Polya Tree Processes

Jianjun Zhang,Chunjuan Qiu,Xianyi Wu

DOI: https://doi.org/10.1016/j.insmatheco.2018.06.007

IF: 2.168

2018-01-01

Insurance Mathematics and Economics

Abstract:In classical models for Bayesian ratemaking, claims are usually assumed to be independent over risks. However, this assumption may be violated because there are situations that could derive possible dependence among the insured individuals. This paper aims to investigate the typical problem of experience ratemaking to account for a special type of dependence that is known as common effects in the literature. Polya tree processes are employed to model the common effects and, by means of an MCMC scheme, the corresponding Bayesian premiums are numerically computed. This provides a useful alternative to the well known results on Bayesian ratemaking with common effects.
Fitting Com-Poisson Mixtures to Bimodal Count Data

Smarajit Bose,Galit Shmueli,Pragya Sur,Paromita Dubey

2013-01-01

Abstract:Bi-modal truncated count distributions are frequently observed in aggregate surveys and ratings when respondents are mixed in their opinion. They also arise in censored count data, where the highest category might create an additional mode. The Poisson distribution is the most common distribution for fitting count data and can be modified to achieve mixtures of truncated Poisson distributions. However, it is suitable only for modeling equi-dispersed distributions and is limited in its ability to capture bi-modality. The Conway-Maxwell-Poisson (CMP) distribution is a two-parameter generalization of the Poisson distribution that allows for over- and under-dispersion. While the CMP is much more flexible, it still cannot capture bi-modality. In this work, we propose a mixture of CMPs for capturing a wide range of truncated count data, which can exhibit unimodal behavior (with equi-, under- or over-dispersion) as well as bimodal behaviour. We present methods for estimating the parameters of a mixture of two CMP distributions using an EM approach. Our approach introduces a special two-step optimization within the M-step to estimate multiple parameters. The methods are illustrated using simulated and real data.
Modeling Bimodal Discrete Data Using Conway-Maxwell-Poisson Mixture Models

pragya sur,galit shmueli,smarajit bose,paromita dubey

DOI: https://doi.org/10.1080/07350015.2014.949343

2015-01-01

Abstract:Bimodal truncated count distributions are frequently observed in aggregate survey data and in user ratings when respondents are mixed in their opinion. They also arise in censored count data, where the highest category might create an additional mode. Modeling bimodal behavior in discrete data is useful for various purposes, from comparing shapes of different samples (or survey questions) to predicting future ratings by new raters. The Poisson distribution is the most common distribution for fitting count data and can be modified to achieve mixtures of truncated Poisson distributions. However, it is suitable only for modeling equidispersed distributions and is limited in its ability to capture bimodality. The Conway-Maxwell-Poisson (CMP) distribution is a two-parameter generalization of the Poisson distribution that allows for over- and underdispersion. In this work, we propose a mixture of CMPs for capturing a wide range of truncated discrete data, which can exhibit unimodal and bimodal behavior. We present methods for estimating the parameters of a mixture of two CMP distributions using an EM approach. Our approach introduces a special two-step optimization within the M step to estimate multiple parameters. We examine computational and theoretical issues. The methods are illustrated for modeling ordered rating data as well as truncated count data, using simulated and real examples.
Robust Bernoulli mixture models for credit portfolio risk

Jonathan Ansari,Eva Lütkebohmert

2024-11-18

Abstract:This paper presents comparison results and establishes risk bounds for credit portfolios within classes of Bernoulli mixture models, assuming conditionally independent defaults that are stochastically increasing with a common risk factor. We provide simple and interpretable conditions for conditional default probabilities that imply a comparison of credit portfolio losses in convex order. In the case of threshold models, the ranking of portfolio losses is based on a pointwise comparison of the underlying copulas. Our setting includes as special case the well-known Gaussian copula model but allows for general tail dependencies, which are crucial for modeling credit portfolio risks. Moreover, our results extend the classical parameterized models, such as the industry models CreditMetrics and KMV Portfolio Manager, to a robust setting where individual parameters or the copula modeling the dependence structure can be ambiguous. A simulation study and a real data example under model uncertainty offer evidence supporting the effectiveness of our approach.

Risk Management,Mathematical Finance
Estimation and hypothesis testing in finite mixture models

Murray Aitkin, Donald B Rubin

1985-09-01

Abstract:Finite mixture models are a useful class of models for application to data. When sample sizes are not large and the number of underlying densities is in question, likelihood ratio tests based on joint maximum likelihood estimation of the mixing parameter, λ, and the parameter of the underlying densities, θ, are problematical. Our approach places a prior distribution on λ and estimates θ by maximizing the likelihood of the data given θ with λ integrated out. Advantages of this approach, computational issues using the EM algorithm and directions for further work are discussed. The technique is applied to two examples.
Mixed Poisson families with real-valued mixing distributions

F. William Townes

2024-07-25

Abstract:Mixed Poisson distributions provide a flexible approach to the analysis of count data with overdispersion, zero inflation, or heavy tails. Since the Poisson mean must be nonnegative, the mixing distribution is typically assumed to have nonnegative support. We show this assumption is unnecessary and real-valued mixing distributions are also possible. Informally, the mixing distribution merely needs to have a light (subexponential) left tail and a small amount of probability mass on negative values. We provide several concrete examples, including the mixed Poisson-extreme stable family, where the mixing distribution has a power law tail.

Probability
Mixed memory markov models: Decomposing complex stochastic processes as mixtures of simpler ones

Lawrence K Saul, Michael I Jordan

1999-10-01

Abstract:We study Markov models whose state spaces arise from the Cartesian product of two or more discrete random variables. We show how to parameterize the transition matrices of these models as a convex combination—or mixture—of simpler dynamical models. The parameters in these models admit a simple probabilistic interpretation and can be fitted iteratively by an Expectation-Maximization (EM) procedure. We derive a set of generalized Baum-Welch updates for factorial hidden Markov models that make use of this parameterization. We also describe a simple iterative procedure for approximately computing the statistics of the hidden states. Throughout, we give examples where mixed memory models provide a useful representation of complex stochastic processes.
Choice of mixture Poisson models based on Extreme value theory

Samuel Valiquette,Frédéric Mortier,Jean Peyhardi,Gwladys Toulemonde

DOI: https://doi.org/10.48550/arXiv.2204.10649

2022-04-22

Abstract:Count data are omnipresent in many applied fields, often with overdispersion due to an excess of zeroes or extreme values. With mixtures of Poisson distributions representing an elegant and appealing modelling strategy, we focus here on the challenging problem of identifying a suitable mixing distribution and study how extreme value theory can be used. We propose an original strategy to select the most appropriate candidate among three categories: Fr{é}chet, Gumbel and pseudo-Gumbel. Such an approach is presented with the aid of a decision tree and evaluated with numerical simulations.

Statistics Theory
EM Estimation for the Poisson-Inverse Gamma Regression Model with Varying Dispersion: An Application to Insurance Ratemaking

George Tzougas

DOI: https://doi.org/10.3390/risks8030097

2020-09-11

Risks

Abstract:This article presents the Poisson-Inverse Gamma regression model with varying dispersion for approximating heavy-tailed and overdispersed claim counts. Our main contribution is that we develop an Expectation-Maximization (EM) type algorithm for maximum likelihood (ML) estimation of the Poisson-Inverse Gamma regression model with varying dispersion. The empirical analysis examines a portfolio of motor insurance data in order to investigate the efficiency of the proposed algorithm. Finally, both the a priori and a posteriori, or Bonus-Malus, premium rates that are determined by the Poisson-Inverse Gamma model are compared to those that result from the classic Negative Binomial Type I and the Poisson-Inverse Gaussian distributions with regression structures for their mean and dispersion parameters.
Statistical inference for a mixture of Markov jump processes

Halina Frydman,Budhi Surya

DOI: https://doi.org/10.48550/arXiv.2103.02755

2022-04-11

Abstract:We estimate a general mixture of Markov jump processes. The key novel feature of the proposed mixture is that the transition intensity matrices of the Markov processes comprising the mixture are entirely unconstrained. The Markov processes are mixed with distributions that depend on the initial state of the mixture process. The new mixture is estimated from its continuously observed realizations using the EM algorithm, which provides the maximum likelihood (ML) estimates of the mixture's parameters. We derive the asymptotic properties of the ML estimators. To obtain estimated standard errors of the ML estimates of the mixture's parameters, an explicit form of the observed Fisher information matrix is derived. In its new form, the information matrix simplifies the conditional expectation of outer product of the complete-data score function in the Louis (1982) general matrix formula for the observed Fisher information matrix. Simulation study verifies the estimates' accuracy and confirms the consistency and asymptotic normality of the estimators. The developed methods are applied to a medical dataset, for which the likelihood ratio test rejects the constrained mixture in favor of the proposed unconstrained one. This application exemplifies the usefulness of a new unconstrained mixture for identification and characterization of homogeneous subpopulations in a heterogeneous population.

Methodology
EM estimation for bivariate mixed poisson INAR(1) claim count regression models with correlated random effects

Zezhun Chen,Angelos Dassios,George Tzougas

DOI: https://doi.org/10.1007/s13385-023-00351-7

2023-06-07

European Actuarial Journal

Abstract:This article considers bivariate mixed Poisson INAR(1) regression models with correlated random effects for modelling correlations of different signs and magnitude among time series of different types of claim counts. This is the first time that the proposed family of INAR(1) models is used in a statistical or actuarial context. For expository purposes, the bivariate mixed Poisson INAR(1) claim count regression models with correlated Lognormal and Gamma random effects paired via a Gaussian copula are presented as competitive alternatives to the classical bivariate Negative Binomial INAR(1) claim count regression model which only allows for positive dependence between the time series of claim count responses. Our main achievement is that we develop novel alternative Expectation-Maximization type algorithms for maximum likelihood estimation of the parameters of the models which are demonstrated to perform satisfactory when the models are fitted to Local Government Property Insurance Fund data from the state of Wisconsin.

Experience Rating with Poisson Mixtures

About the Justification of Experience Rating: Bonus Malus System and a new Poisson Mixture Model

Bayesian estimation and prediction for certain mixtures

Bivariate phase-type distributions for experience rating in disability insurance

Erlang mixture modeling for Poisson process intensities

Mixture Models With a Prior on the Number of Components

mixture Poisson exponential–inverse Gaussian regression model

Bivariate Gamma Mixture of Experts Models for Joint Insurance Claims Modeling

A Mixture Model of Proportional Reversed Hazard Rate

Maximum-penalized-likelihood estimation for independent and Markov-dependent mixture models.

Bayesian Ratemaking with Common Effects Modeled by Mixture of Polya Tree Processes

Fitting Com-Poisson Mixtures to Bimodal Count Data

Modeling Bimodal Discrete Data Using Conway-Maxwell-Poisson Mixture Models

Robust Bernoulli mixture models for credit portfolio risk

Estimation and hypothesis testing in finite mixture models

Mixed Poisson families with real-valued mixing distributions

Mixed memory markov models: Decomposing complex stochastic processes as mixtures of simpler ones

Choice of mixture Poisson models based on Extreme value theory

EM Estimation for the Poisson-Inverse Gamma Regression Model with Varying Dispersion: An Application to Insurance Ratemaking

Statistical inference for a mixture of Markov jump processes

EM estimation for bivariate mixed poisson INAR(1) claim count regression models with correlated random effects