代写 R 1.

1.
STAD70 Assignment 6
Antithetic variables are closely related to stratification with only 2 strata, i.e. with positive & negative values for Z . Their only difference is that for antithetic variables you have the same absolute values for positive and negative variables, whereas with stratification you can have different absolute values for positive and negative variables.
a.
b. c.
For what type of function f does the use of antithetic variables completely eliminate
the variability of the Monte-Carlo estimator, i.e. Var Y   1 Var Yi  Yi   0 AVn 2
whereYf(Z),Yf(Z)andZ~iid N(0,1) iiiii
2n
For the type of function you identified in part a., show that the variance of the
stratified estimator (with 2 equiprobable strata of size n each) is also not 0 and that it is smaller than that from part b., i.e. 0  Var[YStr ]  Var[Y ] .

,
For the type of function you identified in part a., show that the variance of the simpleMonteCarloestimatorisnotnecessarily0,i.e.Var[Y] 1VarYi0.
2.
A European gap option has a strike price K1 and a trigger price K2 . The trigger price determines whether or not the gap option is exercised, while the strike price determines
the amount of the payoff. For example, the European gap call option has payoff given by: ST  K1, if ST  K2 , presented in the following plots:
 0,ifSTK2
Obviously, if the strike price is equal to the trigger price then the gap option is an ordinary option. Note that gap options can result in negative payoff at expiration (so they should not really be called “options”). Moreover, gap options are different from barrier options, in that for barrier options the barrier can be reached prior to expiration, whereas for gap options only the final stock price is compared to the trigger.
a. Estimate the price of a European gap call option with S0  55, K1  50, K2  60, T 1yr, r  2%,   20% and assuming the price follows GBM. Use n 100,000
paths and calculate the standard error of your result. Compare that to the exact
price given by the R function GapOption in the package fExoticOptions.
b. Repeat part a. using antithetic variables, and report the new price estimate and
standard error.

3.
c. d.
Repeat part a. using stratification with m  5,10, & 20 equiprobable strata, and report the new price estimate and standard error.
Would you expect the variance reduction methods in part b. and c. to work better
for in-the-money, at-the-money, or out-of-the-money options? Justify your answer. A forward start option is an option that starts at a specified future date with an
expiration date set even further in the future. For example, consider a forward start call option with start date in 1yr and expiration date in 2yrs: at t  0 you pay the option premium, and at t  1 you receive a European call option with expiration date t  2 and strike price K  S1 (this is called an at-the-money forward start call, because the strike
price is equal to the asset price at t 1, which is unknown at t  0).
a.
Estimate the price of an at-the-money forward start European call option with
S0 70,Tstart 1, Texp 2yr,r2%, 20% andassumingthepricefollowsGBM.
Use n  10, 000 paths and calculate the standard error of your result. You should use the Black-Scholes formula to explicitly calculate the value of the option at Tstart , i.e. you should generate STStart and then calculate its “payoff” as the Black-Scholes price
of a European at-the-money call option with 1 year to maturity (use BSprice function provided in Lecture 8’s R code). Compare that to the exact price given by the R function ForwardStartOption in the package fExoticOptions. Repeat part a. using antithetic variables, and report the new price estimate and standard error.
Repeat part a. using stratification with m  10 & n  1, 000 equiprobable strata, and report the new price estimate and standard error.
Show that the price of the at-the-money forward call option is actually a simple linear function of S0 (Hint: use risk-neutral pricing and write the payoff of the option
at Tstart as the Black-Scholes price of an at-the-money call with maturity at Texp )
An α-quantile option is an option whose payoff is determined by the αth quantile of
b. c.
d.
4.
the underlying price from time 0 until expiration. For example, the median call option
is the
(i.e. α=50%) with strike K has payoff
TT median S   K , where median S 

t
median asset price over [0,T]. Assume that the underlying asset follows GBM with
S0  60, r  5%,   25% , and use simulation to estimate the price of median call with K 70&T 1.
a. Perform simple Monte Carlo using path discretization with t j T
estimate.
b. Repeat part a. using the Black-Scholes price of the European call as a control variate,
and report the new price estimate and standard error.
t
t0  t0
 j T , j  0, , m m
m  50 in order to approximate the continuous median St  t0
by the discrete medianS(tj)j0 .Use n10,000pathsandreportthestandarderrorofthe
m
and

5.
An Asian option is an option whose payoff is determined by the average underlying price over some pre-set period of time, typically the time till expiration. For example, theAsiancalloptionwithstrike K haspayoff A(0,T)K ,where A(0,t) isthe
average asset price over (0,T) . There are two types of averaging: 1T
arithmeticaverage:A(0,T)T 0 S(t)dt
1T geometricaverage:A(0,T)exp  lnS(t)dt
T 0
If the underlying asset follows GBM, then the price of an Asian option with geometric averaging can be calculated explicitly. For Asian options with arithmetic averaging, however, we rely exclusively on numerical methods for pricing.
Assume that the underlying asset follows GBM with S0  70, r  5%,   25% , and use
simulation to estimate the price of an arithmetic Asian call with K  80 & T  1 . a. Perform simple Monte Carlo using path discretization with m  50 in order to
approximateA(0,T)1TS(t)dtbyA(0,T)1m S(ti).Usen10,000paths T0 m j1
and report the standard error of the estimate.
b. Repeat part a. using the price of the geometric Asian call as a control variate, and
report the new price estimate and standard error. For implementing the control variate, you can get the exact price of the geometric Asian option using the R function GeometricAverageRateOption in the package fExoticOptions.
c. Repeat part a. using importance sampling. Change the probability measure so that E[ST ]  K , and experiment with at least two more measure changes of your choice.
Report the new price estimates and standard errors.