<< стр. 2(всего 2)СОДЕРЖАНИЕ
.
2П†0 (2y
в€’ x)

Inverting this relationship between (x, y) and (H, C),

P [H в€€ в€†H, C в€€ в€†C] = в€’2П†0 (2y в€’ x)в€†xв€†y

conп¬Ѓrming that the joint density of (H, C) is given by в€’2П†0 (2y в€’ x) for x < y.
In order to get the joint density of the High and the Close when the drift is
non-zero, we need only multiply by the ratio of the two normal density functions
of the close
fВµ (x)
f0 (x)
SIMULATING BARRIER AND LOOKBACK OPTIONS 295

Figure 5.8: Conп¬Ѓrmation of the joint density of (H, C)

and this gives the more general result in the table below.
The table below summarizes many of our distributional results for a Brown-
ian motion process with drift on the interval [0, 1],

dSt = Вµdt + ПѓdWt , with S0 = O.

Statistic Density Conditions
в€’в€ћ < x < y,
X = C в€’ O,
f (y, x) = в€’2П†0 (2y в€’ x) exp(Вµx в€’ Вµ2 /2) and y > 0, Пѓ = 1
Y =H в€’O
given O
fY |X (y|x) = 2(2y в€’ x)eв€’2y(yв€’x)
Y |X y > x, Пѓ = 1
exp((Пѓ 2 /2)
Z = Y (Y в€’ X) given O, X
exp(Пѓ2 /2)
(L в€’ O)(L в€’ C) given (O, C)
exp(Пѓ2 /2)
(H в€’ O)(H в€’ C) given (O, C)
Hв€’O
drift ОЅ = 0, given O, 2H в€’ O в€’ C
U [0, 1]
2Hв€’Oв€’C
Lв€’O
drift ОЅ = 0, given O, 2L в€’ O в€’ C
U [0, 1]
2Lв€’Oв€’C
Cв€’O
drift ОЅ = 0, given H, O
U [в€’1, 1]
2Hв€’Oв€’C
TABLE 5.1: Some distributional results for High, Close and Low.
296 CHAPTER 5. SIMULATING THE VALUE OF OPTIONS

We now consider brieп¬‚y the case of non-zero drift for a geometric Brownian
motion. Fortunately, all that needs to be changed in the results above is the
marginal distribution of ln(C) since all conditional distributions given the value
of C are the same as in the zero-drift case. Suppose an option has payoп¬Ђ on
maturity П€(C) if an upper barrier at level Oeb , b > 0 is not breached. We have
already seen that to accommodate the п¬Ѓletering eп¬Ђect of this knock-out barrier
we should determine, numerically or by simulation, the expected value

b(b + ln(O/C))
E[П€(C)(1 в€’ exp{в€’2 })]
Пѓ2 T

the expectation conditional (as always) on the value of the open O. The eп¬Ђect
of a knock-out lower barrier at Oeв€’a is essentially the same but with b replaced
by a, namely
a(a + ln(C/O))
E[П€(C)(1 в€’ exp{в€’2 })].
Пѓ2T
In the next section we consider the eп¬Ђect of two barriers, both an upper and a
lower barrier.

One Process, Two barriers.

We have discussed a simple device above for generating jointly the high and the
close or the low and the close of a process given the value of the open. The joint
distribution of H, L, C given the value of O or the distribution of C in the case
of upper and lower barriers is more problematic. Consider a single factor model
and two barriers- an upper and a lower barrier. Note that the high and the
low in any given interval is dependent, but if we simulate a path in relatively
short segments, by п¬Ѓrst generating n increments and then generating the highs
and lows within each increment, then there is an extremely low probability
that the high and low of the process will both lie in the same short increment.
For example for a Brownian motion with the time interval partitioned into 5
equal subintervals, the probability that the high and low both occur in the
SIMULATING BARRIER AND LOOKBACK OPTIONS 297

same increment is less than around 0.011 whatever the drift. If we increase the
number of subintervals to 10, this is around 0.0008. This indicates that provided
we are willing to simulate highs, lows and close in ten subintervals, pretending
that within subintervals the highs and lows are conditionally independent, the
error in our approximation is very small.

An alternative, more computationally intensive, is to diп¬Ђerentiate the inп¬Ѓnite
series expression for the probability P (H В· b, L в‰Ґ a, C = u|O = 0] A п¬Ѓrst step
in this direction is the the following result, obtained from the reп¬‚ection principle
with two barriers.

Theorem 47 For a Brownian motion process

dSt = Вµdt + dWt , S0 = 0

deп¬Ѓned on [0, 1] and for в€’a < u < b,

P (L < в€’a or H > b|C = u)
в€ћ
1X
[П†{2n(a + b) + u} + П†{2n(a + b) в€’ 2a в€’ u}
=
П†(u) n=1

+ П†{в€’2n(b + a) + u} + П†{2n(b + a) + 2a + u}]

where П† is the N (0, 1) probability density function.

Proof. The proof is a well-known application of the reп¬‚ection principal.
It is suп¬ѓcient to prove the result in the case Вµ = 0 since the conditional
distribution of L, H given C does not depend on Вµ (A statistician would say
that C is a suп¬ѓcient statistic for the drift parameter). Denote the following
paths determined by their behaviour on 0 < t < 1. All paths are assumed to
end at C = u.
298 CHAPTER 5. SIMULATING THE VALUE OF OPTIONS

= H > b (path goes above b)
A+1
path goes above b and then falls below в€’a
=
A+2
goes above b then falls below в€’a then rises above b
=
A+3
etc.
L < в€’a
=
Aв€’1
path falls below в€’a then rises above b
=
Aв€’2
falls below в€’a then rises above b then falls below в€’a
=
Aв€’3
etc.
For an arbitrary event A, denote by P (A|u) probability of the event conditional
on C = u. Then according to the reп¬‚ection principal the probability that the
Brownian motion leaves the interval [в€’a, b] is given from an inclusion-exclusion
argument by

P (A+1 |u) в€’ P (A+2 |u) + P (A+3 |u) в€’ В· В· В· (5.30)

+P (Aв€’1 |u) в€’ P (Aв€’2 |u) + P (Aв€’3 |u) В· В· В·

This can be veriп¬Ѓed by considering the paths in Figure 5.9. (It should be noted
here that, as in our application of the reп¬‚ection principle in the one-barrier case,
the reп¬‚ection principle allows us to show that the number of paths in two sets is
the same, and this really only translates to probability in the case of a discrete
sample space, for example a simple random walk that jumps up or down by a
п¬Ѓxed amount in discrete time steps. This result for Brownian motion obtains if
we take a limit over a sequence of simple random walks approaching a Brownian
motion process.)
Note that

П†(2b в€’ u)
P (A+1 |u) =
П†(u)
П†{2n(a + b) + u}
P (A+2n |u) =
П†(u)
П†{2n(a + b) в€’ 2a в€’ u}
P (A+(2nв€’1) |u) =
П†(u)
SIMULATING BARRIER AND LOOKBACK OPTIONS 299

Figure 5.9: The Reп¬‚ection principle with Two Barriers

and
П†(в€’2a в€’ u)
P (Aв€’1 |u) =
П†(u)
П†{в€’2n(b + a) + u}
P (Aв€’2n |u) =
П†(u)
П†{2n(b + a) + 2a + u}
P (Aв€’(2n+1) |u) = .
П†(u)
The result then obtains from substitution in (5.30).

As a consequence of this result we can obtain an expression for P (a < L В·
H < b, u < C < v) (see also Billingsley, (1968), p. 79) for a Brownian motion
on [0, 1] with zero drift:

P (a, b, u, v) = P (a < L В· H < b, u < C < v)
в€ћ
X
О¦[v + 2k(b в€’ a)] в€’ О¦[u + 2k(b в€’ a)]
=
k=в€’в€ћ
в€ћ
X
в€’ О¦[2b в€’ u + 2k(b в€’ a)] в€’ О¦[2b в€’ v + 2k(b в€’ a)]. (5.31)
k=в€’в€ћ
300 CHAPTER 5. SIMULATING THE VALUE OF OPTIONS

where О¦ is the standard normal cumulative distribution function. From (5.31)
we derive the joint density of (L, H, C) by taking the limit P (a, b, u, u + Оґ)/Оґ as
Оґ в†’ 0, and taking partial derivatives with respect to a and b:
в€ћ
X
k 2 П†00 [u + 2k(b в€’ a)] в€’ k(1 + k)П†00 [2b в€’ u + 2k(b в€’ a)]
f (a, b, u) = 4
k=в€’в€ћ
в€ћ
X
k П† [u + 2k(b в€’ a)] в€’ k(1 + k)П†00 [2b в€’ u + 2k(b в€’ a)]
2 00
=4
k=1

+ k2 П†00 [u в€’ 2k(b в€’ a)] + k(1 в€’ k)П†00 [2b в€’ u в€’ 2k(b в€’ a)] (5.32)

for a < u < b.
From this it is easy to see that the conditional cumulative distribution func-
tion of L given C = u, H = b is given by on a В· u В· b (where в€’2П†0 (2b в€’ u)
is the joint p.d.f. of H, C) by
в€‚2
в€‚bв€‚v P (a, b, u, v)| v=u
(5.33)
F (a|b, u) = 1 +
2П†0 (2b в€’ u)
в€ћ
X
в€’1
{в€’kП†0 [u + 2k(b в€’ a)] + (1 + k)П†0 [2b в€’ u + 2k(b в€’ a)]
=0
П† (2b в€’ u)
k=1

+ kП†0 [u в€’ 2k(b в€’ a)] + (1 в€’ k)П†0 [2b в€’ u в€’ 2k(b в€’ a)]}

This allows us to simulate both the high and the low, given the open and the
close by п¬Ѓrst simulating the high and the close using в€’2П†0 (2b в€’ u) as the joint
p.d.f. of (H, C) and then simulating the low by inverse transform from the
cumulative distribution function of the form (5.33).

Survivorship Bias

It is quite common for retrospective studies in п¬Ѓnance, medicine and to be
subject to what is often called вЂњsurvivorship biasвЂќ. This is a bias due to the
fact that only those members of a population that remained in a given class
(for example the survivors) remain in the sampling frame for the duration of
the study. In general, if we ignore the вЂњdrop-outsвЂќ from the study, we do so
SURVIVORSHIP BIAS 301

at risk of introducing substantial bias in our conclusions, and this bias is the
survivorship bias.

Suppose for example we have hired a stable of portfolio managers for a large
pension plan. These managers have a responsibility for a given portfolio over
a period of time during which their performance is essentially under continuous
review and they are subject to one of several possible decisions. If returns below
a given threshhold, they are deemed unsatisfactory and п¬Ѓred or converted to
another line of work. Those with exemplary performance are promoted, usually
to an administrative position with little direct п¬Ѓnancial management. And those
between these two вЂњabsorbingвЂќ barriers are retained. After a period of time,
T, an amibitious graduate of an unnamed Ivey league school working out of
head oп¬ѓce wishes to compare performance of those still employed managing
portfolios. How are should the performance measures reп¬‚ect the п¬Ѓltering of
those with unusually good or unusually bad performance? This is an example
of a process with upper and lower absorbing barriers, and it is quite likely
that the actual value of these barriers diп¬Ђers from one employee to another, for
example the son-in-law of the CEO has a substantially diп¬Ђerent barriers than the
math graduate fresh out of UW. However, let us ignore this diп¬Ђerence, at least
for the present, and concentrate on a diп¬Ђerence that is much harder to ignore
in the real world, the diп¬Ђerence between the volatility parameters of portfolios,
possibly in diп¬Ђerent sectors of the market, controlled by diп¬Ђerent managers.
For example suppose two managers were responsible for funds that began and
ended the year at the same level and had approximately the same value for the
lower barrier as in Table 5.2. For each the value of the volatility parameter
Пѓ was estimated using individual historical volatilities and correlations of the
component investments.
302 CHAPTER 5. SIMULATING THE VALUE OF OPTIONS

Portfolio Open price Close Price Lower Barrier Volatility

56 5
1 40 30 .5
8

56 1
2 40 30 .2
4
Table 5.2

Suppose these portfolios (or their managers) have been selected retrospec-
tively from a list of вЂњsurvivorsвЂќ which is such that the low of the portfolio value
never crossed a barrier at l = Oeв€’a (bankruptcy of fund or termination or
demotion of manager, for example) and the high never crossed an upper barrier
at h = Oeb . However, for the moment let us assume that the upper barrier is so
high that its inп¬‚uence can be neglected, so that the only absorbtion with any
substantial probability is at the lower barrier. We interested in the estimate of
return from the two portfolios, and a preliminary estimate indicates a continu-
ously compounded rate of return from portfolio 1 of R1 = ln(56.625/40) = 35%
and from portfolio two of R2 = ln(56.25/40) = 34%. Is this diп¬Ђerence signiп¬Ѓcant
and are these returns reasonably accurate in view of the survivorship bias?
We assume a geometric Brownian motion for both portfolios,

(5.34)
dSt = ВµSt dt + ПѓSt dWt ,

and deп¬Ѓne O = S(0), C = S(T ),

H = max S(t), L = min S(t)
0tT 0tT

with parameters Вµ, Пѓ possibly diп¬Ђerent.
In this case it is quite easy to determine the expected return or the value of
any performance measure dependent on C conditional on survival, since this is
essentially the same as a problem already discussed, the valuation of a barrier
option. According to (5.27), the probability that a given Brownian motion
process having open 0 and close c strikes a barrier placed at l < min(0, c) is

zl
}
exp{в€’2
Пѓ2 T
SURVIVORSHIP BIAS 303

with
zl = l(l в€’ c).

Converting this statement to the Geometric Brownian motion (5.34), the prob-
ability that a geometric Brownian motion process with open O and close c
breaches a lower barrier at l is

zl
P [L В· l|O, C] = exp{в€’2 }
Пѓ2T

with
zl = ln(O/l) ln(C/l) = a(a + ln(C/O)).

Of course the probability that a particular path with this pair of values (O, C)
is a вЂњsurvivorвЂќ is 1 minus this or

zl
1 в€’ exp{в€’2 }. (5.35)
Пѓ2 T

When we observe the returns or the closing prices C of survivors only, the results
have been п¬Ѓltered with probability (5.35). In other words if the probability
density function of C without any barriers at all is f (c) (in our case this is a
lognormal density with parameters that depend on Вµ and Пѓ) then the density
function of C of the survivors in the presence of a lower barrier is proportional
to

ln(O/l) ln(c/l)
f (c)[1 в€’ exp{в€’2 }]
Пѓ2 T
l 2 ln(O/l) 2a
= f (c)(1 в€’ ( )О» ), with О» = = 2 > 0.
Пѓ2T
c ПѓT

It is interesting to note the eп¬Ђect of this adjustment on the moments of C for
various values of the parameters. For example consider the expected value of C
conditional on survival
Rв€ћ l
cf (c)(1 в€’ ( c )О» )dc
l
E(C|L в‰Ґ l] = R в€ћ l
f (c)(1 в€’ ( c )О» )dc
l
E[CI(C в‰Ґ l)] в€’ lО» E[C 1в€’О» I(C в‰Ґ l)]
(5.36)
=
P [C в‰Ґ l] в€’ lО» E[C в€’О» I(C в‰Ґ l)]
304 CHAPTER 5. SIMULATING THE VALUE OF OPTIONS

and this is easy to evaluate in the case of interest in which C has a lognormal
distribution. In fact the same kind of calculation is used in the development of
the Black-Scholes formula. In our case C = exp(Z) where Z is N (ВµT, Пѓ 2 T )
and so for any p and l > 0, we have from (3.11), using the fact that E(C|O) =
O exp{ВµT + Пѓ2 T /2}, (and assuming O is п¬Ѓxed),
в€љ
1
E[C p I(C > l)] = Op exp{pВµT + p2 Пѓ 2 T /2}О¦( в€љ (a + ВµT ) + Пѓ T p)
ПѓT

To keep things slightly less combersome, let us assume that we observe the
geometric Brownian motion for a period of T = 1. Then (5.36) results in
2 2 2
1 1
OeВµ+Пѓ /2
О¦( Пѓ (a + Вµ) + Пѓ) в€’ Oeв€’aО»+(1в€’О»)Вµ+(1в€’О») Пѓ /2 О¦( Пѓ (a + Вµ) + Пѓ(1 в€’ О»))
1 1
О¦( Пѓ (a + Вµ)) в€’ eв€’О»aв€’О»Вµ+О»2 Пѓ2 /2 О¦( Пѓ (a + Вµ) в€’ ПѓО»)

Let there be no bones about it. At п¬Ѓrst blush this is still a truly ugly and
opaque formula. We can attempt to beautify it by re-expressing it in terms more
like those in the Black-Scholes formula, putting

1 1
(Вµ в€’ a), and d2 (0) = (a + Вµ),
d2 (О») =
Пѓ Пѓ
d1 (О») = d2 (О») + Пѓ, d1 (0) = d2 (0) + Пѓ.

These are analogous to the values of d1, d2 in the Black-Scholes formula in the
case О» = 0. Then
2 2 2
eВµ+Пѓ /2
О¦(d1 (0)) в€’ eв€’О»a+(1в€’О»)Вµ+(1в€’О») Пѓ /2 О¦(d1 (О»))
E[C|L в‰Ґ l] = O (5.37)
.
О¦(d2 (0)) в€’ eв€’О»aв€’О»Вµ+О»2 Пѓ2 /2 О¦(d2 (О»))

What is interesting is how this conditional expectation, the expected close for
the survivors, behaves as a function of the volatility parameter Пѓ. Although this
is a rather complicated looking formula, we can get a simpler picture (Figure
5.10) using a graph with the drift parameter Вµ chosen so that E(C) = 56.25
is held п¬Ѓxed. We assume a = в€’ ln(30/40) (consistent with Table 5.2)and vary
the value of Пѓ over a reasonable range from Пѓ = 0.1 (a very stable investment)
through Пѓ = .8 (a highly volatile investment). In Figure 5.10 notice that for
small volatility, e.g. for Пѓ В· 0.2, the conditional expectation E[C|L в‰Ґ 30]
SURVIVORSHIP BIAS 305

Figure 5.10: E[C|L в‰Ґ 30] for various values of (Вµ, Пѓ) chosen such that E(C) =
56.25.

remains close to its unconditional value E(C) but for Пѓ в‰Ґ 0.3 it increases almost
linearly in Пѓ to around 100 for Пѓ = 0.8. The intuitive reason for this dramatic
increase is quite simple. For large values of Пѓ the process п¬‚uctuates more,
and only those paths with very large values of C have abeen able to avoid
the absorbing barrier at l = 30. Two comparable portfolios with unconditional
return about 40% will show radically diп¬Ђerent apparent returns in the presence
of an absorbing barrier. If Пѓ = 20% then the survivorвЂ™s return will still average
around 40%, but if Пѓ = 0.8, the survivorвЂ™s returns average close to 150%. The
practical implications are compelling. If there is any form of survivorship bias
(as there usually is), no measure of performance should be applied to the returns
from diп¬Ђerent investments, managers, or portfolios without an adjustment for
the risk or volatility.

In the light of this discussion we can return to the comparison of the two
portfolios in Table 5.2. Evidently there is little bias in the estimate of returns
for portfolio 2, since in this case the volatility is small Пѓ = 0.2. However there
is very substantial bias associated with the estimate for portfolio 1, Пѓ = 0.5.
In fact if we repeat the graph of Figure 5.10 assuming that the unconditional
306 CHAPTER 5. SIMULATING THE VALUE OF OPTIONS

Figure 5.11: The Eп¬Ђect of Surivorship bias for a Brownian Motion

return is around 8% we discover that E[C|L в‰Ґ 30] is very close to 56 5 when
8

Пѓ = 0.5 indicating that this is a more reasonable estimator of the performance
of portfolio 1.

For a Brownian motion process it is easy to demonstrate graphically the
nature of the surivorship bias. In Figure 5.11, the points under the graph of
the probability density which are shaded correspond to those which whose low
fell below the absorbing barrier at l = 30. The points in the unshaded region
correspond to the survivors. The expected value of the return conditional on
survival is the mean return (x-cooredinate of the center of mass) of those points
chosen uniformly under the density but above the lower curve, in the region
labelled вЂњsurvivorsвЂќ. Note that if the mean Вµ of the unconditional density
approaches the barrier (here at 30) , this region approaches a narrow band
along the top of the curve and to the right of 30. Similarly if the unconditional
standard deviation or volatility increases, the unshaded region stretches out to
the right in a narrow band and the conditional mean increases.
SURVIVORSHIP BIAS 307

We arrive at the following seemingly paradoxical conclusions which make it
imperative to adjust for survivorship bias The conditional mean, conditional on
survivorship, may increase as the volatility increases even if the unconditional
mean decreases.

Let us return to the problem with both an upper and lower barrier and
consider the distribution of returns conditional on the low never passing a barrier
Oeв€’a and the high never crossing a barrier at Oeb ( representing a fund buyout,
recruitment of manager by competitor or promotion of fund manager to Vice
President).It is common in process control to have an upper and lower barrier
and to intervene if either is crossed, so we might wish to study those processes
for which no intervention was required. Similarly, in a retrospective study we
may only be able to determine the trajectory of a particle which has not left
a given region and been lost to us. Again as an example, we use the following
data on two portfolio managers, both observed conditional on survival, for a
period of one year.

Portfolio Open price Close Price Lower Barrier Upper Barrier Volatility

56 5
1 40 30 100 .5
8

56 1
2 40 30 100 .2
4

If П† denotes the standard normal p.d.f., then the conditional probability
density function of ln(C/O) given that Oeв€’a < L < H < Oeb is proportional
uв€’Вµ
1
to where, as before
Пѓ П†( Пѓ )w(u)

2 2 2 2
w(u) = 1 в€’ eв€’2b(bв€’u)/Пѓ + eв€’2(a+b)(a+bв€’u)/Пѓ в€’ eв€’2a(a+u)/Пѓ + eв€’2(a+b)(a+b+u)/Пѓ в€’ E(W ),
в€’ ln(L)
ln(H) b a
], and
W = I[f rac1( )> ] + I[f rac1( )>
a+b a+b a+b a+b
a = в€’ln(30/40).
b = ln(100/40),
308 CHAPTER 5. SIMULATING THE VALUE OF OPTIONS

The expected return conditional on survival when the drift is Вµ is given by
Z b
uв€’Вµ
1
E(ln(C/O)|30 < L < H < 100) = uw(u)П†( )du.
Пѓ Пѓ
в€’a

where w(u) is the weight function above. Therefore a moment estimator of the
drift for the two portfolios is determined by setting this expected return equal
to the observed return, and solving for Вµi the equation
Z b
u в€’ Вµi
1
uw(u)П†( )du = Ri , i = 1, 2.
Пѓi Пѓi
в€’a

The solution is, for portfolio 1, Вµ1 = 0 and for portfolio 2, Вµ2 = 0.3. Thus the
observed values of C are completely consistent with a drift of 30% per annum
for portfolio 2 and a zero drift for portfolio 1. The bias again very strongly
eп¬Ђects the portfolio with the greater volatility and estimators of drift should
account for this substantial bias. Ignoring the survivorship bias has led in the
mutual funds.

Problems

1. If the values of dj are equally spaced, i.e. if dj = jв€†, j = ..., в€’2, в€’1, 0, 1, ...and
with S0 = 0, ST = C and M = max(S0 , ST ), show that

P [C > u and Cв€’M is even]
в€†
E[H|C = u] = M + в€† .
P [C = u]

2. Let W (t) be a standard Brownian motion on [0, 1] with W0 = 0. Deп¬Ѓne
C = W (1) and H = max{W (t); 0 В· t В· 1}. Show that the joint probabil-
ity density function of (C, H) is given by

f (c, h) = 2П†(c)(2h в€’ c)eв€’2h(hв€’c) , for h > max(0, c)

where П†(c) is the standard normal probability density function.
PROBLEMS 309

3. Use the results of Problem 2 to show that the joint probability density
function of the random variables

Y = exp{в€’(2H в€’ C)2 /2}

and C is a uniform density on the region {(x, y); y < exp(x2 /2)}.

4. Let X(t) be a Brownian motion on [0, 1], i.e. Xt satisп¬Ѓes

dXt = Вµdt + ПѓdWt , and X0 = 0.

Deп¬Ѓne C = X(1) and H = max{X(t); 0 В· t В· 1}. Find the joint proba-
bility density function of (C, H).
310 CHAPTER 5. SIMULATING THE VALUE OF OPTIONS

 << стр. 2(всего 2)СОДЕРЖАНИЕ