We study two continuous-time Stackelberg games between a life insurance buyer and seller over a random time horizon. The buyer invests in a risky asset and purchases life insurance, and she maximizes a mean-variance criterion applied to her wealth at death. The seller chooses the insurance premium rate to maximize its expected wealth at the buyer’s random time of death. We consider two life insurance games: one with term life insurance and the other with whole life insurance—the latter with pre-commitment of the constant investment strategy. In the term life insurance game, the buyer chooses her life insurance death benefit and investment strategy continuously from a time-consistent perspective. We find the buyer’s equilibrium control strategy explicitly, along with her value function, for the term life insurance game by solving the extended Hamilton–Jacobi–Bellman equations. By contrast, in the whole life insurance game, the buyer pre-commits to a constant life insurance death benefit and a constant amount to invest in the risky asset. To solve the whole life insurance problem, we first obtain the buyer’s objective function and then we maximize that objective function over constant controls. Under both models, the seller maximizes its expected wealth at the buyer’s time of death, and we use the resulting optimal life insurance premia to find the Stackelberg equilibria of the two life insurance games. We also analyze the effects of the parameters on the Stackelberg equilibria, and we present some numerical examples to illustrate our results.

We provide explicit small-time formulae for the at-the-money implied volatility, skew, and curvature in a large class of models, including rough volatility models and their multi-factor versions. Our general setup encompasses both European options on a stock and VIX options, thereby providing new insights on their joint calibration. The tools used are essentially based on Malliavin calculus for Gaussian processes. We develop a detailed theoretical and numerical analysis of the two-factor rough Bergomi model and provide insights on the interplay between the different parameters for joint SPX–VIX smile calibration.

We present a closed-form solution to a discounted optimal stopping zero-sum game in a model based on a generalised geometric Brownian motion with coefficients depending on its running maximum and minimum processes. The optimal stopping times forming a Nash equilibrium are shown to be the first times at which the original process hits certain boundaries depending on the running values of the associated maximum and minimum processes. The proof is based on the reduction of the original game to the equivalent coupled free-boundary problem and the solution of the latter problem by means of the smooth-fit and normal-reflection conditions. We show that the optimal stopping boundaries are partially determined as either unique solutions to the appropriate system of arithmetic equations or unique solutions to the appropriate first-order nonlinear ordinary differential equations. The results obtained are related to the valuation of the perpetual lookback game options with floating strikes in the appropriate diffusion-type extension of the Black–Merton–Scholes model.

Complicated option pricing models attract much attention in financial industries, as they produce relatively better accurate values by taking into account more realistic assumptions such as market liquidity, uncertain volatility and so forth. We propose a new hybrid method to accurately explore the behaviour of the nonlinear pricing model in illiquid markets, which is important in financial risk management. Our method is based on the Newton iteration technique and the Fréchet derivative to linearize the model. The linearized equation is then discretized by a differential quadrature method in space and a quadratic trapezoid rule in time. It is observed through computations that the accurate solutions for the model emerge using very few grid points and time elements, compared with the finite difference method in the literature. Furthermore, this method also helps to avoid consideration of the convergence issues of the Newton approach applied to the nonlinear algebraic system containing many unknowns at each time step if an implicit method is used in time discretization. It is important to note that the Fréchet derivative supports to enhance the convergence order of the proposed iterative scheme.

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.

We construct a new stochastic interest rate model with two stochastic factors, by introducing a stochastic long-run equilibrium level into the Vasicek interest rate model which follows another Ornstein–Uhlenbeck process. With the interest rate under the Black–Scholes model being assumed to follow the newly proposed model, a closed-form representation of European option prices is successfully presented, when the analytical characteristic function of the underlying log-price under a forward measure is derived. To assess the model performance, a preliminary empirical study is conducted using S&P 500 index and its options, with the Vasicek model and an alternative two-factor Vasicek model taken as benchmarks.

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 consider the super-replication problem for a class of exotic options known as life-contingent options within the framework of the Black–Scholes market model. The option is allowed to be exercised if the death of the option holder occurs before the expiry date, otherwise there is a compensation payoff at the expiry date. We show that there exists a minimal super-replication portfolio and determine the associated initial investment. We then give a characterisation of when replication of the option is possible. Finally, we give an example of an explicit super-replicating hedge for a simple life-contingent option.

We investigate an optimal stopping problem for the expected value of a discounted payoff on a regime-switching geometric Brownian motion under two constraints on the possible stopping times: only at exogenous random times, and only during a specific regime. The main objectives are to show that an optimal stopping time exists as a threshold type and to derive expressions for the value functions and the optimal threshold. To this end, we solve the corresponding variational inequality and show that its solution coincides with the value functions. Some numerical results are also introduced. Furthermore, we investigate some asymptotic behaviors.

This paper considers the first passage times to constant boundaries and the two-sided exit problem for Lévy process with a characteristic exponent in which at least one of the two jumps having rational Laplace transforms. The joint distribution of the first passage times and undershoot/overshoot are obtained. The processes recover many models that have appeared in the literature such as the compound Poisson risk models, the perturbed compound Poisson risk models, and their dual ones. As applications, we obtain the solutions for popular path-dependent options such as lookback and barrier options in terms of Laplace transforms.

We show that an estimate by de la Peña, Ibragimov, and Jordan for ${\mathbb{E}}(X-c)^+$, with c a constant and X a random variable of which the mean, the variance, and $\mathbb{P}(X \leqslant c)$ are known, implies an estimate by Scarf on the infimum of ${\mathbb{E}}(X \wedge c)$ over the set of positive random variables X with fixed mean and variance. This also shows, as a consequence, that the former estimate implies an estimate by Lo on European option prices.

We formulate a centrally planned portfolio selection problem with the investor and the manager having S-shaped utilities under a recently popular first-loss contract. We solve for the closed-form optimal portfolio, which shows that a first-loss contract can sometimes behave like an option contract. We propose an asymptotic approach to investigate the portfolio. This approach can be adopted to illustrate economic insights, including the fact that the portfolio under a convex contract becomes more conservative when the market state is better. Furthermore, we discover a means of Pareto improvement by simultaneously considering the investor’s utility and increasing the manager’s incentive rate. This is achieved by establishing the collection of Pareto points of a single contract, proving that it is a strictly decreasing and strictly concave frontier, and comparing the Pareto frontiers of different contracts. These results may be helpful for the illustration of risk choices and the design of Pareto-optimal contracts.

This paper examines the issue of derivative pricing within the framework of a fractional stochastic volatility model. We present a deterministic partial differential equation system to derive an approximate expression for the derivative price. The proposed approach allows for the stochastic volatility to be expressed as a composition of deterministic functions of time and a fractional Ornstein–Uhlenbeck process. We apply this method to the European option pricing under the fractional Stein–Stein volatility model, demonstrating its feasibility and reliability through numerical simulations. Our numerical simulations also illustrate the impact of the parameters in the fractional stochastic volatility model on the option price.

In Chen and Liang previous work, a model, together with its well-posedness, was established for credit rating migrations with different upgrade and downgrade thresholds (i.e. a buffer zone, also called dead band in engineering). When positive dividends are introduced, the model in Chen and Liang (SIAM Financ. Math. 12, 941–966, 2021) may not be well-posed. Here, in this paper, a new model is proposed to include the realistic nonzero dividend scenarios. The key feature of the new model is that partial differential equations in Chen and Liang (SIAM Financ. Math. 12, 941–966, 2021) are replaced by variational inequalities, thereby creating a new free boundary, besides the original upgrading and downgrading free boundaries. Well-posedness of the new model, together with a few financially meaningful properties, is established. In particular, it is shown that when time to debt paying deadline is long enough, the underlying dividend paying company is always in high grade rating, that is, only when time to debt paying deadline is within a certain range, there can be seen the phenomenon of credit rating migration.

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 we obtain a duality result for the exponential utility maximization problem where trading is subject to quadratic transaction costs and the investor is required to liquidate her position at the maturity date. As an application of the duality, we treat utility-based hedging in the Bachelier model. For European contingent claims with a quadratic payoff, we compute the optimal trading strategy explicitly.

We define and study properties of implied volatility for American perpetual put options. In particular, we show that if the market prices are derived from a local volatility model with a monotone volatility function, then the corresponding implied volatility is also monotone as a function of the strike price.

Consider the problem of determining the Bayesian credibility mean $E(X_{n+1}|X_1,\cdots, X_n),$ whenever the random claims $X_1,\cdots, X_n,$ given parameter vector $\boldsymbol{\Psi},$ are sampled from the K-component mixture family of distributions, whose members are the union of different families of distributions. This article begins by deriving a recursive formula for such a Bayesian credibility mean. Moreover, under the assumption that using additional information $Z_{i,1},\cdots,Z_{i,m},$ one may probabilistically determine a random claim $X_i$ belongs to a given population (or a distribution), the above recursive formula simplifies to an exact Bayesian credibility mean whenever all components of the mixture distribution belong to the exponential families of distributions. For a situation where a 2-component mixture family of distributions is an appropriate choice for data modelling, using the logistic regression model, it shows that: how one may employ such additional information to derive the Bayesian credibility model, say Logistic Regression Credibility model, for a finite mixture of distributions. A comparison between the Logistic Regression Credibility (LRC) model and its competitor, the Regression Tree Credibility (RTC) model, has been given. More precisely, it shows that under the squared error loss function, it shows the LRC’s risk function dominates the RTC’s risk function at least in an interval which about $0.5.$ Several examples have been given to illustrate the practical application of our findings.

We study a stochastic differential equation with an unbounded drift and general Hölder continuous noise of order $\lambda \in (0,1)$. The corresponding equation turns out to have a unique solution that, depending on a particular shape of the drift, either stays above some continuous function or has continuous upper and lower bounds. Under some mild assumptions on the noise, we prove that the solution has moments of all orders. In addition, we provide its connection to the solution of some Skorokhod reflection problem. As an illustration of our results and motivation for applications, we also suggest two stochastic volatility models which we regard as generalizations of the CIR and CEV processes. We complete the study by providing a numerical scheme for the solution.