Aggregate implements an efficient fast Fourier transform (FFT)-based algorithm to approximate compound probability distributions. Leveraging FFT-based methods offers advantages over recursion and simulation-based approaches, providing speed and accuracy to otherwise time-consuming calculations. Combining user-friendly features and an expressive domain-specific language called DecL, Aggregate enables practitioners and nonprogrammers to work with complex distributions effortlessly. The software verifies the accuracy of its FFT-based numerical approximations by comparing their first three moments to those calculated analytically from the specified frequency and severity. This moment-based validation, combined with carefully chosen default parameters, allows users without in-depth knowledge of the underlying algorithm to be confident in the results. Aggregate supports a wide range of frequency and severity distributions, policy limits and deductibles, and reinsurance structures and has applications in pricing, reserving, risk management, teaching, and research. It is written in Python.

This paper studies a hierarchical risk model where an accident can cause a combination of different types of claims, whose sizes could be dependent. In addition, the frequencies of accidents that cause the different combinations of claims are dependent. We first derive formulas for computing risk measures, such as the Tail Conditional Expectation and Tail Variance of the aggregate losses for a portfolio of businesses. Then, we present formulas for performing the associated capital allocation to different types of claims in the portfolio. The main tool we used is the moment (or size-biased) transform of the multivariate distributions.

This paper presents the AffineMortality R package which performs parameter estimation, goodness-of-fit analysis, simulation, and projection of future mortality rates for a set of affine mortality models for use in pricing and reserving. The computational routines build on the univariate Kalman Filtering approach of Koopman and Durbin ((2000). Journal of Time Series Analysis, 21(3), 281–296.) along other numerical methods to enhance the robustness of the results. This paper provides a discussion of how the package works in order to effectively estimate and project survival curves, and describes the available functions. Illustration of the package for mortality analysis of the US male data set is provided.

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.

This study investigates the benefits and drawbacks of pension plan consolidation by quantifying the impact of mergers of heterogeneous plans on different stakeholders in a unique Canadian implementation of defined benefit plans. Using a comprehensive framework that combines a realistic economic scenario generator, a stochastic mortality model that captures differences among subpopulations, a cost model with economies of scale, and a dynamic asset allocation methodology, we evaluate the combined effect of asset- and liability-side changes on three groups of measures: plan-related risk measures assessing profits from an economic capital perspective, consumption-based metrics to understand the impact on members, and contribution risk measures capturing the risk from the employer’s viewpoint. We apply the framework to a hypothetical and empirically relevant merger and find that consolidation is favorable under most circumstances: the positive impacts of better diversification and economies of scale continue to outweigh the negative effects of heterogeneity even when the merging plans have different mortality expectations, different maturity levels, or modest differences in initial funded ratios.

We study a discrete-time life cycle retirement planning problem for individual workers with four distinct investment options: self-management with dynamic investment (S), self-management with benchmark investment (B), hire-management with flexible allocation ($\text{H}_{1}$), and hire-management with alpha focus ($\text{H}_{2}$). We examine the investment strategies and consumption patterns during the defined contribution fund accumulation period, ending with a life annuity purchase at retirement to finance post-retirement consumption. Based on the calibrated model using US data, we employ numerical dynamic programming technique to optimize worker’s financial decisions. Our analysis reveals that, despite the agency risk, delegated investments can add value to a worker’s lifetime utility, with the $\text{H}_2$ option yielding the best lifetime utility outcome. However, after taking the fund management fee into consideration, we find that both the $\text{H}_1$ and $\text{H}_2$ options may not offer additional value compared to the S option, yet they still surpass the B option in performance.

Risk measurement and econometrics are the two pillars of actuarial science. Unlike econometrics, risk measurement allows taking into account decision-makers’ risk aversion when analyzing the risks. We propose a hybrid model that captures decision-makers’ regression-based approach to study risks, focusing on explanatory variables while paying attention to risk severity. Our model considers different loss functions that quantify the severity of the losses that are provided by the risk manager or the actuary. We present an explicit formula for the regression estimators for the proposed risk-based regression problem and study the proposed results. Finally, we provide a numerical study of the results using data from the insurance industry.

The association between economic variables and the frequency and duration of disability income insurance (DII) claims is well established. Across many jurisdictions, heightened levels of unemployment have been associated with both a higher incidence and a longer duration of DII claims. This motivated us to derive an asset portfolio for which the total asset value moves in line with the level of unemployment, thus, providing a natural match for the DII portfolio liabilities. To achieve this, we develop an economic tracking portfolio where the asset weights in the portfolio are chosen so that the portfolio value changes in a way that reflects, as closely as possible, the level of unemployment. To the best of our knowledge, this is the first paper applying economic tracking portfolios to hedge economic risk in DII. The methodology put forward to establish this asset-liability matching portfolio is illustrated using DII data from the UK between 2004 and 2016. The benefits of our approach for claims reserving in DII portfolios are illustrated using a simulation study.

Many pension plans and private retirement products contain annuity factors, converting the funds at some future time into lifelong income. In general model settings like, for example, the Li-Lee mortality model, analytical values for the annuity factors are not available and one has to rely on numerical techniques. Their computation typically requires nested simulations as they depend on the interest rate level and the mortality tables at the time of retirement. We exploit the flexibility and efficiency of feed-forward neural networks (NNs) to value the annuity factors at the time of retirement. In a numerical study, we compare our deep learning approach to (least-squares) Monte-Carlo, which can be represented as a special case of the NN.

The extended $3/2$ short rate model is a mean reverting model of the short rate which, for suitably chosen parameters, permits a sensible term structure of bond yields and closed-form valuation formulae of zero-coupon bonds and options. This article supplies proofs of the formulae for the expected present values of future cash flows under the real-world probability measure, known as actuarial valuation. Finally, we give formulae for asymptotic levels of bond yields and formulae for bond option prices for the extended $3/2$ model, under particular conditions on its parameters.

We consider a new approach in the definition of two-dimensional heavy-tailed distributions. Specifically, we introduce the classes of two-dimensional long-tailed, of two-dimensional dominatedly varying, and of two-dimensional consistently varying distributions. Next, we define the closure property with respect to two-dimensional convolution and to joint max-sum equivalence in order to study whether they are satisfied by these classes. Further, we examine the joint-tail behavior of two random sums, under generalized tail asymptotic independence. Afterward, we study the closure property under scalar product and two-dimensional product convolution, and by these results, we extended our main result in the case of jointly randomly weighted sums. Our results contained some applications where we establish the asymptotic expression of the ruin probability in a two-dimensional discrete-time risk model.

We design the insurance contract when the insurer faces arson-type risks. We show that the optimal contract must be manipulation-proof. Therefore, it is continuous, has a bounded slope, and satisfies the no-sabotage condition when arson-type actions are free. Any contract that mixes a deductible, coinsurance, and an upper limit is manipulation-proof. A key feature of our models is that we provide a simple, general, and entirely elementary proof of manipulation-proofness that is easily adapted to different settings. We also show that the ability to perform arson-type actions reduces the insured’s welfare as less coverage is offered in equilibrium.

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.

This paper utilizes neural networks (NNs) for cycle detection in the insurance industry. The efficacy of NNs is compared on simulated data to the standard methods used in the underwriting cycles literature. The results show that NN models perform well in detecting cycles even in the presence of outliers and structural breaks. The methodology is applied to a granular data set of prices per risk profile from the Brazilian insurance industry.

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.

We investigate bottom-up risk aggregation applied by insurance companies facing reserve risk from multiple lines of business. Since risk capitals should be calculated in different time horizons and calendar years, depending on the regulatory or reporting regime (Solvency II vs IFRS 17), we study correlations of ultimate losses and correlations of one-year losses in future calendar years in lines of business. We consider a multivariate version of a Hertig’s lognormal model and we derive analytical formulas for the ultimate correlation and the one-year correlations in future calendar years. Our main conclusion is that the correlation coefficients that should be used in a bottom-up aggregation formula depend on the time horizon and the future calendar year where the risk emerges. We investigate analytically and numerically properties of the ultimate and the one-year correlations, their possible values observed in practice, and the impact of misspecified correlations on the diversified risk capital.