Comparisons of aggregate claim numbers and amounts: a study of heterogeneity

In this paper, we investigate how the heterogeneity among occurrence probabilities and claim severities affects the aggregate claim numbers and aggregate claim amount for an insurance portfolio. We show that higher heterogeneity (and dependence) among occurrence probabilities results in both smaller aggregate claim numbers and aggregate claim amount in the sense of the mean residual lifetime order. We also prove that as the heterogeneity among the claims increases, the aggregate claim amount increases in the sense of the usual stochastic order when the vector of occurrence probabilities is left tail weakly stochastic arrangement increasing. These theoretical findings are applied to (i) study ordering properties of convolutions of binomial random variables, (ii) provide upper bounds for the mean residual lifetime functions of the aggregate claim numbers and amount, and (iii) compare stop-loss premiums and risk capital of different insurance portfolios.     

Moment generating functions of compound renewal sums with discounted claims

Léveillé & Garrido (2001a, 2001b) have obtained recursive formulas for the moments of compound renewal sums with discounted claims, which incorporate both, Andersen's (1957) generalization of the classical risk model, where the claim number process is an ordinary renewal process, and Taylor's (1979), where the joint effect of the claims cost inflation and investment income on a compound Poisson risk process is considered.

In this paper, assuming certain regularity conditions, we improve the preceding results by examining more deeply the asymptotic and finite time moment generating functions of the discounted aggregate claims process. Examples are given for claim inter-arrival times and claim severity following phase-type distributions, such as the Erlang case.     

Ruin probabilities in multivariate risk models with periodic common shock

We propose a multidimensional risk model where the common shock affecting all classes of insurance business is arriving according to a non-homogeneous periodic Poisson process. In this multivariate setting, we derive upper bounds of Lundberg-type for the probability that ruin occurs in all classes simultaneously using the martingale approach via piecewise deterministic Markov processes theory. These results are numerically illustrated in a bivariate risk model, where the beta-shape periodic claim intensity function is considered. Under the assumption of dependent heavy-tailed claims, asymptotic bounds for the finite-time ruin probabilities associated to three types of ruin in this multivariate framework are investigated.     

Spatial modelling of claim frequency and claim size in non-life insurance

In this paper, models for claim frequency and average claim size in non-life insurance are considered. Both covariates and spatial random effects are included allowing the modelling of a spatial dependency pattern. We assume a Poisson model for the number of claims, while claim size is modelled using a Gamma distribution. However, in contrast to the usual compound Poisson model, we allow for dependencies between claim size and claim frequency. A fully Bayesian approach is followed, parameters are estimated using Markov Chain Monte Carlo (MCMC). The issue of model comparison is thoroughly addressed. Besides the deviance information criterion and the predictive model choice criterion, we suggest the use of proper scoring rules based on the posterior predictive distribution for comparing models. We give an application to a comprehensive data set from a German car insurance company. The inclusion of spatial effects significantly improves the models for both claim frequency and claim size, and also leads to more accurate predictions of the total claim sizes. Further, we detect significant dependencies between the number of claims and claim size. Both spatial and number of claims effects are interpreted and quantified from an actuarial point of view.     

“The Future of Risk Classification in the Age of Predictive DNA-Based Testing”, Donald C. Chambers,January 1999

Abstract

The correlation among multiple lines of business plays an important role in quantifying the uncertainty of loss reserves for insurance portfolios. To accommodate correlation, most multivariate loss-reserving methods focus on the pairwise association between corresponding cells in multiple run-off triangles. However, such practice usually relies on the independence assumption across accident years and ignores the calendar year effects that could affect all open claims simultaneously and induce dependencies among loss triangles. To address this issue, we study a Bayesian log-normal model in the prediction of outstanding claims for dependent lines of business. In addition to the pairwise correlation, our method allows for an explicit examination of the correlation due to common calendar year effects. Further, different specifications of the calendar year trend are considered to reflect valuation actuaries’ prior knowledge of claim development. In a case study, we analyze an insurance portfolio of personal and commercial auto lines from a major U.S. property-casualty insurer. It is shown that the incorporation of calendar year effects improves model fit significantly, though it contributes substantively to the predictive variability. The availability of the realizations of predicted claims permits us to perform a retrospective test, which suggests that extra prediction uncertainty is indispensable in modern risk management practices.     

Insurance ratemaking using a copula-based multivariate Tweedie model

The Tweedie distribution, featured with a mass probability at zero, is a convenient tool for insurance claims modeling and pure premium determination in general insurance. Motivated by the fact that an insurance policy typically provides multiple types of coverage, we propose a copula-based multivariate Tweedie regression for modeling the semi-continuous claims while accommodating the association among different types. The proposed approach also allows for dispersion modeling, resulting in a multivariate version of the double generalized linear model. We demonstrate the application in insurance ratemaking using a portfolio of policyholders of automobile insurance from the state of Massachusetts in the United States.     

On the discounted aggregate claim costs until ruin in dependent Sparre Andersen risk processes

In this paper, a dependent Sparre Andersen risk process in which the joint density of the interclaim time and the resulting claim severity satisfies the factorization as in Willmot and Woo is considered. We study a generalization of the Gerber–Shiu function (i) whose penalty function further depends on the surplus level immediately after the second last claim before ruin; and (ii) which involves the moments of the discounted aggregate claim costs until ruin. The generalized discounted density with a moment-based component proposed in Cheung plays a key role in deriving recursive defective renewal equations. We pay special attention to the case where the marginal distribution of the interclaim times is Coxian, and the required components in the recursion are obtained. A reverse type of dependency structure, where the claim severities follow a combination of exponentials, is also briefly discussed, and this leads to a nice explicit expression for the expected discounted aggregate claims until ruin. Our results are applied to generate some numerical examples involving (i) the covariance of the time of ruin and the discounted aggregate claims until ruin; and (ii) the expectation, variance and third central moment of the discounted aggregate claims until ruin.     

Extremes on the discounted aggregate claims in a time dependent risk model

This paper presents an extension of the classical compound Poisson risk model for which the inter-claim time and the forthcoming claim amount are no longer independent random variables (rv's). Asymptotic tail probabilities for the discounted aggregate claims are presented when the force of interest is constant and the claim amounts are heavy tail distributed rv's. Furthermore, we derive asymptotic finite time ruin probabilities, as well as asymptotic approximations for some common risk measures associated with the discounted aggregate claims. A simulation study is performed in order to validate the results obtained in the free interest risk model.     

Third cumulant for multivariate aggregate claim models

The third cumulant for the aggregated multivariate claims is considered. A formula is presented for the general case when the aggregating variable is independent of the multivariate claims. Two important special cases are considered. In the first one, multivariate skewed normal claims are considered and aggregated by a Poisson variable. The second case is dealing with multivariate asymmetric generalized Laplace and aggregation is made by a negative binomial variable. Due to the invariance property the latter case can be derived directly, leading to the identity involving the cumulant of the claims and the aggregated claims. There is a well-established relation between asymmetric Laplace motion and negative binomial process that corresponds to the invariance principle of the aggregating claims for the generalized asymmetric Laplace distribution. We explore this relation and provide multivariate continuous time version of the results. It is discussed how these results that deal only with dependence in the claim sizes can be used to obtain a formula for the third cumulant for more complex aggregate models of multivariate claims in which the dependence is also in the aggregating variables.     

Ruin under stochastic dependence between premium and claim arrivals

We investigate, focusing on the ruin probability, an adaptation of the Cramér–Lundberg model for the surplus process of an insurance company, in which, conditionally on their intensities, the two mixed Poisson processes governing the arrival times of the premiums and of the claims respectively, are independent. Such a model exhibits a stochastic dependence between the aggregate premium and claim amount processes. An explicit expression for the ruin probability is obtained when the claim and premium sizes are exponentially distributed.     

On a conjecture related to the ruin probability for nonhomogeneous exponentially distributed claims

Recently, some recursive formulas have been obtained for the ruin probability evaluated at or before claim instants for a surplus process under the assumptions that the claim sizes are independent, nonhomogeneous Erlang distributed, and independent of the inter-claim revenues, which are assumed to be independent, identically distributed, following an arbitrary distribution. Based on numerical examples, a conjecture has also been stated relating the order in which the claims arrive to the magnitude of the corresponding ruin probability. In this paper, we prove this conjecture in the particular case when the claims are all exponentially distributed with different parameters.     

On The Expected Discounted Penalty function for Lévy Risk Processes

Abstract

Dufresne et al. (1991) introduced a general risk model defined as the limit of compound Poisson processes. Such a model is either a compound Poisson process itself or a process with an infinite number of small jumps. Later, in a series of now classical papers, the joint distribution of the time of ruin, the surplus before ruin, and the deficit at ruin was studied (Gerber and Shiu 1997, 1998a, 1998b; Gerber and Landry 1998). These works use the classical and the perturbed risk models and hint that the results can be extended to gamma and inverse Gaussian risk processes.

In this paper we work out this extension to a generalized risk model driven by a nondecreasing Lévy process. Unlike the classical case that models the individual claim size distribution and obtains from it the aggregate claims distribution, here the aggregate claims distribution is known in closed form. It is simply the one-dimensional distribution of a subordinator. Embedded in this wide family of risk models we find the gamma, inverse Gaussian, and generalized inverse Gaussian processes. Expressions for the Gerber-Shiu function are given in some of these special cases, and numerical illustrations are provided.     

The Expected Discounted Penalty at Ruin for a Markov-Modulated Risk Process Perturbed by Diffusion

Abstract

A Markov-modulated risk process perturbed by diffusion is considered in this paper. In the model the frequencies and distributions of the claims and the variances of the Wiener process are influenced by an external Markovian environment process with a finite number of states. This model is motivated by the flexibility in modeling the claim arrival process, allowing that periods with very frequent arrivals and ones with very few arrivals may alternate. Given the initial surplus and the initial environment state, systems of integro-differential equations for the expected discounted penalty functions at ruin caused by a claim and oscillation are established, respectively; a generalized Lundberg’s equation is also obtained. In the two-state model, the expected discounted penalty functions at ruin due to a claim and oscillation are derived when both claim amount distributions are from the rational family. As an illustration, the explicit results are obtained for the ruin probability when claim sizes are exponentially distributed. A numerical example also is given for the case that two classes of claims are Erlang(2) distributed and of a mixture of two exponentials.     

Multivariate Cox Hidden Markov models with an application to operational risk

ABSTRACT

Modeling multivariate time-series aggregate losses is an important actuarial topic that is very challenging due to the fact that losses can be serially dependent with heterogeneous dependence structures across loss types and business lines. In this paper, we investigate a flexible class of multivariate Cox Hidden Markov Models for the joint arrival process of loss events. Some of the nice properties possessed by this class of models, such as closed-form expressions, thinning properties and model versatility are discussed in details. We provide the expectation-maximization (EM) algorithm for efficient model calibration. Applying the proposed model to an operational risk dataset, we demonstrate that the model offers sufficient flexibility to capture most characteristics of the observed loss frequencies. By modeling the log-transformed loss severities through mixture of Erlang distributions, we can model the aggregate losses. Finally, out-of-sample testing shows that the proposed model is adequate to predict short-term future operational risk losses.     

Martingale Approach for Moments of Discounted Aggregate Claims

We examine the Laplace transform of the distribution of the shot noise process using the martingale. Applying the piecewise deterministic Markov processes theory and using the relationship between the shot noise process and the accumulated/discounted aggregate claims process, the Laplace transform of the distribution of the accumulated aggregate claims is obtained. Assuming that the claim arrival process follows the Poisson process and claim sizes are assumed to be exponential and mixture of exponential, we derive the explicit expressions of the actuarial net premiums and variances of the discounted aggregate claims, which are the annuities paid continuously. Numerical examples are also provided based on them.     

Compound Poisson Model with Covariates

Abstract

Pet insurance in North America continues to be a growing industry. Unlike in Europe, where some countries have as much as 50% of the pet population insured, very few pets in North America are insured. Pricing practices in the past have relied on market share objectives more so than on actual experience. Pricing still continues to be performed on this basis with little consideration for actuarial principles and techniques. Developments of mortality and morbidity models to be used in the pricing model and new product development are essential for pet insurance. This paper examines insurance claims as experienced in the Canadian market. The time-to-event data are investigated using the Cox’s proportional hazards model. The claim number follows a nonhomogenous Poisson process with covariates. The claim size random variable is assumed to follow a lognormal distribution. These two models work well for aggregate claims with covariates. The first three central moments of the aggregate claims for one insured animal, as well as for a block of insured animals, are derived. We illustrate the models using data collected over an eight-year period.     

Unallocated loss adjustment expense reserving

In non-life insurance, the provision for outstanding claims (the claims reserve) should include future loss adjustment expenses, i.e. administrative expenses to settle the claims, and therefore we have to estimate the expected Unallocated Loss Adjustment Expenses (ULAE) – expenses that are not attributable to individual claims, such as salaries at the claims handling department. The ULAE reserve has received little attention from European actuaries in the literature, supposedly because of the lack of detailed data for estimation and evaluation. Having good estimation procedures will, however, become even more important with the introduction of the Solvency II regulations, which require unbiased estimation of future cash flows for all expenses. We present a model for ULAE at the individual claim level that includes both fixed and variable costs. This model leads to an estimate of the ULAE reserve at the aggregate (line-of-business) level, as demonstrated in a numerical example from a Swedish non-life insurer.     

Longitudinal modeling of insurance claim counts using jitters

Modeling insurance claim counts is a critical component in the ratemaking process for property and casualty insurance. This article explores the usefulness of copulas to model the number of insurance claims for an individual policyholder within a longitudinal context. To address the limitations of copulas commonly attributed to multivariate discrete data, we adopt a ‘jittering’ method to the claim counts which has the effect of continuitizing the data. Elliptical copulas are proposed to accommodate the intertemporal nature of the ‘jittered’ claim counts and the unobservable subject-specific heterogeneity on the frequency of claims. Observable subject-specific effects are accounted in the model by using available covariate information through a regression model. The predictive distribution together with the corresponding credibility of claim frequency can be derived from the model for ratemaking and risk classification purposes. For empirical illustration, we analyze an unbalanced longitudinal dataset of claim counts observed from a portfolio of automobile insurance policies of a general insurer in Singapore. We further establish the validity of the calibrated copula model, and demonstrate that the copula with ‘jittering’ method outperforms standard count regression models.     

Exact solutions for ruin probability in the presence of an absorbing upper barrier

Abstract

An explicit solution for the probability of ruin in the presence of an absorbing upper barrier was developed by Segerdahl (1970) for the particular case in which both the interoccurrence times between successive claims and the single claim amounts follow an exponential distribution with unit mean. In this paper we show that his method of solution may be extended to produce explicit solutions for two more general types of single claim amount distribution. These are the gamma distribution, denoted γ(a), where a is an integer, and the mixed exponential distribution. Comparisons are drawn between this approach when the upper barrier tends to infinity, and the classical solution for ruin probability in these particular cases given in Cramér (1955).     

Mean and dispersion modelling for policy claims costs

A model for the statistical analysis of the total amount of insurance paid out on a policy is developed and applied. The model simultaneously deals with the number of claims (zero or more) and the amount of each claim. The number of claims is from a Poisson-based discrete distribution. Individual claim sizes are from a continuous right skewed distribution. The resulting distribution of total claim size is a mixed discrete-continuous model, with positive probability of a zero claim. The means and dispersions of the claim frequency and claim size distribution are modeled in terms of risk factors. The model is applied to a car insurance data set.