Pricing Variance Swaps for Stochastic Volatilities with Delay and Jumps

We study the valuation of the variance swaps under stochastic volatility with delay and jumps. In our model, the volatility of the underlying stock price process not only incorporates jumps, which are found to be active empirically, but also exhibits past dependence: the behavior of a stock price right after a given time t depends not only on the situation at t but also on the whole past history of the process S t up to time t as well. The jump part in our model is finally represented by a general version of compound Poisson processes. We provide some analytical closed forms for the expectation of the realized variance for the stochastic volatility with delay and jumps. We also present a lower bound for delay as a measure of risk. As applications of our analytical solutions, a numerical example using S&P60 Canada Index 1998–2002 is then provided to price variance swaps.


Introduction
Variance swaps are forward contracts on future realized variance, the square of the realized volatility, which provide an easy way for investors to gain exposure to the future realized variance of the asset returns instead of directly exposure to the underlying assets.The market for such derivatives develops quickly after the collapse of LCTM in 1998 when the volatilities increased to an unprecedented high level and many investors are now interested in these derivatives to hedge volatility.Recently, several papers address the valuation of variance swaps or other volatility derivatives see 1-6 .The most recent ones are 7, 8 .Zhu and Lian 7 presented a highly efficient approach to price variance swaps with discrete sampling times.They found a closed-form exact solution for the partial differential equation system based on Heston's two-factor stochastic volatility model.The effects of jumps and International Journal of Stochastic Analysis discrete sampling on volatility and variance swaps have been investigated in 8 .For realistic contract specifications and model parameters, they found that the effect of discrete sampling is typically small while the effect of jumps can be significant.
It is known that the assumption of constant volatility in Black-Scholes model 9 is incompatible with some empirical study in the real market.Several stochastic volatility models see 10,11 are developed and able to fit some important characteristics in the market-like skews and smirks.Kazmerchuk et al. 12,13 proposed a stochastic volatility model assuming that the volatility σ t, S t depends on S t S t θ for θ ∈ −τ, 0 , namely, stochastic volatility with delay.In 13 , the Black-Scholes formula for security markets with delayed response was found, and, in 12 , the continuous-time GARCH model for stochastic volatility with delay was proposed and studied.In 14 , the stochastic volatility model with delay to price variance swaps was discussed and some analytical close forms for expectation of the realized continuously sampled variance were found both in stationary regime and in general case.The key features of this model are the following: i continuous-time analogue of discrete-time GARCH model, ii mean reversion, iii containing the same source of randomness as stock price, iv completeness of the market, v incorporating the expectation of log-return.However, there are two ways to develop this model and make it more realistic.One way is to consider multifactor stochastic volatility with delay which is developed in 15 .The other way is to consider stochastic volatility with delay and jumps which is the major work of this paper.
There are various works showing delayed response is an important factor in stock prices.Some statistical studies of stock prices see 16,17 indicate the dependence on past returns.For example, Kind et al. 18 obtained a diffusion approximation result for processes satisfying some equations with past-dependent coefficients, and they applied this result to a model of option pricing, in which the underlying asset price volatility depends on the past evolution to obtain a generalized asymptotic Black-Scholes formula.Hobson and Rogers 19 suggested a new class of nonconstant volatility models, which can be extended to include the aforementioned level-dependent model and share many characteristics with the stochastic volatility model.The volatility is nonconstant and can be regarded as an endogenous factor in the sense that it is defined in terms of the past behavior of the stock price.This is done in such a way that the price and volatility form a multidimensional Markov process.Chang and Youree 20 studied the pricing of a European contingent claim for the B, S -securities markets with a hereditary price structure in the sense that the rate of change of the unit price of the bond account and rate of change of the stock account S depend not only on the current unit price but also on their historical prices.The price dynamics for the bank account and that of the stock account are described by a linear functional differential equation and a linear stochastic functional differential equation, respectively.They show that the rational price for a European contingent claim is independent of the mean growth rate of the stock.Later, Chang and Youree 21 generalized the celebrated Black-Scholes formula to include the B, S -securities market with hereditary price structure.Clearly related to our work is the work by Mohammed et al. 22 devoted to the derivation of a delayed Black-Scholes formula for the B, S -securities market using PDE approach.Hobson and Rogers 19 also observed in their past-dependent model that the resulting implied volatility is U-shaped as a function of strike price.However, they dealt with only a special case where the model can be reduced to a system of SDEs.Unfortunately, not every past-dependent model can be reduced to a system of SDEs, and a more sophisticated approach, as developed in this paper, is needed.
Our model of stochastic volatility exhibits past dependence: the behavior of a stock price right after a given time t depends not only on the situation at t but also on the whole past history of the process S t up to time t.This draws some similarities with fractional Brownian motion models due to a long-range dependence property.Our work is also based on the GARCH 1,1 model see 23 or, more general, and the work of Duan 24 where he showed that it is possible to use the GARCH model as the basis for an internally consistent option pricing model.If we write down the last equation in differential form, we can get the continuous-time GARCH with expectation of log-returns of zero: If we incorporate nonzero expectation of log-return using It ô Lemma for ln S t /S t − τ , then we arrive at our continuous-time GARCH model for stochastic volatility with delay: We should mention that in the work of Kind et al. 18 , a past-dependent model was defined by diffusion approximation.In their model, the volatility depends on the quadratic variation of the process, while our model deals with more general dependence of the volatility on the history of the process over a finite interval.
In addition to delay, jumps are another evidence in the financial market.During the last decade, financial models based on jumps processes have acquired increasing popularity in risk management and option pricing applications.A good reference is 25 , which provides a self-contained overview of the theoretical, numerical, and empirical aspects of using jump processes in financial modeling.Stochastic volatility models with jumps are also included in this book.Some attempts have been made to incorporate jumps in stochastic volatility to price variance and volatility swaps see 26, 27 .The literature has mainly focused on two approaches: 1 time-varying volatility models that allow for market extremes to be outcome of normally distributed shocks that have a randomly changing variance and 2 models that incorporate discontinuous jumps in the asset price.Neither stochastic volatility models nor jump models have alone proven entirely empirically successful.For example, in the time-series literature, the models run into problems explaining large price movements such as the October 1987 crash.Hence, a price jump cannot explain the enormous increase in implied volatility following the crash of 1987.In response to these issues, researchers have proposed models that incorporate both stochastic volatility and jumps components.

International Journal of Stochastic Analysis
Eraker et al. 28 use returns data to investigate the performance of models with jumps in volatility using the class of jump-in-volatility models proposed by Duffie et al. 29 .The results in 28 show that the jump-in-volatility models provide a significant better fit to the returns data.
The key risk factors considered in option pricing models, besides the diffusive price risk of the underlying asset, are stochastic volatility and jumps, both in the asset price and its volatility.Models that include some or all of these factors were developed in 11, 29-32 .The importance of jumps in volatility has become apparent in recent studies, which try to explain the time series properties of both stock and option prices, for example, 28, 33 .In an asset allocation context, the main papers analyzing the impact of jumps are 34-36 .In the presence of jumps markets are incomplete and the analysis far less tractable.Technical issues aside, jumps are important because they represent a significant source of nondiversifiable risk as discussed at length in 37 .
There is currently fairly compelling evidence for jumps in the level of financial prices.The most convincing evidence comes from recent nonparametric work using high-frequency data as in 38, 39 among others.Also, the paper by Todorov and Tauchen 40 conducts a nonparametric analysis of the market volatility dynamics using high-frequency data on the VIX index compiled by the CBOE and the S&P 500 index.The data suggest that stock market volatility is best described as a pure jump process without a continuous component.Their results imply that a plausible model for stochastic volatility is a model of pure-jump type whose driving jumps come from a very active Lévy process.
The jumps in stock market volatility are found to be so active that this discredits many recently proposed stochastic volatility models without jumps see 37, 40 .Another advantage of our stochastic volatility model with delay and jumps is mean-reversion: the volatility is allowed to mean revert.Such models have shown some success in modeling interest rate e.g., 41 .The sharp decline of option implied that spot volatility following the extreme peak caused by the 1987 crash would be indicative of such a model.
In this paper, we incorporate a jump part in the stochastic volatility model with delay proposed by Swishchuk 14 to price variance swaps.We find some analytical closed forms for the expectation of the realized continuously sampled variance for stochastic volatility with delay and jumps.The jump part in our model is finally represented by a general version of compound Poisson processes, and the expectation and the covariance of the jump sizes are assumed to be deterministic functions.We note that after adding jumps, the model still keeps those good features of the previous model such as continuous-time analogue of GARCH model and mean reversion.But it is more realistic and still quick to implement.Besides, we also present a lower bound for delay as a measure of risk.As applications of our analytical solutions, a numerical example using S&P 60 Canada index 1998-2002 is also provided to price variance swaps with delay and jumps.We show that the price differences with and without jumps depend on the jump intensity of the underlying price processes, that is, how you count and estimate the number of jumps.
The rest of the paper is organized as follows.In Section 2, we introduce some basic concepts of variance swaps, and, in Section 3, we recall the pricing model of variance swaps for stochastic volatility with delay presented by Swishchuk 14 .The pricing model and analytical formulae for stochastic volatility with delay and jumps are discussed in Section 4. The main and new results are formulated in Theorems 4.1-4.6 of this section.A lower bound for delay as a measure of risk is presented in Section 5. Finally, we give a numerical example in Section 6 and conclude in Section 7.

Variance Swaps
A variance swap is a forward contract on realized variance, the square of the realized volatility.Its payoff at expiration is equal to where v R σ 2 R is the realized stock variance quoted in annual terms over the life of the contract, K var is the delivery price for variance, and N is the notional amount of the swap in dollars per annualized volatility point squared.The holder of a variance swap at expiration receives N dollars for every point by which the stock's realized variance σ 2 R has exceeded the variance delivery price K var .We note that usually N αI, where α is a converting parameter such as 1 per volatility square and I is a long-short index 1 for long and −1 for short .
The measure of realized variance which will be used is defined at the beginning of the contract.The continuous time realized variance over the life of the contract T is generally in the following form here σ 2 t is a stock variance : and its discrete version is defined as follows: where we neglect 1/n n i 1 log S t i /S t i−1 since we assume that the mean of the returns is of the order 1/n and can be neglected.The scaling by n/T ensures that these quantities are annualized daily if the maturity T is expressed in years days .
The discrete-time variance v n is unbiased variance estimation for σ t .It can be shown that see 4 Valuing a variance forward contract or swap is no different from valuing any other derivative security.The value of a forward contract on future realized variance with strike price K var is the expected present value of the future payoff in the risk-neutral world: where r is the risk-free discount rate corresponding to the expiration date T and E * denotes the expectation under the risk-neutral measure.Since the price of entering a forward contract International Journal of Stochastic Analysis is always zero, which means it costs nothing at the beginning of the contract, and by the riskneutral pricing framework, we find that the fair delivery price of a variance swap K var is just the expectation of the realized variance under the risk-neutral probability measure: and, consequently, the value of a variance swap at some time t during the contract's life given the initial specific delivery price K can be expressed by Thus, for valuating variance swaps we only need to know E * v R , namely, mean value of the underlying variance under the risk-neutral measure.
In this paper, we are interested in valuing variance swaps for security markets when stochastic volatility σ t, S t has delay response, that is, σ t, S t depends on S t S t−τ , τ > 0, and S t is a stock price at time t ∈ 0, T .
In this way, We need to calculate E * v R to get the fair delivery price K var of a variance swap, and the price of it given the specific delivery price also follows.

Stochastic Volatility with Delay and without Jumps
In this section, we recall the model and approach of pricing variance swaps for stochastic volatility with delay and without jumps presented in the paper of Swishchuk 14 .
In our model, we assume that the price of the underlying asset S t follows the following stochastic delay differential equation SDDE : where μ ∈ R is the mean rate of return, the volatility term σ > 0 is a continuous and bounded function, and W t is a Brownian motion on a probability space Ω, F, P with a filtration F t .We also let r > 0 be the risk-free rate of return of the market.Throughout the paper, we denote S t S t − τ , t > 0 and the initial data of S t is defined by S t ϕ t , where ϕ t is a deterministic function with t ∈ −τ, 0 , τ > 0.

Now we can rewrite the equation of S t as the following expression:
dS t μS t dt σ t, S t S t dW t .

3.2
The asset volatility is defined as the solution of the following equation: where the parameters α, γ, τ, and V are positive constants and 0 < α γ < 1.The Brownian motion W t is the same as the one in 3.2 .We suppose that the following conditions are satisfied, C1 σ t, S t satisfies local Lipschitz and growth conditions; Condition C1 guarantees the existence and uniqueness of a solution of 3.1 and 3.2 see 42 .Condition C2 guarantees the existence of It ô integral in 3.3 , and Condition C3 guarantees the existence of risk-neutral measure P * see below .
To price the variance swaps, we need to calculate E * v E * σ 2 , that is, the expectation of the variance under the risk-neutral probability measure.We should find such probability measure that makes the discounted asset price D t e −rt S t a local martingale and then the risk-neutral pricing formula applies.
Let θ t μ − r /σ t, S t be the market price of risk, which is adapted to the filtration F t ; then by Girsanov's theorem for single Brownian motion, we obtain in the following.
1 There is a probability measure P * equivalent to P such that is its Radon-Nikodym density.
2 The discounted asset price D t is a positive local martingale with respect to P * , and is a standard Brownian motion with respect to P * .

International Journal of Stochastic Analysis
Therefore, in the risk-neutral world, the underlying asset price S t follows the process dS t rS t dt σ t, S t S t dW * t , 3.6 and the asset volatility is defined then as follows: where W * t is defined in 3.5 see 14 for more details .
Let us take the expectations under risk-neutral measure P * on both sides of the equation above.Denoting v t E * σ 2 t, S t , we obtain the following deterministic delay differential equation: We note that this equation has stationary particular solution and approximate general solution.
In this way, we have obtained the following result in 14 .
Theorem 3.1.If conditions (C1)-(C3) are satisfied, then the price of a variance swap at time t given delivery price K in this case should be 1 − e −γT Tγ . 3.9

Pricing Model of Variance Swaps for Stochastic Volatility with Delay and Jumps
In the following section, we will derive some analytical closed formulae for the expectation of the realized variance for stochastic volatility with delay and jumps.First of all, we need to define the jumps and add them to the stochastic volatility model with delay.In 3.3 , we find that the volatility is driven by a Brownian motion W t , but it is more realistic to consider that it is driven by a process which consists a jump part.We represent jumps in stochastic volatility by general compound Poisson processes and try to write the stochastic volatility in the following form: Our purpose is to valuate variance swaps when the stochastic volatility satisfies this general equation.In order to get and check the results, we first consider two simple cases which are easier to model and implement but fundamental and still capture some characteristics of the market.In Section 4.1, we discuss the case that the jump size y t is always equal to constant one, that is, the jump part is represented by t t−τ dN s , just simple Poisson processes.In Section 4.2, the jump part is still compound Poisson processes denoted as t t−τ y s dN s but the jump size y t is assumed to be identically independent distributed random variable with mean value ξ and variance η.The general case is discussed in Section 4.3.We can compare our results with the model in 14 and will see that it is a special case of our model after adding jumps in stochastic volatility.Finally, in Section 4.3 we will show that the model for stochastic volatility with delay and jumps keeps those good features of the model in 14 .

Simple Poisson Process Case
We assume that the price of the underlying asset S t follows the following stochastic delay differential equation SDDE : Recall that our purpose is to calculate E * v E * σ 2 t, S t , the expectation of the variance under the risk-neutral measure.Since we assume that there is no Poisson process in the asset price, the change of measure is not different from the model we discussed in Section 3 to make it risk-neutral.
Let θ t μ − r /σ t, S t be the market price of risk, which is adapted to the filtration F t , and, then by Girsanov's theorem for single Brownian motion, we obtain the risk-neutral probability measure which makes the discounted asset price a local martingale and W * t t 0 θ s ds W t is a Brownian motion under this probability measure.
Note that the change of measure does not change the Poisson intensity λ since it is independent of the Brownian motion.
In the risk-neutral world, the volatility can be defined as follows: Now let us take the expectation under risk-neutral probability P * on both sides of the equation.Note that the Brownian motion and the Poisson process are independent.Letting v t E * σ 2 t, S t , we obtain the following deterministic delay differential equation:

4.6
From this equation, if the intensity of the Poisson process λ 0, then it is the same as 3.8 , the case of stochastic volatility with delay and without jumps.Notice that 4.6 has a stationary solution Hence, the expectation of the realized variance or say the fair delivery price K var of variance swap for stochastic volatility with delay and Poisson jump in stationary regime under risk-neutral measure P * is equal to and the price P of a variance swap at time t given delivery price K in this case should be 4.9 There is no way to write a solution in explicit form for arbitrarily given initial data.But we can understand an approximate behavior of solutions of 4.6 by looking at its eigenvalues.Let us substitute v t X Ce ρt into 4.6 , where X is defined in 4.7 .Then, the characteristic equation for ρ is which is equivalent to when ρ / 0 The only solution to this equation is ρ ≈ −γ, assuming that γ is sufficiently small.Then the behavior of any solution is stable near X, and for large values of t.The approximation has the order O t , which follows from the solution of 4.10 .

International Journal of Stochastic Analysis
Hence, the expectation of the realized variance or say the fair delivery price K var of variance swap for stochastic volatility with delay and Poisson jump in general case under risk-neutral measure P * is equal to where C is given by 4.15 .
Of course, 4.16 can also be written as and the price P of a variance swap at time t given delivery price K in this case should be where X is given by 4.7 .Summarizing, we have the following result.
where X is given by 4.7 .

Compound Poisson Process Case
In the section we will consider the jumps represented by a compound Poisson process, and since it allows the jumps size to be a random number but not always one in Poisson process, the model is more realistic.Our approach in the last section can be easily used in compound Poisson process case.
In the risk-neutral world, the volatility can be defined as follows:

4.23
Hence, the expectation of the realized variance or say the fair delivery price K var of variance swap for stochastic volatility with delay and compound Poisson jump in stationary regime under riskneutral measure P * is equal to and the price P of a variance swap at time t given delivery price K in this case should be

4.25
In general case, we substitute v t X Ce ρt in 4.21 where X is defined in 4.23 .Then the characteristic equation for ρ is which is also the same as that in the last section.

International Journal of Stochastic Analysis
Therefore, the only solution to this equation is ρ ≈ −γ, and by the same method, we have, Hence, the expectation of the realized variance or say the fair delivery price K var of variance swap for stochastic volatility with delay and compound Poisson jump in general case under riskneutral measure P * is equal to where C is given by 4.28 .
Of course, 4.29 can also be written as and the price P of a variance swap at time t given delivery price K in this case should be where X is given by 4.23 .
Remark 4.3.It is interesting to see that when τ 0, which means there is no delay in the model, we have that International Journal of Stochastic Analysis 15

More General Case
In the previous section, we assume that the mean value and variance of the jump size y t in the compound Poisson process are constants.Now we consider a more general case in which they are deterministic functions.The approach used in this section is different from the previous ones, which is a more general method and can be applied to derive the same formulae in the previous simple cases.
In the risk-neutral world, the volatility still satisfies the following equation:

4.38
Now take the expectation under risk-neutral probability, we have that

4.40
We After taking the first derivative of this equation, we obtain This is a second-order delay differential equation with constant coefficients, and so Laplace transform can be applied to find its solution with initial condition v t σ t, S t , t ∈ −τ, 0 , which is already known see 44 or 45 .
Let us denote the Laplace transform of a function f t as and do the Laplace transform for 4.41 : By change of variable and the property of Laplace transform, 4.43 yields

4.44
The characteristic function of 4.41 is Therefore, Applying the inverse transform see 45 , we have that

International Journal of Stochastic Analysis
By the initial condition, Hence, the expectation of the realized variance or say the fair delivery price K var of variance swap for stochastic volatility with delay and compound Poisson jump under risk-neutral measure P * can be obtained by Remark 4.7.From expression 4.49 , we note that only a double integral needs to be solved to get the expectation of the realized variance.If the functions A t , B t , and C s, t are in simple forms, say any integrable functions, the integral can even be solved explicitly.If we only consider the simple case in Sections 4.1 or 4.2, the final analytical formulae only need algebraic computation and no numerical approximation is needed.Thus, we can say that the stochastic volatility model with delay and jumps is easy to implement and time saving.

Delay as a Measure of Risk
Volatility is the standard deviation of the change in value of a financial instrument with specific time horizon.It is often used to quantify the risk of the instrument over that time period.The higher the volatility is, the riskier the security is.The variance is a square of volatility and is also a measure of risk of a financial instrument.In this section we look at the estimation for variance through the delay to minimize the risk.The compound Poisson case in Section 4.2 is a special case of the model we discussed in Section 4.3, but it is more basic and also shows some important characteristics.Therefore, in this section, we focus on this case.By Section 4.2, we obtain 1 − e −γT γT .
To reduce the risk we need to take into account the following relationship with respect to the delay τ which follows from 5.1 :

5.3
In this case, there is a way to control the delay.

Numerical Example
In this section, we apply the analytical solutions from Section 4.2 to price the variance swap of the S&P 60 Canada index for five years January 1998-February 2002 see 46 .
In the end of February 2002, we wanted to price the fixed leg of a variance swap based on the S&P 60 Canada index.The statistics on log-returns of S&P 60 Canada index for 5 years January 1997-February 2002 is presented in Table 1.
From the histogram of the S&P 60 Canada index log-returns on a 5-year historical period 1,300 observations from January 1998 to February 2002 leptokurtosis may be seen in the histogram.If we take a look at the graph of the S&P 60 Canada index log-returns on a 5-year historical period we may see volatility clustering in the returns series.These facts indicate the conditional heteroscedasticity.
There are several parameters which need to be estimated from the data, the jump intensity λ, the mean value ξ, and the variance η of the jump size.We use the following method to detect the jumps.If the difference between the spot log-return and the mean of our log-return series is larger than triple of the standard deviation, then we say that a jump occurs at that time point see 47 .Clewlow and Strickland 47 estimated jumps in asset price.In our case, the volatility is a trading instrument and can be considered as an asset.In this way, we can apply here their methodology.We count the number of the jumps, denoted as N, and the estimation of jump intensity λ N/1300.The jump size is defined as the difference between the log-returns at a jump time point and the previous time point, and the sample mean and variance of these data are unbiased estimation of ξ and η.
A GARCH 1,1 regression is applied to the series, and the results are obtained as in  Repeating this approach for a series of delays up to 5 days and series of maturities up to 100 years, we compare our results with the no jump model and obtain the following plot see the appendix, Figures 1 and 2 of S&P 60 Canada index variance swap.Figure 3

An Important Remark with Respect to Our Model with Jumps
We note that we provided back testing in 48 to show that estimated variance swap prices are close to realized prices for stochastic volatility without jumps.Figure 1 shows that, in presence of jumps, the estimated delivery prices are close for models with and without jumps the order of estimation is 10 −4 .This marginal difference between the delivery prices E * v with and without jumps is because our λ is very small since we calculated the λ 0.0115 with respect to the number of jumps, 15, divided by the number of samples, 1300.We count one jump if the price change is bigger than 3 times the standard deviation.For example,   without jumps, E * v 0.0001124, and, with jumps λ 0.0115 , E * v 0.00011365 our case .In this way, the relative change with respect to no jumps is 1.07%.If we increase the jump intensity to, for example, λ 0.1, then E * v 0.0001233 see 4.29 for calculation , and the relative change with respect to no jumps is 9.7%.For λ 0.2 the relative change becomes 19.4%, and for λ 0.5 the relative change would be 48.6%!All in all, in our case, with λ 0.0115, even 1% change cannot be ignored in the market, because our value for E * v is just for one unit of expected variance.But if one buys many contracts and times a notional amount N, then 1% change is not that small in value anymore.

Conclusion
In this paper we studied stochastic volatility model with delay and jumps to price variance swaps.We applied a general approach to derive the analytical close forms for expectation of the realized continuously sampled variance for stochastic volatility with delay and jumps.The jump part in our model is represented by a general version of compound Poisson processes.The key features of the model are the following: i continuous-time analogue of discrete-time GARCH model, ii mean reversion, iii containing the same source of randomness as stock price, iv completeness of the market, v incorporating the expectation of log-return, and vi incorporating the jumps in volatility.The model is also easy to implement and time saving.Besides, we presented a lower bound for delay as a measure of risk.From the numerical example, we found that after adding jumps in volatility, the expectation of the realized variance is higher than the one without jumps for variance swaps.It is easy to explain it since the existence of jumps means the market is more risky which asks for higher cost of variance swaps.
For further study, we may add jumps in the spot price of the underlying asset which is also an important characteristic found in real markets.Besides, the jump part can be in more complicated form, for example, t t−τ σ s, S s dN s .We may also consider applying our methods to price volatility swaps which make more sense for investors to hedge volatility.Also, one of the referees suggested considering a continuum of delay times, not only just one delay time, and an associated "memory function" to weight their contributions which would make more sense for introduced model.However, all these topics are very interesting but are beyond the scope of the present paper and will be considered in our future research work.

Appendix Figures
Figure 1 depicts the dependence of delivery price on maturity for S&P 60 Canada index 1998-2002 for stochastic volatility SV with and without jumps.As we can see the delivery price for SV with jumps is higher than that for SV without jumps because the security in this case is riskier.
Figure 2 depicts the dependence of delivery price on delay for S&P 60 Canada index 1998-2002 for SV with and without jumps.Again, the price is obviously higher.We note that, if the delay is zero, then the price approaches the constant value when maturity T is large i.e., see formula 4.36 .If the delay is not zero, then the delivery price is proportional to the length of delay.
Figure 3 depicts the dependence of delivery price on jump intensity for S&P 60 Canada index 1998-2002 for SV with jumps.The higher the intensity of jumps is, the higher the delivery price is.
Figure 4 depicts the dependence of delivery price on delay and jump intensity for S&P 60 Canada index 1998-2002 .This figure is a combination of Figures 2 and 3 for SV with jumps and delay.
Figure 5 depicts the dependence of delivery price on delay and maturity for S&P 60 Canada index 1998-2002 .This figure is a combination of Figures 1 and 2 for SV with jumps and delay.
Figure 6 depicts the dependence of delivery price on jump intensity and maturity for S&P 60 Canada index 1998-2002 .This figure is a combination of Figures 1 and 3 for SV with jumps and delay.

Theorem 4 . 1 .
Consider stock price satisfying 4.2 with stochastic volatility in 4.3 .Under conditions (C1)-(C3) (Section 3) the price of the variance swap at time t given delivery price K is C s, u − C s, u − τ ds du and G 0 −τ 0 −τ C s, u ds du.Taking into account 4.36 , 4.37 , and 4.39 , equation 4.35 becomes

Figure 1 :Figure 2 :
Figure 1: Dependence of delivery price on maturity S&P 60 Canada index .

Figure 3 :
Figure 3: Dependence of delivery price on jump intensity S&P 60 Canada index .

Figure 4 :
Figure 4: Dependence of delivery price on delay and jump intensity S&P 60 Canada index .
u r i t y Dependence of delivery price on delay and maturity D el a y

Figure 5 :
Figure 5: Dependence of delivery price on delay and maturity S&P 60 Canada index .

Figure 6 :
Figure 6: Dependence of delivery price on jump intensity and maturity S&P 60 Canada index .
Brownian motion, N t is a Poisson process with intensity λ, and y t is the jump size at time t.We assume that E y t A t , E y s y t C s, t , s < t, and E y 2 t B t C t, t and A t , B t , and C s, t are all deterministic functions.
Our first step is still taking the expectation under risk-neutral probability P * on both sides of the equation.Note that the Brownian motion and the compound Poisson process are independent.Letting v t E * σ 2 t, S t , we obtain the following equation: * t is a Brownian motion, N t is a Poisson process with intensity λ, and y t is the jump size at time t which is identically independent normally distributed random variable.We assume that the mean of y t is ξ and the variance of y t is η.Note that the Poisson intensity λ and the jump size y t do not change in risk-neutral world since they are independent of the Brownian motion.
, where A t , B t , and C s, t are all deterministic functions.Note that the change of measure does not change the Poisson intensity λ and the distribution of jump size y t , since they are independent of the Brownian motion.Let v t E * σ 2 t, S t and take the expectation under risk-neutral probability P * on both sides of 4.34 .Noting that the Brownian motion and the Poisson process are independent, we obtain the following equation: * t is a Brownian motion, N t is a Poisson process with intensity λ, and y t is the jump size at time t.We assume that E y t A t , E y s y t C s, t , s < t, and E y 2 To get the expectation of the realized variance in the risk-neutral world E * v , we have to find a solution to 4.40 , a nonhomogeneous integrodifferential equation with delay.

Table 1 :
Statistics on log-returns of S&P 60 Canada index.

Table 2 .
This table allows to generate different input variables to the volatility swap model.
r 0.02 and τ 1 day .Parameter v may be found from the expression V C/ 1−α−β , where C 2.58×10 −6 is defined in

Table 2 .
Thus, V 0.00020991 and see the appendix depicts the dependence of delivery price of variance swaps on jump intensity.Figures4, 5, and 6 see the appendix depict the dependence of delivery price of variance swaps on two variables of maturity, delay, and jump intensity.