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.

In the traditional multidimensional credibility models developed by Jewell ((1973) Operations Research Center, pp. 73–77.), the estimation of the hypothetical mean vector involves complex matrix manipulations, which can be challenging to implement in practice. Additionally, the estimation of hyperparameters becomes even more difficult in high-dimensional risk variable scenarios. To address these issues, this paper proposes a new multidimensional credibility model based on the conditional joint distribution function for predicting future premiums. First, we develop an estimator of the joint distribution function of a vector of claims using linear combinations of indicator functions based on past observations. By minimizing the integral of the expected quadratic distance function between the proposed estimator and the true joint distribution function, we obtain the optimal linear Bayesian estimator of the joint distribution function. Using the plug-in method, we obtain an explicit formula for the multidimensional credibility estimator of the hypothetical mean vector. In contrast to the traditional multidimensional credibility approach, our newly proposed estimator does not involve a matrix as the credibility factor, but rather a scalar. This scalar is composed of both population information and sample information, and it still maintains the essential property of increasingness with respect to the sample size. Furthermore, the new estimator based on the joint distribution function can be naturally extended and applied to estimate the process covariance matrix and risk premiums under various premium principles. We further illustrate the performance of the new estimator by comparing it with the traditional multidimensional credibility model using bivariate exponential-gamma and multivariate normal distributions. Finally, we present two real examples to demonstrate the findings of our study.

This paper considers variable annuity (VA) contracts embedded with guaranteed minimum accumulation benefit (GMAB) riders when policyholder’s proceeds are taxed upon early surrender or maturity. These contracts promise the return of the premium paid by the policyholder, or a higher rolled-up value, at the end of the investment period. A partial differential equation valuation framework which exploits the numerical method of lines is used to determine fair fees that render the policyholder and insurer breakeven. Two taxation regimes are considered: one where capital gains are allowed to offset losses and a second where gains do not offset losses. Most insurance providers highlight the tax-deferred features of VA contracts. We show that the regime under which the insured is taxed significantly impacts prices. If losses are allowed to offset gains then this enhances the market, increasing the policyholder’s willingness to participate in the market compared to the case when losses are not allowed to offset gains. With fair fees from the policyholder’s perspective, we show that the net profit is generally positive for insurance companies offering the contract as a naked option without any hedge. We also show how investment policy, as reflected in the Sharpe ratio, impacts and interacts with policyholder persistency.

Modeling multivariate dependence in high dimensions is challenging, with popular solutions constructing multivariate copula as a composition of lower dimensional copulas. Pair-copula constructions do so by using bivariate linking copulas, but their parametrization, in size, being quadratic in the dimension, is not quite parsimonious. Besides, the number of regular vines grows super-exponentially with the dimension. One parsimonious solution is factor copulas, and in particular, the one-factor copula is touted for its simplicity – with the number of parameters linear in the dimension – while being able to cater to asymmetric non-linear dependence in the tails. In this paper, we add nuance to this claim from the point of view of a popular measure of multivariate tail dependence, the tail dependence matrix (TDM). We focus on the one-factor copula model with the linking copula belonging to the BB1 family, pointing out later the applicability of our results to a wider class of linking copulas. For this model, we derive tail dependence coefficients and study their basic properties as functions of the parameters of the linking copulas. Based on this, we study the representativeness of the class of TDMs supported by this model with respect to the class of all possible TDMs. We establish that since the parametrization is linear in the dimension, it is no surprise that the relative volume is zero for dimensions greater than three, and hence, by necessity, we present a novel manner of evaluating the representativeness that has a combinatorial flavor. We formulate the problem of finding the best representative one-factor BB1 model given a target TDM and suggest an implementation along with a simulation study of its performance across dimensions. Finally, we illustrate the results of the paper by modeling rainfall data, which is relevant in the context of weather-related insurance.

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.

Consider a company whose business carries the potential for investment losses and is additionally vulnerable to exogenous shocks. The unpredictability of the shocks makes it challenging for both the company and the regulator to accurately assess their impact, potentially leading to an underestimation of solvency capital when employing traditional approaches. In this paper, we utilize a stylized model to conduct an extreme value analysis of the tail risk of the company under a Fréchet-type and a Gumbel-type shock. Our main results explicitly demonstrate the different roles of investment risk and shock risk in driving large losses. Furthermore, we derive asymptotic estimates for the value at risk and expected shortfall of the total loss. Numerical studies are conducted to examine the accuracy of the obtained estimates.

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.

Reinsurers may default when they have to pay large claims to insurers but are unable to fulfill their obligations due to various reasons such as catastrophic events, underwriting losses, inadequate capitalization, or financial mismanagement. This paper studies the problem of optimal reinsurance design from the perspectives of both the insurer and reinsurer when the insurer faces the potential default risk of the reinsurer. If the insurer aims to minimize the convex distortion risk measure of his retained loss, we prove the optimality of a stop-loss treaty when the promised ceded loss function is charged by the expected value premium principle and the reinsurer offers partial recovery in the event of default. For any fixed premium loading set by the reinsurer, we then derive the explicit expressions of optimal deductible levels for three special distortion functions, including the TVaR, Gini, and PH transform distortion functions. Under these three explicit distortion risk measures adopted by the insurer, we seek the optimal safety loading for the reinsurer by maximizing her net profit where the reserve capital is determined by the TVaR measure and the cost is governed by the expectation. This procedure ultimately leads to the Bowley solution between the insurer and the reinsurer. We provide several numerical examples to illustrate the theoretical findings. Sensitivity analyses demonstrate how different settings of default probability, recovery rate, and safety loading affect the optimal deductible values. Simulation studies are also implemented to analyze the effects induced by the default probability and recovery rate on the Bowley solution.

This paper proposes a theoretical insurance model to explain well-documented loss underreporting and to study how strategic underreporting affects insurance demand. We consider a utility-maximizing insured who purchases a deductible insurance contract and follows a barrier strategy to decide whether she should report a loss. The insurer adopts a bonus-malus system with two rate classes, and the insured will move to or stay in the more expensive class if she reports a loss. First, we fix the insurance contract (deductibles) and obtain the equilibrium reporting strategy in semi-closed form. A key result is that the equilibrium barriers in both rate classes are strictly greater than the corresponding deductibles, provided that the insured economically prefers the less expensive rate class, thereby offering a theoretical explanation to underreporting. Second, we study an optimal deductible insurance problem in which the insured strategically underreports losses to maximize her utility. We find that the equilibrium deductibles are strictly positive, suggesting that full insurance, often assumed in related literature, is not optimal. Moreover, in equilibrium, the insured underreports a positive amount of her loss. Finally, we examine how underreporting affects the insurer’s expected profit.

The theory of utility is a well-known method of constructing insurance premiums (see e.g., Newton et al. (1986) Actuarial Mathematics. Itasca, Illinois: The Society of Actuaries.). Furman and Zitikis ((2008) Insurance: Mathematics and Economics, 42, 459–465.) proposed an alternative method using the mean value of a weighted random variable. According to this approach, for various choices of weighting, popular premiums such as net premium, modified variance premium, Esscher premium, and Kamps premium are obtained. On the other hand, some premiums cannot be obtained with this method, such as the premium of the exponential principle. In this paper, we provide a complementary theory by introducing a family of unimodal weighted distributions for which the mode is a premium principle.