A quantification of the financial implications of the design of a funded, collective defined contribution (CDC) pension scheme is presented and illustrated. It is done through an attribution analysis, which allows the importance of various elements of CDC scheme design to be determined. The model of a CDC scheme analysed is based lightly on the first CDC scheme set to be approved in the UK. In the CDC scheme analysed, contributions are fixed and the initial benefit accrued by each contribution is fixed. Once accrued, benefits are subsequently adjusted annually in response to changes in assumptions and returns. An attribution of the benefit payments shows that this design gives higher benefits to the first generations and lower benefits to the last generations, for a scheme which starts with no members. The contributions paid also affect the balance of benefits paid between generations. Too high a contribution is to the advantage of the first generations. Too low a contribution is in the interests of the later generations. The conclusion, within the simple model considered, is that a constant benefit accrual is an important design choice. Its financial consequences across all generations should be carefully analysed, if it is intended to be implemented. Additionally, contributions should be reviewed regularly in such a CDC scheme, to ensure that cross-subsidies are not borne excessively by particular generations.

A detailed analysis of management and performance fees for asset managers and investment funds is undertaken. While fund fees are considered as a cost of capital for investors, the structuring of such fee mechanisms in a fund can also influence a fund manager’s decisions and investment strategy, thereby also influencing the investment performance of the investors funds. The study undertaken will allow for an assessment of the effect of fee structures and the potential for asymmetric incentives to arise that may promote adverse risk-taking behaviours by the fund manager, to the detriment of the investor or retiree who places a portion of their retirement savings into such a managed fund with such fee structures. As such, understanding the mechanism of fee charging as well as pricing the fees correctly is vital. An exploration of the application of actuarial distortion pricing methods for complete and incomplete market valuation is performed on a variety of path-dependent option-like performance fee structures for various funds in the European and American markets. Furthermore, several scenario analysis and sensitivity studies are undertaken. The class of Net Asset Value models adopted are Lévy processes, and the pricing is performed via Monte Carlo techniques.

This paper considers the pricing of long-term options on assets such as housing, where either government intervention or the economic nature of the asset limits large falls in prices. The observed asset price is modelled by a geometric Brownian motion (“the notional price”) reflected at a lower barrier. The resulting observed price has standard dynamics but with localised intervention at the barrier, which allows arbitrage with interim losses; this is funded by the government’s unlimited powers of intervention, and its exploitation is subject to credit constraints. Despite the lack of an equivalent martingale measure for the observed price, options on this price can be expressed as compound options on the arbitrage-free notional price, to which standard risk-neutral arguments can be applied. Because option deltas tend to zero when the observed price approaches the barrier, hedging with the observed price gives the same results as hedging with the notional price and so exactly replicates option payoffs. Hedging schemes are not unique, with the cheapest scheme for any derivative being the one which best exploits the interventions at the barrier. The price of a put is clear: direct replication has a lower initial cost than synthetic replication, and the replication portfolio always has positive value. The price of a call is ambiguous: synthetic replication has a lower initial cost than direct replication, but the replication portfolio may give interim losses. So the preferred replication strategy (and hence price) of a call depends on what margin payments need to be made on these losses.

Determining accurate capital requirements is a central activity across the life insurance industry. This is computationally challenging and often involves the acceptance of proxy errors that directly impact capital requirements. Within simulation-based capital models, where proxies are being used, capital estimates are approximations that contain both statistical and proxy errors. Here, we show how basic error analysis combined with targeted exact computation can entirely eliminate proxy errors from the capital estimate. Consideration of the possible ordering of losses, combined with knowledge of their error bounds, identifies an important subset of scenarios. When these scenarios are calculated exactly, the resulting capital estimate can be made devoid of proxy errors. Advances in the handling of proxy errors improve the accuracy of capital requirements.

This paper provides a stochastic model, consistent with Solvency II and the Delegated Regulation, to quantify the capital requirement for demographic risk. In particular, we present a framework that models idiosyncratic and trend risks exploiting a risk theory approach in which results are obtained analytically. We apply the model to non-participating policies and quantify the Solvency Capital Requirement for the aforementioned risks in different time horizons.

We derive formulae for the moments of the time of ruin in both ordinary and modified Sparre Andersen risk models without specifying either the inter-claim time distribution or the individual claim amount distribution. We illustrate the application of our results in the special case of exponentially distributed claims, as well as for the following ordinary models: the classical risk model, phase-type(2) risk models, and the Erlang( $\mathscr{n}$ ) risk model. We also show how the key quantities for modified models can be found.

The financing of long-term care and the planning of care capacity are of increasing interest due to demographic changes and the ageing population in many countries. Since many care-intensive conditions begin to manifest at higher ages, a better understanding and assessment of the expected costs, required infrastructure, and number of qualified personnel are essential. To evaluate the overall burden of institutional care, we derive a model based on the duration of stay in dependence and the intensity of help provided to elderly individuals. This article aims to model both aspects using novel longitudinal data from nursing homes in the canton of Geneva in Switzerland. Our data contain comprehensive health and care information, including medical diagnoses, levels of dependence, and physical and psychological impairments on 21,758 individuals. We build an accelerated failure time model to study the influence of selected factors on the duration of care and a beta regression model to describe the intensity of care. We show that apart from age and gender, the duration of stay before death is mainly affected by the underlying diseases and the number of different diagnoses. Simultaneously, care intensity is driven by the individual level of dependence and specific limitations. Using both evaluations, we approximate the overall care severity for individual profiles. Our study sheds light on the relevant medical, physical, and psychological health indicators that need to be accounted for, not only by care providers but also by policy-makers and insurers.

We propose a new estimator for the ultimate prediction uncertainty within the famous Mack’s distribution-free chain-ladder model, which can be proved to be unbiased (conditionally given the first triangle column) under some additional technical assumptions. A peculiar behaviour of the unbiased estimator is given by its possible negativity. This is a drawback which might be worth trading off for the unbiasedness property, since there is empirical evidence that the likelihood of a negative realisation is extremely low. This offers an alternative to the well-known Mack and BBMW formulas since the latters can be proved to be biased. However, we also show that this novel estimator, as well as the Mack and BBMW formulas, can (with non-negligible probability) materially fail to estimate the true uncertainty.

In recent papers, Bonus-Malus Scales (BMS) estimated using data have been considered as an alternative to longitudinal data and hierarchical data approaches to model the dependence between different contracts for the same insured. Those papers, however, did not discuss in detail how to construct and understand BMS models, and many of the BMS’s basic properties were not discussed. The first objective of this paper is to correct this situation by explaining the logic behind BMS models and by describing those properties. More particularly, we will explain how BMS models are linked with simple count regression models that have covariates associated with the past claims experience. This study could help actuaries to understand how and why they should use BMS models for experience rating. The second objective of this paper is to create artificial past claims history for each insured. This is done by combining recent panel data theory with BMS models. We show that this addition significantly improves the prediction capacity of the BMS and provides a temporary solution for insurers who do not have enough historical data. We apply the BMS model to real data from a major Canadian insurance company. Results are analysed deeply to identify specific aspects of the BMS model.

This paper proposes a shift in the valuation and production of long-term annuities, away from the classical risk-neutral methodology towards a methodology using the real-world probability measure. The proposed production method is applied to three examples of annuity products, one having annual payments linked to a mortality index and the savings account and the others having annual payments linked to a mortality index and an equity index with a guarantee that is linked to the same mortality index and the savings account. Out-of-sample hedge simulations demonstrate the effectiveness of the proposed less-expensive production method. In contrast to classical risk-neutral production, which revolves around the savings account as reference unit, the long-term best-performing portfolio, the numéraire portfolio of the equity market, is employed as the fundamental reference unit in the production of the annuity. The numéraire portfolio is the strictly positive, tradable portfolio that when used as denominator or benchmark makes all benchmarked non-negative portfolios supermartingales. Under real-world valuation, the initial benchmarked value of a benchmarked contingent claim equals its real-world conditional expectation. The proposed real-world valuation and production can lead to significantly lower values of long-term annuities and their less-expensive production than suggested by the risk-neutral approach.

The Covid-19 pandemic caused an alarming mortality stress. The evidence shows that a significant proportion of people who die from Covid-19 are in a frail state. According to this consideration, we assume that the mortality shocks are related to a group of the individuals with some co-morbidities at Covid-19 diagnosis. In other words, the mortality shocks present a specific characterisation, which consists of a causal connection with pre-existing conditions, and the phenomenon could be described as a mortality acceleration. In this paper, an Accelerated Mortality Model is proposed in order to capture the different effects on mortality that depend on the evolution of the pandemic and the presence of co-morbidities at diagnosis. Furthermore, we assess the impact of Covid-19 mortality acceleration on a set of traditional life insurance contracts. We observe that, although mortality acceleration by Covid-19 affects more markedly the elderly and unhealthy sub-populations, it could be considered as a temporary shock with a limited impact on the life insurance market.

There has been a growing interest among pension plan sponsors in envisioning how the mortality experience of their active and deferred members may turn out to be if a pandemic similar to the COVID-19 occurs in the future. To address their needs, we propose in this paper a stochastic model for simulating future mortality scenarios with COVID-alike effects. The proposed model encompasses three parameter levels. The first level includes parameters that capture the long-term pattern of mortality, whereas the second level contains parameters that gauge the excess age-specific mortality due to COVID-19. Parameters in the first and second levels are estimated using penalised quasi-likelihood maximisation method which was proposed for generalised linear mixed models. Finally, the third level includes parameters that draw on expert opinions concerning, for example, how likely a COVID-alike pandemic will occur in the future. We illustrate our proposed model with data from the United States and a range of expert opinions.

The COVID-19 pandemic interrupts the relatively steady trend of improving longevity observed in many countries over the last decades. We claim that this needs to be addressed explicitly in many mortality modelling applications, for example, in the life insurance industry. To support this position, we provide a descriptive analysis of the mortality development of several countries up to and including the year 2020. Furthermore, we perform an empirical and theoretical investigation of the impact a mortality jump has on the parameters, forecasts and implied present values of the popular Lee–Carter mortality model. We find that COVID-19 has resulted in substantial mortality shocks in many countries. We show that such shocks have a large impact on point and interval forecasts of death rates and, consequently, on the valuation of mortality-related insurance products. We obtain similar findings under the Cairns–Blake–Dowd mortality model, which demonstrates that the effects caused by COVID-19 show up in a variety of models. Finally, we provide an overview of approaches to handle extreme mortality events such as the COVID-19 pandemic in mortality modelling.

Estimating tail risk measures for portfolios of complex variable annuities is an important enterprise risk management task which usually requires nested simulation. In the nested simulation, the outer simulation stage involves projecting scenarios of key risk factors under the real-world measure, while the inner simulations are used to value pay-offs under guarantees of varying complexity, under a risk-neutral measure. In this paper, we propose and analyse an efficient simulation approach that dynamically allocates the inner simulations to the specific outer scenarios that are most likely to generate larger losses. These scenarios are identified using a proxy calculation that is used only to rank the outer scenarios, not to estimate the tail risk measure directly. As the proxy ranking will not generally provide a perfect match to the true ranking of outer scenarios, we calculate a measure based on the concomitant of order statistics to test whether further tail scenarios are required to ensure, with given confidence, that the true tail scenarios are captured. This procedure, which we call the dynamic importance allocated nested simulation approach, automatically adjusts for the relationship between the proxy calculations and the true valuations and also signals when the proxy is not sufficiently accurate.

We examine the impact of asset allocation and contribution rates on the risk of defined benefit (DB) pension schemes, using both a run-off and a shorter 3-year time horizon. Using the 3-year horizon, which is typically preferred by regulators, a high bond allocation reduces the spread of the distribution of surplus. However, this result is reversed when examined on a run-off basis. Furthermore, under both the 3-year horizon and the run-off, the higher bond allocation reduces the median level of surplus. Pressure on the affordability of DB schemes has led to widespread implementation of the so-called de-risking strategies, such as moving away from predominantly equity investments to greater bond investments. If the incentives produced by shorter term risk assessments are contributing to this shift, they might be harming the long-term financial health of the schemes. Contribution rates have relatively lower impact on the risk.

Understanding core statistical properties and data features in mortality data are fundamental to the development of machine learning methods for demographic and actuarial applications of mortality projection. The study of statistical features in such data forms the basis for classification, regression and forecasting tasks. In particular, the understanding of key statistical structure in such data can aid in improving accuracy in undertaking mortality projection and forecasting when constructing life tables. The ability to accurately forecast mortality is a critical aspect for the study of demography, life insurance product design and pricing, pension planning and insurance-based decision risk management. Though many stylised facts of mortality data have been discussed in the literature, we provide evidence for a novel statistical feature that is pervasive in mortality data at a national level that is as yet unexplored. In this regard, we demonstrate in this work a strong evidence for the existence of long memory features in mortality data, and second that such long memory structures display multifractality as a statistical feature that can act as a discriminator of mortality dynamics by age, gender and country. To achieve this, we first outline the way in which we choose to represent the persistence of long memory from an estimator perspective. We make a natural link between a class of long memory features and an attribute of stochastic processes based on fractional Brownian motion. This allows us to use well established estimators for the Hurst exponent to then robustly and accurately study the long memory features of mortality data. We then introduce to mortality analysis the notion from data science known as multifractality. This allows us to study the long memory persistence features of mortality data on different timescales. We demonstrate its accuracy for sample sizes commensurate with national-level age term structure historical mortality records. A series of synthetic studies as well a comprehensive analysis of real mortality death count data are studied in order to demonstrate the pervasiveness of long memory structures in mortality data, both mono-fractal and multifractal functional features are verified to be present as stylised facts of national-level mortality data for most countries and most age groups by gender. We conclude by demonstrating how such features can be used in kernel clustering and mortality model forecasting to improve these actuarial applications.

Raw mortality data often exhibit irregular patterns due to randomness. Graduation refers to the act of smoothing crude mortality rates. In this paper, we propose a flexible and robust methodology for graduating mortality rates using adaptive P-splines. Since the observed data at high ages are often sparse and unreliable, we use an exponentially increasing penalty. We use mortality data of England and Wales and model male and female mortality rates jointly by means of penalties, achieving borrowing of information between the two sexes.

The classical credibility theory is a cornerstone of experience rating, especially in the field of property and casualty insurance. An obstacle to putting the credibility theory into practice is the conversion of available prior information into a precise choice of crucial hyperparameters. In most real-world applications, the information necessary to justify a precise choice is lacking, so we propose an imprecise credibility estimator that honestly acknowledges the imprecision in the hyperparameter specification. This results in an interval estimator that is doubly robust in the sense that it retains the credibility estimator’s freedom from model specification and fast asymptotic concentration, while simultaneously being insensitive to prior hyperparameter specification.

In most industrialised countries, one of the major societal challenges is the demographic change coming along with the ageing of the population. The increasing life expectancy observed over the last decades underlines the importance to find ways to appropriately cover the financial needs of the elderly. A particular issue arises in the area of health, where sufficient care must be provided to a growing number of dependent elderly in need of long-term care (LTC) services. In many markets, the offering of life insurance products incorporating care options and LTC insurance products is generally scarce. In our research, we therefore examine a life annuity product with an embedded care option potentially providing additional financial support to dependent persons. To evaluate the care option, we determine the minimum price that the annuity provider requires and the policyholder’s willingness to pay for the care option. For the latter, we employ individual utility functions taking account of the policyholder’s condition. We base our numerical study on recently developed transition probability data from Switzerland. Our findings give new and realistic insights into the nature and the utility of life annuity products proposing an embedded care option for tackling the financing of LTC needs.

Initial insurance losses are often reported with a textual description of the claim. The claims manager must determine the adequate case reserve for each known claim. In this paper, we present a framework for predicting the amount of loss given a textual description of the claim using a large number of words found in the descriptions. Prior work has focused on classifying insurance claims based on keywords selected by a human expert, whereas in this paper the focus is on loss amount prediction with automatic word selection. In order to transform words into numeric vectors, we use word cosine similarities and word embedding matrices. When we consider all unique words found in the training dataset and impose a generalised additive model to the resulting explanatory variables, the resulting design matrix is high dimensional. For this reason, we use a group lasso penalty to reduce the number of coefficients in the model. The scalable, analytical framework proposed provides for a parsimonious and interpretable model. Finally, we discuss the implications of the analysis, including how the framework may be used by an insurance company and how the interpretation of the covariates can lead to significant policy change. The code can be found in the TAGAM R package (github.com/scottmanski/TAGAM).