﻿ Estimation of Stochastic Volatility with a Compensated Poisson Jump Using Quadratic Variation

Applied Mathematics
Vol.08 No.07(2017), Article ID:77969,14 pages
10.4236/am.2017.87077

Estimation of Stochastic Volatility with a Compensated Poisson Jump Using Quadratic Variation

Perpetual Saah Andam1, Joseph Ackora-Prah1, Sure Mataramvura2

1Department of Mathematics, Kwame Nkrumah University of Science and Technology (KNUST), Kumasi, Ghana

2Department of Actuarial Science, University of Cape Town, Cape Town, South Africa

Received: June 9, 2017; Accepted: July 24, 2017; Published: July 27, 2017

ABSTRACT

The degree of variation of trading prices with respect to time is volatility- measured by the standard deviation of returns. We present the estimation of stochastic volatility from the stochastic differential equation for evenly spaced data. We indicate that, the price process is driven by a semi-martingale and the data are evenly spaced. The results of Malliavin and Mancino [1] are extended by adding a compensated poisson jump that uses a quadratic variation to calculate volatility. The volatility is computed from a daily data without assuming its functional form. Our result is well suited for financial market applications and in particular the analysis of high frequency data for the computation of volatility.

Keywords:

Stochastic Volatility, Compensated Poisson Jump, Quadratic Variation

1. Introduction

Volatility is the downward and upward movements of the market. Naturally stock prices attract volatility. Economists argue that when demand is high, price goes up and producers are willing to sell more because a lot of people want to buy the same goods. Also higher supply leads to price fall and producers are willing to supply less because a lot of people do not want to pay more for goods that can be easily found in the market. The biggest driver of volatility is a drop in the market which later goes up for a while and comes down again. Mostly, when the economy is in a very precarious state, investors intervene to help solve the crisis which causes the stock market to revive and recover from its decline [2] . High volatility brings panic to investors and this fear has an impact on the financial market. Volatility of a financial instrument measures the uncertainty of the returns. It is an important parameter in pricing of an asset, portfolio management, investment analysis and risk management. Volatility influences cash flow from selling asset at a specific future date. Higher volatility means a greater chance of a fall in the price of an asset. Investors find it difficult to determine the movement of the prices of stock because of the stochastic nature of volatility which makes it difficult to predict stock prices. In 1973, Fischer Black and Myron Scholes published a paper called “The Pricing of Options and Corporate Liabilities’’ [3] in which volatility was assumed to be constant. In reality, real stock returns have higher kurtosis; there is a high chance of large price change. Empirical evidence has shown that, volatility assumed constant is restrictive, therefore estimation of volatility is the first step to adjust the Black-Scholes model to real data. Volatility cannot be directly observed and must therefore be estimated from historical data.

In most recent financial econometrics literature, volatility is regarded as a high degree of persistence and time-varying [4] . The development of new methods for estimating volatility still play active role in research. There is much evidence of volatility affecting financial assets and so estimating the ex-post volatility is an important tool for financial economics research. The common way to estimate volatility is by summing the squared returns. Although methods for volatility estimate work under some assumptions, it runs into a problem when the data exhibit some kinds of jumps in real life. Mostly since volatility cannot be observed, most researchers fit it in parametric econometrics models like Generalized Autoregressive Conditional Heteroscedasticity (GARCH) and its family [5] . They also compute implied volatilities using Black-Scholes for option pricing or compute the absolute returns or ex-post squared [6] . Daily data available on prices of asset which are speculative allows quadratic variation method to be used to measure the activities of returns in financial market. Most researchers have advocated for quadratic variation as a nonparametric method to estimate volatility of an asset. Most researches done on quadratic variation use tick-by-tick intra daily data, weekly or monthly data [6] [7] [8] . Quadratic variation controls micro-structure effect, it is consistent and feasible [9] and the error of realized volatility of its moment and asymptotic distribution has been studied [7] . For Pure Jump and mixed Jump-Diffusion Processes, it is approxi- mately free of measurement error under general conditions [6] . Estimation of instantaneous volatility or spot volatility for high frequency data using quadratic variation is being researched recently and has been tested on some simulations [10] [11] and it has numerical derivative involved in it [12] [13] [14] . Quadratic variation have been used to estimate the instantaneous volatility of an asset price without jump diffusion process [1] . A proposed extension of a rolling sample variance with a continuous record asymptotic analysis of the use of quadratic variation for the estimation of the asset returns has been done [15] . The inte- grated volatility was treated as a continuous time stochastic process using a high frequency data and established the link between integrated volatility estimators and spot volatilities [15] .

Unlike many other researches that has been done on quadratic variation, there has not been a consideration of the addition of compensated poisson jump to the stochastic differential equation model. In view of this, we are motivated to propose a quadratic variation as a nonparametric method for estimating instan- taneous volatility from a stochastic differential equation model by including compensated poisson jump as an extension. This jump process has an intensity which is able to capture the steepness and skewness in volatility smiles for short-dated options. The presence of the jump will add more reality to financial market applications.

Section 2 reviews some mathematical preliminaries used in estimating vola- tility.

Section 3 presents the estimation of stochastic volatility using quadratic va- riation.

Section 4 presents the analysis of the results using data.

Section 5 concludes the work.

2. Mathematical Preliminaries

In this section, we review some mathematical theories which are used in esti- mating volatility.

Definition 1 A real valued process X defined on the filtered probability space $\left(\Omega ,\mathcal{F},{\mathcal{F}}_{t},ℙ\right),t\ge 0$ is called a semi-martingale if it can be decomposed in the form of

$X\left(t\right)=M\left(t\right)+A\left(t\right)$

where $M$ is a local martingale and $A$ is a càdlàg adapted process of locally bounded variation [16] .

Definition 2 Suppose $\left(M,d\right)$ is a metric space, and $E\subseteq ℝ$ , a function $f:E\to M$ is a càdlàg function if $f$ is right-continuous with left limits. That is for every $t\in E$ ,

• the left limit $f\left({t}_{-}\right):=li{m}_{s↑t}f\left(s\right)$ exists and

• the right limit $f\left({t}_{+}\right):=li{m}_{s↓t}f\left(s\right)$ exists and equals $f\left(t\right)$ [17] .

Definition 3 A poisson process is a stochastic process ${N}_{t},t\in {ℝ}_{+}$ which has a jump size of 1 only and has a constant path between two jumps. At time, $t$ , the value of ${N}_{t}$ is

${N}_{t}=\underset{j=1}{\overset{\infty }{\sum }}\text{ }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)$

where

${I}_{\left[{T}_{j},\infty \right)}\left(t\right)=\left(\begin{array}{l}1\text{ }\text{if}\text{\hspace{0.17em}}t\ge {T}_{j}\hfill \\ 0\text{ }\text{if}\text{\hspace{0.17em}}0\le t<{T}_{j}.\hfill \end{array}$

$j\ge 1$ and ${T}_{j}$ is the increasing family of jump times of ${N}_{t}:{\mathrm{lim}}_{j\to \infty }{T}_{j}=+\infty$ where

${N}_{t}$ satisfies:

• Independent increments: $\forall \left\{0\le {t}_{0}<{t}_{1}<\cdots <{t}_{n},n\ge 1\right\}$ the random variables ${N}_{{t}_{1}}-{N}_{{t}_{0}},\cdots ,{N}_{{t}_{n}}-{N}_{{t}_{n-1}}$ are independent.

• Stationary increments: ${N}_{t+h}-{N}_{s+h}$ has the same distribution as ${N}_{t}-{N}_{s}$ , $\forall h>0$ and $0\le s\le t.$ [18]

Definition 4 The compensated poisson process ${M}_{t}={N}_{t}-\lambda t$ , $t\in {ℝ}_{+}$ is a martingale with respect to its own filtration ${\mathcal{F}}_{t}$ [18] .

Definition 5 Let ${V}_{t}$ be a real-valued stochastic process defined on a pro- bability space $\left(\Omega ,\mathcal{F},ℙ\right)$ and with time $t$ that ranges over non-negative real numbers then the pth variation is defined as,

${V}_{t}=\underset{‖\Pi ‖\to 0}{\mathrm{lim}}\underset{k=1}{\overset{n}{\sum }}{|{V}_{{t}_{k}}-{V}_{{t}_{k-1}}|}^{p}$

where $\Pi \in \left[0,t\right]$ , and $‖\Pi ‖$ is the norm of the partition $0={t}_{0}<{t}_{1}<{t}_{2}<\cdots <{t}_{n}=t$ such that we have $‖\Pi ‖=\mathrm{max}\left\{\left({t}_{i}-{t}_{i-1}\right),\forall i=1,\cdots ,n\right\}$ if the above sum converges [19] .

This is true under certain conditions for example, $p=1$ defines the first variation or total variation process, for $p=2$ , the pth variation equals the quadratic variation if the sum converges. Also, it is a bounded variation if and only if for $p=1,{V}_{t}<+\infty$ .

For a generalized Itô processes,

${X}_{t}={X}_{0}+{\int }_{0}^{t}\text{ }\text{ }{\sigma }_{s}\text{d}{B}_{s}+{\int }_{0}^{t}\text{ }\text{ }{\mu }_{s}\text{d}s,$

where $B$ is a Brownian motion, its quadratic variation is given by

${\left[X\right]}_{t}={\int }_{0}^{t}\text{ }\text{ }{\sigma }_{s}^{2}\text{d}s$ [20] .

The quadratic variation of a compensated poisson process ${M}_{t}={N}_{t}-\lambda t$ is

${\left[M\right]}_{t}=\underset{s\le t}{\sum }\text{ }\text{ }\Delta {M}_{s}^{2}={N}_{t}$ [21] .

3. Stochastic Volatility Estimation

A poisson process added to the stochastic differential equation with Brownian motion dependent embodies discontinuities in the returns on stock. This helps fit a better market data with regards to the reflection of the reality in the stock market and also accurately calculate volatility because it incorporates the effect of uncertainty over the jump size [22] . When there is a jump-diffusion process, it can describe the stock prices more accurately at the expense of making the market incomplete because jumps in stock prices cannot be hedged in traded securities. When the market is incomplete then the payoffs of the option cannot be replicated which will make it difficult to price the option [22] . In order to describe the stock prices more accurately, a compensated poisson jump is in- corporated in a stochastic differential equation.

Suppose that the price process $p\left(t\right)$ of an asset follows the Itô process of a semi-martingale form:

$\text{d}p\left(t\right)=\alpha \left(t,B\right)\text{d}t+\sigma \left(t,B\right)\text{d}B\left(t\right)+\text{d}M\left(t\right),$ (1)

where $M\left(t\right)=N\left(t\right)-\lambda t$ and $N\left(t\right)$ is a poisson process with intensity $\lambda$ , $\alpha$ is the drift, $\sigma$ is the volatility, time is represented by $t$ and the Brownian motion is $B$ then $\alpha$ and $\sigma$ are adapted to a random process satisfying the conditions in (QA) below;

$E\left[{\int }_{0}^{T}{\left(\alpha \left(t\right)\right)}^{2}\text{d}t\right]<\infty$

(QA)

$E\left[{\int }_{0}^{T}{\left(\sigma \left(t\right)\right)}^{2}\text{d}t\right]<\infty$

where $\alpha \left(t\right)=\alpha \left(t,B\right),\sigma \left(t\right)=\sigma \left(t,B\right)$ and it’s bounded by $|\alpha |+|\sigma |\le c$ , where $c$ is a constant.

The time interval can be split into subsequent equal intervals which contain a certain number of observations. The volatility can be expressed in terms of its mean quadratic variation. From equation (1), $\alpha =0$ implies an efficient market [23] . The drift term makes no contribution to the quadratic variation process [24] . The theorem below is the main results of this paper which exhibits the estimation of volatility with a compensated poisson jump.

Lemma 6 For the time interval $\left[t,t+\epsilon \right]$ , where $\epsilon$ denotes the interval between the two consecutive points, the volatility is constant in the interval and its value of observed quotations of the volatility is the mean quadratic variation:

${\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right].$

Proof.

$\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right]=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\frac{1}{\epsilon }\mathbb{E}\left[{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s|{\mathcal{F}}_{t}\right].$

Assuming $\sigma$ is càdlàg adapted process, if ${\mathcal{F}}_{t}$ is generated by the full observation of the price process, then using Fatou’s Lemma we have,

$\begin{array}{l}\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\frac{1}{\epsilon }\mathbb{E}\left[{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s|{\mathcal{F}}_{t}\right]\\ =\mathbb{E}\left[\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right]=\mathbb{E}\left[\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\frac{\frac{\text{d}}{\text{d}\epsilon }{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\frac{\text{d}}{\text{d}\epsilon }\left(\epsilon \right)}|{\mathcal{F}}_{t}\right]\\ =\mathbb{E}\left[\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\frac{\text{d}}{\text{d}\epsilon }{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s|{\mathcal{F}}_{t}\right]=\mathbb{E}\left[\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}{\sigma }^{2}\left(t+\epsilon \right)|{\mathcal{F}}_{t}\right]\\ =\mathbb{E}\left[{\sigma }^{2}\left(t\right)|{\mathcal{F}}_{t}\right]={\sigma }^{2}\left(t\right).\end{array}$

Hence,

$\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right]={\sigma }^{2}\left(t\right).$

The expression in Lemma 6 shows the volatility process integrated over a given interval $\left[t,t+\epsilon \right]$ . This expressions gives the general case for obtaining ${\sigma }^{2}\left(t\right)$ . On the other hand, ${\sigma }^{2}\left(t\right)$ can be estimated from a given model, for instance equation (1), as shown in the theorem below.

If volatility is estimated, it will help identify the type of financial model the stock price follows; whether it is geometric Brownian motion or Ornstein- Uhlenbeck process and many more for the purpose of option pricing. This will also help investors to evaluate the rate of their investment risks. The theorem below is the volatility estimate based on Equation (1).

Theorem 7 Suppose that the price process, $p$ , is given as $\text{d}p\left(t\right)=\alpha \left(t,B\right)\text{d}t+\sigma \left(t,B\right)\text{d}B\left(t\right)+\text{d}M\left(t\right)$ , then its volatility, $Vol\left(p\right)\left(t\right)$ is;

$Vol\left(p\right)\left(t\right)\equiv {\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}-\Delta \left(\underset{j=1}{\overset{\infty }{\sum }}\text{ }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)\right)}{\epsilon }|{\mathcal{F}}_{t}\right]$

where $\mathbb{E}\left[\cdot |{\mathcal{F}}_{t}\right]$ denotes the conditional expectation operator with respect to the s-field ${\mathcal{F}}_{t}$ .

Proof. Given a price process $\text{d}p\left(t\right)=\alpha \left(t,B\right)\text{d}t+\sigma \left(t,B\right)\text{d}B\left(t\right)+\text{d}M\left(t\right)$ , we estimate the volatility for evenly spaced time intervals using quadratic variation as follows:

$\text{d}p\left(t\right)=\alpha \left(t,B\right)\text{d}t+\sigma \left(t,B\right)\text{d}B\left(t\right)+\text{d}M\left(t\right)$

${\left(\text{d}p\left(t\right)\right)}^{2}={\left(\alpha \left(t,B\right)\text{d}t+\sigma \left(t,B\right)\text{d}B\left(t\right)+\text{d}M\left(t\right)\right)}^{2}.$

Since $\left(\text{d}B\left(t\right)\text{d}t,{\left(\text{d}t\right)}^{2}\right)\to 0$ and ${\left(\text{d}B\left(t\right)\right)}^{2}\to \text{d}t$ [25] , and $M\left(t\right)$ is a marti- ngale, it implies that $\left[\text{d}M\left(t\right),\text{d}t\right]=0,\left[\text{d}M\left(t\right),\text{d}B\left(t\right)\right]=0,$ and we obtain,

${\left(\text{d}p\left(t\right)\right)}^{2}={\sigma }^{2}\left(t\right)\text{d}t+{\left(\text{d}M\left(t\right)\right)}^{2}$

${\left(\text{d}p\left(t\right)\right)}^{2}-{\left(\text{d}M\left(t\right)\right)}^{2}={\sigma }^{2}\left(t\right)\text{d}t.$ (2)

If $\left(p,M\right):\left[t,t+\epsilon \right]\to ℝ$ , then integrating Equation (2) with respect to the interval gives:

${\int }_{t}^{t+\epsilon }{\left(\text{d}p\left(s\right)\right)}^{2}-{\int }_{t}^{t+\epsilon }{\left(\text{d}M\left(s\right)\right)}^{2}={\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s$

${\int }_{t}^{t+\epsilon }\text{d}p\left(s\right)\cdot \text{d}p\left(s\right)-{\int }_{t}^{t+\epsilon }\text{d}M\left(s\right)\cdot \text{d}M\left(s\right)={\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s.$

With respect to $p$ and $M$ , we have their quadratic variation as $\text{d}p\left(s\right)\cdot \text{d}p\left(s\right)=\text{d}{\left[p,p\right]}_{\left[s\right]}=\text{d}{\left[p\right]}_{\left[s\right]}$ and $\text{d}M\left(s\right)\cdot \text{d}M\left(s\right)=\text{d}{\left[M,M\right]}_{\left[s\right]}=\text{d}{\left[M\right]}_{\left[s\right]}$ , this implies that,

$\text{d}{\left[M\right]}_{\left[s\right]}=\text{d}\left[N-\lambda s,N-\lambda s\right]=\text{d}{\left({N}_{s}-\lambda s\right)}^{2}=\text{d}{N}_{s}.$

Then we have,

${\int }_{t}^{t+\epsilon }\text{d}{\left[p\right]}_{\left[s\right]}-{\int }_{t}^{t+\epsilon }\text{d}{\left[M\right]}_{\left[s\right]}={\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s$

${\int }_{t}^{t+\epsilon }\text{d}{\left[p\right]}_{\left[s\right]}-{\int }_{t}^{t+\epsilon }\text{d}{N}_{s}={\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s$

${\left[p\right]}_{\left[t+\epsilon ,t\right]}-{N}_{\left[t+\epsilon ,t\right]}={\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s$

${\left({p}_{t+\epsilon }-{p}_{t}\right)}^{2}-\left({N}_{t+\epsilon }-{N}_{t}\right)={\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s.$

Dividing both sides by $\epsilon$ , we have,

$\frac{{\left({p}_{t+\epsilon }-{p}_{t}\right)}^{2}}{\epsilon }-\frac{{N}_{t+\epsilon }-{N}_{t}}{\epsilon }=\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }.$

Taking the conditional expectation of both sides we have,

$\mathbb{E}\left[\left(\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}}{\epsilon }-\frac{\left({N}_{t+\epsilon }-{N}_{t}\right)}{\epsilon }\right)|{\mathcal{F}}_{t}\right]=\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right]$

$⇒\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}-\left({N}_{t+\epsilon }-{N}_{t}\right)}{\epsilon }|{\mathcal{F}}_{t}\right]=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right].$

${N}_{t+\epsilon }-{N}_{t}$ can be written as $\Delta {N}_{t}$ , but from Definition (3),

${N}_{t}={\sum }_{j=1}^{\infty }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)$ , then it follows that,

$\begin{array}{l}\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}-\Delta \left(\underset{j=1}{\overset{\infty }{\sum }}\text{ }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)\right)}{\epsilon }|{\mathcal{F}}_{t}\right]\\ =\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right].\end{array}$

From Lemma 6, we have,

${\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right].$

Therefore the volatility of a price process which has an evenly spaced interval is:

$Vol\left(p\right)\left(t\right)\equiv {\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}-\Delta \left(\underset{j=1}{\overset{\infty }{\sum }}\text{ }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)\right)}{\epsilon }|{\mathcal{F}}_{t}\right].$

This operator $\left(\mathbb{E}\left[\cdot |{\mathcal{F}}_{t}\right]\right)$ is generated by the full observation of the price process until time $t$ where ${\mathcal{F}}_{t}$ is the observation of single path of the market evolution within a given period.

This theorem is a link to Lemma 6, that is:

$\begin{array}{c}Vol\left(p\right)\left(t\right)\equiv {\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right]\\ =\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}-\Delta \left(\underset{j=1}{\overset{\infty }{\sum }}\text{ }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)\right)}{\epsilon }|{\mathcal{F}}_{t}\right].\end{array}$

If the price process is independent of its past information, then $Vol\left(p\right)\left(t\right)\equiv {\sigma }^{2}\left(t\right)$ will be,

${\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}-\Delta \left(\underset{j=1}{\overset{\infty }{\sum }}\text{ }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)\right)}{\epsilon }|{\mathcal{F}}_{t}\right]=0$

since the martingale part of the price process will vanish and the jump is esti- mated from the price process and hence the volatility will be constant.

From Theorem 7, it can be seen that the compensated poisson jump has a negative effect on the price process. Also, Theorem 7 requires high frequency data since it is the time evolution of the price of an asset in its semi-martingale form with a compensated poisson jump. Theorem 7 can be further simplified to

$Vol\left(p\right)\left(t\right)\equiv {\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}}{\epsilon }|{\mathcal{F}}_{t}\right]-\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\frac{\Delta \left(\underset{j=1}{\overset{\infty }{\sum }}\text{ }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)\right)}{\epsilon }$

but it cannot be determined numerically, because $\mathbb{E}\left[\cdot |{\mathcal{F}}_{t}\right]$ cannot be deduced from the observations but it can be determined when the distribution of the price process is known. For example, if the distribution of the price process follows a normal distribution then, we can apply antithetic variate Monte Carlo simulation to it.

Antithetic variate is defined as;

$\begin{array}{l}\frac{1}{\epsilon }{\int }_{0}^{t}{\left(p\left(s+\epsilon \right)-p\left(s\right)\right)}^{2}\\ =\frac{1}{n\epsilon }\underset{s=1}{\overset{n}{\sum }}\left(\frac{{\left(p\left(s+\epsilon \right)-p\left(s\right)\right)}^{2}+\stackrel{¯}{{\left(p\left(s+\epsilon \right)-p\left(s\right)\right)}^{2}}}{2}\right)\end{array}$

where $\stackrel{¯}{{\left(p\left(s+\epsilon \right)-p\left(s\right)\right)}^{2}}={\left(p\left(-s+\epsilon \right)-p\left(-s\right)\right)}^{2}$ , but $p\left(s+\epsilon \right)=p\left({s}_{1}+\epsilon ,{s}_{2}+\epsilon ,\cdots ,{s}_{n}+\epsilon \right)$ and $p\left(-s+\epsilon \right)=p\left(-{s}_{1}+\epsilon ,-{s}_{2}+\epsilon ,\cdots ,-{s}_{n}+\epsilon \right)=p\left(1-{s}_{1}+\epsilon ,1-{s}_{2}+\epsilon ,\cdots ,1-{s}_{n}+\epsilon \right)$ .

${s}_{1},{s}_{2},\cdots ,{s}_{n}$ are independent random numbers and are uniformly distributed on $\left(0,1\right)$ , $1-{s}_{1},1-{s}_{2},\cdots ,1-{s}_{n}$ are also uniformly distributed on $\left(0,1\right)$ and it makes ${\left(p\left(s+\epsilon \right)-p\left(s\right)\right)}^{2}$ has the same distribution as ${\left(p\left(-s+\epsilon \right)-p\left(-s\right)\right)}^{2}$ since $1-{s}_{1},1-{s}_{2},\cdots ,1-{s}_{n}$ are negatively correlated with ${s}_{1},{s}_{2},\cdots ,{s}_{n}$ .

Then $\mathbb{E}\left[\cdot |{\mathcal{F}}_{t}\right]$ can be computed as:

$\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}}{\epsilon }|{\mathcal{F}}_{t}\right]=\frac{1}{\epsilon }{\int }_{0}^{t}{\left(p\left(s+\epsilon \right)-p\left(s\right)\right)}^{2}.$ (3)

The algorithm below estimates $\mathbb{E}\left[\cdot |{\mathcal{F}}_{t}\right]$ using antithetic variate Monte Carlo simulation.

When simulating an evenly spaced data, there is a need to have a benchmark in order to assess the good prediction of the model proposed. To achieve this, we use the true volatility as the benchmark calculated as:

$\text{truevolatility}\left({\sigma }^{2}\right)=\underset{t=1}{\overset{n}{\sum }}{\left(p\left(t+1\right)-p\left(t\right)\right)}^{2}.$ (4)

Equation (4) assists in assessing the performance of different estimators used. In this work, the estimators used are; the true volatility, the instantaneous volatility with a compensated poisson jump and the instantaneous volatility proposed by Malliavin and Mancino [1] .

Malliavin and Mancino [1] proposed;

$Vol\left(p\right)\left(t\right)\equiv {\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\int }_{t}^{t+\epsilon }{\sigma }^{2}\left(s\right)\text{d}s}{\epsilon }|{\mathcal{F}}_{t}\right]$ (5)

for calculating the instantaneous volatility without a jump. In order to apply data to Equation (5), we use antithetic variate Monte Carlo simulation to evaluate it.

From the Algorithm above, 1 to 8 is used to estimate the true volatility using Equation (4). The algorithm numbered 1 to 13 is used to generate the antithetic variate Monte Carlo simulation using Equation (5) and the numbered algorithm from 1 to 19 is used to compute the instantaneous volatility with a compensated poisson jump using Theorem 7.

4. Empirical Analysis of Theorem Using Data

Stock prices indicate the strength and performance of the company. Any data can be used but with reference to this analysis we use data from AngloGold Ashanti Ghana (AGA). A daily high frequency stock prices data from AGA containing the closing prices from 16 April 2008 to 17 June 2016 was plotted with time (in years) as shown in Figure 1.

Closer to the later part of 2015 to 2016, the graph exhibits some sort of a jump. At the top of bear market, returns are earned than the bottom of a bull market. Increase in stock prices means investors are buying stocks rather than selling it and the fall in the price of stocks means investors are selling their stocks rather than buying. Also, an increase in the stock price of a company means there is a profit making and a risk reduction of layoff of workers in the company. When there is an increase in stock price, good returns are generated. When the stock price of a company falls continuously then there is a higher chance of a takeover by another company, a merge or even close down of the company.

The volatility is determined from the stock price of AngloGold Ashanti. It is used to determine the performance and strength of the stock price. It also help generate the investors return made from the stock price. The graphs below show the volatility of the stock prices obtained from AngloGold Ashanti.

The graph in Figure 2 indicates the plot of the true volatility using Equation (4). The graph in Figure 3 reveals the plot of the instantaneous volatility without a jump by Malliavin and Mancino (2009) [1] using Equation (5). The graph in Figure 4 shows the plot of the instantaneous volatility with a compensated poisson jump using Theorem 7 with the Algorithm stated above. The stock price is assumed to be normally distributed then antithetic variate Monte Carlo simulation is applied to it. The higher the volatility, the greater the returns and the higher the risk associated with the type of investment. Comparing the graphs in the stock prices as in Figure 1 and the volatility graphs in Figures 2-4, it shows that from 2009 to the mid of 2010, the stock prices increased sharply and its corresponding volatility was high which implies that investors were buying stocks rather than selling it. Also when there was a drastic decrease in stock prices from 2011 to 2014, the volatility also decreased sharply. Closer to 2015, the stock prices shot up again which made it highly volatile and might be

Figure 1. Plot of stock prices from AGA.

Figure 2. Plot of the True volatility.

Figure 3. Plot of Instantaneous volatility without a jump Malliavin and Mancino (2009) [1] .

attributed to the merge that took place between the then Ashanti Goldfields Corporation limited and AngloGold to form AngloGold Ashanti [26] . The trend of volatility helps the investor to re-balance its weights on portfolio of stocks and

Figure 4. Plot of Instantaneous volatility with a compensated poisson jump.

also help the investor knows the performance and strength of the company.

Comparison of Results

The results obtained by Malliavin and Mancino [1] was

$Vol\left(p\right)\left(t\right)\equiv {\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}}{\epsilon }|{\mathcal{F}}_{t}\right]$

from the model $\text{d}p\left(t\right)=\alpha \left(t,B\right)\text{d}t+\sigma \left(t,B\right)\text{d}B\left(t\right),$ and when compensated poisson jump was added to it we obtained

$Vol\left(p\right)\left(t\right)\equiv {\sigma }^{2}\left(t\right)=\underset{\epsilon \to {0}^{+}}{\mathrm{lim}}\mathbb{E}\left[\frac{{\left(p\left(t+\epsilon \right)-p\left(t\right)\right)}^{2}-\Delta \left(\underset{j=1}{\overset{\infty }{\sum }}\text{ }\text{ }{I}_{\left[{T}_{j},\infty \right)}\left(t\right)\right)}{\epsilon }|{\mathcal{F}}_{t}\right]$

From this, it can be seen that the compensated poisson jump has an effect on the price process. Also, error analysis of the estimates were done using the Root Mean Square Error(RMSE) calculated as;

$\text{RMSE}=\sqrt{\frac{1}{N}\underset{i=1}{\overset{N}{\sum }}{\left(\text{estimated}-\text{truevolatility}\right)}^{2}}.$

The RMSE of the instantaneous volatility proposed by Malliavin and Mancino [1] was 3.0210 and the RMSE of the instantaneous volatility with a compensated poisson jump was 3.0184. This means that, the introduction of compensated poisson jump reduced the error margin and is a good predictive model than the instantaneous volatility proposed by Malliavin and Mancino [1] . This means that the stochastic differential equation with a compensated poisson jump is a better representative model for stock prices than stochastic differential equation without a compensated poisson jump.

5. Conclusion

We have established the theoretical basis for estimating stochastic volatility with the presence of a compensated poisson jump for univariate settings for evenly spaced observations using quadratic variation. The conditional expectation, $\mathbb{E}\left[\cdot |{\mathcal{F}}_{t}\right]$ was calculated using antithetic variate Monte Carlo simulation of the price process for implementation of data to be possible. We analyzed this using data from AGA, Ghana. Also, the RMSE of the instantaneous volatility with a compensated poisson jump gave a smaller value as compared to the RMSE of the instantaneous volatility proposed by Malliavin and Mancino [1] . It is therefore advisable for financial and economic analyst to employ jump processes when using such data for forecasting else; there will be error incorporated. This will help minimize the risk (in terms of loss) associated in investing in assets that are highly volatile like gold. It will be interesting to investigate the effect of other types of jumps on evenly spaced data.

Acknowledgements

This work was carried out with financial support from the government of Canada’s International Development Research Centre (IDRC), and within the framework of AIMS Research for Africa Project and also our gratitude to the Department of Mathematics, KNUST, University of Cape Town and African Institute for Mathematical Science, Ghana.

Cite this paper

Andam, P.S., Ackora-Prah, J. and Mataramvura, S. (2017) Estimation of Stochastic Volatility with a Compensated Poisson Jump Using Quadratic Variation. Applied Mathematics, 8, 987-1000. https://doi.org/10.4236/am.2017.87077

References

1. 1. Malliavin, P. and Mancino, M.E. (2009) A Fourier Transform Method for Nonparametric Estimation of Multivariate Volatility. Annals of Statistics, 37, 1983-2010. https://doi.org/10.1214/08-AOS633

2. 2. Ibbotson, R.G. (2011) Why Does Market Volatility Matter? Yale School of Management, Yale School of Management.

3. 3. Black, F. and Scholes, M. (1973) Pricing of Options and Corporate Liabilities. Journal of Political Economy, 81, 637-654. https://doi.org/10.1086/260062

4. 4. Andersen, T.G., Bollerslev, T. and Diebold, F.X. (2002) Parametric and Nonparametric Volatility Measurement. NBER Technical Working Paper, 279.

5. 5. Barucci, E. and Reno, R. (2002) On Measuring Volatility and the GARCH Forecasting Performance. Journal of International Financial Markets, Institutions and Money, 12, 183-200.

6. 6. Andersen, T.G., Bollerslev, T., Diebold, F.X. and Labys, P. (2001) The Distribution of Realized Exchange Rate Volatility. Journal of the American Statistical Association, 96, 42-55.

7. 7. Barndorff-Nielsen, O.E. and Shephard, N. (2002) Econometric Analysis of Realized Volatility and Its Use in Estimating Stochastic Volatility Models. Royal Statistical Society, 64, 253-280.

8. 8. Zhang, L., Mykland, P.A. and Ait-Sahalia, Y. (2005) A Tale of Two Time Scales: Determining Integrated Volatility with Noisy High-Frequency Data. American Statistical Association, Theory and Methods, 100, 1394-1411.

9. 9. Large, J. (2007) Estimating Quadratic Variation When Quoted Prices Change by a Constant Increment. University of Oxford, Oxford.

10. 10. Alvarez, A., Panloup, F. and Savy, N. (2010) Estimation of the Instantaneous Volatility. American Mathematical Society. http://www.arxiv.org/abs/0812.3538v4

11. 11. Zu, Y. and Boswijk, H.P. (2013) Estimating Spot Volatility with High-Frequency Financial Data. Discussion Paper Series.

12. 12. Comte, F. and Renault, E. (1998) Long Memory in Continuous Time Stochastic Volatility. Mathematical Finance, 8, 4 291-323.

13. 13. Foster, D.P. and Nelson, D.B. (1996) Continuous Record Asymptotics for Rolling Sample Variance Estimators. Econometrica, 64, 139-174. https://doi.org/10.2307/2171927

14. 14. Mykland, P.A. and Zhang, L. (2006) ANOVA for Diffusions and Ito Processes. The Annals of Statistics, 34, 1931-1963.

15. 15. Andreou, E. and Ghysels, E. (2002) Rolling-Sample Volatility Estimators: Some New Theoretical, Simulation, and Empirical Results. Journal of Business and Economic Statistics, 20, 363-376. http://www.jstor.org/stable/1392123 https://doi.org/10.1198/073500102288618504

16. 16. He, S., Wang, J. and Yan, J. (1992) Semimartingale Theory and Stochastic Calculus. Science Press, CRC Press Inc.

17. 17. Billingsley, P. (1995) Probability and Measure. John Wiley and Sons, Hoboken.

18. 18. Privault, N. (2016) Notes on Stochastic Finance: Stochastic Calculus for Jump Processes. Technical Report, Nanyang Technological University, Singapore. http://www.ntu.edu.sg/home/nprivault/indext.html

19. 19. Karandikar, R.L. and Rao, B.V. (2007) On Quadratic Variation of Martingales. Proceedings of Indian Academy of Sciences, 124, 3 457-469.

20. 20. Protter, P.E. (2004) Stochastic Integration and Differential Equations. In: Stochastic Modeling and Applied Probability, 2nd Edition, Springer-Verlag, Berlin, Heidelberg.

21. 21. Jeanblanc, M. (2007) Jump Processes. Universite d’Evry Val d’Essonne, Cimpa School.

22. 22. Navas, J.F. (2003) Calculation of Volatility in a Jump-Diffusion Model. The Journal of Derivatives, 11, 66-72. https://doi.org/10.3905/jod.2003.319217

23. 23. Kanatani, T. (2004) High Frequency Data and Realized Volatility. Ph.D. Thesis, Graduate School of Economics, Kyoto University, Kyoto.

24. 24. Calderon, C.P. (2009) Lecture III: Review of Classic Quadratic Variation Results and Relevance to Statistical Inference in Finance. PASI, Rice University/Numerica Corporation.

25. 25. Oksendal, B. (2003) Stochastic Differential Equations: An Introduction with Applications. 6th Edition, Springer-Verlag, New York. https://doi.org/10.1007/978-3-642-14394-6

26. 26. AngloGold Ashanti (2011) History. Johannesburg, South Africa. http://www.anglogoldashanti.com/en/About-Us/history/Pages/default.aspx