Interfaces with Other Disciplines
A data-driven framework for consistent financial valuation and risk measurement

https://doi.org/10.1016/j.ejor.2020.07.011Get rights and content

Highlights

  • A consistent framework for option pricing and risk-measurment.

  • Applicable to case of limited option quotes.

  • General approach, requires only an observable underlying time series.

  • Fast and accurate computational procedures for practical application.

Abstract

In this paper, we propose a general data-driven framework that unifies the valuation and risk measurement of financial derivatives, which is especially useful in markets with thinly-traded derivatives. We first extract the empirical characteristic function from market-observable time series for the underlying asset prices, and then utilize Fourier techniques to obtain the physical nonparametric density and cumulative distribution function for the log-returns process, based on which we compute risk measures. Then we risk-neutralize the nonparametric density and distribution functions to model-independently valuate a variety of financial derivatives, including path-independent European options and path-dependent exotic contracts. By estimating the state-price density explicitly, and utilizing a convenient basis representation, we are able to greatly simplify the pricing of exotic options all within a consistent model-free framework. Numerical examples, and an empirical example using real market data (Brent crude oil prices) illustrate the accuracy and versatility of the proposed method in handling pricing and risk management of multiple financial contracts based solely on observable time series data.

Introduction

The problem of decision making is central to the field of operations research, where problems with embedded “optionality” play a central role. Examples include capital budgeting of R&D investment projects (Pennings & Lint, 1997), adoption of a new technology (Bobtcheff & Villeneuve, 2010), and the investment project value creation (Borgonovo, Gatti, & Peccati, 2010). To date, the operations research literature studying derivatives markets has focused mainly on parametric options pricing. With the well-accepted belief that “all models are wrong”, there have been investigations on the potential model risk in financial valuation (see Coqueret & Tavin, 2016), as well as a growing interest in the development of nonparametric/model-independent or data-driven techniques. This paper contributes a consistent data-driven framework for the valuation of derivative securities and their risk measurement.

With the rapid development of high frequency trading, the financial market has generated volumes of trades and price quote data suitable for and in need of analysis. A practical question that confronts investors, traders or fund managers is how to fairly determine the price and assess the risk of financial assets and their associated derivatives given the historical times series data of the asset, while maintaining a minimal use of parametric model assumptions. This motivates us to explore the possibility of designing a generic and consistent financial valuation and risk management framework that is fully data-driven, and is well-suited for the case of limited trading in derivatives. The goal of this research is to create a bridge that allows data to be directly translated into useful information that can be utilized for both valuation and risk management purposes. What we propose in this paper is to directly let the data meet the valuation problem without resorting to parametric models, i.e. “let the data directly speak”.

Previously, there have been several successful and popular parametric and nonparametric methods developed and employed for the valuation and risk management of financial derivatives. On the nonparametric side, there are significant developments in nonparametric density estimation through kernel density estimators (see Botev, Grotowski, Kroese, 2010, Figueroa-López, Li, 2020). Although most of the literatures lie in the realm of statistics, the methods proposed can be applied equally well to financial applications. For example, nonparametric density estimators have been applied to estimate transition densities of asset prices (Aït-Sahalia & Park, 2016), the state price densities (Song, Xiu, 2016, Zhang, Brooks, King, 2009), and risk measures such as value-at-risk and expected shortfall (Martins-Filho, Yao, Torero, 2018, Patton, Ziegel, Chen, 2019, Wang, Zhao, 2016). For a review of some of the earlier nonparametric literature pertaining to financial applications, please refer to Fan (2005) and references therein.

Another strand of nonparametric literature focuses on recovering risk-neutral densities from options prices data. Most literature stems from the celebrated (Breeden & Litzenberger, 1978) result, which states that the risk-neutral density is essentially the second derivative of the European call option price with respect to the strike, up to a risk-free discount factor, see for example (Jackwerth & Rubinstein, 1996). It is worthwhile to point out an alternative approach proposed in Rompolis and Tzavalis (2008), which first invokes the Gram-Charlier expansions, and then utilizes the Bakshi and Madan (2000) spanning result. The above referenced work (and related works) assumes that the market already prices a broad spectrum of vanilla options. In contrast, we focus on the other common case when the underlying asset is observable, but the corresponding options market is not available or not mature. An example is the cryptocurrency market, e.g. the Bitcoin market, which still awaits the introduction of options written on Bitcoins. This is also the case when we face with a thinly-traded derivative or an over-the-counter (OTC) derivative whose trading data is obscure or hard to obtain. Examples in this category arise in the diverse commodity markets for energy, agriculture, and metals (see Li & Linetsky, 2014), where the majority of contracts have little or no daily trading volume.1 Hence we are unable to recover risk-neutral densities from option prices, but are instead given properties of the underlying asset process alone. The research question we aim to address is how to develop a consistent and arbitrage-free option valuation framework given only data of the underlying asset.

Naturally there are pros and cons associated with various parametric and nonparametric methods. Parametric models are widely used in industry for option pricing, since they are convenient and possess desirable statistical properties when the underlying model’s assumptions hold. In the literature on parametric models for options pricing, popular stochastic models include stochastic volatility models (see Bernard, Cui, & McLeish, 2017) and exponential Lévy models (see Lian, Zhu, Elliott, & Cui, 2017). Various financial contracts with path-dependent payoffs are considered, such as the arithmetic Asian options (see Cui, Lee, & Liu, 2018b). However, misspecification of the true underlying model results in inconsistent estimators. Thus it is preferable to develop an alternative approach that minimizes the impacts of choices of models. There is a strand of literature that incorporates model uncertainty in the valuation, however, it is still not fully nonparametric, or data-driven. This research is in line with recent efforts in a data-driven approach to simulation-based model calibration (Peng, Fu, Heidergott, & Lam, 2019), which aim to reduce or eliminate model-misspecification errors. Instead of carrying out inference of input parameters, our focus here is on directly extracting information from the output-level data (i.e. historical time series of asset prices) to estimate prices of financial contracts and compute risk measures within a unified framework.

Our approach is nonparametric in nature, but it inherits some nice properties akin to parametric methods in that no-arbitrage constraints are satisfied by design. By utilizing the continuous nonparametric density estimator designed in this paper, and pricing through integration with respect to this density, the smoothness of the pricing function is guaranteed similar to that with a parametric approach. Thus our approach is fully data-driven while at the same time it also preserves the smoothness of the pricing function. This distinguishes our method from other nonparametric methods. There exists literature that combines parametric and nonparametric features2. A paper closely related is Barone-Adesi et al. (2008), in which they propose a parametric asymmetric GARCH model, while they estimate the innovation distribution in a nonparametric way. There are at least two differences from our approach. First, their basis is a GARCH model assumption, while we assume an underlying Lévy model. Second, they use Monte Carlo simulation to price derivatives, while our method allows for analytical Fourier inversion valuations. A detailed comparison between the two approaches is left to future research.

Besides valuation and pricing, accurate and prompt risk measurements are also important in everyday financial practice. For example, as required by the Dodd-Frank act3 or Basel III4, financial institutions have to report their financial positions and control the risk exposures to illiquid assets. Similarly the solvency requirements (e.g. Solvency II Directive5) in Europe imposes that insurance and reinsurance companies have to hold certain amount of capital (aka solvency capital requirements (SCR)) for risk planning. The risk measures serve this reporting purpose, and two of the most popular ones are the value at risk (VaR), and the expected shortfall (ES), which is also named as the conditional value at risk (CVaR). There have been some recent literatures on the nonparametric estimation of risk measures (see Chen, 2007) and also portfolio optimization problems with nonparametric risk measures as constraints (see Cui, Sun, Zhu, Jiang, & Li, 2018a). In a recent paper (Cai & Wang, 2008), they consider the application of the double kernel local linear estimator and the weighted Nadaraya-Watson estimator, and construct a distribution estimator that is continuous and differentiable. In contrast, the method proposed in this paper allows us to obtain risk estimates in a consistent manner within the same framework as valuation.

The approach taken in this paper centers around the interplay between observable financial time series under the physical measure, and the need for risk management under the physical measure and valuation of financial derivatives under the risk-neutral measure. Instead of formulating a parametric model such as a diffusion model, or adopting the risk-neutralization of a well-known times series model (such as the Heston & Nandi, 2000 GARCH model), we propose a new alternative approach, which is model-independent, and fully data-driven. The novelty of our method lies in that we carry out the risk-neutralization directly at the empirical distribution function level. Thus we are able to utilize time series to valuate different types of options contracts under the risk-neutral measure, including both path-independent and path-dependent options.

The contributions of the paper are three-fold:

  • 1.

    We propose a fully data-driven approach for financial valuation and risk management, which is based on the empirical characteristic function and the novel use of Fourier techniques. The input is historical time series of the financial asset, while the output are risk measures, and option prices for various types of payoffs, ranging from vanilla European call options to path-dependent Bermudan and Asian options.

  • 2.

    Our proposed approach is accurate and efficient, when certain conditions are met for the data generating process (DGP). We test the method using both known parametric models and real market data DGP, and demonstrate that our method yields consistently accurate computation of risk measures and valuation of a variety of financial contracts. We also carry out backtests of our method, which yield highly-accurate results.

  • 3.

    We propose a general data-driven framework which enables consistent no arbitrage option pricing and risk calculation without sacrificing accuracy. In general, market practitioners use the risk-neutral measure to price, and calculate VaR and ES under the physical measure. Our proposed approach uses a consistent estimator for both, and can be treated as a unified calibration method for financial valuation and risk management, which is especially useful for opaque OTC markets or markets with thinly-traded derivatives.

The rest of the paper is organized as follows: Section 2 describes the main methodology behind the nonparametric option pricing, where we propose in Section 2.1.1 the density estimator based on B-splines. Section 3 discusses the model-free option pricing method base on the Fourier method, i.e. the projection (PROJ) method. in particular, two risk-neutralization methods are discussed and compared with each other. Section 4 considers the application to risk measurement. More specifically, we consider the Value at Risk (VaR) in Section 4.1 and the Expected Shortfall (ES) in Section 4.2. Section 5 presents the numerical experiments and an empirical example, which illustrate the accuracy and efficiency of our proposed data-driven approach. We also carry out backtests in this section. Section 6 concludes the paper with discussions of potential future research directions.

Section snippets

Nonparametric Estimation

Assume that the market is equipped with a complete probability space (Ω,F,Ft,P), where P is the physical or real-world probability measure. In the following discussions, we assume that there exists a risk-neutral probability measure P*, which will serve as the pricing measure, and that there exists a constant rate of interest rR. Denote the risk-neutral probability density function of the asset as p*( · ). The price of a T-maturity European call option with terminal payoff (STK)+:=max{STK,0}

Model-Free Option Pricing

In the field of option pricing, parametric models are abundant in mature and liquid markets, and model parameters are calibrated (tuned) by minimizing the discrepancy between model predictions and market-quoted option prices. Alternatively, nonparametric and semi-parametric methods have been proposed to mitigate the dependence on model assumptions and alleviate the potential for “model misspecification error” (see Aït-Sahalia, 1996, Gagliardini, Gourieroux, Renault, 2011, Gagliardini,

Risk Measurement

While the risk-neutral measure is used to price financial derivatives, the physical (market) measure is generally used when simulating losses and calculating risk measures. In this section, we provide formulas which allow for the calculation of the risky asset’s risk measures based on the pricing measure, up to a shift which accounts for the risk neutrality during pricing. That is, we use f˜Da(x) to measure the risk, and f˘Da(x)=f˜Da(xθ) to price derivatives, where the shift θ ensures the

Numerical Experiments

In this section, we carry out several numerical experiments demonstrating the accuracy of the proposed general nonparametric approach in applications to vanilla and exotic options pricing and the estimation of financial risk measures.12 Unless otherwise specified, prices are computed using the mean-correcting EMM of Section 3.1.2. In Section 5.4, we make a comparison to the Esscher

Conclusions and Future Research

In this paper, we have introduced a new consistent approach to carry out both risk-neutral valuation (e.g. options valuation) and physical or market risk measurement (e.g. VaR and expected shortfall). The approach is particularly useful for markets where the corresponding derivative market is thinly-traded or opaque. Numerical and empirical examples demonstrate the effectiveness of the proposed approach, including simulated data and historical backtests. Potential future research directions are

References (93)

  • Z. Cui et al.

    A general framework for discretely sampled realized variance derivatives in stochastic volatility models with jumps

    European J. Operational Research

    (2017)
  • Z. Cui et al.

    A general framework for time-changed Markov processes and applications

    European J. Operational Research

    (2019)
  • Z. Cui et al.

    Single-transform formulas for pricing Asian options in a general approximation framework under Markov processes

    European Journal of Operational Research

    (2018)
  • G. Fusai et al.

    Pricing discretely monitored Asian options under Lévy processes

    J. of Banking & Finance

    (2008)
  • P. Gagliardini et al.

    Semi-parametric estimation of American option prices

    J. of Econometrics

    (2013)
  • J. Kirkby et al.

    A unified approach to Bermudan and barrier options under stochastic volatility models with jumps

    J. of Economic Dynamics and Control

    (2017)
  • A. Leitao et al.

    On the data-driven COS method

    Applied Math. and Computation

    (2018)
  • G. Lian et al.

    Semi-analytical valuation for discrete barrier options under time-dependent Lévy processes

    Journal of Banking & Finance

    (2017)
  • R. Merton

    Option pricing when underlying stock returns are discontinuous

    J. of Financial Econ.

    (1976)
  • L. Ortiz-Gracia et al.

    Efficient VaR and expected shortfall computations for nonlinear portfolios within the delta-gamma approach

    Applied Math. and Computation

    (2014)
  • A.J. Patton et al.

    Dynamic semiparametric models for expected shortfall (and value-at-risk)

    J. of Econometrics

    (2019)
  • E. Pennings et al.

    The option value of advanced R & D

    European J. Operational Research

    (1997)
  • C. Phelan et al.

    Fluctuation identities with continuous monitoring and their application to the pricing of barrier options

    European J. Operational Research

    (2018)
  • Z. Song et al.

    A tale of two option markets: Pricing kernels and volatility risk

    J. of Econometrics

    (2016)
  • C.-S. Wang et al.

    Conditional Value-at-Risk: Semiparametric estimation and inference

    J. of Econometrics

    (2016)
  • D. Xiu

    Hermite polynomial based expansion of European option prices

    J. of Econometrics

    (2014)
  • X. Zhang et al.

    A Bayesian approach to bandwidth selection for multivariate kernel regression with an application to state-price density estimation

    J. of Econometrics

    (2009)
  • Z. Zhao et al.

    Nonparametric inference of discretely sampled stable Lévy processes

    J. of Econometrics

    (2009)
  • Y. Aït-Sahalia

    Nonparametric pricing of interest rate derivative securities

    Econometrica

    (1996)
  • Y. Aït-Sahalia

    Maximum likelihood estimation of discretely sampled diffusions: A closed-form approximation approach

    Econometrica

    (2002)
  • Y. Aït-Sahalia et al.

    Nonparametric estimation of state-price densities implicit in financial asset prices

    The J. of Finance

    (1998)
  • D. Applebaum

    Lévy Processes and Stochastic Calculus

    (2009)
  • L. Ballotta et al.

    Estimation of multivariate asset models with jumps

    J. of Financial and Quantitative Analysis

    (2019)
  • R.W. Banz et al.

    Prices for state-contingent claims: Some estimates and applications

    J. of Business

    (1978)
  • G. Barone-Adesi et al.

    A GARCH option pricing model with filtered historical simulation

    Review of Financial Studies

    (2008)
  • J. Berkowitz

    Testing density forecasts with applications to risk management

    J. Business and Econ. Stat.

    (2001)
  • C. Bernard et al.

    On the martingale property in stochastic volatility models based on time-homogeneous diffusions

    Mathematical Finance

    (2017)
  • Z.I. Botev et al.

    Kernel density estimation via diffusion

    Annals of Statistics

    (2010)
  • S. Boyarchenko et al.

    Non-Gaussian Merton-Black-Scholes Theory

    (2002)
  • D.T. Breeden et al.

    Prices of state-contingent claims implicit in option prices

    J. of Business

    (1978)
  • S. Campbell

    A review of backtesting and backtesting procedures

    J. Risk

    (2007)
  • A. Černỳ et al.

    An improved convolution algorithm for discretely sampled Asian options

    Quantitative Finance

    (2011)
  • R. Chatterjee et al.

    An efficient and stable method for short maturity Asian options

    J. of Futures Markets

    (2018)
  • S.X. Chen

    Nonparametric estimation of expected shortfall

    J. of Financial Econometrics

    (2007)
  • O. Christensen

    An Introduction to Frames and Riesz Bases

    (2003)
  • P. Christoffersen

    Evaluating interval forecasts

    Int. Econ. Review

    (1998)
  • Cited by (21)

    • Nonparametric density estimation with nonuniform B-spline bases

      2024, Journal of Computational and Applied Mathematics
    • A review of the operations literature on real options in energy

      2023, European Journal of Operational Research
      Citation Excerpt :

      Nonetheless, valuation and hedging of risky energy cash flows will continue to face the challenges posed by market incompleteness (Swindle, 2016), thus fueling the growth of practically relevant models and methodology. On one side, model error considerations will push modelers toward data driven approaches (see, e.g., Bertsimas & Dunn, 2019, Cui, Kirkby, & Nguyen, 2021). On the other side, past data may become increasingly irrelevant to a future that may look very different from the past, e.g., because sudden weather related changes will be much more frequent than they have been observed so far.

    • Smiles & smirks: Volatility and leverage by jumps

      2022, European Journal of Operational Research
      Citation Excerpt :

      Further research could also be carried out to develop expansions of the shape characteristics of the implied volatility surface generated by the models proposed in this paper on the line of Aït-Sahalia, Li, & Li (2020) and Aït-Sahalia et al. (2021). An equally interesting topic would be the extension of the model-free approach to calibration of Cui, Kirkby, & Nguyen (2021) to allow for time changes as advocated in this paper. Other relevant directions would be the analysis of the performance of the general framework proposed in this paper against VIX derivatives and other variance products in the spirit of Fouque & Saporito (2018) and Kaeck & Seeger (2020), and in terms of joint calibration of the smile of S&P options and VIX options (see Guyon, 2020, for further details).

    • Nonparametric density estimation and bandwidth selection with B-spline bases: A novel Galerkin method

      2021, Computational Statistics and Data Analysis
      Citation Excerpt :

      Related spline-based techniques have also been studied, such as logsplines (Kooperberg and Stone, 1991, 1992, 2004; Koo, 1996), smoothing splines (Gu and Qiu, 1993; Gu, 1993), penalized B-splines (P-splines) (Eilers and Marx, 1996), and shape-constrained splines (Papp and Alizadeh, 2014). Local basis expansion methods are especially useful in contexts where intensive numerical computations are conducted after estimation, such as in risk-management (Leitao and Ortiz-Gracia, 2020; Cui et al., 2021a,b) or forward problems in uncertainty quantification (Colombo et al., 2018; Ditkowski et al., 2020; Rahman, 2020). For example, in Cui et al. (2021a), the authors derive closed-form expressions for various portfolio risk measures, such as value at risk (VaR) and expected shortfall, due to the tractability1 of the B-spline basis.

    • Pricing equity-linked death benefits by complex Fourier series expansion in a regime-switching jump diffusion model

      2021, Applied Mathematics and Computation
      Citation Excerpt :

      To get around this problem, Ruijter et al. [22] propose spectral filters to accelerate the convergence rate. See also Cui et al. [14,23,24] and Shi and Zhang [25]. Accordingly, the parameters are given in Table 3, which are also considered in Cui et al. [26].

    View all citing articles on Scopus
    View full text