Robustly Hedging Variable Annuities With Guarantees Under Jump and Volatility Risks

Recent variable annuities offer participation in the equity market and attractive protection against downside movements. Accurately quantifying this additional equity market risk and robustly hedging options embedded in the guarantees of variable annuities are new challenges for insurance companies. Due to sensitivities of the benefits to tails of the account value distribution, a simple Black–Scholes model is inadequate in preventing excessive liabilities. A model which realistically describes the real world price dynamics over a long time horizon is essential for the risk management of the variable annuities. In this article, both jump risk and volatility risk are considered for risk management of lookback options embedded in guarantees with a ratchet feature. We evaluate relative performances of delta hedging and dynamic discrete risk minimization hedging strategies. Using the underlying as the hedging instrument, we show that, under a Black–Scholes model, local risk minimization hedging can be significantly better than delta hedging. In addition, we compare risk minimization hedging using the underlying with that of using standard options. We demonstrate that, under a Merton's jump diffusion model, hedging using standard options is superior to hedging using the underlying in terms of the risk reduction. Finally, we consider a market model for volatility risks in which the at‐the‐money implied volatility is a state variable. We compute risk minimization hedging by modeling at‐the‐money Black–Scholes implied volatility explicitly; the hedging effectiveness is evaluated, however, under a joint model for the underlying price and implied volatility. Our computational results suggest that, when implied volatility risk is suitably modeled, risk minimization hedging using standard options, compared to hedging using the underlying, can potentially be more effective in risk reduction under both jump and volatility risks.     

Implied integrated variance and hedging

If the volatility is stochastic, stock price returns and European option prices depend on the time average of the variance, i.e. the integrated variance, not on the path of the volatility. Applying a Bayesian statistical approach, we compute a forward-looking estimate of this variance, an option-implied integrated variance. Simultaneously, we obtain estimates of the correlation coefficient between stock price and volatility shocks, and of the parameters of the volatility process. Due to the convexity of the Black–Scholes formula with respect to the volatility, pricing and hedging with Black–Scholes-type formulas and the implied volatility often lead to inaccuracies if the volatility is stochastic. Theoretically, this problem can be avoided by using Hull–White-type option pricing and hedging formulas and the integrated variance. We use the implied integrated variance and Hull–White-type formulas to hedge European options and certain volatility derivatives.     

An Empirical Comparison of Two Stochastic Volatility Models using Indian Market Data

We conduct an empirical comparison of hedging strategies for two different stochastic volatility models proposed in the literature. One is an asymptotic expansion approach and the other is the risk-minimizing approach applied to a Markov-switched geometric Brownian motion. We also compare these with the Black–Scholes delta hedging strategies using historical and implied volatilities. The derivatives we consider are European call options on the NIFTY index of the Indian National Stock Exchange. We compare a few cases with profit and loss data from a trading desk. We find that for the cases that we analyzed, by far the better results are obtained for the Markov-switched geometric Brownian motion.     

Cross-sectional tests of deterministic volatility functions

We study the cross-sectional performance of option pricing models in which the volatility of the underlying stock is a deterministic function of the stock price and time. For each date in our sample of FTSE 100 index option prices, we fit an implied binomial tree to the panel of all European style options with different strike prices and maturities and then examine how well this model prices a corresponding panel of American style options. We find that the implied binomial tree model performs no better than an ad-hoc procedure of smoothing Black–Scholes implied volatilities across strike prices and maturities. Our cross-sectional results complement the time-series findings of Dumas et al. [J. Finance 53 (1998) 2059].     

Generative Bayesian neural network model for risk-neutral pricing of American index options

Financial models with stochastic volatility or jumps play a critical role as alternative option pricing models for the classical Black–Scholes model, which have the ability to fit different market volatility structures. Recently, machine learning models have elicited considerable attention from researchers because of their improved prediction accuracy in pricing financial derivatives. We propose a generative Bayesian learning model that incorporates a prior reflecting a risk-neutral pricing structure to provide fair prices for the deep ITM and the deep OTM options that are rarely traded. We conduct a comprehensive empirical study to compare classical financial option models with machine learning models in terms of model estimation and prediction using S&P 100 American put options from 2003 to 2012. Results indicate that machine learning models demonstrate better prediction performance than the classical financial option models. Especially, we observe that the generative Bayesian neural network model demonstrates the best overall prediction performance.     

A revised option pricing formula with the underlying being banned from short selling

An important issue in derivative pricing that hasn't been explored much until very recently is the impact of short selling to the price of an option. This paper extends a recent publication in this area to the case in which a ban of short selling of the underlying alone is somewhat less ‘effective’ than the extreme case discussed by Guo and Zhu [Equal risk pricing under convex trading constraints. J. Econ. Dyn. Control, 2017, 76, 136–151]. The case presented here is closer to reality, in which the effect of a ban on the underlying of an option alone may quite often be ‘diluted’ due to market interactions of the underlying asset with other correlated assets. Under a new assumption that there exists at least a correlated asset in the market, which is allowed to be short sold and thus can be used by traders for hedging purposes even though short selling of the underlying itself is banned, a new closed-form equal-risk pricing formula for European options is successfully derived. The new formula contains two distinguishable advantages; (a) it does not induce any significantly extra burden in terms of numerically computing option values, compared with the effort involved in using the Black–Scholes formula, which is still popularly used in finance industry today; (b) it remains simple and elegant as only one additional parameter beyond the Black–Scholes formula is introduced, to reflect the dilution effect to the ban as a result of market interactions.     

Implied volatility and skewness surface

The Homoscedastic Gamma (HG) model characterizes the distribution of returns by its mean, variance and an independent skewness parameter. The HG model preserves the parsimony and the closed form of the Black–Scholes–Merton (BSM) while introducing the implied volatility (IV) and skewness surface. Varying the skewness parameter of the HG model can restore the symmetry of IV curves. Practitioner’s variants of the HG model improve pricing (in-sample and out-of-sample) and hedging performances relative to practitioners’ BSM models, with as many or less parameters. The pattern of improvements in Delta-Hedged gains across strike prices accord with predictions from the HG model. These results imply that expanding around the Gaussian density does not offer sufficient flexibility to match the skewness implicit in options. Consistent with the model, we also find that conditioning on implied skewness increases the predictive power of the volatility spread for excess returns.     

Are classical option pricing models consistent with observed option second-order moments? Evidence from high-frequency data

As a means of validating an option pricing model, we compare the ex-post intra-day realized variance of options with the realized variance of the associated underlying asset that would be implied using assumptions as in the Black and Scholes (BS) model, the Heston, and the Bates model. Based on data for the S&P 500 index, we find that the BS model is strongly directionally biased due to the presence of stochastic volatility. The Heston model reduces the mismatch in realized variance between the two markets, but deviations are still significant. With the exception of short-dated options, we achieve best approximations after controlling for the presence of jumps in the underlying dynamics. Finally, we provide evidence that, although heavily biased, the realized variance based on the BS model contains relevant predictive information that can be exploited when option high-frequency data is not available.     

Option pricing based on hybrid GARCH-type models with improved ensemble empirical mode decomposition

The exploration of option pricing is of great significance to risk management and investments. One important challenge to existing research is how to describe the underlying asset price process and fluctuation features accurately. Considering the benefits of ensemble empirical mode decomposition (EEMD) in depicting the fluctuation features of financial time series, we construct an option pricing model based on the new hybrid generalized autoregressive conditional heteroskedastic (hybrid GARCH)-type functions with improved EEMD by decomposing the original return series into the high frequency, low frequency and trend terms. Using the locally risk-neutral valuation relationship (LRNVR), we obtain an equivalent martingale measure and option prices with different maturities based on Monte Carlo simulations. The empirical results indicate that this novel model can substantially capture volatility features and it performs much better than the M-GARCH and Black–Scholes models. In particular, the decomposition is consistently helpful in reducing option pricing errors, thereby proving the innovativeness and effectiveness of the hybrid GARCH option pricing model.     

A new closed-form formula for pricing European options under a skew Brownian motion

In this paper, we present a new pricing formula based on a modified Black–Scholes (B-S) model with the standard Brownian motion being replaced by a particular process constructed with a special type of skew Brownian motions. Although Corns and Satchell [2007. “Skew Brownian Motion and Pricing European Options.” The European Journal of Finance 13 (6): 523–544] have worked on this model, the results they obtained are incorrect. In this paper, not only do we identify precisely where the errors in Although Corns and Satchell [2007. “Skew Brownian Motion and Pricing European Options”. The European Journal of Finance 13 (6): 523–544] are, we also present a new closed-form pricing formula based on a newly proposed equivalent martingale measure, called ‘endogenous risk neutral measure’, by which only endogenous risks should and can be fully hedged. The newly derived option pricing formula takes the B-S formula as a special case and it does not induce any significant additional burden in terms of numerically computing option values, compared with the effort involved in computing the B-S formula.     

Model risk of the implied GARCH-normal model

Model risk causes significant losses in financial derivative pricing and hedging. Investors may undertake relatively risky investments due to insufficient hedging or overpaying implied by flawed models. The GARCH model with normal innovations (GARCH-normal) has been adopted to depict the dynamics of the returns in many applications. The implied GARCH-normal model is the one minimizing the mean square error between the market option values and the GARCH-normal option prices. In this study, we investigate the model risk of the implied GARCH-normal model fitted to conditional leptokurtic returns, an important feature of financial data. The risk-neutral GARCH model with conditional leptokurtic innovations is derived by the extended Girsanov principle. The option prices and hedging positions of the conditional leptokurtic GARCH models are obtained by extending the dynamic semiparametric approach of Huang and Guo [Statist. Sin., 2009, 19, 1037–1054]. In the simulation study we find significant model risk of the implied GARCH-normal model in pricing and hedging barrier and lookback options when the underlying dynamics follow a GARCH-t model.     

Hedging jump risk,expected returns and risk premia in jump-diffusion economies

This article presents a pure exchange economy that extends Rubinstein [Bell J. Econ. Manage. Sci., 1976, 7, 407–425] to show how the jump-diffusion option pricing model of Black and Scholes [J. Political Econ., 1973, 81, 637–654] and Merton [J. Financ. Econ., 1976, 4, 125–144] evolves in gamma jumping economies. From empirical analysis and theoretical study, both the aggregate consumption and the stock price are unknown in determining jumping times. By using the pricing kernel, we determine both the aggregate consumption jump time and the stock price jump time from the equilibrium interest rate and CCAPM (Consumption Capital Asset Pricing Model). Our general jump-diffusion option pricing model gives an explicit formula for how the jump process and the jump times alter the pricing. This innovation with predictable jump times enhances our analysis of the expected stock return in equilibrium and of hedging jump risks for jump-diffusion economies.     

An approximate distribution of delta-hedging errors in a jump-diffusion model with discrete trading and transaction costs

We introduce a jump-diffusion model for asset returns with jumps drawn from a mixture of normal distributions and show that this model adequately fits the historical data of the S&P500 index. We consider a delta-hedging strategy (DHS) for vanilla options under the diffusion model (DM) and the proposed jump-diffusion model (JDM), assuming discrete trading intervals and transaction costs, and derive an approximation for the probability density function (PDF) of the profit-and-loss (P&L) of the DHS under both models. We find that, under the log-normal model of Black–Scholes–Merton, the actual PDF of the P&L can be well approximated by the chi-squared distribution with specific parameters. We derive an approximation for the P&L volatility in the DM and JDM. We show that, under both DM and JDM, the expected loss due to transaction costs is inversely proportional to the square root of the hedging frequency. We apply mean–variance analysis to find the optimal hedging frequency given the hedger's risk tolerance. Since under the JDM it is impossible to reduce the P&L volatility by increasing the hedging frequency, we consider an alternative hedging strategy, following which the P&L volatility can be reduced by increasing the hedging frequency.     

Using the short-lived arbitrage model to compute minimum variance hedge ratios: application to indices,stocks and commodities

The short-lived arbitrage model has been shown to significantly improve in-sample option pricing fit relative to the Black–Scholes model. Motivated by this model, we imply both volatility and virtual interest rates to adjust minimum variance hedge ratios. Using several error metrics, we find that the hedging model significantly outperforms the traditional delta hedge and a current benchmark hedge based on the practitioner Black–Scholes model. Our applications include hedges of index options, individual stock options and commodity futures options. Hedges on gold and silver are especially sensitive to virtual interest rates.     

A hyperbolic diffusion model for stock prices

In the present paper we consider a model for stock prices which is a generalization of the model behind the Black–Scholes formula for pricing European call options. We model the log-price as a deterministic linear trend plus a diffusion process with drift zero and with a diffusion coefficient (volatility) which depends in a particular way on the instantaneous stock price. It is shown that the model possesses a number of properties encountered in empirical studies of stock prices. In particular the distribution of the adjusted log-price is hyperbolic rather than normal. The model is rather successfully fitted to two different stock price data sets. Finally, the question of option pricing based on our model is discussed and comparison to the Black–Scholes formula is made. The paper also introduces a simple general way of constructing a zero-drift diffusion with a given marginal distribution, by which other models that are potentially useful in mathematical finance can be developed.     

The pricing of basket-spread options

Since the pioneering paper of Black and Scholes was published in 1973, enormous research effort has been spent on finding a multi-asset variant of their closed-form option pricing formula. In this paper, we generalize the Kirk [Managing Energy Price Risk, 1995] approximate formula for pricing a two-asset spread option to the case of a multi-asset basket-spread option. All the advantageous properties of being simple, accurate and efficient are preserved. As the final formula retains the same functional form as the Black–Scholes formula, all the basket-spread option Greeks are also derived in closed form. Numerical examples demonstrate that the pricing and hedging errors are in general less than 1% relative to the benchmark results obtained by numerical integration or Monte Carlo simulation with 10 million paths. An implicit correction method is further applied to reduce the pricing errors by factors of up to 100. The correction is governed by an unknown parameter, whose optimal value is found by solving a non-linear equation. Owing to its simplicity, the computing time for simultaneous pricing and hedging of basket-spread option with 10 underlying assets or less is kept below 1 ms. When compared against the existing approximation methods, the proposed basket-spread option formula coupled with the implicit correction turns out to be one of the most robust and accurate methods.     

Hedging errors with Leland's option model in the presence of transaction costs

Nonzero transaction costs invalidate the Black–Scholes [1973. Journal of Political Economy 81, 637–654] arbitrage argument based on continuous trading. Leland [1985. Journal of Finance 40, 1283–1301] developed a hedging strategy which modifies the Black–Scholes hedging strategy with a volatility adjusted by the length of the rebalance interval and the rate of the proportional transaction cost. Kabanov and Safarian [1997. Finance and Stochastics 1, 239–250] calculated the limiting hedging error of the Leland strategy and pointed out that it is nonzero for the approximate pricing of an European call option, in contradiction to Leland's claim. As a further contribution, we first identify the mathematical flaw in the argument of Leland's claim and then quantify the expected percentage of hedging losses in terms of the hedging frequency and the level of the option strike price.     

Discrete time hedging with liquidity risk

We study a discrete time hedging and pricing problem in a market with liquidity costs. Using Leland’s discrete time replication scheme [Leland, H.E., 1985. Journal of Finance, 1283–1301], we consider a discrete time version of the Black–Scholes model and a delta hedging strategy. We derive a partial differential equation for the option price in the presence of liquidity costs and develop a modified option hedging strategy which depends on the size of the parameter for liquidity risk. We also discuss an analytic method of solving the pricing equation using a series solution.     

A PDE method for estimation of implied volatility

In this paper it is proved that the Black–Scholes implied volatility satisfies a second order non-linear partial differential equation. The obtained PDE is then used to construct an algorithm for fast and accurate polynomial approximation for Black–Scholes implied volatility that improves on the existing numerical schemes from literature, both in speed and parallelizability. We also show that the method is applicable to other problems, such as approximation of implied Bachelier volatility.