Improved Variance Reduced Monte-Carlo Simulation of in-the-Money Options

Pricing derivatives with Monte-Carlo simulations involve standard errors that typically decrease at a rate proportional to N 1 where N is the sample size. Several approaches have been discussed to reduce the empirical variance for a given sample size. This article analyzes the joint application of the put-call-parity approach and importance sampling to variance reduced option pricing. For this purpose, we examine non-path-dependent and path-dependent options. For European options, we observe dramatic variance reduction, especially for in-the-money options. Also for arithmetic Asian options, a significant variance reduction is achieved.


Introduction
Monte-Carlo simulations are frequently used to estimate prices of financial products for which no analytical formulae exist, e.g., several path-dependent options.A difficulty related to this approach is that the empirical variance of estimators decreases slowly, typically at a rate 1 N ∝ where N is the sample size.Several techniques have been discussed to reduce the empirical variance for a given sample size N [1].Importance sampling turns out to be a particularly effective variance reduction technique [2]- [5].
An alternative approach to variance reduction at least for in-the-money options is the application of the put-call-parity: instead of simulating an in-the-money call price, the corresponding out-of-the-money put price can be simulated.Then, the call price can be calculated from the put-call-parity yielding a variance reduced estimator [6] [7].This analysis will investigate how the joint application of the put-call-parity and importance sampling can lead to significant variance reduction in the simulation of Monte-Carlo estimators.The article is organized as follows: in Section 2 the method of importance sampling is introduced.Section 3 explains how put-call-parities can be employed for variance reduction.Next, the combination of the two approaches is introduced in Section 4. Numerical results of the achieved variance reductions are presented in Section 5 for both non-path-dependent and path-dependent options.A conclusion follows in Section 6.

Importance Sampling
Foundations of Importance Sampling This section introduces the concept of importance sampling, which serves to reduce the empirical variance of estimators.
The expectation value of a non-negative function of a random variable X with probability density p is calculated as An estimator of α with N i.i.d.realizations 1 , , N X X  of the random variable X is ( ) With any other probability density p′ , α can be rewritten as The fraction p p′ is usually called "Radon-Nikodým derivative".An unbiased estimator of this expression is Its variance vanishes if we choose To obtain a a new probability density this product must be normalized (note that 0 h ≥ by definition): For this purpose, the integral must be calculated, which is the original problem to be solved in Equation (1).However, an approximation of the optimal density can already lead to significant variance reduction [1].
Monte-Carlo simulation of the Feynman-Kac formula Option prices can be calculated as the expectation value of the derivative's discounted pay-off function.In the following, a variance reduced estimator of the resulting Feynman-Kac formula will be presented.Note that the Feynman-Kac formula also covers the case of path-dependent options, i.e. the case where h does not only depend on the terminal value of the underlying's price T S but also on other values A. Müller t S with 0 t T < < .Instead of writing t S in Equation ( 9), one could write t X and define t X as a bivariate random process whose first element is the underlying price t S itself and the second element is the integral of t S with respect to time.For simplicity, a constant risk-free interest rate r is assumed.For the purpose of numerical calculations, the Feynman-Kac formula can be expressed as using a finite-dimensional approximation on a grid k t k τ τ = + ∆ with Euler densities ( ) and n discretization steps (for details see Appendix B in [5]).We set n T τ = and consequently n T S S = .
In accordance with Equation ( 2), a simple Monte-Carlo estimator of this expression is Importance sampling by adding an additional drift term In 2014, Singer presented a variance reduced Monte-Carlo estimator for Equation ( 9) which relies on importance sampling by adding an additional drift term to the stochastic differential equation where ( ) W t a Wiener process.Here we just present the result, while the lenghty derivation can be found in [5].Melchior and Öttinger have introduced a similar approach in 1995 [2].
The optimal importance sampling density opt p p ′ = is chosen as indicated in Equation ( 7).It can be shown that this density is equivalent to a new stochastic differential equation with additional drift term of the form In the Black-Scholes model with f rS = and g S σ = , one obtains ( ) using the abbreviation for the option price elasticity.
Here, C is the Feynman-Kac formula from Equation (9).Again, the problem of Equations ( 7) and ( 8) materializes: To describe the optimal stochastic differential equation, knowledge of the unknown quantity C is required.
However, variance reduction can also be achieved by approximating C. For this purpose, in our analysis we use the Black-Scholes formula for European options [8].This formula is the best possible choice when simulating European options in the Black-Scholes model, as it yields the analytically exact value.But also for other models of the underlying price (e.g., the CEV model or models with stochastic volatility) or for other options (e.g., Asian options) it represents a useful approximation.
In the univariate Black-Scholes case that is considered in this article, the Radon-Nikodým derivative for one step on the discretized grid yields Ultimately, for the variance reduced Monte-Carlo estimator with sample size N, n discretization steps and sampling from opt p one obtains As mentioned, ik  can be calculated from the Black-Scholes formula [8].Depending on whether h represents the pay-off function of a call option or a put option, one obtains [9] ( ) ( )

Variance Reduction with Put-Call-Parities for European and Arithmetic Asian Options
European options Simple non-arbitrage arguments require that for a European call with pay-off function ( ) the following put-call-parity holds true where 0 C is the initial call price at 0 t = and 0 P the put price [10]: the following relation is required to avoid arbitrage [7]: Variance reduction Reider (1994) [6] suggested to apply put-call-parities to variance reduced importance sampling.In-the-money call/put options can be estimated more efficiently by first estimating the corresponding put/call option with the same parameters and subsequently calculating the required call/put price from the put-call-parity.Conversely, if the approach is applied to out-of-the-money options, the empirical variance is generally increased.
Put-call-parities allow to split a stochastic quantity into deterministic components and a residual stochastic quantity.The deterministic components' standard errors vanish by definition.If the residual stochastic quantity is reduced in absolute size compared to the stochastic quantity of interest, the residual stochastic quantity's standard error will generally tend to be lower as well, all else being equal.According to Gaussian propagation of uncertainty, variance reduction is achieved [7].

Joint Application of the Put-Call-Parity and Importance Sampling
In a previous article, it is suggested that a combined application of the put-call-parity and importance sampling might be particularly attractive [7].The author discusses that importance sampling in many cases is especially attractive for out-of-the-money options.Therefore, the valuation of an in-the-money call/put option could be conducted more efficiently by pricing the corresponding out-of-the-money put/call option with importance sampling and then calculating the desired option price from the put-call-parity.More specifically, the following procedure to price call options is examined numerically in the next section: 1) To obtain a benchmark against which variance reductions can be analyzed, a simple Monte-Carlo simulation of the estimator (11) is conducted.Both the estimator and its empirical variance are calculated.
2) To simulate the option price of a call option with a given set of parameters, the corresponding put option with the same set of parameters is simulated using the variance reduced Monte-Carlo estimator (17). is calculated using the Black-Scholes formula for put options as described in Equation (18).
3) The call price is calculated from the put-call-parity (20) resp.( 23).The empirical variance of this estimated call price is the same as the empirical variance of the previously simulated put price.This is the case as all other terms in Equation (20) resp.( 23) have zero variance.
4) The variance reduction ratio ("VarRatio") is calculated as the ratio between the empirical variance of the benchmark estimator and the empirical variance of the variance reduced estimator calculated employing importance sampling and the put-call-parity.

Numerical Results
In this section, we apply the introduced combined approach of importance sampling and employing the putcall-parity to simulate variance reduced Monte-Carlo estimators.To make results comparable, we present two examples that have served previously in the analysis of other approaches discussed in literature [11].More specifically, we simulate prices of European call options and of arithmetic Asian call options and evaluate achieved variance reductions.
According to the Black-Scholes model, the underlying t S is simulated as follows: ( ) ( ) We use a discretized time-grid k t k τ τ = + ∆ with n equidistant steps.For benchmark simulations without importance sampling 0 k =  follows for all k.For simulations involving importance sampling we choose k  as described in Section 2. European options First, we consider the example of a European call option with pay-off function ( ) Column "PCP-IS method" of Table 1 presents the results of the approach described in section 4. For comparison, results of two other approaches previously discussed in literature are presented as well in the columns "GHS method" and "ZLG method".The first method was presented by Glasserman, Heidelberger and Shahabuddin in 1999 [4] and the second by Zhao, Liu and Gu in 2013 [11].The numerical results for these two methods are taken from [11].
The approach presented in this article achieves higher variance reductions for all considered parameter constellations compared to the other two approaches.Especially for (deep-)in-the-money options ( 30 K = ), we observe dramatic variance reductions of several orders of magnitude.For in-the-money options, the approach turns out to be more effective for low-volatility conditions.
Two further aspects should be noted: First, in the case of out-of-the money options a direct application of the introduced importance sampling procedure would be more appropriate than the joint application employing the put-call-parity.E.g., when foregoing the put-call-parity for 55 K = , we obtain variance reduction ratios of 59.3 .This is not surprising, as the put-call-parity approach increases empirical variance for out-of-the-money options [7].
Second, the combination of the put-call-parity approach with importance sampling outperforms the two standalone approaches (put-call-parity and importance sampling) for in-the-money options.The standalone importance sampling approach in the case of deep-in-the money calls only leads to limited variance reduction.For the standalone put-call-parity approach applied to the low-volatility case 0.1 σ = , we observe that not even one of the trajectories.To approximate the option price elasticity  , the Black-Scholes formula for European options with corresponding parameters has been used.Again, we truncate  .However, in the case of arithmetic Asian options we observe best results for a truncation of its absolute value to the interval [ ] 1; 2.5 .In the same structure as before, Table 2 shows that we observe the best results for in-the-money options.For low-volatility conditions, we achieve variance reductions of several orders of magnitude.The presented approach outperforms the other methods for in-the-money options.For out-of-the-money options, we achieve no or low variance reduction by employing the combined approach.Again, better results could be achieved by foregoing the put-call-parity.
For European options, we achieved higher variance reduction than for arithmetic Asian options.This is not surprising, as we calculated the option price elasticity  from the Black-Scholes formula for European options in both examples.Thus, we approximated the optimal importance sampling density more accurately for European options and consequently a higher variance reduction should be expected.More accurate estimations of the option price elasticity for arithmetic Asian options might lead to increased variance reduction ratios.
Again, the numerical results for the GHS method and the ZLG method are taken from [11].

Conclusions
We conducted variance reduced Monte-Carlo simulations for European and arithmetic Asian options.For inthe-money options, we achieve significant variance reduction by combining importance sampling with the putcall-parity approach.The presented approach delivers best results in low-volatility environments.For European options, it outperforms previously published results for all parameter constellations.Also for arithmetic Asian options, we obtain higher variance reduction ratios for in-the-money options.This analysis does not consider the combination with other variance reduction techniques, e.g., antithetic variates or stratification, which could lead to further increased variance reduction ratios.
Arithmetic Asian options Also for arithmetic Asian call options with the pay-off function reasons of numerical stability, we truncate the absolute value of the option price elasticity  to the interval simulated yields a positive contribution to the Monte-Carlo put option estimator.The result is an empirical variance of zero.Consequently, no reasonable variance reduction can be determined.Arithmetic Asian options The second example examines an arithmetic Asian option with pay-off function

Table 1 .
Estimated variance reduction ratios for European call option.

Table 2 .
Estimated variance reduction ratios for arithmetic Asian call option.