Abstract The Nelson–Siegel model is widely used in fixed income markets to produce yield curve dynamics. The multiple time-dependent parameter model conveniently addresses the level, slope, and curvature dynamics of the yield curves. In this study, we present a novel state-space functional regression model that incorporates a dynamic Nelson–Siegel (DNS) model and functional regression formulations applied to a multi-economy setting. This framework offers distinct advantages in explaining the relative spreads in yields between a reference economy and a response economy. To address the inherent challenges of model calibration, a kernel principal component analysis is employed to transform the representation of functional regression into a finite-dimensional, tractable estimation problem. A comprehensive empirical analysis is conducted to assess the efficacy of the functional regression approach, including an in-sample performance comparison with the DNS model. We conducted the stress testing analysis of the yield curves’ term structure within a dual economy framework. The bond ladder portfolio was examined through a case study focused on spread modeling using historical data for US Treasury and UK bonds.

Abstract Credibility theory provides a fundamental framework in actuarial science for estimating policyholder premiums by blending individual claims experience with overall portfolio data. Bühlmann and Bühlmann–Straub credibility models are widely used because, in the Bayesian hierarchical setting, they are the best linear Bayes estimators, minimizing the Bayes risk (expected squared error loss) within the class of linear estimators given the experience data for a particular risk class. To improve estimation accuracy, quadratic credibility models incorporate higher-order terms, capturing more information about the underlying risk structure. This study develops a robust quadratic credibility (RQC) framework that integrates second-order polynomial adjustments of robustly transformed ground-up loss data, such as winsorized moments, to improve stability in the presence of extreme claims or heavy-tailed distributions. Extending semi-linear credibility, RQC maintains interpretability while enhancing statistical efficiency. We establish its asymptotic properties, derive closed-form expressions for the RQC premium, and demonstrate its superior performance in reducing mean square error (MSE). We additionally derive semi-linear credibility structural parameters using winsorized data, further strengthening the robustness of credibility estimation. Analytical comparisons and empirical applications highlight RQC’s ability to capture claim heterogeneity, offering a more reliable and equitable approach to premium estimation. This research advances credibility theory by introducing a refined methodology that balances efficiency, robustness, and practical applicability across diverse insurance settings.

Abstract The matrixdist R package provides a comprehensive suite of tools for the statistical analysis of matrix distributions, including phase-type, inhomogeneous phase-type, discrete phase-type, and related multivariate distributions. This paper introduces the package and its key features, including the estimation of these distributions and their extensions through expectation-maximization algorithms, as well as the implementation of regression through the proportional intensities and mixture-of-experts models. Additionally, the paper provides an overview of the theoretical background, discusses the algorithms and methods implemented in the package, and offers practical examples to illustrate the application of matrixdist in real-world actuarial problems. The matrixdist R package aims to provide researchers and practitioners a wide set of tools for analyzing and modeling complex data using matrix distributions.

Abstract Healthcare costs tend to increase with age. In particular, in the case of illness, the last year before death can be an exceptionally costly period as the need for healthcare increases. Using a novel private insurance dataset containing over one million records of claims submitted by individuals to their health insurance providers during the last year of life, our research seeks to shed light on the costs before death in Switzerland. Our work documents how spending patterns change with proximity to dying. We use machine learning algorithms to identify and quantify the key effects that drive a person’s spending during this critical period. Our findings provide a more profound understanding of the costs associated with hospitalization before death, the role of age, and the variation in costs based on the services, including care services, which individuals require.

Abstract Cause-of-death mortality forecasting, a key topic in public health and actuarial science, is a challenging task due to the difficulty of modeling that accounts for dependencies among causes of death. While several cause-of-death mortality models have been proposed to address this difficulty, little attention has been paid to improving their predictive performance. Recently, purely data-driven approaches using tensor decomposition methods have been introduced to cause-of-death mortality modeling, demonstrating strong out-of-sample predictive performance compared to existing models. However, these methods have difficulties in the interpretability of multi-rank tensor components to achieve strong predictive performance. In response, we propose a novel tensor-based cause-of-death mortality model by replacing the tensor decomposition with a convolutional autoencoder with a one-dimensional latent layer that provides a Lee-Carter-like time-series factor; the model also provides the age sensitivity of cause-specific log mortality to the time-series factor. Due to the representational capability of the neural network, our model achieves better predictive performance compared to the existing tensor decomposition-based models, despite the simplified latent layer for model interpretability.

Abstract Understanding and forecasting mortality by cause is an essential branch of actuarial science, with wide-ranging implications for decision-makers in public policy and industry. To accurately capture trends in cause-specific mortality, it is critical to consider dependencies between causes of death and produce forecasts by age and cause coherent with aggregate mortality forecasts. One way to achieve these aims is to model cause-specific deaths using compositional data analysis (CODA), treating the density of deaths by age and cause as a set of dependent, nonnegative values that sum to one. A major drawback of standard CODA methods is the challenge of zero values, which frequently occur in cause-of-death mortality modeling. Thus, we propose using a compositional power transformation, the $\alpha$ -transformation, to model cause-specific life-table death counts. The $\alpha$ -transformation offers a statistically rigorous approach to handling zero value subgroups in CODA compared to ad hoc techniques: adding an arbitrarily small amount. We illustrate the $\alpha$ -transformation in England and Wales and US death counts by cause from the Human Cause-of-Death database, for cardiovascular-related causes of death. The results demonstrate the $\alpha$ -transformation improves forecast accuracy of cause-specific life-table death counts compared with log-ratio-based CODA transformations. The forecasts suggest declines in the proportions of deaths from major cardiovascular causes (myocardial infarction and other ischemic heart diseases).

Abstract In this paper, we provide a theoretical framework justifying the existence of a correlation risk premium in a market with two traded assets. We prove that risk-neutral dependence can differ substantially from real-world dependence by characterizing the set of risk-neutral martingale measures. This implies that implied correlation can be significantly different with the realized correlation. Depending on the choice of the market regarding the pricing measure, implied correlation can be high or low. We label the difference between risk-neutral and real-world correlation the “correlation gap” and make the connection with correlation risk premium. We show how dispersion trading can be used to exploit this correlation gap and demonstrate how there can exist a negative correlation risk premium in the financial market.

Abstract Generalized Poisson (GP) distribution was introduced in Consul & Jain ((1973). Technometrics, 15(4), 791–799.). Since then it has found various applications in actuarial science and other areas. In this paper, we focus on the distributional properties of GP and its related distributions. In particular, we study the distributional properties of distributions in the $\mathcal{H}$ family, which includes GP and generalized negative binomial distributions as special cases. We demonstrate that the moment and size-biased transformations of distributions within the $\mathcal{H}$ family remain in the same family, which significantly extends the results presented in Ambagaspitiya & Balakrishnan ((1994). ASTINBulletin: the Journal of the IAA, 24(2), 255–263.) and Ambagaspitiya ((1995). Insurance Mathematics and Economics, 2(16), 107–127.). Such findings enable us to provide recursive formulas for evaluating risk measures, such as Value-at-Risk and conditional tail expectation of the compound GP distributions. In addition, we show that the risk measures can be calculated by making use of transform methods, such as fast Fourier transform. In fact, the transformation method showed a remarkable time advantage over the recursive method. We numerically compare the risk measures of the compound sums when the primary distributions are Poisson and GP. The results illustrate the model risk for the loss frequency distribution.

The publisher would like to apologise for an error introduced during the production process of the above article (Macdonald, 2024). The page references to four sources from the Scandinavian Actuarial Journal were incorrect. This erratum contains the correct references:

Abstract This paper introduces gemact, a Python package for actuarial modeling based on the collective risk model. The library supports applications to risk costing and risk transfer, loss aggregation, and loss reserving. We add new probability distributions to those available in scipy, including the (a, b, 0) and (a, b, 1) discrete distributions, copulas of the Archimedean family, the Gaussian, the Student t and the Fundamental copulas. We provide an implementation of the AEP algorithm for calculating the cumulative distribution function of the sum of dependent, nonnegative random variables, given their dependency structure specified with a copula. The theoretical framework is introduced at the beginning of each section to give the reader with a sufficient understanding of the underlying actuarial models.