Next Article in Journal
News-Driven Expectations and Volatility Clustering
Next Article in Special Issue
Editorial Statement for Mathematical Finance
Previous Article in Journal
Acknowledgement to Reviewers of Journal of Risk and Financial Management in 2019
Previous Article in Special Issue
CVaR Regression Based on the Relation between CVaR and Mixed-Quantile Quadrangles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Equity Option Pricing with Systematic and Idiosyncratic Volatility and Jump Risks

Business School, Nanjing Normal University, Nanjing 210023, China
J. Risk Financial Manag. 2020, 13(1), 16; https://doi.org/10.3390/jrfm13010016
Submission received: 21 November 2019 / Revised: 10 January 2020 / Accepted: 12 January 2020 / Published: 17 January 2020
(This article belongs to the Special Issue Mathematical Finance with Applications)

Abstract

:
Recently, a large number of empirical studies indicated that individual equity options exhibit a strong factor structure. In this paper, the importance of systematic and idiosyncratic volatility and jump risks on individual equity option pricing is analyzed. First, we propose a new factor structure model for pricing the individual equity options with stochastic volatility and jumps, which takes into account four types of risks, i.e., the systematic diffusion, the idiosyncratic diffusion, the systematic jump, and the idiosyncratic jump. Second, we derive the closed-form solutions for the prices of both the market index and individual equity options by utilizing the Fourier inversion. Finally, empirical studies are carried out to show the superiority of our model based on the S&P 500 index and the stock of Apple Inc. on options. The out-of-sample pricing performance of our proposed model outperforms the other three benchmark models especially for short term and deep out-of-the-money options.
JEL Classification:
G13

1. Introduction

Most of the existing literature studies on option pricing are for index options, and there are very few about equity options. One approach to modeling equity options is to employ the state-of-the-art model in the index option literature, a stochastic volatility model with jumps (see, for example, Bates 1996, 2000; Bakshi et al. 1997; Duffie et al. 2000; Eraker et al. 2003; Broadie et al. 2007; Christoffersen et al. 2012; Andersen et al. 2015; Bardgett et al. 2019), but to ignore any underlying factor structure.
In Bakshi and Kapadia (2003a), the research results indicated that the volatility risk premium is negative in index options by examining the statistical properties of delta hedged option portfolios, i.e., a portfolio of a long call option position hedged by a short position in the stock. On the one hand, stock returns have a significant market component; the emergence of market volatility risk premiums is bound to have an impact on individual equity option pricing. On the other hand, from the economic point of view, the risk neutral distributions of individual equities are systematically different from the market index. Thus, it is necessary to explore how volatility risk is priced in individual equity options, which also can produce additional insights into the pricing structure of individual equity options (see Bakshi et al. 2003). As is well known, the beta of a stock represents the sensitivity of the risk of the individual equity with respect to the systematic risk of the market and is very useful for portfolio construction in the capital asset pricing model. Therefore, under the assumption that stock returns include a market component and an idiosyncratic component, Bakshi and Kapadia (2003b) developed a factor model for equity option valuation and investigated the pricing of market volatility risk in individual equity options. Their empirical results showed that volatility risk premiums in equity options are smaller than in index options.
Afterwards, Fouque and Kollman (2011) proposed a continuous-time capital asset pricing model (CAPM) where the dynamics of the market index have a stochastic volatility driven by a fast mean reverting process. Moreover, they derived the analytical approximation pricing formulas for both the market index and individual equity call options using a singular perturbation method. Meanwhile, a calibration method for the beta parameter was also presented based on the estimated model parameters of both the market index and individual equity option prices. Subsequently, Fouque and Tashman (2012) extended the constant beta-parameter factor model of Fouque and Kollman (2011) by considering a piecewise-linear relationship between the individual asset and the market index and proposed a regime switching factor model for the pricing market index and individual equity options. Supposing that stock return is linearly related to market index return in terms of the beta parameter, Carr and Madan (2012) developed a factor model for individual equity option pricing under a purely discontinuous Lévy process via fast Fourier transform, in which the variance gamma process for the dynamics of both the market index and stock was taken as an example for illustration. By supposing a continuous-time CAPM with Lévy processes, Wong et al. (2012) also derived analytical solutions to the index and equity options and explored the corresponding static hedging with index futures. Christoffersen et al. (2018) empirically studied the equity volatility levels, skews, and term structures by using equity option prices and principal component analysis. The results indicated that the equity options had a strong factor structure, and then, they developed an equity option pricing model with a CAPM factor structure and stochastic volatility, which allowed for mean reverting stochastic volatility for the dynamics of both the market factor and individual equity.
Recently, Xiao and Zhou (2018) proposed a GARCH-jump model for individual stock returns that took into account four types of risks: the systematic and idiosyncratic jumps and the systematic and idiosyncratic diffusive volatility. By using a dataset consisting of the S&P 500 index and 15 individual stock prices, their empirical results indicated that idiosyncratic jumps were a key determinant of expected stock return.1 Instead of using only stock returns, Kapadia and Zekhnini (2019) used both stock and option data to decompose the four risk premiums associated with systematic and idiosyncratic diffusive and jump risks and also documented that idiosyncratic jumps are important determinants of the mean returns of a stock from both an ex post and ex ante perspective.
Motivated by the above mentioned insights, we propose to price individual equity options in stochastic volatility jump-diffusion models with a market factor structure, which can be seen as a generalized version of Christoffersen et al. (2018). Specifically, in our proposed model, the individual equity prices are driven by the market factor, as well as an idiosyncratic component that also has stochastic volatility and jump. Due to the model belonging to the affine class, we derive the closed-form solutions for the prices of both the market index and individual equity options by utilizing the Fourier inversion. In addition, we provide the empirical results to test the pricing performance of the proposed factor model based on the S&P 500 index and the stock of Apple Inc. (AAPL) on options. Toward this end, we empirically compare the pricing performance of the proposed model with those of the other three classical two factor stochastic volatility models being taken as benchmark models. Empirical results presented here confirm that the equity option pricing model considering systematic and idiosyncratic volatility and jump risks may offer a good competitor to the models of Bates (2000), Christoffersen et al. (2009), or Christoffersen et al. (2018) for some other option markets.
The remainder of the paper proceeds as follows. In Section 2, we present a novel factor model for equity option valuation and derive the corresponding closed-form solutions. In Section 3, empirical studies are carried out to show the pricing performance of our proposed model. Finally, some conclusions are stated in Section 4.

2. Equity Option Valuation

In this section, we introduce a general class of stochastic volatility models with jumps for the dynamics of both the market factor and individual equity prices and derive closed form solutions to the prices of the European equity call and put options.

2.1. Model Description

Consider a filtered probability space ( Ω , F , Q ) with information filtration { F t } 0 t T satisfying the usual conditions (increase, right-continuous, and augmented), where Q is a risk neutral measure. We model an equity market consisting of N firms with a single market factor, I t (usually approximated by a market index in practice). The individual stock prices are denoted by S t i , for i = 1 , 2 , , N . For the sake of convenience, we ignore the superscript i, and denote ( S t ) t 0 the pricing process of an individual stock. Investors also have access to a risk free bond that pays a return rate of r. To start, the market factor I t evolves under a risk neutral measure Q as:
d I t I t = r d t + V I , t d W 1 , t I + R ( e y 1 ) N ˜ y ( d t , d y ) ,
d V I , t = κ I ( θ I V I , t ) d t + σ I V I , t d W 2 , t I ,
where I t stands for the value of I t before a possible jump occurs, y R = R { 0 } , V I , t is the variance of market factor, θ I denotes the long run variance, κ I captures the mean reversion speed of V I , t to θ I , σ I measures the volatility of volatility, 2 κ I θ I σ I 2 to ensure that the process V I , t remains strictly positive2, W 1 , t I and W 2 , t I are correlated standard Brownian motions, i.e., the innovations to the market return and volatility are correlated with correlation coefficient ρ I , Cov d W 1 , t I , d W 2 , t I = ρ I d t , and N ˜ y ( d t , d y ) = N y ( d t , d y ) ν y ( d y ) d t is a compensated jump measure, where N y ( d t , d y ) is the jump measure and the Lévy kernel (or density) ν y ( d y ) satisfies R min ( 1 , y 2 ) ν y ( d y ) < .
Furthermore, we separate the effects of the market factor on individual equities’ returns into two types of risks: the systematic diffusive volatility and jump. More specifically, the diffusive random variation of individual equities’ returns is dependent on the Brownian motion that drives market returns through the coefficient β d i f f . In addition, the discontinuous movements in the market return can also trigger jumps in individual equities’ returns through the coefficient β j u m p . Therefore, the individual equity prices are driven by the market factor, as well as an idiosyncratic component that also has stochastic volatility and jump, whose process under a risk neutral measure Q follows:3
d S t S t = r d t + β d i f f V I , t d W 1 , t I Systematic   diffusive + R ( e β j u m p y 1 ) N ˜ y ( d t , d y ) Systematic   jump + V S , t d W 1 , t S Idiosyncratic   diffusive + R ( e ξ 1 ) N ˜ ξ ( d t , d ξ ) Idiosyncratic   jump ,
d V S , t = κ S ( θ S V S , t ) d t + σ S V S , t d W 2 , t S ,
where S t stands for the value of S t before a possible jump occurs, ξ R = R { 0 } , V S , t is the idiosyncratic variance of individual equity, θ S denotes the long run idiosyncratic variance, κ S captures the mean reversion speed of V S , t to θ S , σ S measures the volatility of idiosyncratic variance, 2 κ S θ S σ S 2 to ensure that the process V S , t remains strictly positive4, W 1 , t S and W 2 , t S are correlated standard Brownian motions, i.e., the innovations to the idiosyncratic return and volatility are correlated with correlation coefficient ρ S , Cov d W 1 , t S , d W 2 , t S = ρ S d t , but they are independent of Brownian motions in the market factor, i.e., Cov d W i , t S , d W j , t I = 0 for i , j = 1 , 2 , and N ˜ ξ ( d t , d y ) = N ξ ( d t , d ξ ) ν ξ ( d ξ ) d t is a compensated jump measure, where N ξ ( d t , d ξ ) is the jump measure and the Lévy kernel (or density) ν ξ ( d ξ ) satisfies R min ( 1 , ξ 2 ) ν ξ ( d ξ ) < .

2.2. Characteristic Function

In order to be able to derive the pricing formulas for the European call and put equity options, we are particularly interested in the characteristic function of the logarithm asset price. Given the dynamics of the underlying asset price under the Q measure, we consider the conditional characteristic function of log-asset price X T = ln S T given the market information up to time t, which is denoted by φ ( x , υ 1 , υ 2 , t , T ; ϕ ) :
φ ( x , υ 1 , υ 2 , t , T ; ϕ ) = E Q e i ϕ X T | X t = x , V I , t = υ 1 , V S , t = υ 2 E t Q e i ϕ X T ,
where E t Q [ · ] denotes the condition expectation under the Q measure, t T , and i = 1 .
Lemma 1.
Suppose that the market factor I t and individual equity price S t are driven by Equations (1) and (3), respectively. Then, the conditional characteristic function of log-asset price X T = ln S T is given by:
φ ( x , υ 1 , υ 2 , t , T ; ϕ ) = exp A ( τ ) x + B ( τ ) υ 1 + C ( τ ) υ 2 + D ( τ ) ,
where:
A ( τ ) = i ϕ , B ( τ ) = κ I i ϕ β d i f f σ I ρ I d 1 σ I 2 1 e d 1 τ 1 g 1 e d 1 τ , C ( τ ) = κ S i ϕ σ S ρ S d 2 σ S 2 1 e d 2 τ 1 g 2 e d 2 τ , D ( τ ) = i ϕ r + R e i ϕ β j u m p y 1 i ϕ e β j u m p y 1 ν y ( d y ) I 1 + R e i ϕ ξ 1 i ϕ e ξ 1 ν ξ ( d ξ ) I 2 τ + κ I θ I σ I 2 κ I i ϕ β d i f f σ I ρ I d 1 τ 2 ln 1 g 1 e d 1 τ 1 g 1 + κ S θ S σ S 2 κ S i ϕ σ S ρ S d 2 τ 2 ln 1 g 2 e d 2 τ 1 g 2 , g 1 = κ I i ϕ β d i f f σ I ρ I d 1 κ I i ϕ β d i f f σ I ρ I + d 1 , g 2 = κ S i ϕ σ S ρ S d 2 κ S i ϕ σ S ρ S + d 2 , d 1 = i ϕ β d i f f σ I ρ I κ I 2 + β d i f f 2 σ I 2 ( i ϕ + ϕ 2 ) , d 2 = i ϕ σ S ρ S κ S 2 + σ S 2 ( i ϕ + ϕ 2 ) ,
and τ = T t .
Proof. 
To obtain the conditional characteristic function of log-asset price X T = ln S T , we first take the following transformation by using the Itô lemma for Equation (3):
d ln S t = r 1 2 β d i f f 2 V I , t 1 2 V S , t R e β j u m p y 1 ν y ( d y ) R e ξ 1 ν ξ ( d ξ ) d t + β d i f f V I , t d W 1 , t I + β j u m p R y N y ( d t , d y ) + V S , t d W 1 , t S + R ξ N ξ ( d t , d ξ ) .
The Feynman–Kac formula states that φ ( x , υ 1 , υ 2 , t , T ; ϕ ) is governed by the following partial integro-differential equation (PIDE):
φ τ = r 1 2 β d i f f 2 V I , t 1 2 V S , t R e β j u m p y 1 ν y ( d y ) R e ξ 1 ν ξ ( d ξ ) φ x + 1 2 β d i f f 2 υ 1 + υ 2 2 φ x 2 + κ I ( θ I υ 1 ) φ υ 1 + 1 2 σ I 2 υ 1 2 φ υ 1 2 + κ S ( θ S υ 2 ) φ υ 2 + 1 2 σ S 2 υ 2 2 φ υ 2 2 + β d i f f σ I ρ I υ 1 2 φ x υ 1 + σ S ρ S υ 2 2 φ x υ 2 + R φ ( x + β j u m p y , υ 1 , υ 2 , t , T ; ϕ ) φ ( x , υ 1 , υ 2 , t , T ; ϕ ) ν y ( d y ) + R φ ( x + ξ , υ 1 , υ 2 , t , T ; ϕ ) φ ( x , υ 1 , υ 2 , t , T ; ϕ ) ν ξ ( d ξ ) , φ ( x , υ 1 , υ 2 , t , T ; ϕ ) | t = T = e i ϕ X T .
Due to the affine structure of our model, we postulate φ ( x , υ 1 , υ 2 , t , T ; ϕ ) admitting the form of (6). Substituting Equation (6) into the above PIDE (7) gives the following system of ordinary differential equations (ODEs) for A ( τ ) , B ( τ ) , C ( τ ) , and D ( τ ) :
A ( τ ) τ = 0 , B ( τ ) τ = 1 2 σ I 2 B 2 ( τ ) + β d i f f σ I ρ I A ( τ ) κ I B ( τ ) 1 2 β d i f f 2 A ( τ ) A 2 ( τ ) , C ( τ ) τ = 1 2 σ S 2 C 2 ( τ ) + σ S ρ S A ( τ ) κ S C ( τ ) 1 2 A ( τ ) A 2 ( τ ) , D ( τ ) τ = r A ( τ ) + κ I θ I B ( τ ) + κ S θ S C ( τ ) + R e A ( τ ) β j u m p y 1 A ( τ ) e β j u m p y 1 ν y ( d y ) + R e A ( τ ) ξ 1 A ( τ ) e ξ 1 ν ξ ( d ξ ) ,
where the boundary conditions are given as A ( 0 ) = i ϕ and B ( 0 ) = C ( 0 ) = D ( 0 ) = 0 .
By solving the above ODEs, we can obtain the characteristic function (6).  □
Lemma 2.
Suppose that the market factor I t is driven by Equation (1). Then, the conditional characteristic function of log-market factor Z T = ln I T is given by:
ψ ( z , υ 1 , t , T ; ϕ ) = E Q e i ϕ Z T | Z t = z , V I , t = υ 1 = exp A ˜ ( τ ) z + B ˜ ( τ ) υ 1 + D ˜ ( τ ) ,
where:
A ˜ ( τ ) = i ϕ , B ˜ ( τ ) = κ I i ϕ σ I ρ I d σ I 2 1 e d τ 1 g e d τ , D ˜ ( τ ) = i ϕ r + R e i ϕ y 1 i ϕ e y 1 ν y ( d y ) I 3 τ + κ I θ I σ I 2 κ I i ϕ σ I ρ I d τ 2 ln 1 g e d τ 1 g , g = κ I i ϕ σ I ρ I d κ I i ϕ σ I ρ I + d , d = i ϕ σ I ρ I κ I 2 + σ I 2 ( i ϕ + ϕ 2 ) ,
and τ = T t .
Proof. 
Similar to the proof of Lemma 1, we can easily verify the above results.  □

2.3. Valuation of the European Index and Equity Options

Once the characteristic function is found, it is straightforward to calculate the prices of European options by using Fourier inversion. Let C ( S t , T , K ) and P ( S t , T , K ) be the prices of the European equity call and put options at time t with strike price K and maturity T under the risk neutral measure Q , respectively. Then, these option prices are determined by:
C ( S t , T , K ) = e r τ E t Q max ( S T K , 0 )
and:
P ( S t , T , K ) = e r τ E t Q max ( K S T , 0 )
where τ = T t is the time to maturity.
Theorem 1.
Suppose that the market factor I t and the individual equity price S t are driven by Equations (1) and (3), respectively. Then, the prices of the European equity call and put options with strike price K and maturity τ = T t are given by:
C ( S t , T , K ) = S t Π 1 S t , T , K ; β d i f f , β j u m p K e r τ Π 2 S t , T , K ; β d i f f , β j u m p
and:
P ( S t , T , K ) = K e r τ 1 Π 2 S t , T , K ; β d i f f , β j u m p S t 1 Π 1 S t , T , K ; β d i f f , β j u m p
where the risk neutral probability distribution functions Π 1 and Π 2 are defined by:
Π 1 S t , T , K ; β d i f f , β j u m p = 1 2 + e r τ π S t 0 + e i ϕ ln K φ ( x , υ 1 , υ 2 , t , T ; ϕ i ) i ϕ d ϕ
and:
Π 2 S t , T , K ; β d i f f , β j u m p = 1 2 + 1 π 0 + e i ϕ ln K φ ( x , υ 1 , υ 2 , t , T ; ϕ ) i ϕ d ϕ ,
where φ ( x , υ 1 , υ 2 , t , T ; ϕ ) is the conditional characteristic function of ln S T , which can be seen in Equation (6), and [ · ] indicates the real part of a complex number.
Proof. 
In order to get the pricing formulas of the European equity call and put options, let us first introduce a change of measure from Q to Q ˜ by the following Radon–Nikodym derivative:
d Q ˜ d Q = e r ( T t ) S T S t .
We denote the conditional characteristic function of X T = ln S T under the Q ˜ measure by φ ˜ ( x , υ 1 , υ 2 , t , T ; ϕ ) . Then, φ ˜ ( x , υ 1 , υ 2 , t , T ; ϕ ) can be expressed as:
φ ˜ ( x , υ 1 , υ 2 , t , T ; ϕ ) = E t Q ˜ e i ϕ X T = E t Q e r ( T t ) S T S t e i ϕ X T = e r ( T t ) x E t Q e i ( ϕ i ) X T = e r ( T t ) x φ ( x , υ 1 , υ 2 , t , T ; ϕ i ) .
Thus, the price of a European equity call option C ( S t , T , K ) can be calculated by utilizing φ ( x , υ 1 , υ 2 , t , T ; ϕ ) and φ ˜ ( x , υ 1 , υ 2 , t , T ; ϕ ) :
C ( S t , T , K ) = e r τ E t Q max ( S T K , 0 ) = e r τ E t Q S T 1 { S T K } K e r τ E t Q 1 { S T K } = S t E t Q ˜ 1 { S T K } K e r τ E t Q 1 { S T K } = S t E t Q ˜ 1 { X T ln K } K e r τ E t Q 1 { X T ln K } = S t Π 1 S t , T , K ; β d i f f , β j u m p K e r τ Π 2 S t , T , K ; β d i f f , β j u m p .
Once the conditional characteristic function φ ( x , υ 1 , υ 2 , t , T ; ϕ ) is obtained, we can easily calculate the probability distribution functions Π 1 S t , T , K ; β d i f f , β j u m p and Π 2 S t , T , K ; β d i f f , β j u m p according to the Lévy inversion formula:
Π 1 S t , T , K ; β d i f f , β j u m p = 1 2 + 1 π 0 + e i ϕ ln K φ ˜ ( x , υ 1 , υ 2 , t , T ; ϕ ) i ϕ d ϕ
and:
Π 2 S t , T , K ; β d i f f , β j u m p = 1 2 + 1 π 0 + e i ϕ ln K φ ( x , υ 1 , υ 2 , t , T ; ϕ ) i ϕ d ϕ ,
A similar approach can be used to derive the pricing formula for the European equity put option.  □
In a similar way, we also can present the pricing formulas for the European index call and put options.
Theorem 2.
Suppose that the market factor I t is driven by Equation (1). Then, the time t prices of the European index call and put options with strike price K and maturity τ = T t are given by:
C ( I t , T , K ) = I t Π ˜ 1 I t , T , K K e r τ Π ˜ 2 I t , T , K
and:
P ( I t , T , K ) = K e r τ 1 Π ˜ 2 I t , T , K I t 1 Π ˜ 1 I t , T , K
where the risk neutral probability distribution functions Π 1 and Π 2 are defined by:
Π ˜ 1 I t , T , K = 1 2 + e r τ π I t 0 + e i ϕ ln K ψ ( z , υ 1 , t , T ; ϕ i ) i ϕ d ϕ
and:
Π ˜ 2 I t , T , K = 1 2 + 1 π 0 + e i ϕ ln K ψ ( z , υ 1 , t , T ; ϕ ) i ϕ d ϕ ,
where ψ ( z , υ 1 , t , T ; ϕ ) is the conditional characteristic function of ln I T , which can be seen in Equation (8).

3. Empirical Studies

In this section, we empirically compare the pricing performance of our proposed model with those of the classical two factor stochastic volatility models, such as Bates (2000) (two variance SVmodel with price jumps, 2-SVJmodel), Christoffersen et al. (2009) (two-variance SV model, 2-SV model), and Christoffersen et al. (2018) (two-variance SV model with a single market factor, 2-FSVmodel), being taken as benchmark models.

3.1. Data Description

We used the S&P 500 index (SPX) to proxy for the market factor and AAPL as the individual equity. We employed the delayed market quotes on arbitrary date 8 May 2019, which was the last date available at the time of writing, as the in-sample data to calibrate the risk neutral parameters, and those on 9 May 2019 were used for the out-of-sample test. We used mid-quotes to represent the option prices. To eliminate the sample noise in raw option data, we adopted some filtering rules commonly used within the related literature: (i) we omitted those options with fewer than seven days and more than 365 days to maturity; (ii) all observations with zero trading volume were discarded; (iii) all options with implied volatility equal to zero and larger than 1.0 were discarded. In addition, for the convenience of the empirical analysis in the following, we only considered the sample data of the index call options and individual equity call options with the same expiration date. Thus, we focused only on ten maturities slices, namely on the maturities of 24 May 2019, 31 May 2019, 7 June 2019, 14 June 2019, 21 June 2019, 19 July 2019, 16 August 2019, 20 September 2019, 18 October 2019, and 17 January 2020.
After these filters, we had a total of 401 observations for the S&P 500 index call option on 8 May 2019. The individual equity option sample contained 233 call options on 8 May 2019 and 264 call options on 9 May 2019, respectively. Due to the life of an option being usually less than one year, we chose the three month U.S. Treasury Bill Rate to substitute for the risk free interest rate. All of the data were downloaded from the Chicago Board Options Exchange (http://www.cboe.com/).

3.2. Parameter Estimation

Our proposed model allowed a general distribution for jump components of the market factor and individual equity price and thus could be easily introduced to the special cases such that the jump components follow the compound Poisson process of Merton (1976) and Kou (2002), etc. For different types of Lévy kernels, different forms of our model can be presented. In order to keep consistent with Bates (2000) for comparative analysis, in the following, we assumed that the jump components of the dynamics for the market factor and individual equity followed compound Poisson processes and the jump magnitude was drawn from the log-normal distribution of Merton (1976). Thus, the Lévy kernels for the market factor and individual equity, respectively, are given by:
ν y ( d y ) = λ I 1 2 π δ I 2 exp y μ I 2 2 δ I 2 d y
and:
ν ξ ( d ξ ) = λ S 1 2 π δ S 2 exp ξ μ S 2 2 δ S 2 d ξ ,
where λ j , for j = I , S , denotes the jump intensity, μ j is the mean of the jump size, and δ j is the variance of the jump size. Then, the integrals I i , for i = 1 , 2 , 3 , in Lemmas 1 and 2 can be calculated as follows:
I 1 = λ I e i ϕ β j u m p μ I 1 2 ϕ 2 β j u m p 2 δ I 2 1 i ϕ e β j u m p μ I + 1 2 β j u m p 2 δ I 2 1 , I 2 = λ S e i ϕ μ S 1 2 ϕ 2 δ S 2 1 i ϕ e μ S + 1 2 δ S 2 1 ,
and
I 3 = λ I e i ϕ μ I 1 2 ϕ 2 δ I 2 1 i ϕ e μ I + 1 2 δ I 2 1 .
Based on Theorems 1 and 2, we employed a two step calibration procedure (see, for example, Wong et al. 2012; Christoffersen et al. 2018) to estimate the model parameters. First, we calibrated the market index dynamic Θ I based on the S&P 500 index option price alone. Second, we used the equity option price to calibrate the individual equity dynamic Θ S conditional on estimates of Θ I . Consider the situation in which an investor wants to hedge his or her equity position with index options and hedging horizon T. For brevity, we further suppose that the investor observes index option prices and equity option prices both with maturity T, the same as hedging horizon. Specifically, the dataset contains M t index option prices C ( I t , T , K i ) , for i = 1 , 2 , , M t , and N t equity option prices C ( S t , T , K j ) , for j = 1 , 2 , , N t .
In the calibration process, the risk neutral model parameters were backed out by minimizing a loss function capturing the fit between the theoretical model and market prices. We employed the root mean squared errors (RMSE) as the objective function. The first step calibrated the risk neutral parameters for the index process, which are calibrated by:
RMSE ( I ) = arg min Θ I 1 M t i = 1 M t C i , m a r k e t ( I t , T , K i ) C i , m o d e l Θ I ( I t , T , K i ) 2 ,
where C i , m a r k e t ( I t , T , K i ) is the market price of the index call option contract from the sample and C i , m o d e l Θ I ( I t , T , K i ) represents the model price calculated using Equation (15) and the vector of model input parameters Θ I .
The second calibrated the beta and the parameters for the idiosyncratic risk:
RMSE ( S ) = arg min Θ S 1 N t j = 1 N t C j , m a r k e t ( S t , T , K j ) C j , m o d e l Θ S ( S t , T , K j ) 2 ,
where C j , m a r k e t ( S t , T , K j ) is the market price of the equity call option contract from the sample and C j , m o d e l Θ S ( S t , T , K j ) represents the model price calculated using Equation (13) and the vector of model input parameters Θ S .
On the basis of the above calibration method, Table 1 presents the risk neutral parameter estimates across various model specifications. Note that the values of the diffusive beta β d i f f and jump beta β j u m p for our proposed model were 0.3891 and 0.8429, respectively. The corresponding value of β d i f f for the 2-FSV model was 0.2457. Obviously, both our proposed model and the 2-FSV model showed that AAPL tended to have a relatively low exposure to diffusive market movements. However, the jump exposure coefficient β j u m p = 0.8429 indicated that the AAPL had a strong exposure to market jumps, which meant that the factor structure of the jumps was much stronger than the one of the diffusive movements. The reason for this result may be related to the sample data we selected. If we can get more sample data in the future, we will do an in-depth analysis. Moreover, we also can see that the values of correlation ρ were strongly negative for four models, capturing the so-called leverage effect both in the index and individual equity.

3.3. Pricing Performance

In this subsection, we present the empirical results for the calibrated models. In order to investigate the impacts of the systematic and idiosyncratic volatility and jump risks on equity option pricing, we took the 2-FSV, 2-SV, and 2-SVJ models as benchmark models to evaluate the pricing performance of our proposed model.
Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 exhibit the predicted prices of the four model specifications and market prices listed on 9 May 2019, with 11, 16, 21, 26, 31, 51, 71, 96, 116, and 181 trading days to expiry, respectively. Here, the predicted prices (out-of-sample pricing) were calculated by the in-sample calibration parameters reported in Table 1. One can clearly observe from the left panels of Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 that the option prices obtained by theoretical models were generally closer to the market prices for different strike prices. To further investigate the pricing performance of the four models, the right panels of Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 show the relative price differences (relative errors) between the theoretical model prices and market prices.5 For simplicity, we refer to a call option as deep out-of-the-money (DOTM) if S / K 0.90 ; out-of-the-money (OTM) if 0.90 < S / K 0.97 ; at-the-money (ATM) if 0.97 < S / K 1.03 ; in-the-money (ITM) if 0.97 < S / K 1.10 ; and deep in-the-money (ITM) if 1.10 < S / K . Moreover, we considered options less than 60 days to expiration as short term; options with 60–120 days to expiration as medium term; and options larger than 120 days to expiration as long term. For the options with 11, 16, 21, 26, 31, and 51 trading days to expiry, the relative pricing errors produced by our proposed model were all significantly lower than those of 2-FSV, 2-SV, and 2-SVJ models in the case of DOTM options, while the relative errors of all models were slightly higher.
It is also worth noting that the pricing performance of the stochastic model with jump behavior was much better than that of the model without jump in the case of deep out-of-money. For the options with 71, 96, 116, and 181 trading days to expiry, we did not find the same conclusions as the above short term options. In conclusion, the pricing performance of equity option valuation model considering market and idiosyncratic volatility and jump risks was significantly improved for short term and DOTM options.
To summarize the model calibration results, we also adopted the RMSE as a measure of the goodness of fit. Table 2 reports the out-of-sample pricing errors for the four models across different maturities. Note from Table 2 that our proposed model generally outperformed the other three models in terms of out-of-sample pricing errors. In fact, the same was true for in-sample, whose pricing errors were generally lower than those of the out-of-sample. We will not repeat them here. To measure the extent to which a model was better or worse than another, we defined the improvement rate as the relative differences between the pricing errors from the benchmark model and our proposed model, i.e.,
Improvement   rate = RMSE b e n c h m a r k RMSE o u r RMSE b e n c h m a r k × 100 %
where RMSE o u r and RMSE b e n c h m a r k denote the RMSE implied by our model and benchmark model, respectively. A positive (or negative) value of improvement rate meant that our model yielded lower (or higher) pricing errors than benchmark model, implying that the pricing performance of the former was better (or worse) than that of the latter by a percentage of that value.
From the last column of Table 2, we can see that our model was superior to the 2-SVJ model across different maturities, which meant that it was necessary to consider the market factor structure in equity option pricing. From the third last column of Table 2, the improvement rate indicated that our model slightly outperformed the 2-FSV model in terms of short term options, but was worse than that of both medium and long term. In spite of this, our empirical study presented here could at least illustrate that the equity option pricing model considering systematic and idiosyncratic volatility and jump risks may offer a good competitor of the models of Bates (2000), Christoffersen et al. (2009), or Christoffersen et al. (2018) for some other equity option markets.

4. Conclusions

In Christoffersen et al. (2018), the issues of the equity volatility levels, skews, and term structures were investigated by using equity option prices and the principal component analysis method. Their empirical results indicated that the equity options had a strong factor structure, and then, they developed an equity option pricing model with a CAPM factor structure and stochastic volatility. In addition, jumps in stock returns of individual firms were triggered by either systematic events or idiosyncratic shocks. Some recent studies indicated that idiosyncratic jumps were a key important determinant of expected stock; see, for example, Xiao and Zhou (2018), Kapadia and Zekhnini (2019) and Bégin et al. (2020).
Motivated by these insights, we developed a novel model for pricing individual equity options that incorporated a market factor structure, which could be seen as a generalized version of the work by Christoffersen et al. (2018). Specifically, in our model, the individual equity prices were driven by the market factor, as well as an idiosyncratic component that also had stochastic volatility and jump. Due to our model belonging to the affine class, we derived the closed-form solutions for the prices of both the market index and individual equity options by utilizing the Fourier inversion. In addition, we provided the empirical results to test the pricing performance of our proposed factor model based on the S&P 500 index and the AAPL stock on options. Toward this end, we empirically compared the pricing performance of our proposed model with those of the other three classical two factor stochastic volatility models being taken as benchmark models. The out-of-sample pricing performance of equity option valuation model considering market and idiosyncratic volatility and jump risks was significantly improved for short term and DOTM options. In conclusion, the empirical results presented here at least confirmed that the equity option pricing model considering systematic and idiosyncratic volatility and jump risks may offer as good competitor of the models of Bates (2000), Christoffersen et al. (2009), or Christoffersen et al. (2018) for some other option markets.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 71901124) and the Natural Science Foundation of Jiangsu Province (Grant No. BK20190695).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Andersen, Torben G., Nicola Fusari, and Viktor Todorov. 2015. The risk premia embedded in index options. Journal of Financial Economics 117: 558–84. [Google Scholar] [CrossRef] [Green Version]
  2. Bakshi, Gurdip, Charles Cao, and Zhiwu Chen. 1997. Empirical performance of alternative option pricing models. Journal of Finance 52: 2003–49. [Google Scholar] [CrossRef]
  3. Bakshi, Gurdip, Nikunj Kapadia, and Dilip Madan. 2003. Stock return characteristics, skew laws, and the differential pricing of individual equity options. Review of Financial Studies 16: 101–43. [Google Scholar] [CrossRef]
  4. Bakshi, Gurdip, and Nikunj Kapadia. 2003a. Delta-hedged gains and the negative market volatility risk premium. Review of Financial Studies 16: 527–66. [Google Scholar] [CrossRef]
  5. Bakshi, Gurdip, and Nikunj Kapadia. 2003b. Volatility risk premiums embedded in individual equity options: Some new insights. Journal of Derivatives 11: 45–54. [Google Scholar] [CrossRef]
  6. Bardgett, Chris, Elise Gourier, and Markus Leippold. 2019. Inferring volatility dynamics and risk premia from the S&P 500 and VIX markets. Journal of Financial Economics 131: 593–618. [Google Scholar]
  7. Bates, David S. 1996. Jumps and stochastic volatility: Exchange rate processes implicit in Deutsche mark options. Review of Financial Studies 9: 69–107. [Google Scholar] [CrossRef]
  8. Bates, David S. 2000. Post-’87 crash fears in the S&P 500 futures option market. Journal of Econometrics 94: 181–238. [Google Scholar]
  9. Broadie, Mark, Mikhail Chernov, and Michael Johannes. 2007. Model specification and risk premia: Evidence from futures options. Journal of Finance 62: 1453–90. [Google Scholar] [CrossRef] [Green Version]
  10. Bégin, Jean-François, Christian Dorion, and Geneviève Gauthier. 2020. Idiosyncratic jump risk matters: Evidence from equity returns and options. Review of Financial Studies 33: 155–211. [Google Scholar] [CrossRef]
  11. Carr, Peter, and Dilip B. Madan. 2012. Factor models for option pricing. Asia-Pacific Financial Markets 19: 319–29. [Google Scholar] [CrossRef]
  12. Cheang, Gerald H. L., Carl Chiarella, and Andrew Ziogas. 2013. The representation of American options prices under stochastic volatility and jump-diffusion dynamics. Quantitative Finance 13: 241–53. [Google Scholar] [CrossRef]
  13. Cheang, Gerald H. L., and Len Patrick Dominic M. Garces. 2019. Representation of exchange option prices under stochastic volatility jump-diffusion dynamics. Quantitative Finance. [Google Scholar] [CrossRef] [Green Version]
  14. Christoffersen, Peter, Kris Jacobs, and Chayawat Ornthanalai. 2012. Dynamic jump intensities and risk premiums: Evidence from S&P 500 returns and options. Journal of Financial Economics 106: 447–72. [Google Scholar]
  15. Christoffersen, Peter, Mathieu Fournier, and Kris Jacobs. 2018. The factor structure in equity options. Review of Financial Studies 31: 595–637. [Google Scholar] [CrossRef] [Green Version]
  16. Christoffersen, Peter, Steven Heston, and Kris Jacobs. 2009. The shape and term structure of the index option smirk: Why multifactor stochastic volatility models work so well. Management Science 55: 1914–32. [Google Scholar] [CrossRef] [Green Version]
  17. Duffie, Darrell, Jun Pan, and Kenneth Singleton. 2000. Transform analysis and asset pricing for affine jump diffusions. Econometrica 68: 1343–76. [Google Scholar] [CrossRef] [Green Version]
  18. Eraker, Biørn, Miichael Johannes, and Nicholas Polson. 2003. The Impact of Jumps in Volatility and Returns. Journal of Finance 58: 1269–1300. [Google Scholar] [CrossRef]
  19. Fouque, Jean-Pierre, and Adam P. Tashman. 2012. Option pricing under a stressed-beta model. Annals of Finance 8: 183–203. [Google Scholar] [CrossRef] [Green Version]
  20. Fouque, Jean-Pierre, and Eli Kollman. 2011. Calibration of stock betas from skews of implied volatilities. Applied Mathematical Finance 18: 119–37. [Google Scholar] [CrossRef]
  21. Kapadia, Nishad, and Morad Zekhnini. 2019. Do idiosyncratic jumps matter? Journal of Financial Economics 131: 666–92. [Google Scholar] [CrossRef]
  22. Kou, Steven G. 2002. A jump-diffusion model for option pricing. Management Science 48: 1086–101. [Google Scholar] [CrossRef] [Green Version]
  23. Merton, Robert C. 1976. Option pricing when underlying stock returns are discontinuous. Journal of Financial Economics 3: 125–44. [Google Scholar] [CrossRef] [Green Version]
  24. Wong, Hoi Ying, Edwin Kwan Hung Cheung, and Shiu Fung Wong. 2012. Lévy betas: Static hedging with index futures. Journal of Futures Markets 32: 1034–59. [Google Scholar] [CrossRef]
  25. Xiao, Xiao, and Chen Zhou. 2018. The decomposition of jump risks in individual stock returns. Journal of Empirical Finance 47: 207–28. [Google Scholar] [CrossRef]
1.
In fact, the work of Xiao and Zhou (2018) is a complement to the recent studies that disentangle the four types of risks in equity premiums, such as Bégin et al. (2020), who developed a GARCH-jump model in which an individual firm’s systematic and idiosyncratic risk have both a Gaussian diffusive and a jump component. Their empirical results showed that normal diffusive and jump risks have drastically different effects on the expected return of individual stocks by using 20 years of returns and options on the S&P 500 and 260 stocks.
2.
One can refer to Assumption 2.1 of Cheang et al. (2013) and Cheang and Garces (2019) for a more detailed explanation.
3.
Obviously, our proposed model for the dynamics of the market factor and individual equity prices is an extension of Christoffersen et al. (2018). In fact, our model also can be regarded as a further generalization of Cheang et al. (2013) and Cheang and Garces (2019) by taking into account the factor structure.
4.
One can refer to the Assumption 2.1 of Cheang et al. (2013) and Cheang and Garces (2019) for a more detailed explanation.
5.
The relative error is defined by | C m o d e l C m a r k e t | C m a r k e t × 100 % , where C m o d e l and C m a r k e t denote the theoretical model option prices and the real market prices, respectively.
Figure 1. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 24 May 2019.
Figure 1. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 24 May 2019.
Jrfm 13 00016 g001
Figure 2. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 31 May 2019.
Figure 2. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 31 May 2019.
Jrfm 13 00016 g002
Figure 3. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 7 June 2019.
Figure 3. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 7 June 2019.
Jrfm 13 00016 g003
Figure 4. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 14 June 2019.
Figure 4. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 14 June 2019.
Jrfm 13 00016 g004
Figure 5. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 21 June 2019.
Figure 5. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 21 June 2019.
Jrfm 13 00016 g005
Figure 6. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 19 July 2019.
Figure 6. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 19 July 2019.
Jrfm 13 00016 g006
Figure 7. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 16 August 2019.
Figure 7. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 16 August 2019.
Jrfm 13 00016 g007
Figure 8. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 20 September 2019.
Figure 8. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 20 September 2019.
Jrfm 13 00016 g008
Figure 9. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 18 October 2019.
Figure 9. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 18 October 2019.
Jrfm 13 00016 g009
Figure 10. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 17 January 2019.
Figure 10. The comparison of predicted prices of four model specifications and market prices on 9 May 2019, with maturity T = 17 January 2019.
Jrfm 13 00016 g010
Table 1. Estimated parameters. Note: This table shows the average of the estimated parameters obtained by minimizing the root mean squared pricing errors between the market price and the model price for each option on 8 May 2019. Standard errors are reported in parentheses.
Table 1. Estimated parameters. Note: This table shows the average of the estimated parameters obtained by minimizing the root mean squared pricing errors between the market price and the model price for each option on 8 May 2019. Standard errors are reported in parentheses.
Parameters Our 2-FSV 2-SV 2-SVJ
SPX AAPL SPX AAPL AAPL AAPL
V I , 0 / V 1 , 0 0.0133 0.0119 0.0239 0.0181
(0.0000) (0.0000) (0.0002) (0.0001)
V S , 0 / V 2 , 0 0.0470 0.0514 0.0197 0.0176
(0.0000) (0.0000) (0.0002) (0.0002)
κ I / κ 1 0.2496 0.2929 0.3489 0.4064
(0.0212) (0.0148) (0.0118) (0.0311)
κ S / κ 2 0.2454 0.1504 0.4131 0.4108
(0.0288) (0.0797) (0.0729) (0.0171)
θ I / θ 1 0.2820 0.3066 0.3314 0.2817
(0.0181) (0.0317) (0.0534) (0.0348)
θ S / θ 2 0.2303 0.3683 0.2447 0.3415
(0.0190) (0.0590) (0.0365) (0.0423)
σ I / σ 1 0.3472 0.3932 0.1615 0.1898
(0.0127) (0.0137) (0.0081) (0.0106)
σ S / σ 2 0.1496 0.1640 0.2206 0.1970
(0.0056) (0.0135) (0.0386) (0.0059)
λ I 0.0450
(0.0017)
λ S 0.3413 0.3065
(0.2463) (0.1194)
μ I 0.1657
(0.0599)
μ S 0.0889 0.0333
(0.0391) (0.0042)
δ I 0.0850
(0.0113)
δ S 0.0679 0.0534
(0.0078) (0.0013)
β d i f f 0.3891 0.2457
(0.0381) (0.0983)
β j u m p 0.8429
(0.8091)
ρ I / ρ 1 −0.9290 −0.8498 −0.9222 −0.7445
(0.0063) (0.0080) (0.0096) (0.0297)
ρ S / ρ 2 −0.9926 −0.8938 −0.7673 −0.7817
(0.0001) (0.0469) (0.1632) (0.0549)
Table 2. Out-of-sample pricing errors. Note: This table shows the out-of-sample pricing errors across different maturities. Pricing errors are reported as the root mean squared errors (RMSE) of option prices for four models.
Table 2. Out-of-sample pricing errors. Note: This table shows the out-of-sample pricing errors across different maturities. Pricing errors are reported as the root mean squared errors (RMSE) of option prices for four models.
RMSEOur2-FSV2-SV2-SVJImprovement Rate
MaturityOur vs. 2-FSVOur vs. 2-SVOur vs. 2-SVJ
24 May 20190.25730.25740.25960.27070.0373%0.8803%4.9568%
31 May 20190.25070.25080.25640.26520.0392%2.2499%5.4846%
7 June 20190.23430.23470.25270.24740.1764%7.2947%5.3044%
14 June 20190.19920.20410.22610.20992.4278%11.9155%5.0858%
21 June 20190.18240.18270.18730.19160.1399%2.5963%4.7934%
19 July 20190.32560.33010.33260.33831.3434%2.0948%3.7368%
16 August 20190.28560.28350.28790.2922−0.7573%0.7946%2.2384%
20 September 20190.31770.31590.31620.3222−0.5932%-0.4851%1.4002%
18 October 20190.11850.11800.12150.1272−0.4458%2.4886%6.8593%
17 January 20200.48820.48820.48930.4943−0.0071%0.2182%1.2201%

Share and Cite

MDPI and ACS Style

Li, Z. Equity Option Pricing with Systematic and Idiosyncratic Volatility and Jump Risks. J. Risk Financial Manag. 2020, 13, 16. https://doi.org/10.3390/jrfm13010016

AMA Style

Li Z. Equity Option Pricing with Systematic and Idiosyncratic Volatility and Jump Risks. Journal of Risk and Financial Management. 2020; 13(1):16. https://doi.org/10.3390/jrfm13010016

Chicago/Turabian Style

Li, Zhe. 2020. "Equity Option Pricing with Systematic and Idiosyncratic Volatility and Jump Risks" Journal of Risk and Financial Management 13, no. 1: 16. https://doi.org/10.3390/jrfm13010016

APA Style

Li, Z. (2020). Equity Option Pricing with Systematic and Idiosyncratic Volatility and Jump Risks. Journal of Risk and Financial Management, 13(1), 16. https://doi.org/10.3390/jrfm13010016

Article Metrics

Back to TopTop