We propose a new neighbouring prediction model for mortality forecasting. For each mortality rate at age x in year t, mx,t, we construct an image of neighbourhood mortality data around mx,t, that is, Ꜫmx,t (x1, x2, s), which includes mortality information for ages in [x-x1, x+x2], lagging k years (1 ≤ k ≤ s). Combined with the deep learning model – convolutional neural network, this framework is able to capture the intricate nonlinear structure in the mortality data: the neighbourhood effect, which can go beyond the directions of period, age, and cohort as in classic mortality models. By performing an extensive empirical analysis on all the 41 countries and regions in the Human Mortality Database, we find that the proposed models achieve superior forecasting performance. This framework can be further enhanced to capture the patterns and interactions between multiple populations.

We construct a double common factor model for projecting the mortality of a population using as a reference the minimum death rate at each age among a large number of countries. In particular, the female and male minimum death rates, described as best-performance or best-practice rates, are first modelled by a common factor model structure with both common and sex-specific parameters. The differences between the death rates of the population under study and the best-performance rates are then modelled by another common factor model structure. An important result of using our proposed model is that the projected death rates of the population being considered are coherent with the projected best-performance rates in the long term, the latter of which serves as a very useful reference for the projection based on the collective experience of multiple countries. Our out-of-sample analysis shows that the new model has potential to outperform some conventional approaches in mortality projection.

It has become standard practice in the non-life insurance industry to employ generalized linear models (GLMs) for insurance pricing. However, these GLMs traditionally work only with a priori characteristics of policyholders, while nowadays we increasingly have a posteriori information of individual customers available across multiple product categories. In this paper, we therefore develop a framework to capture this a posteriori information over several product lines using a dynamic claim score. More specifically, we extend the bonus-malus-panel model of Boucher and Inoussa (2014) and Boucher and Pigeon (2018) to include claim scores from other product categories and to allow for nonlinear effects of these scores. The application of the proposed multi-product framework to a Dutch property and casualty insurance portfolio shows that customers’ individual claims experience can have a significant impact on the risk classification. Moreover, it indicates that considerably more profits can be gained by accounting for their multi-product claims experience.

A machine learning approach to zero-inflated Poisson (ZIP) regression is introduced to address common difficulty arising from imbalanced financial data. The suggested ZIP can be interpreted as an adaptive weight adjustment procedure that removes the need for post-modeling re-calibration and results in a substantial enhancement of predictive accuracy. Notwithstanding the increased complexity due to the expanded parameter set, we utilize a cyclic coordinate descent optimization to implement the ZIP regression, with adjustments made to address saddle points. We also study how various approaches alleviate the potential drawbacks of incomplete exposures in insurance applications. The procedure is tested on real-life data. We demonstrate a significant improvement in performance relative to other popular alternatives, which justifies our modeling techniques.

Using telematics technology, insurers are able to capture a wide range of data to better decode driver behavior, such as distance traveled and how drivers brake, accelerate, or make turns. Such additional information also helps insurers improve risk assessments for usage-based insurance, a recent industry innovation. In this article, we explore the integration of telematics information into a classification model to determine driver heterogeneity. For motor insurance during a policy year, we typically observe a large proportion of drivers with zero accidents, a lower proportion with exactly one accident, and a far lower proportion with two or more accidents. We here introduce a cost-sensitive multi-class adaptive boosting (AdaBoost) algorithm we call SAMME.C2 to handle such class imbalances. We calibrate the algorithm using empirical data collected from a telematics program in Canada and demonstrate an improved assessment of driving behavior using telematics compared with traditional risk variables. Using suitable performance metrics, we show that our algorithm outperforms other learning models designed to handle class imbalances.

This paper provides a method to assess the risk relief deriving from a foreign expansion by a life insurance company. We build a parsimonious continuous-time model for longevity risk that captures the dependence across different ages in domestic versus foreign populations. We calibrate the model to portray the case of a UK annuity portfolio expanding internationally toward Italian policyholders. The longevity risk diversification benefits of an international expansion are sizable, in particular when interest rates are low. The benefits are judged based on traditional measures, such as the Risk Margin or volatility reduction, and on a novel measure, the Diversification Index.

Even though the trend in mortality improvements has experienced several permanent changes in the past, the uncertainty regarding future mortality trends is often left unmodeled when pricing longevity-linked securities. In this paper, we present a stochastic modeling framework for the valuation of longevity-linked securities which explicitly considers the risk of random future changes in the long-term mortality trend. We construct a set of meaningful probability distortions which imply equivalent risk-adjusted pricing measures under which the basic model structure is preserved. Inspired by risk-based capital requirements for (re)insurers, we also establish a cost-of-capital pricing approach which then serves as the appropriate reference framework for finding a reasonable range for the market price of longevity risk. In a numerical application, we demonstrate that our model produces plausible risk loadings and show that a greater proportion of the risk loading is allocated to longer maturities when the risk of random future mortality trend changes is adequately modeled.

Catastrophic loss data are known to be heavy-tailed. Practitioners then need models that are able to capture both tail and modal parts of claim data. To this purpose, a new parametric family of loss distributions is proposed as a gamma mixture of the generalized log-Moyal distribution from Bhati and Ravi (2018), termed the generalized log-Moyal gamma (GLMGA) distribution. While the GLMGA distribution is a special case of the GB2 distribution, we show that this simpler model is effective in regression modeling of large and modal loss data. Regression modeling and applications to risk measurement are illustrated using a detailed analysis of a Chinese earthquake loss data set, comparing with the results of competing models from the literature. To this end, we discuss the probabilistic characteristics of the GLMGA and statistical estimation of the parameters through maximum likelihood. Further illustrations of the applicability of the new class of distributions are provided with the fire claim data set reported in Cummins et al. (1990) and a Norwegian fire losses data set discussed recently in Bhati and Ravi (2018).

Catastrophe insurance markets fail to provide sufficient protections against natural catastrophes, whereas they have the capacity to absorb the losses. In this paper, we assume the catastrophic risks are dependent and extremely heavy-tailed, and insurers have limited liability to cover losses up to a certain amount. We provide a comprehensive study to show that the diversification in the catastrophe insurance markets can be transited from suboptimal to preferred by increasing the number of insurers in the market. This highlights the importance of coordination among insurers and the government intervention in encouraging insurers to participate in the catastrophe insurance market to exploit risk sharing. Simulation studies are provided to illuminate the key findings of our results.

Target date funds (TDFs) are becoming increasingly popular investment choices among investors with long-term prospects. Examples include members of superannuation funds seeking to save for retirement at a given age. TDFs provide efficient risk exposures to a diversified range of asset classes that dynamically match the risk profile of the investment payoff as the investors age. This is often achieved by making increasingly conservative asset allocations over time as the retirement date approaches. Such dynamically evolving allocation strategies for TDFs are often referred to as glide paths. We propose a systematic approach to the design of optimal TDF glide paths implied by retirement dates and risk preferences and construct the corresponding dynamic asset allocation strategy that delivers the optimal payoffs at minimal costs. The TDF strategies we propose are dynamic portfolios consisting of units of the growth-optimal portfolio (GP) and the risk-free asset. Here, the GP is often approximated by a well-diversified index of multiple risky assets. We backtest the TDF strategies with the historical returns of the S&P500 total return index serving as the GP approximation.

This article proposes a complex economic scenario generator that nests versions of well-known actuarial frameworks. The generator estimation relies on the Bayesian paradigm and accounts for both model and parameter uncertainty via Markov chain Monte Carlo methods. So, to the question is less more?, we answer maybe, but it depends on your criteria. From an in-sample fit perspective, on the one hand, a complex economic scenario generator seems better. From the conservatism, forecasting and coverage perspectives, on the other hand, the situation is less clear: having more complex models for the short rate, term structure and stock index returns is clearly beneficial. However, that is not the case for inflation and the dividend yield.

The number of people who receive a stable income for life from a closed pooled annuity fund is studied. Income stability is defined as keeping the income within a specified tolerance of the initial income in a fixed proportion of future scenarios. The focus is on quantifying the effect of the number of members, which drives the level of idiosyncratic longevity risk in the fund, on the income stability. To do this, investment returns are held constant, and systematic longevity risk is omitted. An analytical expression that closely approximates the number of fund members who receive a stable income is derived and is seen to be independent of the mortality model. An application of the result is to calculate the length of time for which the pooled annuity fund can provide the desired level of income stability.

Variable annuity (VA) policies are typically issued on mutual funds invested in both fixed income and equity asset classes. However, due to the lack of specialized models to represent the dynamics of fixed income fund returns, the literature has primarily focused on studying long-term investment guarantees on single-asset equity funds. This article develops a mixed bond and equity fund model in which the fund return is linked to movements of the yield curve. Theoretical motivation for our proposed specification is provided through an analogy with a portfolio of rolling horizon bonds. Moreover, basis risk between the portfolio return and its risk drivers is naturally incorporated into our framework. Numerical results show that the fit of our model to Canadian VA data is adequate. Finally, the valuation of VAs is illustrated and it is found that the prevailing interest rate environment can have a substantial impact on guarantee costs.

In the context of the Solvency II directive, the operation of an internal risk model is a possible way for risk assessment and for the determination of the solvency capital requirement of an insurance company in the European Union. A Monte Carlo procedure is customary to generate a model output. To be compliant with the directive, validation of the internal risk model is conducted on the basis of the model output. For this purpose, we suggest a new test for checking whether there is a significant change in the modeled solvency capital requirement. Asymptotic properties of the test statistic are investigated and a bootstrap approximation is justified. A simulation study investigates the performance of the test in the finite sample case and confirms the theoretical results. The internal risk model and the application of the test is illustrated in a simplified example. The method has more general usage for inference of a broad class of law-invariant and coherent risk measures on the basis of a paired sample.

The primary objective of this scholarly work is to develop two estimation procedures – maximum likelihood estimator (MLE) and method of trimmed moments (MTM) – for the mean and variance of lognormal insurance payment severity data sets affected by different loss control mechanism, for example, truncation (due to deductibles), censoring (due to policy limits), and scaling (due to coinsurance proportions), in insurance and financial industries. Maximum likelihood estimating equations for both payment-per-payment and payment-per-loss data sets are derived which can be solved readily by any existing iterative numerical methods. The asymptotic distributions of those estimators are established via Fisher information matrices. Further, with a goal of balancing efficiency and robustness and to remove point masses at certain data points, we develop a dynamic MTM estimation procedures for lognormal claim severity models for the above-mentioned transformed data scenarios. The asymptotic distributional properties and the comparison with the corresponding MLEs of those MTM estimators are established along with extensive simulation studies. Purely for illustrative purpose, numerical examples for 1500 US indemnity losses are provided which illustrate the practical performance of the established results in this paper.

Modeling policyholders’ lapse behaviors is important to a life insurer, since lapses affect pricing, reserving, profitability, liquidity, risk management, and the solvency of the insurer. In this paper, we apply two machine learning methods to lapse modeling. Then, we evaluate the performance of these two methods along with two popular statistical methods by means of statistical accuracy and profitability measure. Moreover, we adopt an innovative point of view on the lapse prediction problem that comes from churn management. We transform the classification problem into a regression question and then perform optimization, which is new to lapse risk management. We apply the aforementioned four methods to a large real-world insurance dataset. The results show that Extreme Gradient Boosting (XGBoost) and support vector machine outperform logistic regression (LR) and classification and regression tree with respect to statistic accuracy, while LR performs as well as XGBoost in terms of retention gains. This highlights the importance of a proper validation metric when comparing different methods. The optimization after the transformation brings out significant and consistent increases in economic gains. Therefore, the insurer should conduct optimization on its economic objective to achieve optimal lapse management.

This paper investigates a high-dimensional vector-autoregressive (VAR) model in mortality modeling and forecasting. We propose an extension of the sparse VAR (SVAR) model fitted on the log-mortality improvements, which we name “spatially penalized smoothed VAR” (SSVAR). By adaptively penalizing the coefficients based on the distances between ages, SSVAR not only allows a flexible data-driven sparsity structure of the coefficient matrix but simultaneously ensures interpretable coefficients including cohort effects. Moreover, by incorporating the smoothness penalties, divergence in forecast mortality rates of neighboring ages is largely reduced, compared with the existing SVAR model. A novel estimation approach that uses the accelerated proximal gradient algorithm is proposed to solve SSVAR efficiently. Similarly, we propose estimating the precision matrix of the residuals using a spatially penalized graphical Lasso to further study the dependency structure of the residuals. Using the UK and France population data, we demonstrate that the SSVAR model consistently outperforms the famous Lee–Carter, Hyndman–Ullah, and two VAR-type models in forecasting accuracy. Finally, we discuss the extension of the SSVAR model to multi-population mortality forecasting with an illustrative example that demonstrates its superiority in forecasting over existing approaches.

In various applications of heavy-tail modelling, the assumed Pareto behaviour is tempered ultimately in the range of the largest data. In insurance applications, claim payments are influenced by claim management and claims may, for instance, be subject to a higher level of inspection at highest damage levels leading to weaker tails than apparent from modal claims. Generalizing earlier results of Meerschaert et al. (2012) and Raschke (2020), in this paper we consider tempering of a Pareto-type distribution with a general Weibull distribution in a peaks-over-threshold approach. This requires to modulate the tempering parameters as a function of the chosen threshold. Modelling such a tempering effect is important in order to avoid overestimation of risk measures such as the value-at-risk at high quantiles. We use a pseudo maximum likelihood approach to estimate the model parameters and consider the estimation of extreme quantiles. We derive basic asymptotic results for the estimators, give illustrations with simulation experiments and apply the developed techniques to fire and liability insurance data, providing insight into the relevance of the tempering component in heavy-tail modelling.

Target benefit (TB) plans that incorporate intergenerational risk sharing have been demonstrated to be welfare improving over the long term. However, there has been little discussion of the short-term benefits for members in a defined benefit (DB) plan that is transitioning to TB. In this paper, we adopt a two-step approach that is designed to ensure the long-term sustainability of the new plan, without unduly sacrificing the benefit security of current retirees. We propose a cohort-based transition plan for reducing intergenerational inequity. Our study is based on simulations using an economic scenario generator with some theoretical results under simplified settings.

Assessing conditional tail risk at very high or low levels is of great interest in numerous applications. Due to data sparsity in high tails, the widely used quantile regression method can suffer from high variability at the tails, especially for heavy-tailed distributions. As an alternative to quantile regression, expectile regression, which relies on the minimization of the asymmetric l2-norm and is more sensitive to the magnitudes of extreme losses than quantile regression, is considered. In this article, we develop a new estimation method for high conditional tail risk by first estimating the intermediate conditional expectiles in regression framework, and then estimating the underlying tail index via weighted combinations of the top order conditional expectiles. The resulting conditional tail index estimators are then used as the basis for extrapolating these intermediate conditional expectiles to high tails based on reasonable assumptions on tail behaviors. Finally, we use these high conditional tail expectiles to estimate alternative risk measures such as the Value at Risk (VaR) and Expected Shortfall (ES), both in high tails. The asymptotic properties of the proposed estimators are investigated. Simulation studies and real data analysis show that the proposed method outperforms alternative approaches.