This paper investigates the precise large deviations of the net loss process in a two-dimensional risk model with consistently varying tails and dependence structures, and gives some asymptotic formulas which hold uniformly for all x varying in t-intervals. The study is among the initial efforts to analyze potential risk via large deviation results for the net loss process of the two-dimensional risk model, and can provide a novel insight to assess the operation risk in a long run by fully considering the premium income factors of the insurance company.

We consider a robust optimal investment–reinsurance problem to minimize the goal-reaching probability that the value of the wealth process reaches a low barrier before a high goal for an ambiguity-averse insurer. The insurer invests its surplus in a constrained financial market, where the proportion of borrowed amount to the current wealth level is no more than a given constant, and short-selling is prohibited. We assume that the insurer purchases per-claim reinsurance to transfer its risk exposure to a reinsurer whose premium is computed according to the mean–variance premium principle. Using the stochastic control approach based on the Hamilton–Jacobi–Bellman equation, we derive robust optimal investment–reinsurance strategies and the corresponding value functions. We conclude that the behavior of borrowing typically occurs with a lower wealth level. Finally, numerical examples are given to illustrate our results.

This paper studies a bi-dimensional compound risk model with quasi-asymptotically independent and consistently varying-tailed random numbers of claims and establishes an asymptotic formula for the finite-time sum-ruin probability. Additionally, some results related to tail probabilities of random sums are presented, which are of significant interest in their own right. Some numerical studies are carried out to check the accuracy of the asymptotic formula.

In this paper, we explore the applications of Tail Variance (TV) as a measure of tail riskiness and the confidence level of using Tail Conditional Expectation (TCE)-based risk capital. While TCE measures the expected loss of a risk that exceeds a certain threshold, TV measures the variability of risk along its tails. We first derive analytical formulas of TV and TCE for a large variety of probability distributions. These formulas are useful instruments for relevant research works on tail risk measures. We then propose a distribution-free approach utilizing TV to estimate the lower bounds of the confidence level of using TCE-based risk capital. In doing so, we introduce sharpened conditional probability inequalities, which halve the bounds of conventional Markov and Cantelli inequalities. Such an approach is easy to implement. We further investigate the characterization of tail risks by TV through an exploration of TV’s asymptotics. A distribution-free limit formula is derived for the asymptotics of TV. To further investigate the asymptotic properties, we consider two broad distribution families defined on tails, namely, the polynomial-tailed distributions and the exponential-tailed distributions. The two distribution families are found to exhibit an asymptotic equivalence between TV and the reciprocal square of the hazard rate. We also establish asymptotic relationships between TCE and VaR for the two families. Our asymptotic analysis contributes to the existing research by unifying the asymptotic expressions and the convergence rate of TV for Student-t distributions, exponential distributions, and normal distributions, which complements the discussion on the convergence rate of univariate cases in [28]. To show the usefulness of our results, we present two case studies based on real data from the industry. We first show how to use conditional inequalities to assess the confidence of using TCE-based risk capital for different types of insurance businesses. Then, for financial data, we provide alternative evidence for the relationship between the data frequency and the tail categorization by the asymptotics of TV.

A new risk measure, the Lambda Value-at-Risk (VaR), was proposed from a theoretical point of view as a generalization of the ordinary VaR in the literature. Motivated by the recent developments in risk sharing problems for the VaR and other risk measures, we study the optimization of risk sharing for the Lambda VaR. Explicit formulas of the inf-convolution and sum-optimal allocations are obtained with respect to the left Lambda VaRs, the right Lambda VaRs, or a mixed collection of the left and right Lambda VaRs. The inf-convolution of Lambda VaRs constrained to comonotonic allocations is investigated. Explicit formula for worst-case Lambda VaRs under model uncertainty induced by likelihood ratios is also given.

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.

Multivariate regular variation is a key concept that has been applied in finance, insurance, and risk management. This paper proposes a new dependence assumption via a framework of multivariate regular variation. Under the condition that financial and insurance risks satisfy our assumption, we conduct asymptotic analyses for multidimensional ruin probabilities in the discrete-time and continuous-time cases. Also, we present a two-dimensional numerical example satisfying our assumption, through which we show the accuracy of the asymptotic result for the discrete-time multidimensional insurance risk model.

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.

Traditional techniques for calculating outstanding claim liabilities such as the chain-ladder are notoriously at risk of being distorted by outliers in past claims data. Unfortunately, the literature in robust methods of reserving is scant, with notable exceptions such as Verdonck & Debruyne (2011, Insurance: Mathematics and Economics, 48, 85–98) and Verdonck & Van Wouwe (2011, Insurance: Mathematics and Economics, 49, 188–193). In this paper, we put forward two alternative robust bivariate chain-ladder techniques to extend the approach of Verdonck & Van Wouwe (2011, Insurance: Mathematics and Economics, 49, 188–193). The first technique is based on Adjusted Outlyingness (Hubert & Van der Veeken, 2008. Journal of Chemometrics, 22, 235–246) and explicitly incorporates skewness into the analysis while providing a unique measure of outlyingness for each observation. The second technique is based on bagdistance (Hubert et al., 2016. Statistics: Methodology, 1–23) which is derived from the bagplot; however; it is able to provide a unique measure of outlyingness and a means to adjust outlying observations based on this measure.

Furthermore, we extend our robust bivariate chain-ladder approach to an N-dimensional framework. The implementation of the methods, especially beyond bivariate, is not trivial. This is illustrated on a trivariate data set from Australian general insurers and results under the different outlier detection and treatment mechanisms are compared.

In this paper, the pricing of equity warrants under a class of fractional Brownian motion models is investigated numerically. By establishing a new nonlinear partial differential equation (PDE) system governing the price in terms of the observable stock price, we solve the pricing system effectively by a robust implicit-explicit numerical method. This is fundamentally different from the documented methods, which first solve the price with respect to the firm value analytically, by assuming that the volatility of the firm is constant, and then compute the price with respect to the stock price and estimate the firm volatility numerically. It is shown that the proposed method is stable in the maximum-norm sense. Furthermore, a sharp theoretical error estimate for the current method is provided, which is also verified numerically. Numerical examples suggest that the current method is efficient and can produce results that are, overall, closer to real market prices than other existing approaches. A great advantage of the current method is that it can be extended easily to price equity warrants under other complicated models.

We establish new results on the strictly stationary solution to an iterated function system. When the driving sequence is stationary and ergodic, though not independent, the strictly stationary solution may admit no moment but we show an exponential control of the trajectories. We exploit these results to prove, under mild conditions, the consistency of the quasi-maximum likelihood estimator of GARCH(p,q) models with non-independent innovations.

In this paper, we consider a financial or insurance system with a finite number of individual risks described by real-valued random variables. We focus on two kinds of risk measures, referred to as the tail moment (TM) and the tail central moment (TCM), which are defined as the conditional moment and conditional central moment of some individual risk in the event of system crisis. The first-order TM and the second-order TCM coincide with the popular risk measures called the marginal expected shortfall and the tail variance, respectively. We derive asymptotic expressions for the TM and TCM with any positive integer orders, when the individual risks are pairwise asymptotically independent and have distributions from certain classes that contain both light-tailed and heavy-tailed distributions. The formulas obtained possess concise forms unrelated to dependence structures, and hence enable us to estimate the TM and TCM efficiently. To demonstrate the wide application of our results, we revisit some issues related to premium principles and optimal capital allocation from the asymptotic point of view. We also give a numerical study on the relative errors of the asymptotic results obtained, under some specific scenarios when there are two individual risks in the system. The corresponding asymptotic properties of the degenerate univariate versions of the TM and TCM are discussed separately in an appendix at the end of the paper.

The survival energy model (SEM) is a recently introduced novel approach to mortality prediction, which offers a cohort-wise distribution function of the time of death as the first hitting time of a “survival energy” diffusion process to zero. In this study, we propose a novel SEM that can serve as a suitable candidate in the family of prediction models. We also proposed a method to improve the prediction in an earlier work. We further examine the practical advantages of SEM over existing mortality models.

Risk measurements are clearly central to risk management, in particular for banks, (re)insurance companies, and investment funds. The question of the appropriateness of risk measures for evaluating the risk of financial institutions has been heavily debated, especially after the financial crisis of 2008/2009. Another concern for financial institutions is the pro-cyclicality of risk measurements. In this paper, we extend existing work on the pro-cyclicality of the Value-at-Risk to its main competitors, Expected Shortfall, and Expectile: We compare the pro-cyclicality of historical quantile-based risk estimation, taking into account the market state. To characterise the latter, we propose various estimators of the realised volatility. Considering the family of augmented GARCH(p, q) processes (containing well-known GARCH models and iid models, as special cases), we prove that the strength of pro-cyclicality depends on the three factors: the choice of risk measure and its estimators, the realised volatility estimator and the model considered, but, no matter the choices, the pro-cyclicality is always present. We complement this theoretical analysis by performing simulation studies in the iid case and developing a case study on real data.

This article examines the impact of the largest claims reinsurance treaties on loss reserve of the ceding company. The largest claims reinsurance, known as LCR, and ECOMOR reinsurance treaties are considered to be the two most appropriate reinsurance treaties for large or catastrophe claims. Then, it studies the impact of such treaties on loss reserves. Through a simulation study, it shown that, under a more general situation, the LCR treaty can be a more efficient (in some sense, see below) treaty than the ECOMOR treaty for the ceding company.

We propose a generalized Cramér–Lundberg model of the risk theory of non-life insurance and study its ruin probability. Our model is an extension of that of Dubey (1977) to the case of multiple insureds, where the counting process is a mixed Poisson process and the continuously varying premium rate is determined by a Bayesian rule on the number of claims. We use two proofs to show that, for each fixed value of the safety loading, the ruin probability is the same as that of the classical Cramér–Lundberg model and does not depend on either the distribution of the mixing variable of the driving mixed Poisson process or the number of claim contracts.

This paper considers a variant of the classical Cramér–Lundberg model that is particularly appropriate in the credit context, with the distinguishing feature that it corresponds to a finite number of obligors. The focus is on computing the ruin probability, i.e. the probability that the initial reserve, increased by the interest received from the obligors and decreased by the losses due to defaults, drops below zero. As well as an exact analysis (in terms of transforms) of this ruin probability, an asymptotic analysis is performed, including an efficient importance-sampling-based simulation approach.

The base model is extended in multiple dimensions: (i) we consider a model in which there may, in addition, be losses that do not correspond to defaults, (ii) then we analyze a model in which the individual obligors are coupled via a regime switching mechanism, (iii) then we extend the model so that between the losses the reserve process behaves as a Brownian motion rather than a deterministic drift, and (iv) we finally consider a set-up with multiple groups of statistically identical obligors.

We consider the optimal portfolio and consumption problem for a jump-diffusion process with regime switching. Under the criterion of maximizing the expected discounted total utility of consumption, two methods, namely, the dynamic programming principle and the stochastic maximum principle, are used to obtain the optimal result for the general objective function, which is the solution to a system of partial differential equations. Furthermore, we investigate the power utility as a specific example and analyse the existence and uniqueness of the optimal solution. Under the constraints of no-short-selling and nonnegative consumption, closed-form expressions for the optimal strategy and the value function are derived. Besides, some comparisons between the optimal results for the jump-diffusion model and the pure diffusion model are carried out. Finally, we discuss our optimal results in some special cases.

Rough volatility is a well-established statistical stylized fact of financial assets. This property has led to the design and analysis of various new rough stochastic volatility models. However, most of these developments have been carried out in the mono-asset case. In this work, we show that some specific multivariate rough volatility models arise naturally from microstructural properties of the joint dynamics of asset prices. To do so, we use Hawkes processes to build microscopic models that accurately reproduce high-frequency cross-asset interactions and investigate their long-term scaling limits. We emphasize the relevance of our approach by providing insights on the role of microscopic features such as momentum and mean-reversion in the multidimensional price formation process. In particular, we recover classical properties of high-dimensional stock correlation matrices.