To evaluate a large portfolio of variable annuity (VA) contracts, many insurance companies rely on Monte Carlo simulation, which is computationally intensive. To address this computational challenge, machine learning techniques have been adopted in recent years to estimate the fair market values (FMVs) of a large number of contracts. It is shown that bootstrapped aggregation (bagging), one of the most popular machine learning algorithms, performs well in valuing VA contracts using related attributes. In this article, we highlight the presence of prediction bias of bagging and use the bias-corrected (BC) bagging approach to reduce the bias and thus improve the predictive performance. Experimental results demonstrate the effectiveness of BC bagging as compared with bagging, boosting, and model points in terms of prediction accuracy.

Regression modelling involving heavy-tailed response distributions, which have heavier tails than the exponential distribution, has become increasingly popular in many insurance settings including non-life insurance. Mixed Exponential models can be considered as a natural choice for the distribution of heavy-tailed claim sizes since their tails are not exponentially bounded. This paper is concerned with introducing a general family of mixed Exponential regression models with varying dispersion which can efficiently capture the tail behaviour of losses. Our main achievement is that we present an Expectation-Maximization (EM)-type algorithm which can facilitate maximum likelihood (ML) estimation for our class of mixed Exponential models which allows for regression specifications for both the mean and dispersion parameters. Finally, a real data application based on motor insurance data is given to illustrate the versatility of the proposed EM-type algorithm.

The transition from defined benefit to defined contribution (DC) pension schemes has increased the interest in target annuitization funds that aim to fund a minimum level of retirement income. Prior literature has studied the optimal investment strategies for DC funds that provide minimum guarantees, but far less attention has been given to portfolio insurance strategies for DC pension funds focusing on retirement income targets. We evaluate the performance of option-based and constant proportion portfolio insurance strategies for a DC fund that targets a minimum level of inflation-protected annuity income at retirement. We show how the portfolio allocation to an equity fund varies depending on the member’s age upon joining the fund, displaying a downward trend through time for members joining the fund before ages in the mid-30s. We demonstrate how both portfolio insurance strategies provide strong protection against downside equity risk in financing a minimum level of retirement income. The option-based strategy generally leads to higher accumulated savings at retirement, whereas the constant proportion strategy provides better downside risk protection robust to equity market jumps/volatilities.

The existence of long memory in mortality data improves the understandings of features of mortality data and provides a new approach for establishing mortality models. The findings of long-memory phenomena in mortality data motivate us to develop new mortality models by extending the Lee–Carter (LC) model to death counts and incorporating long-memory model structure. Furthermore, there are no identification issues arising in the proposed model class. Hence, the constraints which cause many computational issues in LC models are removed. The models are applied to analyse mortality death count data sets from three different countries divided according to genders. Bayesian inference with various selection criteria is applied to perform the model parameter estimation and mortality rate forecasting. Results show that multivariate long-memory mortality model with long-memory cohort effect model outperforms multivariate extended LC cohort model in both in-sample fitting and out-sample forecast. Increasing the accuracy of forecasting of mortality rates and improving the projection of life expectancy is an important consideration for insurance companies and governments since misleading predictions may result in insufficient funds for retirement and pension plans.

A variable annuity (VA) is an equity-linked annuity that provides investment guarantees to its policyholder and its contributions are normally invested in multiple underlying assets (e.g., mutual funds), which exposes VA liability to significant market risks. Hedging the market risks is therefore crucial in risk managing a VA portfolio as the VA guarantees are long-dated liabilities that may span decades. In order to hedge the VA liability, the issuing insurance company would need to construct a hedging portfolio consisting of the underlying assets whose positions are often determined by the liability Greeks such as partial dollar Deltas. Usually, these quantities are calculated via nested simulation approach. For insurance companies that manage large VA portfolios (e.g., 100k+ policies), calculating those quantities is extremely time-consuming or even prohibitive due to the complexity of the guarantee payoffs and the stochastic-on-stochastic nature of the nested simulation algorithm. In this paper, we extend the surrogate model-assisted nest simulation approach in Lin and Yang [(2020) Insurance: Mathematics and Economics, 91, 85–103] to efficiently calculate the total VA liability and the partial dollar Deltas for large VA portfolios with multiple underlying assets. In our proposed algorithm, the nested simulation is run using small sets of selected representative policies and representative outer loops. As a result, the computing time is substantially reduced. The computational advantage of the proposed algorithm and the importance of dynamic hedging are further illustrated through a profit and loss (P&L) analysis for a large synthetic VA portfolio. Moreover, the robustness of the performance of the proposed algorithm is tested with multiple simulation runs. Numerical results show that the proposed algorithm is able to accurately approximate different quantities of interest and the performance is robust with respect to different sets of parameter inputs. Finally, we show how our approach could be extended to potentially incorporate stochastic interest rates and estimate other Greeks such as Rho.

A generalization of range-value-at-risk (RVaR) and tail-value-at-risk (TVaR) for d-dimensional distribution functions is introduced. Properties of these new risk measures are studied and illustrated. We provide special cases, applications and a comparison with traditional univariate and multivariate versions of the TVaR and RVaR.

This paper provides a toolbox for the credibility analysis of frequency risks, with allowance for the seniority of claims and of risk exposure. We use Poisson models with dynamic and second-order stationary random effects that ensure nonnegative credibilities per period. We specify classes of autocovariance functions that are compatible with positive random effects and that entail nonnegative credibilities regardless of the risk exposure. Random effects with nonnegative generalized partial autocorrelations are shown to imply nonnegative credibilities. This holds for ARFIMA(0, d, 0) models. The AR(p) time series that ensure nonnegative credibilities are specified from their precision matrices. The compatibility of these semiparametric models with log-Gaussian random effects is verified. Gaussian sequences with ARFIMA(0, d, 0) specifications, which are then exponentiated entrywise, provide positive random effects that also imply nonnegative credibilities. Dynamic random effects applied to Poisson distributions are retained as products of two uncorrelated and positive components: the first is time-invariant, whereas the autocovariance function of the second vanishes at infinity and ensures nonnegative credibilities. The limit credibility is related to the three levels for the length of the memory in the random effects. The limit credibility is less than one in the short memory case, and a formula is provided.

This paper studies bilateral risk sharing under no aggregate uncertainty, where one agent has Expected-Utility preferences and the other agent has Rank-dependent utility preferences with a general probability distortion function. We impose exogenous constraints on the risk exposure for both agents, and we allow for any type or level of belief heterogeneity. We show that Pareto-optimal risk-sharing contracts can be obtained via a constrained utility maximization under a participation constraint of the other agent. This allows us to give an explicit characterization of optimal risk-sharing contracts. In particular, we show that an optimal risk-sharing contract contains allocations that are monotone functions of the likelihood ratio, where the latter is obtained from Lebesgue’s Decomposition Theorem.

Interest rate is one of the main risks for the liability of the variable annuity (VA) due to its long maturity. However, most existing studies on the risk measures of the VA assume a constant interest rate. In this paper, we propose an efficient two-dimensional willow tree method to compute the liability distribution of the VA with the joint dynamics of the mutual fund and interest rate. The risk measures can then be computed by the backward induction on the tree structure. We also analyze the sensitivity and impact on the risk measures with regard to the market model parameters, contract attributes, and monetary policy changes. It illustrates that the liability of the VA is determined by the long-term interest rate whose increment leads to a decrease in the liability. The positive correlation between the interest rate and mutual fund generates a fat-tailed liability distribution. Moreover, the monetary policy change has a bigger impact on the long-term VAs than the short-term contracts.

This paper presents analytical representations for an optimal insurance contract under distortion risk measure and in the presence of model uncertainty. We incorporate ambiguity aversion and distortion risk measure through the model of Robert and Therond [(2014) ASTIN Bulletin: The Journal of the IAA, 44(2), 277–302.] as per the framework of Klibanoff et al. [(2005) A smooth model of decision making under ambiguity. Econometrica, 73(6), 1849–1892.]. Explicit optimal insurance indemnity functions are derived when the decision maker (DM) applies Value-at-Risk as risk measure and is ambiguous about the loss distribution. Our results show that: (1) under model uncertainty, ambiguity aversion results in a distorted probability distribution over the set of possible models with a bias in favor of the model which yields a larger risk; (2) a more ambiguity-averse DM would demand more insurance coverage; (3) for a given budget, uncertainties about the loss distribution result in higher risk level for the DM.

Modeling taxation of Variable Annuities has been frequently neglected, but accounting for it can significantly improve the explanation of the withdrawal dynamics and lead to a better modeling of the financial cost of these insurance products. The importance of including a model for taxation has first been observed by Moenig and Bauer (2016) while considering a Guaranteed Minimum Withdrawal Benefit (GMWB) Variable Annuity. In particular, they consider the simple Black–Scholes dynamics to describe the underlying security. Nevertheless, GMWB are long-term products, and thus accounting for stochastic interest rate has relevant effects on both the financial evaluation and the policyholder behavior, as observed by Goudenège et al. (2018). In this paper, we investigate the outcomes of these two elements together on GMWB evaluation. To this aim, we develop a numerical framework which allows one to efficiently compute the fair value of a policy. Numerical results show that accounting for both taxation and stochastic interest rate has a determinant impact on the withdrawal strategy and on the cost of GMWB contracts. In addition, it can explain why these products are so popular with people looking for a protected form of investment for retirement.

We study a weighted comonotonic risk-sharing problem among multiple agents with distortion risk measures under heterogeneous beliefs. The explicit forms of optimal allocations are obtained, which are Pareto-optimal. A necessary and sufficient condition is given to ensure the uniqueness of the optimal allocation, and sufficient conditions are given to obtain an optimal allocation of the form of excess of loss or full insurance. The optimal allocation may satisfy individual rationality depending on the choice of the weight. When the distortion risk measure is value at risk or tail value at risk, an optimal allocation is generally of the excess-of-loss form. The numerical examples suggest that a risk is more likely to be shared among agents with heterogeneous beliefs, and the introduction of the weight enables us to prioritize some agents as part of a group sharing a risk.

We present a method for constructing and interpreting weighted premium principles. The method is based on modifying the underlying risk distribution in such a way that the risk-adjusted expected value (or premium) is greater than the expected value of some conveniently chosen function of claims, which defines the insurer’s perception of the risk. Under some assumptions on the function of claims, the method produces distortion premium principles. We provide several examples under different assumptions on the claim arrival process and different functions of claims, including record claims and kth record claims.

In this study, we propose new risk measures from a regulator’s perspective on the regulatory capital requirements. The proposed risk measures possess many desired properties, including monotonicity, translation-invariance, positive homogeneity, subadditivity, nonnegative loading, and stop-loss order preserving. The new risk measures not only generalize the existing, well-known risk measures in the literature, including the Dutch, tail value-at-risk (TVaR), and expectile measures, but also provide new approaches to generate feasible and practical coherent risk measures. As examples of the new risk measures, TVaR-type generalized expectiles are investigated in detail. In particular, we present the dual and Kusuoka representations of the TVaR-type generalized expectiles and discuss their robustness with respect to the Wasserstein distance.

We consider the conditional mean risk allocation for an insurance pool, as defined by Denuit and Dhaene (2012). Precisely, we study the asymptotic behavior of the respective relative contributions of the participants as the total loss of the pool tends to infinity. The numerical illustration in Denuit (2019) suggests that the application of the conditional mean risk sharing rule may produce a linear sharing in the tail of the total loss distribution. This paper studies the validity of this empirical finding in the class of compound Panjer–Katz sums consisting of compound Binomial, compound Poisson, and compound Negative Binomial sums with either Gamma or Pareto severities. It is demonstrated that such a behavior does not hold in general since one term may dominate the other ones conditional of large total loss.