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 endogenous volatility approach to pricing and hedging call options with transaction costs

Standard delta hedging fails to exactly replicate a European call option in the presence of transaction costs. We study a pricing and hedging model similar to the delta hedging strategy with an endogenous volatility parameter for the calculation of delta over time. The endogenous volatility depends on both the transaction costs and the option strike prices. The optimal hedging volatility is calculated using the criterion of minimizing the weighted upside and downside replication errors. The endogenous volatility model with equal weights on the up and down replication errors yields an option premium close to the Leland [J. Finance, 1985 Leland, HE. 1985. Option pricing and replication with transaction costs. J. Finance, 40: 1283–1301. [Crossref], [Web of Science ®] , [Google Scholar], 40, 1283–1301] heuristic approach. The model with weights being the probabilities of the option's moneyness provides option prices closest to the actual prices. Option prices from the model are identical to the Black–Scholes option prices when transaction costs are zero. Data on S&P 500 index cash options from January to June 2008 illustrate the model.     

A novel Monte Carlo approach to hybrid local volatility models

We present in a Monte Carlo simulation framework, a novel approach for the evaluation of hybrid local volatility [Risk, 1994, 7, 18–20], [Int. J. Theor. Appl. Finance, 1998, 1, 61–110] models. In particular, we consider the stochastic local volatility model—see e.g. Lipton et al. [Quant. Finance, 2014, 14, 1899–1922], Piterbarg [Risk, 2007, April, 84–89], Tataru and Fisher [Quantitative Development Group, Bloomberg Version 1, 2010], Lipton [Risk, 2002, 15, 61–66]—and the local volatility model incorporating stochastic interest rates—see e.g. Atlan [ArXiV preprint math/0604316, 2006], Piterbarg [Risk, 2006, 19, 66–71], Deelstra and Rayée [Appl. Math. Finance, 2012, 1–23], Ren et al. [Risk, 2007, 20, 138–143]. For both model classes a particular (conditional) expectation needs to be evaluated which cannot be extracted from the market and is expensive to compute. We establish accurate and ‘cheap to evaluate’ approximations for the expectations by means of the stochastic collocation method [SIAM J. Numer. Anal., 2007, 45, 1005–1034], [SIAM J. Sci. Comput., 2005, 27, 1118–1139], [Math. Models Methods Appl. Sci., 2012, 22, 1–33], [SIAM J. Numer. Anal., 2008, 46, 2309–2345], [J. Biomech. Eng., 2011, 133, 031001], which was recently applied in the financial context [Available at SSRN 2529691, 2014], [J. Comput. Finance, 2016, 20, 1–19], combined with standard regression techniques. Monte Carlo pricing experiments confirm that our method is highly accurate and fast.     

Pricing variance and volatility swaps in a stochastic volatility model with regime switching: discrete observations case

This study presents a set of closed-form exact solutions for pricing discretely sampled variance swaps and volatility swaps, based on the Heston stochastic volatility model with regime switching. In comparison with all the previous studies in the literature, this research, which obtains closed-form exact solutions for variance and volatility swaps with discrete sampling times, serves several purposes. (1) It verifies the degree of validity of Elliott et al.'s [Appl. Math. Finance, 2007, 14(1), 41–62] continuous-sampling-time approximation for variance and volatility swaps of relatively short sampling periods. (2) It examines the effect of ignoring regime switching on pricing variance and volatility swaps. (3) It contributes to bridging the gap between Zhu and Lian's [Math. Finance, 2011, 21(2), 233–256] approach and Elliott et al.'s framework. (4) Finally, it presents a semi-Monte-Carlo simulation for the pricing of other important realized variance based derivatives.     

Options on realized variance by transform methods: a non-affine stochastic volatility model

In this paper we study the pricing and hedging of options on realized variance in the 3/2 non-affine stochastic volatility model by developing efficient transform-based pricing methods. This non-affine model gives prices of options on realized variance that allow upward-sloping implied volatility of variance smiles. Heston's model [Rev. Financial Stud., 1993, 6, 327–343], the benchmark affine stochastic volatility model, leads to downward-sloping volatility of variance smiles—in disagreement with variance markets in practice. Using control variates, we propose a robust method to express the Laplace transform of the variance call function in terms of the Laplace transform of the realized variance. The proposed method works in any model where the Laplace transform of realized variance is available in closed form. Additionally, we apply a new numerical Laplace inversion algorithm that gives fast and accurate prices for options on realized variance, simultaneously at a sequence of variance strikes. The method is also used to derive hedge ratios for options on variance with respect to variance swaps.     

Turbocharging Monte Carlo pricing for the rough Bergomi model

The rough Bergomi model, introduced by Bayer et al. [Quant. Finance, 2016, 16(6), 887–904], is one of the recent rough volatility models that are consistent with the stylised fact of implied volatility surfaces being essentially time-invariant, and are able to capture the term structure of skew observed in equity markets. In the absence of analytical European option pricing methods for the model, we focus on reducing the runtime-adjusted variance of Monte Carlo implied volatilities, thereby contributing to the model’s calibration by simulation. We employ a novel composition of variance reduction methods, immediately applicable to any conditionally log-normal stochastic volatility model. Assuming one targets implied volatility estimates with a given degree of confidence, thus calibration RMSE, the results we demonstrate equate to significant runtime reductions—roughly 20 times on average, across different correlation regimes.     

The representation of American options prices under stochastic volatility and jump-diffusion dynamics

This paper considers the problem of pricing American options when the dynamics of the underlying are driven by both stochastic volatility following a square-root process as used by Heston [Rev. Financial Stud., 1993, 6, 327–343], and by a Poisson jump process as introduced by Merton [J. Financial Econ., 1976, 3, 125–144]. Probability arguments are invoked to find a representation of the solution in terms of expectations over the joint distribution of the underlying process. A combination of Fourier transform in the log stock price and Laplace transform in the volatility is then applied to find the transition probability density function of the underlying process. It turns out that the price is given by an integral dependent upon the early exercise surface, for which a corresponding integral equation is obtained. The solution generalizes in an intuitive way the structure of the solution to the corresponding European option pricing problem obtained by Scott [Math. Finance, 1997, 7(4), 413–426], but here in the case of a call option and constant interest rates.     

Optimal portfolio positioning within generalized Johnson distributions

Many empirical studies have shown that financial asset returns do not always exhibit Gaussian distributions, for example hedge fund returns. The introduction of the family of Johnson distributions allows a better fit to empirical financial data. Additionally, this class can be extended to a quite general family of distributions by considering all possible regular transformations of the standard Gaussian distribution. In this framework, we consider the portfolio optimal positioning problem, which has been first addressed by Brennan and Solanki [J. Financial Quant. Anal., 1981, 16, 279–300], Leland [J. Finance, 1980, 35, 581–594] and further developed by Carr and Madan [Quant. Finance, 2001, 1, 9–37] and Prigent [Generalized option based portfolio insurance. Working Paper, THEMA, University of Cergy-Pontoise, 2006]. As a by-product, we introduce the notion of Johnson stochastic processes. We determine and analyse the optimal portfolio for log return having Johnson distributions. The solution is characterized for arbitrary utility functions and illustrated in particular for a CRRA utility. Our findings show how the profiles of financial structured products must be selected when taking account of non Gaussian log-returns.     

On accurate and provably efficient GARCH option pricing algorithms

The GARCH model has been very successful in capturing the serial correlation of asset return volatilities. As a result, applying the model to options pricing attracts a lot of attention. However, previous tree-based GARCH option pricing algorithms suffer from exponential running time, a cut-off maturity, inaccuracy, or some combination thereof. Specifically, this paper proves that the popular trinomial-tree option pricing algorithms of Ritchken and Trevor (Ritchken, P. and Trevor, R., Pricing options under generalized GARCH and stochastic volatility processes. J. Finance, , 54(1), 377–402.) and Cakici and Topyan (Cakici, N. and Topyan, K., The GARCH option pricing model: a lattice approach. J. Comput. Finance, , 3(4), 71–85.) explode exponentially when the number of partitions per day, n, exceeds a threshold determined by the GARCH parameters. Furthermore, when explosion happens, the tree cannot grow beyond a certain maturity date, making it unable to price derivatives with a longer maturity. As a result, the algorithms must be limited to using small n, which may have accuracy problems. The paper presents an alternative trinomial-tree GARCH option pricing algorithm. This algorithm provably does not have the short-maturity problem. Furthermore, the tree-size growth is guaranteed to be quadratic if n is less than a threshold easily determined by the model parameters. This level of efficiency makes the proposed algorithm practical. The surprising finding for the first time places a tree-based GARCH option pricing algorithm in the same complexity class as binomial trees under the Black–Scholes model. Extensive numerical evaluation is conducted to confirm the analytical results and the numerical accuracy of the proposed algorithm. Of independent interest is a simple and efficient technique to calculate the transition probabilities of a multinomial tree using generating functions.     

Pricing and hedging guaranteed minimum withdrawal benefits under a general Lévy framework using the COS method

This paper extends the Fourier-cosine (COS) method to the pricing and hedging of variable annuities embedded with guaranteed minimum withdrawal benefit (GMWB) riders. The COS method facilitates efficient computation of prices and hedge ratios of the GMWB riders when the underlying fund dynamics evolve under the influence of the general class of Lévy processes. Formulae are derived to value the contract at each withdrawal date using a backward recursive dynamic programming algorithm. Numerical comparisons are performed with results presented in Bacinello et al. [Scand. Actuar. J., 2014, 1–20], and Luo and Shevchenko [Int. J. Financ. Eng., 2014, 2, 1–24], to confirm the accuracy of the method. The efficiency of the proposed method is assessed by making comparisons with the approach presented in Bacinello et al. [op. cit.]. We find that the COS method presents highly accurate results with notably fast computational times. The valuation framework forms the basis for GMWB hedging. A local risk minimisation approach to hedging intra-withdrawal date risks is developed. A variety of risk measures are considered for minimisation in the general Lévy framework. While the second moment and variance have been considered in existing literature, we show that the Value-at-Risk (VaR) may also be of interest as a risk measure to minimise risk in variable annuities portfolios.     

The Volatility and Price Sensitivities of Managerial Stock Option Portfolios and Corporate Hedging

We use estimates of the Black–Scholes sensitivity of managers' stock option portfolios to stock return volatility and the sensitivity of managers' stock and stock option portfolios to stock price to test the relationship between managers' risk preferences and hedging activities. We find that as the sensitivity of managers' stock and stock option portfolios to stock price increases, firms tend to hedge more. However, as the sensitivity of managers' stock option portfolios to stock return volatility increases, firms tend to hedge less.     

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.     

Forecasting Credit Spread Volatility: Evidence from the Japanese Eurobond Market

Different approaches to forecasting the volatility associated with the credit spreads on Yen Eurobonds are investigated. The actual volatility, historical volatility and estimated conditional volatility on spreads derived from a regression-based model with a GARCH and ARMA specification are compared within an adaptation of Black’s (J. Finance, 31, 1976, 361–367) option-pricing model. Surprisingly, the regression forecast over a medium forecasting horizon suggests that historic volatility provides the better forecast. The implications of these results for volatility forecasting and credit spread modelling are also discussed.JEL Classification: C32; G15     

Are tightened trading rules always bad? Evidence from the Chinese index futures market

This paper investigates the impact of tightened trading rules on the market efficiency and price discovery function of the Chinese stock index futures in 2015. The market efficiency and the price discovery of Chinese stock index futures do not deteriorate after these rule changes. Using variance ratio and spectral shape tests, we find that the Chinese index futures market becomes even more efficient after the tightened rules came into effect. Furthermore, by employing Schwarz and Szakmary [J. Futures Markets, 1994, 14(2), 147–167] and Hasbrouck [J. Finance, 1995, 50(4), 1175–1199] price discovery measures, we find that the price discovery function, to some extent, becomes better. This finding is consistent with Stein [J. Finance, 2009, 64(4), 1517–1548], who documents that regulations on leverage can be helpful in a bad market state, and Zhu [Rev. Financ. Stud., 2014, 27(3), 747–789.], who finds that price discovery can be improved with reduced liquidity. It also suggests that the new rules may effectively regulate the manipulation behaviour of the Chinese stock index futures market during a bad market state, and then positively affect its market efficiency and price discovery function.     

Time-varying long-run mean of commodity prices and the modeling of futures term structures

The exploration of the mean-reversion of commodity prices is important for inventory management, inflation forecasting and contingent claim pricing. Bessembinder et al. [J. Finance, 1995, 50, 361–375] document the mean-reversion of commodity spot prices using futures term structure data; however, mean-reversion to a constant level is rejected in nearly all studies using historical spot price time series. This indicates that the spot prices revert to a stochastic long-run mean. Recognizing this, I propose a reduced-form model with the stochastic long-run mean as a separate factor. This model fits the futures dynamics better than do classical models such as the Gibson–Schwartz [J. Finance, 1990, 45, 959–976] model and the Casassus–Collin-Dufresne [J. Finance, 2005, 60, 2283–2331] model with a constant interest rate. An application for option pricing is also presented in this paper.     

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.     

A comparison principle between rough and non-rough Heston models—with applications to the volatility surface

We present a number of related comparison results, which allow one to compare moment explosion times, moment generating functions and critical moments between rough and non-rough Heston models of stochastic volatility. All results are based on a comparison principle for certain non-linear Volterra integral equations. Our upper bound for the moment explosion time is different from the bound introduced by Gerhold, Gerstenecker and Pinter [Moment explosions in the rough Heston model. Decisions in Economics and Finance, 2019, 42, 575–608] and tighter for typical parameter values. The results can be directly transferred to a comparison principle for the asymptotic slope of implied variance between rough and non-rough Heston models. This principle shows that the ratio of implied variance slopes in the rough versus non-rough Heston model increases at least with power-law behavior for small maturities.     

Dynamic factor long memory volatility

In this paper, we develop a long memory orthogonal factor (LMOF) multivariate volatility model for forecasting the covariance matrix of financial asset returns. We evaluate the LMOF model using the volatility timing framework of Fleming et al. [J. Finance, 2001, 56, 329–352] and compare its performance with that of both a static investment strategy based on the unconditional covariance matrix and a range of dynamic investment strategies based on existing short memory and long memory multivariate conditional volatility models. We show that investors should be willing to pay to switch from the static strategy to a dynamic volatility timing strategy and that, among the dynamic strategies, the LMOF model consistently produces forecasts of the covariance matrix that are economically more useful than those produced by the other multivariate conditional volatility models, both short memory and long memory. Moreover, we show that combining long memory volatility with the factor structure yields better results than employing either long memory volatility or the factor structure alone. The factor structure also significantly reduces transaction costs, thus increasing the feasibility of dynamic volatility timing strategies in practice. Our results are robust to estimation error in expected returns, the choice of risk aversion coefficient, the estimation window length and sub-period analysis.     

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.     

An improved convolution algorithm for discretely sampled Asian options

We suggest an improved FFT pricing algorithm for discretely sampled Asian options with general independently distributed returns in the underlying. Our work complements the studies of Carverhill and Clewlow [Risk, 1990, 3(4), 25–29], Benhamou [J. Comput. Finance, 2002, 6(1), 49–68], and Fusai and Meucci [J. Bank. Finance, 2008, 32(10), 2076–2088], and, if we restrict our attention only to log-normally distributed returns, also Ve?e? [Risk, 2002, 15(6), 113–116]. While the existing convolution algorithms compute the density of the underlying state variable by moving forward on a suitably defined state space grid, our new algorithm uses backward price convolution, which resembles classical lattice pricing algorithms. For the first time in the literature we provide an analytical upper bound for the pricing error caused by the truncation of the state space grid and by the curtailment of the integration range. We highlight the benefits of the new scheme and benchmark its performance against existing finite difference, Monte Carlo, and forward density convolution algorithms.