High-cardinality categorical features are pervasive in actuarial data (e.g., occupation in commercial property insurance). Standard categorical encoding methods like one-hot encoding are inadequate in these settings.

In this work, we present a novel Generalised Linear Mixed Model Neural Network (“GLMMNet”) approach to the modelling of high-cardinality categorical features. The GLMMNet integrates a generalised linear mixed model in a deep learning framework, offering the predictive power of neural networks and the transparency of random effects estimates, the latter of which cannot be obtained from the entity embedding models. Further, its flexibility to deal with any distribution in the exponential dispersion (ED) family makes it widely applicable to many actuarial contexts and beyond. In order to facilitate the application of GLMMNet to large datasets, we use variational inference to estimate its parameters—both traditional mean field and versions utilising textual information underlying the high-cardinality categorical features.

We illustrate and compare the GLMMNet against existing approaches in a range of simulation experiments as well as in a real-life insurance case study. A notable feature for both our simulation experiment and the real-life case study is a comparatively low signal-to-noise ratio, which is a feature common in actuarial applications. We find that the GLMMNet often outperforms or at least performs comparably with an entity-embedded neural network in these settings, while providing the additional benefit of transparency, which is particularly valuable in practical applications.

Importantly, while our model was motivated by actuarial applications, it can have wider applicability. The GLMMNet would suit any applications that involve high-cardinality categorical variables and where the response cannot be sufficiently modelled by a Gaussian distribution, especially where the inherent noisiness of the data is relatively high.

We present novel cross-sectional and longitudinal claim count models for vehicle insurance built upon the combinedd actuarial neural network (CANN) framework proposed by Wüthrich and Merz. The CANN approach combines a classical actuarial model, such as a generalized linear model, with a neural network. This blending of models results in a two-component model comprising a classical regression model and a neural network part. The CANN model leverages the strengths of both components, providing a solid foundation and interpretability from the classical model while harnessing the flexibility and capacity to capture intricate relationships and interactions offered by the neural network. In our proposed models, we use well-known log-linear claim count regression models for the classical regression part and a multilayer perceptron (MLP) for the neural network part. The MLP part is used to process telematics car driving data given as a vector characterizing the driving behavior of each insured driver. In addition to the Poisson and negative binomial distributions for cross-sectional data, we propose a procedure for training our CANN model with a multivariate negative binomial specification. By doing so, we introduce a longitudinal model that accounts for the dependence between contracts from the same insured. Our results reveal that the CANN models exhibit superior performance compared to log-linear models that rely on manually engineered telematics features.

While driver telematics has gained attention for risk classification in auto insurance, scarcity of observations with telematics features has been problematic, which could be owing to either privacy concerns or favorable selection compared to the data points with traditional features.

To handle this issue, we apply a data integration technique based on calibration weights for usage-based insurance with multiple sources of data. It is shown that the proposed framework can efficiently integrate traditional data and telematics data and can also deal with possible favorable selection issues related to telematics data availability. Our findings are supported by a simulation study and empirical analysis in a synthetic telematics dataset.

Unstructured data are a promising new source of information that insurance companies may use to understand their risk portfolio better and improve the customer experience. However, these novel data sources are difficult to incorporate into existing ratemaking frameworks due to the size and format of the unstructured data. This paper proposes a framework to use street view imagery within a generalized linear model. To do so, we use representation learning to extract an embedding vector containing useful information from the image. This embedding is dense and low dimensional, making it appropriate to use within existing ratemaking models. We find that there is useful information included in street view imagery to predict the frequency of claims for certain types of perils. This model can be used as in a ratemaking framework but also opens the door to future empirical research on attempting to extract which characteristics within the image leads to increased or decreased predicted claim frequencies. Throughout, we discuss the practical difficulties (technical and social) of using this type of data for insurance pricing.

The distribution-free chain ladder of Mack justified the use of the chain ladder predictor and enabled Mack to derive an estimator of conditional mean squared error of prediction for the chain ladder predictor. Classical insurance loss models, that is of compound Poisson type, are not consistent with Mack’s distribution-free chain ladder. However, for a sequence of compound Poisson loss models indexed by exposure (e.g., number of contracts), we show that the chain ladder predictor and Mack’s estimator of conditional mean squared error of prediction can be derived by considering large exposure asymptotics. Hence, quantifying chain ladder prediction uncertainty can be done with Mack’s estimator without relying on the validity of the model assumptions of the distribution-free chain ladder.

We develop a flexible Gaussian process (GP) framework for learning the covariance structure of Age- and Year-specific mortality surfaces. Utilizing the additive and multiplicative structure of GP kernels, we design a genetic programming algorithm to search for the most expressive kernel for a given population. Our compositional search builds off the Age–Period–Cohort (APC) paradigm to construct a covariance prior best matching the spatio-temporal dynamics of a mortality dataset. We apply the resulting genetic algorithm (GA) on synthetic case studies to validate the ability of the GA to recover APC structure and on real-life national-level datasets from the Human Mortality Database. Our machine learning-based analysis provides novel insight into the presence/absence of Cohort effects in different populations and into the relative smoothness of mortality surfaces along the Age and Year dimensions. Our modeling work is done with the PyTorch libraries in Python and provides an in-depth investigation of employing GA to aid in compositional kernel search for GP surrogates.

With recent epidemics such as COVID-19, H1N1 and SARS causing devastating financial loss to the economy, it is important that insurance companies plan for financial costs of epidemics. This article proposes a new methodology for epidemic and insurance modelling by combining the existing deterministic compartmental models and the Markov multiple state models to facilitate actuarial computations to design new health insurance plans that cover epidemics. Our method is inspired by the seminal paper (Feng and Garrido (2011) North American Actuarial Journal, 15, 112–136.) of Feng and Garrido and complements the work of Hillairet and Lopez et al. in Hillairet and Lopez ((2021) Scandinavian Actuarial Journal, 2021(8), 671–694.) and Hillairet et al. ((2022) Insurance: Mathematics and Economics, 107, 88–101.) In this work, we use the deterministic SIR model and the Eyam epidemic data set to provide numerical illustrations for our method.

We propose a flexible lattice model to evaluate the fair value of insurance contracts embedding both financial and actuarial risk factors. Flexibility relies on the ability of the model to manage different specifications of the correlated processes governing interest rate, mortality, and fund dynamics, thus allowing the insurer to make the most appropriate choices. The model is also able to handle additional guarantees like a surrender opportunity for which explicit formulae are not available being it similar to an American derivative. The model discretizes mortality and interest rate dynamics through two different binomial lattices and then combines them into a bivariate tree characterized by the presence of four branches for each node. The probability of each branch is defined to replicate the correlation affecting the two processes. The bivariate model is useful to compute the value of survival zero coupon bond. When adding another source of risk, such as the fund dynamics for evaluating fund-linked insurance products, we model it through a bivariate tree that captures the influence of the interest rate on its drift term. Then, the mortality risk is embedded by defining a trivariate tree presenting eight branches emanating from each node with probabilities defined in order to capture the correlations of the processes. Extensive numerical experiments assess the model accuracy by considering some stylized policies, but the model application is not limited to them being it able to manage different contract specifications.

Motivated by insurance applications, we propose a new approach for the validation of real-world economic scenarios. This approach is based on the statistical test developed by Chevyrev and Oberhauser ((2022) Journal of Machine Learning Research, 23(176), 1–42.) and relies on the notions of signature and maximum mean distance. This test allows to check whether two samples of stochastic processes paths come from the same distribution. Our contribution is to apply this test to a variety of stochastic processes exhibiting different pathwise properties (Hölder regularity, autocorrelation, and regime switches) and which are relevant for the modelling of stock prices and stock volatility as well as of inflation in view of actuarial applications.

In this paper, we explore how to design the optimal insurance contracts when the insured faces insurable, counterparty, and additive background risk simultaneously. The target is to minimize the mean-variance of the insured’s loss. By utilizing the calculus of variations, an implicit characterization of the optimal ceded loss function is given. An explicit structure of the optimal ceded loss function is also provided by making full use of its implicit characterization. We further derive a much simpler solution when these three kinds of risk have some special dependence structures. Finally, we give a numerical example to illustrate our results.