The implementation of hedging strategies for variable annuity products requires the calculation of market risk sensitivities (or “Greeks”). The complex, path-dependent nature of these products means that these sensitivities are typically estimated by Monte Carlo methods. Standard market practice is to use a “bump and revalue” method in which sensitivities are approximated by finite differences. As well as requiring multiple valuations of the product, this approach is often unreliable for higher-order Greeks, such as gamma, and alternative pathwise (PW) and likelihood-ratio estimators should be preferred. This paper considers a stylized guaranteed minimum withdrawal benefit product in which the reference equity index follows a Heston stochastic volatility model in a stochastic interest rate environment. The complete set of first-order sensitivities with respect to index value, volatility and interest rate and the most important second-order sensitivities are calculated using PW, likelihood-ratio and mixed methods. It is observed that the PW method delivers the best estimates of first-order sensitivities while mixed estimation methods deliver considerably more accurate estimates of second-order sensitivities; moreover there are significant computational gains involved in using PW and mixed estimators rather than simple BnR estimators when many Greeks have to be calculated.

Age and period are the most widely used parameters for forecasting mortality rates. Empirical mortality rate differences in multiple populations often show strong geometric patterns on the two-dimensional age–period plane. The idea of this paper is to take these geometric patterns as the starting point for the development of forecasts. A parametric approach is presented and discussed which uses simple techniques from spatial statistics. The proposed model is statistically parsimonious and yields forecasts that are consistent with the historical data and coherent for multiple populations.

We analyze optimal consumption in the life cycle model by introducing life and pension insurance contracts. The model contains a credit market with biometric risk, and market risk via risky securities. This idealized framework enables us to clarify important aspects of life insurance and pension contracts. We find optimal pension plans and life insurance contracts where the benefits are state dependent. We compare these solutions both to the ones of standard actuarial theory, and to policies offered in practice. Implications of this include what role the insurance industry may play to improve welfare. The relationship between substitution of consumption and risk aversion is highlighted in the presence of a consumption puzzle. One problem related portfolio choice is discussed the horizon problem. Finally, we present some comments on longevity risk and cohort risk.

We develop quantile functions from regression models in order to derive risk margin and to evaluate capital in non-life insurance applications. By utilizing the entire range of conditional quantile functions, especially higher quantile levels, we detail how quantile regression is capable of providing an accurate estimation of risk margin and an overview of implied capital based on the historical volatility of a general insurers loss portfolio. Two modeling frameworks are considered based around parametric and non-parametric regression models which we develop specifically in this insurance setting. In the parametric framework, quantile functions are derived using several distributions including the flexible generalized beta (GB2) distribution family, asymmetric Laplace (AL) distribution and power-Pareto (PP) distribution. In these parametric model based quantile regressions, we detail two basic formulations. The first involves embedding the quantile regression loss function from the nonparameteric setting into the argument of the kernel of a parametric data likelihood model, this is well known to naturally lead to the AL parametric model case. The second formulation we utilize in the parametric setting adopts an alternative quantile regression formulation in which we assume a structural expression for the regression trend and volatility functions which act to modify a base quantile function in order to produce the conditional data quantile function. This second approach allows a range of flexible parametric models to be considered with different tail behaviors. We demonstrate how to perform estimation of the resulting parametric models under a Bayesian regression framework. To achieve this, we design Markov chain Monte Carlo (MCMC) sampling strategies for the resulting Bayesian posterior quantile regression models. In the non-parametric framework, we construct quantile functions by minimizing an asymmetrically weighted loss function and estimate the parameters under the AL proxy distribution to resemble the minimization process. This quantile regression model is contrasted to the parametric AL mean regression model and both are expressed as a scale mixture of uniform distributions to facilitate efficient implementation. The models are extended to adopt dynamic mean, variance and skewness and applied to analyze two real loss reserve data sets to perform inference and discuss interesting features of quantile regression for risk margin calculations.

Pricing approaches for long-tail quota shares are often based on the chain ladder method. Apart from IBNR calculation, common pricing methods require volume measures for accident years in the observation period, and for the quotation period. In practice, in most cases restated premiums are used as the volume measures. The prediction error of the chain ladder method is an important part of the prediction uncertainty of these pricing approaches. There are, however, two sources of uncertainty that are not addressed by the chain ladder model: the stochastic volatility of the claims in the first development year; and the restatement uncertainty, the risk that the restated premium is not a good volume measure. We extend Mack's chain ladder model to cover these two sources of uncertainty, and calculate the mean-squared error of chain ladder pricing approaches with arbitrary weights for the accident years in the observation period. Then we focus on the problem of finding optimal weights for the accident years. First, we assume that the parameters for restatement uncertainty are given, and provide recursion formulas to calculate approximately-optimal weights. Second, we describe a maximum likelihood approach that can be used to estimate the restatement uncertainty.

Various types of structures that enable a group of individuals to pool their mortality risk have been proposed in the literature. Collectively, the structures are called pooled annuity funds. Since the pooled annuity funds propose different methods of pooling mortality risk, we investigate the connections between them and find that they are genuinely different for a finite heterogeneous membership profile. We discuss the importance of actuarial fairness, defined as the expected benefits equalling the contributions for each member, in the context of pooling mortality risk and comment on whether actuarial unfairness can be seen as solidarity between members. We show that, with a finite number of members in the fund, the group self-annuitization scheme is not actuarially fair: some members subsidize the other members. The implication is that the members who are subsidizing the others may obtain a higher expected benefit by joining a fund with a more favorable membership profile. However, we find that the subsidies are financially significant only for very small or highly heterogeneous membership profiles.

An interesting class of stochastic claims reserving methods is given by the models with conditionally independent loss increments (CILI), where the incremental losses are conditionally independent given a risk parameter Θi,j depending on both the accident year i and the development year j. The Bühlmann–Straub credibility reserving (BSCR) model is a particular case of a CILI model where the risk parameter is only depending on i. We consider CILI models with additive diagonal risk (ADR), where the risk parameter is given by the sum of two components, one depending on the accident year i and the other depending on the calendar year t = i + j. The model can be viewed as an extension of the BSCR model including random diagonal effects, which are often declared to be important in loss reserving but rarely are specifically modeled. We show that the ADR model is tractable in closed form, providing credibility formulae for the reserve and the mean square error of prediction (MSEP). We also derive unbiased estimators for the variance parameters which extend the classical Bühlmann–Straub estimators. The results are illustrated by a numerical example and the estimators are tested by simulation. We find that the inclusion of random diagonal effects can be significant for the reserve estimates and, especially, for the evaluation of the MSEP. The paper is written with the purpose of illustrating the role of stochastic diagonal effects. To isolate these effects, we assume that the development pattern is given. In particular, our MSEP values do not include the uncertainty due to the estimation of the development pattern.

This paper addresses the financing of public pensions in a stochastic environment. Traditionally, funded and unfunded pension schemes have been viewed as opposite solutions for the first pillar of public pensions. However, more recently countries as Sweden and Poland have explored mixed solutions that combine pay-as-you-go (PAYG) with funding mechanisms. The aims of this paper are to examine the rationality of such a combination using portfolio theory arguments and to find the optimal split of the contributions between the two systems. We first introduce the classical deterministic model leading to the well-known Samuelson–Aaron rule according to which diversification is never optimal. We then introduce different stochastic models in which the main processes (wage growth, population growth, financial rate of return) are random. In particular, we obtain conditions on parameters to justify diversification and explicit optimal sharing between PAYG and funding. We also introduce the possibility of investing in several financial assets and explore the impact of introducing systematic longevity risk.

The literature on Bayesian chain ladder models is surveyed. Both Mack and cross-classified forms of the chain ladder are considered. Both cases are examined in the context of error terms distributed according to a member of the exponential dispersion family. Tweedie and over-dispersed Poisson errors follow as special cases. Bayesian cross-classified chain ladder models may randomise row, column or diagonal parameters. Column and diagonal randomisation has been largely absent from the literature until recently. The present paper allows randomisation of row and column parameters. The Bayes estimator, the linear Bayes (credibility) estimator, and the MAP estimator are shown to be identical in the Mack case, and in the cross-classified case provided that the error terms are Tweedie distributed. In the Mack case the variance structure is generalised considerably from the existing literature. In the cross-classified case the model structure differs somewhat from the existing literature, and a comparison is made between the two. MAP estimators for the cross-classified case are often given by implicit equations that require numerical solution. Recursive formulas are given for these in the general case of error terms from the exponential dispersion family. The connection between the cross-classified case and Bornhuetter-Ferguson prediction is explored.

Providers of motor extended warranties with limited usage often face difficulty evaluating the impact of usage limits on warranty price because of incomplete usage data. To address this problem, this paper employs a non-parametric interval-censored survival model of time to accumulate a specific usage. This is used to develop an estimator of the probability that a provider is on risk at a specific time in service. The resulting pricing model is applied to a truck extended warranty case study. The case study demonstrates that interval-censored survival models are ideal for use in pricing motor extended warranties with limited usage cover. The results also suggest that employing a usage rate distribution to forecast the number of vehicles on risk can be misleading, especially on an extended warranty with a relatively high usage limit.

One of the most critical problems in property/casualty insurance is to determine an appropriate reserve for incurred but unpaid losses. These provisions generally comprise most of the liabilities of a non-life insurance company. The global provisions are often determined under an assumption of independence between the lines of business. Recently, Shi and Frees (2011) proposed to put dependence between lines of business with a copula that captures dependence between two cells of two different runoff triangles. In this paper, we propose to generalize this model in two steps. First, by using an idea proposed by Barnett and Zehnwirth (1998), we will suppose a dependence between all the observations that belong to the same calendar year (CY) for each line of business. Thereafter, we will then suppose another dependence structure that links the CYs of different lines of business. This model is done by using hierarchical Archimedean copulas. We show that the model provides more flexibility than existing models, and offers a better, more realistic and more intuitive interpretation of the dependence between the lines of business. For illustration, the model is applied to a dataset from a major US property-casualty insurer, where a bootstrap method is proposed to estimate the distribution of the reserve.

In this article, we consider the problem of computing the expected discounted value of a death benefit, e.g. in Gerber et al. (2012, 2013), in a regime-switching economy. Contrary to their proposed discounted density approach, we adopt the Laplace transform to value the contingent options. By this alternative approach, closed-form expressions for the Laplace transforms of the values of various contingent options, such as call/put options, lookback options, barrier options, dynamic fund protection and the dynamic withdrawal benefits, have been obtained. The value of each contingent option can then be recovered by the numerical Laplace inversion algorithm, and this efficient approach is documented by several numerical illustrations. The strength of our methodology becomes apparent when we tackle the valuations of exotic contingent options in the cases when (1) the contracts have a finite expiry date; (2) when the time-until-death variable is uniformly distributed in accordance with De Moivre's law.

The dual risk model assumes that the surplus of a company decreases at a constant rate over time and grows by means of upward jumps, which occur at random times and sizes. It is said to have applications to companies with economical activities involved in research and development. This model is dual to the well-known Cramér-Lundberg risk model with applications to insurance. Most existing results on the study of the dual model assume that the random waiting times between consecutive gains follow an exponential distribution, as in the classical Cramér-Lundberg risk model. We generalize to other compound renewal risk models where such waiting times are Erlang(n) distributed. Using the roots of the fundamental and the generalized Lundberg's equations, we get expressions for the ruin probability and the Laplace transform of the time of ruin for an arbitrary single gain distribution. Furthermore, we compute expected discounted dividends, as well as higher moments, when the individual common gains follow a Phase-Type, PH(m), distribution. We also perform illustrations working some examples for some particular gain distributions and obtain numerical results.

We study an asset allocation problem for a defined-contribution (DC) pension scheme in its accumulation phase. We assume that the amount contributed to the pension fund by a pension plan member is coupled with the salary income which fluctuates randomly over time and contains both a hedgeable and non-hedgeable risk component. We consider an economy in which macroeconomic risks are existent. We assume that the economy can be in one of I states (regimes) and switches randomly between those states. The state of the economy affects the dynamics of the tradeable risky asset and the contribution process (the salary income of a pension plan member). To model the switching behavior of the economy we use a counting process with stochastic intensities. We find the investment strategy which maximizes the expected exponential utility of the discounted excess wealth over a target payment, e.g. a target lifetime annuity.

In this paper, we compare the error in several approximation methods for the cumulative aggregate claim distribution customarily used in the collective model of insurance theory. In this model, it is usually supposed that a portfolio is at risk for a time period of length t. The occurrences of the claims are governed by a Poisson process of intensity μ so that the number of claims in [0,t] is a Poisson random variable with parameter λ = μ t. Each single claim is an independent replication of the random variable X, representing the claim severity. The aggregate claim or total claim amount process in [0,t] is represented by the random sum of N independent replications of X, whose cumulative distribution function (cdf) is the object of study. Due to its computational complexity, several approximation methods for this cdf have been proposed. In this paper, we consider 15 approximations put forward in the literature that only use information on the lower order moments of the involved distributions. For each approximation, we consider the difference between the true distribution and the approximating one and we propose to use expansions of this difference related to Edgeworth series to measure their accuracy as λ = μ t diverges to infinity. Using these expansions, several statements concerning the quality of approximations for the distribution of the aggregate claim process can find theoretical support. Other statements can be disproved on the same grounds. Finally, we investigate numerically the accuracy of the proposed formulas.

We discuss some properties of a class of multivariate mixed Erlang distributions with different scale parameters and describes various distributional properties related to applications in insurance risk theory. Some representations involving scale mixtures, generalized Esscher transformations, higher-order equilibrium distributions, and residual lifetime distributions are derived. These results allows for the study of stop-loss moments, premium calculation, and the risk allocation problem. Finally, some results concerning minimum and maximum variables are derived and applied to pricing joint life and last survivor policies.

In this paper, a new methodology based on the use of the inverse of the circular tangent function that allows us to add a scale parameter (say α) to an initial survival function is presented. The latter survival function is determined as limiting case when α tends to zero. By choosing as parent the classical Pareto survival function, the Pareto ArcTan (PAT) distribution is obtained. After providing a comprehensive analysis of its statistical properties, theoretical results with reference to insurance are illustrated. Its performance is compared, by means of the well-known Norwegian fire insurance data, with other existing heavy-tailed distributions in the literature such as Pareto, Stoppa, Shifted Lognormal, Inverse Gamma and Fréchet distributions.

In this paper, we present recursive formulae for the ruin probability at or before a certain claim arrival instant for some particular continuous time risk model. The claim number process underlying this risk model is a renewal process with either Erlang or a mixture of exponentials inter-claim times (ICTs). The claim sizes (CSs) are independent and distributed in Erlang's family, i.e., they can have different parameters, which yields a non-homogeneous risk process. We present the corresponding recursive algorithm used to evaluate the above mentioned ruin probability and we illustrate it on several numerical examples in which we vary the model's parameters to assess the impact of the non-homogeneity on the resulting ruin probability.

We demonstrate both analytically and numerically that the existing methods for measuring tail dependence in copulas may sometimes underestimate the extent of extreme co-movements of dependent risks and, therefore, may not always comply with the new paradigm of prudent risk management. This phenomenon holds in the context of both symmetric and asymmetric copulas with and without singularities. As a remedy, we introduce a notion of paths of maximal (tail) dependence and utilize the notion to propose several new indices of tail dependence. The suggested new indices are conservative, conform with the basic concepts of modern quantitative risk management, and are capable of differentiating between distinct risky positions in situations when the existing indices fail to do so.

In this paper we consider an extension to the aggregation of the FGM mixed Erlang risks, proposed by Cossette et al. (2013 Insurance: Mathematics and Economics, 52, 560–572), in which we introduce the Sarmanov distribution to model the dependence structure. For our framework, we demonstrate that the aggregated risk belongs to the class of Erlang mixtures. Following results from S. C. K. Lee and X. S. Lin (2010 North American Actuarial Journal, 14(1) 107–130), G. E. Willmot and X. S. Lin (2011 Applied Stochastic Models in Business and Industry, 27(1) 8–22), analytical expressions of the contribution of each individual risk to the economic capital for the entire portfolio are derived under both the TVaR and the covariance capital allocation principle. By analysing the commonly used dependence measures, we also show that the dependence structure is wide and flexible. Numerical examples and simulation studies illustrate the tractability of our approach.