American-style options in jump-diffusion models: estimation and evaluation

We propose dynamic programming coupled with finite elements for valuing American-style options under Gaussian and double exponential jumps à la Merton [J. Financ. Econ., 1976, 3, 125–144] and Kou [Manage. Sci., 2002, 48, 1086–1101], and we provide a proof of uniform convergence. Our numerical experiments confirm this convergence result and show the efficiency of the proposed methodology. We also address the estimation problem and report an empirical investigation based on Home Depot. Jump-diffusion models outperform their pure-diffusion counterparts.     

Risk-managed industry momentum and momentum crashes

This paper investigates Barroso and Santa-Clara’s [J. Financ. Econ., 2008, 116, 111–120] risk-managed momentum strategy in an industry momentum setting. We investigate several traditional momentum strategies including that recently proposed by Novy-Marx [J. Financ. Econ., 2012, 103, 429–453]. We moreover examine the impact of different variance forecast horizons on average pay-offs and also Daniel and Moskowitz’s [J. Financ. Econ., 2016, 122, 221–247] optionality effects. Our results show in general that neither plain industry momentum strategies nor the risk-managed industry momentum strategies are subject to optionality effects, implying that these strategies have no time-varying beta. Moreover, the benefits of risk management are robust across volatility estimators, momentum strategies and subsamples. Finally, the ‘echo effect’ in industries is not robust in subsamples as the strategy works only during the most recent subsample.     

Real-time market microstructure analysis: online transaction cost analysis

Motivated by the practical challenge in monitoring the performance of a large number of algorithmic trading orders, this paper provides a methodology that leads to automatic discovery of causes that lie behind poor trading performance. It also gives theoretical foundations to a generic framework for real-time trading analysis. The common acronym for investigating the causes of bad and good performance of trading is transaction cost analysis Rosenthal [Performance Metrics for Algorithmic Traders, 2009]). Automated algorithms take care of most of the traded flows on electronic markets (more than 70% in the US, 45% in Europe and 35% in Japan in 2012). Academic literature provides different ways to formalize these algorithms and show how optimal they can be from a mean-variance (like in Almgren and Chriss [J. Risk, 2000, 3(2), 5–39]), a stochastic control (e.g. Guéant et al. [Math. Financ. Econ., 2013, 7(4), 477–507]), an impulse control (see Bouchard et al. [SIAM J. Financ. Math., 2011, 2(1), 404–438]) or a statistical learning (as used in Laruelle et al. [Math. Financ. Econ., 2013, 7(3), 359–403]) viewpoint. This paper is agnostic about the way the algorithm has been built and provides a theoretical formalism to identify in real-time the market conditions that influenced its efficiency or inefficiency. For a given set of characteristics describing the market context, selected by a practitioner, we first show how a set of additional derived explanatory factors, called anomaly detectors, can be created for each market order (following for instance Cristianini and Shawe-Taylor [An Introduction to Support Vector Machines and Other Kernel-based Learning Methods, 2000]). We then will present an online methodology to quantify how this extended set of factors, at any given time, predicts (i.e. have influence, in the sense of predictive power or information defined in Basseville and Nikiforov [Detection of Abrupt Changes: Theory and Application, 1993], Shannon [Bell Syst. Tech. J., 1948, 27, 379–423] and Alkoot and Kittler [Pattern Recogn. Lett., 1999, 20(11), 1361–1369]) which of the orders are underperforming while calculating the predictive power of this explanatory factor set. Armed with this information, which we call influence analysis, we intend to empower the order monitoring user to take appropriate action on any affected orders by re-calibrating the trading algorithms working the order through new parameters, pausing their execution or taking over more direct trading control. Also we intend that use of this method can be taken advantage of to automatically adjust their trading action in the post trade analysis of algorithms.     

Weighing asset pricing factors: a least squares model averaging approach

Empirical evidence has demonstrated that certain factors in asset pricing models are more important than others for explaining specific portfolio returns. We propose a technique that evaluates the factors included in popular linear asset pricing models. Our method has the advantage of simultaneously ranking the relative importance of those pricing factors through comparing their model weights. As an empirical verification, we apply our method to portfolios formed following Fama and French [A five-factor asset pricing model. J. Financ. Econ., 2015, 116, 1–22] and demonstrate that models accommodated to our factor rankings do improve their explanatory power in both in-sample and out-of-sample analyses.     

The survival probability of the SABR model: asymptotics and application

The stochastic-alpha-beta-rho (SABR) model is widely used by practitioners in interest rate and foreign exchange markets. The probability of hitting zero sheds light on the arbitrage-free small strike implied volatility of the SABR model (see, e.g. De Marco et al. [SIAM J. Financ. Math., 2017, 8(1), 709–737], Gulisashvili [Int. J. Theor. Appl. Financ., 2015, 18, 1550013], Gulisashvili et al. [Mass at zero in the uncorrelated SABR modeland implied volatility asymptotics, 2016b]), and the survival probability is also closely related to binary knock-out options. Besides, the study of the survival probability is mathematically challenging. This paper provides novel asymptotic formulas for the survival probability of the SABR model as well as error estimates. The formulas give the probability that the forward price does not hit a nonnegative lower boundary before a fixed time horizon.     

A computational definition of financial randomness

The aim of this study is to present an efficient and easy framework for the application of the Least Squares Monte Carlo methodology to the pricing of gas or power facilities as detailed in Boogert and de Jong [J. Derivatives, 2008, 15, 81–91]. As mentioned in the seminal paper by Longstaff and Schwartz [Rev. Financ. Stud. 2001, 113–147], the convergence of the Least Squares Monte Carlo algorithm depends on the convergence of the optimization combined with the convergence of the pure Monte Carlo method. In the context of the energy facilities, the optimization is more complex and its convergence is of fundamental importance in particular for the computation of sensitivities and optimal dispatched quantities. To our knowledge, an extensive study of the convergence, and hence of the reliability of the algorithm, has not been performed yet, in our opinion this is because the apparent infeasibility and complexity uses a very high number of simulations. We present then an easy way to simulate random trajectories by means of diffusion bridges in contrast to Dutt and Welke [J. Derivatives, 2008, 15 (4), 29–47] that considers time-reversal Itô diffusions and subordinated processes. In particular, we show that in the case of Cox-Ingersoll-Ross and Heston models, the bridge approach has the advantage to produce exact simulations even for non-Gaussian processes, in contrast to the time-reversal approach. Our methodology permits performing a backward dynamic programming strategy based on a huge number of simulations without storing the whole simulated trajectory. Generally, in the valuation of energy facilities, one is also interested in the forward recursion. We then design backward and forward recursion algorithms such that one can produce the same random trajectories by the use of multiple independent random streams without storing data at intermediate time steps. Finally, we show the advantages of our methodology for the valuation of virtual hydro power plants and gas storages.     

A Note on Gaussian Estimation of the CKLS and CIR Models with Feedback Effects for Japan

In this note we extend the Gaussian estimation of two factor CKLS and CIR models recently considered in Nowman, K. B. (2001, Gaussian estimation and forecasting of multi-factor term structure models with an application to Japan and the United Kingdom, Asia Pacif. Financ. Markets 8, 23–34) to include feedback effects in the conditional mean as was originally formulated in general continuous time models by Bergstrom, A. R. (1966, Non-recursive models as discrete approximations to systems of stochastic differential equations, Econometrica 34, 173–182) with constant volatility. We use the exact discrete model of Bergstrom, A. R. (1966, Non-recursive models as discrete approximations to systems of stochastic differential equations, Econometrica 34, 173–182) to estimate the parameters which was first used by Brennan, M. J. and Schwartz, E. S. (1979, A continuous time approach to the pricing of bonds, J. Bank. Financ. 3, 133–155) to estimate their two factor interest model but incorporating the assumption of Nowman, K. B. (1997, Gaussian estimation of single-factor continuous time models of the term structure of interest rates, J. Financ. 52, 1695–1706; 2001, Gaussian estimation and forecasting of multi-factor term structure models with an application to Japan and the United Kingdom, Asia Pacif. Financ. Markets 8, 23–34). An application to monthly Japanese Euro currency rates indicates some evidence of feedback from the 1-year rate to the 1-month rate in both the CKLS and CIR models. We also find a low level-volatility effect supporting Nowman, K. B. (2001, Gaussian estimation and forecasting of multi-factor term structure models with an application to Japan and the United Kingdom, Asia Pacif. Financ. Markets 8, 23–34).     

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.     

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.     

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.     

Pricing via recursive quantization in stochastic volatility models

We provide the first recursive quantization-based approach for pricing options in the presence of stochastic volatility. This method can be applied to any model for which an Euler scheme is available for the underlying price process and it allows one to price vanillas, as well as exotics, thanks to the knowledge of the transition probabilities for the discretized stock process. We apply the methodology to some celebrated stochastic volatility models, including the Stein and Stein [Rev. Financ. Stud. 1991, (4), 727–752] model and the SABR model introduced in Hagan et al. [Wilmott Mag., 2002, 84–108]. A numerical exercise shows that the pricing of vanillas turns out to be accurate; in addition, when applied to some exotics like equity-volatility options, the quantization-based method overperforms by far the Monte Carlo simulation.     

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 American options written on two underlying assets

This paper extends the integral transform approach of McKean [Ind. Manage. Rev., 1965, 6, 32–39] and Chiarella and Ziogas [J. Econ. Dyn. Control, 2005, 29, 229–263] to the pricing of American options written on more than one underlying asset under the Black and Scholes [J. Polit. Econ., 1973, 81, 637–659] framework. A bivariate transition density function of the two underlying stochastic processes is derived by solving the associated backward Kolmogorov partial differential equation. Fourier transform techniques are used to transform the partial differential equation to a corresponding ordinary differential equation whose solution can be readily found by using the integrating factor method. An integral expression of the American option written on any two assets is then obtained by applying Duhamel’s principle. A numerical algorithm for calculating American spread call option prices is given as an example, with the corresponding early exercise boundaries approximated by linear functions. Numerical results are presented and comparisons made with other alternative approaches.     

A risk-sensitive stochastic control approach to an optimal investment problem with partial information

We consider an infinite time horizon optimal investment problem where an investor tries to maximize the probability of beating a given index. From a mathematical viewpoint, this is a large deviation probability control problem. As shown by Pham (in Syst. Control Lett. 49: 295–309, 2003; Financ. Stoch. 7: 169–195, 2003), its dual problem can be regarded as an ergodic risk-sensitive stochastic control problem. We discuss the partial information counterpart of Pham (in Syst. Control Lett. 49: 295–309, 2003; Financ. Stoch. 7: 169–195, 2003). The optimal strategy and the value function for the dual problem are constructed by using the solution of an algebraic Riccati equation. This equation is the limit equation of a time inhomogeneous Riccati equation derived from a finite time horizon problem with partial information. As a result, we obtain explicit representations of the value function and the optimal strategy for the problem. Furthermore we compare the optimal strategies and the value functions in both full and partial information cases.

---

Electronic Supplementary Material Supplementary material is available for this article at      

Learning minimum variance discrete hedging directly from the market

Option hedging is a critical risk management problem in finance. In the Black–Scholes model, it has been recognized that computing a hedging position from the sensitivity of the calibrated model option value function is inadequate in minimizing variance of the option hedge risk, as it fails to capture the model parameter dependence on the underlying price (see e.g. Coleman et al., J. Risk, 2001, 5(6), 63–89; Hull and White, J. Bank. Finance, 2017, 82, 180–190). In this paper, we demonstrate that this issue can exist generally when determining hedging position from the sensitivity of the option function, either calibrated from a parametric model from current option prices or estimated nonparametricaly from historical option prices. Consequently, the sensitivity of the estimated model option function typically does not minimize variance of the hedge risk, even instantaneously. We propose a data-driven approach to directly learn a hedging function from the market data by minimizing variance of the local hedge risk. Using the S&P 500 index daily option data for more than a decade ending in August 2015, we show that the proposed method outperforms the parametric minimum variance hedging method proposed in Hull and White [J. Bank. Finance, 2017, 82, 180–190], as well as minimum variance hedging corrective techniques based on stochastic volatility or local volatility models. Furthermore, we show that the proposed approach achieves significant gain over the implied BS delta hedging for weekly and monthly hedging.     

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.     

A super-replication theorem in Kabanov’s model of transaction costs

We prove a general version of the super-replication theorem, which applies to Kabanov’s model of foreign exchange markets under proportional transaction costs. The market is described by a matrix-valued càdlàg bid-ask process evolving in continuous time. We propose a new definition of admissible portfolio processes as predictable (not necessarily right- or left- continuous) processes of finite variation related to the bid-ask process by economically meaningful relations. Under the assumption of existence of a strictly consistent price system (SCPS), we prove a closedness property for the set of attainable vector-valued contingent claims. We then obtain the super-replication theorem as a consequence of that property, thus generalizing to possibly discontinuous bid-ask processes analogous results obtained by Kabanov (Financ. Stoch. 3, 237–248, 1999), Kabanov and Last (Math. Financ. 12, 63–70, 2002) and Kabanov and Stricker (Advances in Finance and Stochastics: Essays in Honour of Dieter Sondermann, pp 125–136, 2002). Rásonyi’s counter-example (Lecture Notes in Mathematics 1832, 394–398, 2003) served as an important motivation for our approach.     

Fast estimation of true bounds on Bermudan option prices under jump-diffusion processes

Fast pricing of American-style options has been a difficult problem since it was first introduced to the financial markets in 1970s, especially when the underlying stocks’ prices follow some jump-diffusion processes. In this paper, we extend the ‘true martingale algorithm’ proposed by Belomestny et al. [Math. Finance, 2009, 19, 53–71] for the pure-diffusion models to the jump-diffusion models, to fast compute true tight upper bounds on the Bermudan option price in a non-nested simulation manner. By exploiting the martingale representation theorem on the optimal dual martingale driven by jump-diffusion processes, we are able to explore the unique structure of the optimal dual martingale and construct an approximation that preserves the martingale property. The resulting upper bound estimator avoids the nested Monte Carlo simulation suffered by the original primal–dual algorithm, therefore significantly improving the computational efficiency. Theoretical analysis is provided to guarantee the quality of the martingale approximation. Numerical experiments are conducted to verify the efficiency of our algorithm.     

Maximizing survival,growth and goal reaching under borrowing constraints

In this paper, we consider the survival, growth and goal reaching maximization problems treated in Browne [Math. Oper. Res., 1997, 22, 468–493] and solve them in a market constrained due to borrowing prohibition. To solve the problems, we first construct an auxiliary market introduced by Cvitanic and Karatzas [Ann. Appl. Probab., 1992, 2, 767–818] and then apply the dynamic programming approach. Via our solutions, an alternative approach is introduced in order to solve the problems defined under an auxiliary market.