We propose a new inference strategy for general population mortality tables based on annual population and death estimates, completed by monthly birth counts. We rely on a deterministic population dynamics model and establish formulas that link the death rates to be estimated with the observables at hand. The inference algorithm takes the form of a recursive and implicit scheme for computing death rate estimates. This paper demonstrates both theoretically and numerically the efficiency of using additional monthly birth counts for appropriately computing annual mortality tables. As a main result, the improved mortality estimators show better features, including the fact that previous anomalies in the form of isolated cohort effects disappear, which confirms from a mathematical perspective the previous contributions by Richards, Cairns et al., and Boumezoued.

Wavelet theory is known to be a powerful tool for compressing and processing time series or images. It consists in projecting a signal on an orthonormal basis of functions that are chosen in order to provide a sparse representation of the data. The first part of this article focuses on smoothing mortality curves by wavelets shrinkage. A chi-square test and a penalized likelihood approach are applied to determine the optimal degree of smoothing. The second part of this article is devoted to mortality forecasting. Wavelet coefficients exhibit clear trends for the Belgian population from 1965 to 2015, they are easy to forecast resulting in predicted future mortality rates. The wavelet-based approach is then compared with some popular actuarial models of Lee–Carter type estimated fitted to Belgian, UK, and US populations. The wavelet model outperforms all of them.

In insurance analytics, textual descriptions of claims are often discarded, because traditional empirical analyses require numeric descriptor variables. This paper demonstrates how textual data can be easily used in insurance analytics. Using the concept of word similarities, we illustrate how to extract variables from text and incorporate them into claims analyses using standard generalized linear model or generalized additive regression model. This procedure is applied to the Wisconsin Local Government Property Insurance Fund (LGPIF) data, in order to demonstrate how insurance claims management and risk mitigation procedures can be improved. We illustrate two applications. First, we show how the claims classification problem can be solved using textual information. Second, we analyze the relationship between risk metrics and the probability of large losses. We obtain good results for both applications, where short textual descriptions of insurance claims are used for the extraction of features.

When modelling subnational mortality rates, we should consider three features: (1) how to incorporate any possible correlation among subpopulations to potentially improve forecast accuracy through multi-population joint modelling; (2) how to reconcile subnational mortality forecasts so that they aggregate adequately across various levels of a group structure; (3) among the forecast reconciliation methods, how to combine their forecasts to achieve improved forecast accuracy. To address these issues, we introduce an extension of grouped univariate functional time-series method. We first consider a multivariate functional time-series method to jointly forecast multiple related series. We then evaluate the impact and benefit of using forecast combinations among the forecast reconciliation methods. Using the Japanese regional age-specific mortality rates, we investigate 1–15-step-ahead point and interval forecast accuracies of our proposed extension and make recommendations.

We present an actuarial claims reserving technique that takes into account both claim counts and claim amounts. Separate (overdispersed) Poisson models for the claim counts and the claim amounts are combined by a joint embedding into a neural network architecture. As starting point of the neural network calibration, we use exactly these two separate (overdispersed) Poisson models. Such a nested model can be interpreted as a boosting machine. It allows us for joint modeling and mutual learning of claim counts and claim amounts beyond the two individual (overdispersed) Poisson models.

Financial products are priced using risk-neutral expectations justified by hedging portfolios that (as accurate as possible) match the product’s payoff. In insurance, premium calculations are based on a real-world best-estimate value plus a risk premium. The insurance risk premium is typically reduced by pooling of (in the best case) independent contracts. As hybrid life insurance contracts depend on both financial and insurance risks, their valuation requires a hybrid valuation principle that combines the two concepts of financial and actuarial valuation. The aim of this paper is to present a novel three-step projection algorithm to valuate hybrid contracts by decomposing their payoff in three parts: a financial, hedgeable part, a diversifiable actuarial part, and a residual part that is neither hedgeable nor diversifiable. The first two parts of the resulting premium are directly linked to their corresponding hedging and diversification strategies, respectively. The method allows for a separate treatment of unsystematic, diversifiable mortality risk and systematic, aggregate mortality risk related to, for example, epidemics or population-wide improvements in life expectancy. We illustrate our method in the case of CAT bonds and a pure endowment insurance contract with profit and compare the three-step method to alternative valuation operators suggested in the literature.

Variable annuities are products offered by pension funds and life offices that provide periodic future payments to the investor and often have ancillary benefits that guarantee survival benefits or sums insured on death. This paper extends the benchmark approach to value and hedge long-dated variable annuities using a combination of cash, bonds and equities under a variety of market models, allowing for dependence between financial and insurance markets. Under a simplified case of independence, the results show that when the discounted index is modelled as a time-transformed squared Bessel process, less-expensive valuation and reserving is achieved regardless of the short rate model or the mortality model.

Marine is the oldest type of insurance coverage. Nevertheless, unlike cargo and hull covers, marine liability is a rather young line of business with claims that can have very heavy and long tails. For reinsurers, the accumulation of losses from an event insured by various Protection and Indemnity clubs is an additional source for very large claims in the portfolio. In this paper, we first describe some recent developments of the marine liability market and then statistically analyze a data set of large losses for this line of business in a detailed manner both in terms of frequency and severity, including censoring techniques and tests for stationarity over time. We further formalize and examine an optimization problem that occurs for reinsurers participating in XL on XL coverages in this line of business and give illustrations of its solution.

Bonus-malus systems (BMSs) are widely used in actuarial sciences. These systems are applied by insurance companies to distinguish the policyholders by their risks. The most known application of BMS is in automobile third-party liability insurance. In BMS, there are several classes, and the premium of a policyholder depends on the class he/she is assigned to. The classification of policyholders over the periods of the insurance depends on the transition rules. In general, optimization of these systems involves the calculation of an appropriate premium scale considering the number of classes and transition rules as external parameters. Usually, the stationary distribution is used in the optimization process. In this article, we present a mixed integer linear programming (MILP) formulation for determining the premium scale and the transition rules. We present two versions of the model, one with the calculation of stationary probabilities and another with the consideration of multiple periods of the insurance. Furthermore, numerical examples will also be given to demonstrate that the MILP technique is suitable for handling existing BMSs.

The generalized linear model (GLM) is a statistical model which has been widely used in actuarial practices, especially for insurance ratemaking. Due to the inherent longitudinality of property and casualty insurance claim datasets, there have been some trials of incorporating unobserved heterogeneity of each policyholder from the repeated observations. To achieve this goal, random effects models have been proposed, but theoretical discussions of the methods to test the presence of random effects in GLM framework are still scarce. In this article, the concept of Bregman divergence is explored, which has some good properties for statistical modeling and can be connected to diverse model selection diagnostics as in Goh and Dey [(2014) Journal of Multivariate Analysis, 124, 371–383]. We can apply model diagnostics derived from the Bregman divergence for testing robustness of a chosen prior by the modeler to possible misspecification of prior distribution both on the naive model, which assumes that random effects follow a point mass distribution as its prior distribution, and the proposed model, which assumes a continuous prior density of random effects. This approach provides insurance companies a concrete framework for testing the presence of nonconstant random effects in both claim frequency and severity and furthermore appropriate hierarchical model which can explain both observed and unobserved heterogeneity of the policyholders for insurance ratemaking. Both models are calibrated using a claim dataset from the Wisconsin Local Government Property Insurance Fund which includes both observed claim counts and amounts from a portfolio of policyholders.

We determine the optimal asset allocation to bonds and stocks using an annually recalculated virtual annuity (ARVA) spending rule for DC pension plan decumulation. Our objective function minimizes downside withdrawal variability for a given fixed value of total expected withdrawals. The optimal asset allocation is found using optimal stochastic control methods. We formulate the strategy as a solution to a Hamilton–Jacobi–Bellman (HJB) Partial Integro Differential Equation (PIDE). We impose realistic constraints on the controls (no-shorting, no-leverage, discrete rebalancing) and solve the HJB PIDEs numerically. Compared to a fixed-weight strategy which has the same expected total withdrawals, the optimal strategy has a much smaller average allocation to stocks and tends to de-risk rapidly over time. This conclusion holds in the case of a parametric model based on historical data and also in a bootstrapped market based on the historical data.

Tontines, retirement products constructed in such a way that the longevity risk is shared in a pool of policyholders, have recently gained vast attention from researchers and practitioners. Typically, these products are cheaper than annuities, but do not provide stable payments to policyholders. This raises the question whether, from the policyholders' viewpoint, the advantages of annuities and tontines can be combined to form a retirement plan which is cheaper than an annuity, but provides a less volatile retirement income than a tontine. In this article, we analyze and compare three approaches of combining annuities and tontines in an expected utility framework: the previously introduced “tonuity”, a product very similar to the tonuity which we call “antine” and a portfolio consisting of an annuity and a tontine. We show that the payoffs of a tonuity and an antine can be replicated by a portfolio consisting of an annuity and a tontine. Consequently, policyholders achieve higher expected utility levels when choosing the portfolio over the novel retirement products tonuity and antine. Further, we derive conditions on the premium loadings of annuities and tontines indicating when the optimal portfolio is investing a positive amount in both annuity and tontine, and when the optimal portfolio turns out to be a pure annuity or a pure tontine.

Several diagonal-based tail dependence indices have been suggested in the literature to quantify tail dependence. They have well-developed statistical inference theories but tend to underestimate tail dependence. For those problems when assessing the maximal strength of dependence is important (e.g., co-movements of financial instruments), the maximal tail dependence index was introduced, but it has so far lacked empirical estimators and statistical inference results, thus hindering its practical use. In the present paper, we suggest an empirical estimator for the index, explore its statistical properties, and illustrate its performance on simulated data.

In this paper, we consider the risk–return trade-off for variable annuities in a Black–Scholes setting. Our analysis is based on a novel explicit allocation of initial wealth over the payments at various horizons. We investigate the relationship between the optimal consumption problem and the design of variable annuities by deriving the optimal so-called assumed interest rate for an investor with constant relative risk aversion preferences. We investigate the utility loss due to deviations from this. Finally, we show analytically how habit-formation-type smoothing of financial market shocks over the remaining lifetime leads to smaller year-to-year volatility in pension payouts, but to increases in the longer-term volatility.

This paper studies deep learning approaches to find optimal reinsurance and dividend strategies for insurance companies. Due to the randomness of the financial ruin time to terminate the control processes, a Markov chain approximation-based iterative deep learning algorithm is developed to study this type of infinite-horizon optimal control problems. The optimal controls are approximated as deep neural networks in both cases of regular and singular types of dividend strategies. The framework of Markov chain approximation plays a key role in building the iterative equations and initialization of the algorithm. We implement our method to classic dividend and reinsurance problems and compare the learning results with existing analytical solutions. The feasibility of our method for complicated problems has been demonstrated by applying to an optimal dividend, reinsurance and investment problem under a high-dimensional diffusive model with jumps and regime switching.

The class of distortion riskmetrics is defined through signed Choquet integrals, and it includes many classic risk measures, deviation measures, and other functionals in the literature of finance and actuarial science. We obtain characterization, finiteness, convexity, and continuity results on general model spaces, extending various results in the existing literature on distortion risk measures and signed Choquet integrals. This paper offers a comprehensive toolkit of theoretical results on distortion riskmetrics which are ready for use in applications.

In this paper, we first derive closed-form formulas for mortality-interest durations and convexities of the prices of life insurance and annuity products with respect to an instantaneously proportional change and an instantaneously parallel movement, respectively, in μ* (the force of mortality-interest), the addition of μ (the force of mortality) and δ (the force of interest). We then build several mortality-interest duration and convexity matching strategies to determine the weights of whole life insurance and deferred whole life annuity products in a portfolio and evaluate the value at risk and the hedge effectiveness of the weighted portfolio surplus at time zero. Numerical illustrations show that using the mortality-interest duration and convexity matching strategies with respect to an instantaneously proportional change in μ* can more effectively hedge the longevity risk and interest rate risk embedded in the deferred whole life annuity products than using the mortality-only duration and convexity matching strategies with respect to an instantaneously proportional shift or an instantaneously constant movement in μ only.

We study the relation between one-year premium risk and ultimate premium risk. In practice, the one-year risk is sometimes related to the ultimate risk by using a so-called emergence pattern formula which postulates a linear relation between both risks. We define the true emergence pattern of the ultimate loss for the one-year premium risk based on a conditional distribution of the ultimate loss derived from a multivariate distribution of the claims development process. We investigate three models commonly used in claims reserving and prove that the true emergence pattern formulas are different from the linear emergence pattern formula used in practice. We show that the one-year risk, when measured by VaR, can be under and overestimated if the linear emergence pattern formula is applied. We present two modifications of the linear emergence pattern formula. These modifications allow us to go beyond the claims development models investigated in the first part and work with an arbitrary distribution of the ultimate loss.

A generalised property exposure rating framework is presented here to address two issues arising in the standard approach to exposure rating, especially in the context of direct insurance and facultative reinsurance (D&F) property pricing:

1. (a) What to do when the main assumption of exposure rating, scalability – that is, that the probability of a given damage ratio does not depend on the maximum possible loss (MPL) but only on the type of property – breaks down.

2. (b) How to take account of the uncertainty around the MPL, that is, the fact that the MPL (unlike the insured value, IV) is an informed estimate rather than a contractual feature.

The first difficulty is addressed by making the exposure rating framework more flexible by introducing exposure curves that are a mixture of scale-dependent and scale-independent losses, and where the weight given to the two components is a function of the MPL. This allows to model the impact on the expected losses of changes in the underlying deductibles, a classic problem in D&F pricing normally solved with the help of deductible impact tables. The second difficulty is addressed by working out the mathematical implications of having a finite probability of exceeding the MPL and extending the exposure curve up to the maximum of MPL and IV. A practical application of this framework – in which the scale-independent and scale-dependent losses are identified with attritional and large losses respectively – is described. A discussion of how this generalised framework can be calibrated based on actual data is included, and implementation code for the framework is made available.

We determine the optimal robust strategy of an individual who seeks to maximize the (penalized) probability of reaching a bequest goal when she is uncertain about the drift of the risky asset and her hazard rate of mortality. We assume the individual can invest in a Black–Scholes market. We solve two optimization problems with ambiguity. The first is to maximize the penalized probability of reaching a bequest goal without life insurance in the market. In the second problem, in addition to investing in the financial market, the individual is allowed to purchase term life insurance to help her reach her bequest goal. As the individual becomes more ambiguity averse concerning the drift of the risky asset, she becomes more conservative with her investment strategy. Also, as she becomes more ambiguity averse about her hazard rate of mortality, numerical work indicates she is more likely to buy life insurance when the ambiguity towards the return of the risky asset is not too large.