• No results found

Convergence of Option Rewards

N/A
N/A
Protected

Academic year: 2021

Share "Convergence of Option Rewards"

Copied!
38
0
0

Loading.... (view fulltext now)

Full text

(1)

Mälardalen University Press Dissertations No. 89

CONVERGENCE OF OPTION REWARDS

  Robin Lundgren 2010                

(2)

This work was funded by the Graduate School in Mathematics and Computing.

Convergence of Option Rewards c

2010 Robin Lundgren

Typeset by the author in LATEX documentation system.

ISSN: 1651-4238

ISBN: 978-91-86135-84-3

Printed by Arkitektkopia, V¨aster˚as, Sweden Copyright © Robin Lundgren, 2010

ISSN 1651-4238

ISBN 978-91-86135-84-3

(3)

Mälardalen University Press Dissertations No. 89

CONVERGENCE OF OPTION REWARDS

Robin Lundgren

Akademisk avhandling

som för avläggande av filosofie doktorsexamen i matematik/tillämpad matematik vid Akademin för utbildning, kultur och kommunikation kommer att offentligen försvaras fredagen den 12 november, 2010, 13.15 i Kappa, Högskoleplan 1, Mälardalens högskola,

Västerås.

Fakultetsopponent: Professor Johan Tysk, Uppsala universitet, Matematiska institutionen

(4)

Abstract

This thesis consists of an introduction and five articles devoted to optimal stopping problems of American type options. In article A, we get general convergence results for the American option rewards for multivariate Markov price processes. These results are used to prove convergence of tree approximations presented in papers A, B, C and E. In article B, we study the problem of optimal reselling for European options. The problem can be transformed to the problem of exercising an American option with two underlying assets. An approximate binomial-trinomial tree algorithm for the reselling model is constructed. In article C, we continue our study of optimal reselling of European options and give the complete solution of the approximation problem. In the article D, we consider general knockout options of American type. A Monte-Carlo method is used to study structure of optimal stopping domains generated by combinations of different pay-off functions and knockout domains. In article E the American option with knock out domains is considered. In order to show convergence of the reward functional the problem is reformulated in such a way that the convergence results in paper A can be applied.

.

ISSN 1651-4238

(5)

v

Abstract

This thesis consists of an introduction and five articles devoted to optimal stopping problems of American type options.

In article A, we get general convergence results for the American option rewards for multivariate Markov price processes. These results are used to prove convergence of tree approximations presented in papers A, B, C and E. In article B, we study the problem of optimal reselling for European options. The problem can be transformed to the problem of exercising an American option with two underlying assets. An approximate binomial-trinomial tree algorithm for the reselling model is constructed. In article C, we continue our study of optimal reselling of European options and give the complete solution of the approximation problem. In the article D, we consider general knockout options of American type. A Monte-Carlo method is used to study structure of optimal stopping domains generated by combinations of different pay-off functions and knockout domains. In article E the American option with knock out domains is considered. In order to show convergence of the reward functional the problem is reformulated in such a way that the convergence results in paper A can be applied.

(6)

vi

Acknowledgements

First of all I would like to thank my supervisor Dmitrii Silvestrov for all his guidance and support during my work, first as a student and now as a PhD student. I also would like to thank my assistant supervisor Anatoliy Malyarenko for interesting discussions and help with LATEX, and my second

assistant supervisor Kimmo Eriksson for support of the project. I want to thank my co-author Alexander Kukush at Kyiv National Taras Shevchenko University, who has given me comments and input regarding my research.

For contributing to my activities, I want to thank the Graduate School of Mathematics and Computing (FMB), ”Knut & Alice Wallenbergs Stiftelses resefond”, ”Magnusson stipendiet” at the Royal Swedish Academy of Science, SVeFUM ”Stiftelsen f¨or Vetenskaplig Forskning och Utbildning i Matematik” and ”Riksbankens Jubileumsfond”.

I want to thank all my colleagues at the Division of Applied Mathematics at M¨alardalen University, especially Richard Bonner, Lars-G¨oran Larsson, Hillevi Gavel, Eija Landquist, Mona Karlsson and Kristina Konpan.

On a more personal note I would like to thank my friends Nils-Hassan Quttineh, Oskar Schyberg and Fredrik Stenberg for interesting discussions and good support.

I also would like to thank my parents and brother, for all their love and support, especially my father for his will to introduce me to more general mathematics in my early years and help in my first years of mathematical studies.

Finally, I would like to thank the most important person in my life, my wife Maryna for all her love and support.

September 2010 Robin Lundgren

(7)

CONTENTS vii

Contents

Abstract v

Acknowledgements vi

List of papers viii

Communications ix

1 Introduction 1

2 Optimal Stopping 1

2.1 Optimal stopping in mathematical finance . . . 4

2.2 American knock out options . . . 8

2.3 Reselling of European options . . . 10

3 Convergence of Option Rewards 12 3.1 Convergence of option rewards for multivariate price processes . . . 13

3.2 Convergence in the reselling and related models . . . 15

4 Summary of the papers 18 4.1 Paper A . . . 18 4.2 Paper B . . . 19 4.3 Paper C . . . 19 4.4 Paper D . . . 20 4.5 Paper E . . . 20 5 Sammanfattning p˚a svenska 21 References 22

(8)

viii

List of papers

A: Lundgren, R., Silvestrov, D. (2009), Convergence of option rewards for multivariate price processes, Research Report 2009:10 Department of Mathematics, Stockholm University, Sweden, 53 pages, ISSN: 1650-0377.

B: Lundgren, R., Silvestrov D.S., Kukush, A.G. (2008), Reselling of op-tions and convergence of option rewards, Journal of Numerical and Applied Mathematics, 1(96), 149–172.

C: Lundgren, R., Silvestrov, D. (2011), Optimal stopping and reselling of European options, In: Rykov, V., Balakrishan, N., Nikulin, M. (Eds.) Mathematical and Statistical Models and Methods in Reliability. Birkh¨auser: Boston, 378-394.

D: Lundgren R. (2011), Simulation studies of stopping domains for Amer-ican knock out options, Journal of Statistical Planning and Inference, (to appear).

E: Lundgren, R. (2010), Convergence of American knock out options in discrete time, Research Report 2010-1, School of Education, Culture, and Communication, Division of Applied Mathematics, M¨alardalen University, Sweden, 18 pages.

Other publications

• Lundgren, R., Silvestrov, D. (2009), Convergence and approximation of option rewards for multivariate price processes, Research Report 2009-1, School of Education, Culture, and Communication, Division of Applied Mathematics, M¨alardalen University, 46 pages.

• Lundgren, R. (2007), Optimal stopping domains for discrete time knock out American options, In: Skiadas, C.H. (Eds.) Recent Advances in Stochastic Modeling and Data Analysis. World Scientific, 613–620. • Lundgren, R. (2007), Structure of optimal stopping domains for

Amer-ican options with knock out domains, Theory of Stochastic Processes, 13(29), no. 4, 98–129. (also published as: Lundgren R. (2007), Struc-ture of optimal stopping domains for American options with knock out domains, Research Report 2007-6, Department of Mathematics and Physics, M¨alardalen University, 38 pages.)

(9)

CONTENTS ix

Communications

Parts of this thesis has been presented in communications given at the fol-lowing international conferences:

1: XII International Conference on Applied Stochastic Models and Data Analysis (ASMDA 2007), Chania, Crete, Greece (2007).

2: X International Summer School ”Insurance and Finance: Science, Practice and Education”, Foros, Crimea, Ukraine (2007).

3: International School Finance, Insurance, and Energy Markets – Sustainable Development, V¨aster˚as, Sweden (2008).

4: XI International Summer School ”Insurance and Finance: Science, Practice and Education”, Foros, Crimea, Ukraine (2008).

5: International Workshop on Applied Probability (IWAP2008), Compi`egne, France (2008).

6: Symposium ”Optimal Stopping with Applications” Turku, Finland, (2009).

7: 6-th St Petersburg Workshop on Simulation, St. Petersburg, Russia, (2009).

8: International Workshop on Applied Probability (IWAP2010), Madrid, Spain, (2010).

(10)
(11)

2 OPTIMAL STOPPING 1

1

Introduction

This thesis primarily deals with convergence of American option rewards. When evaluating options of American type explicit solutions are often not available. The second best alternative available is to construct effective ap-proximate algorithms in discrete time to find the price. To be able to guar-antee that the approximation is correct we need convergence results. The convergence results obtained in this thesis gives conditions for payoff func-tions, transition probabilities and compactness conditions that needs to be satisfied.

Here we will give an introduction to the optimal stopping problem in mathematical finance, the convergence problem and a summary of the papers.

2

Optimal Stopping

In the optimal stopping problem we consider an underlying stochastic process S(t), t ≥ 0 of Markovian type, also called a price process. The optimal stopping problem is the problem of finding a stopping time τ∗that maximizes

some reward functional over the set Mmax,0,T of all Markov stopping times

0 ≤ τ ≤ T (i.e. such that the event {τ ≤ t} depends only on trajectory S(u), u ≤ t for t ≥ 0), and to find the corresponding maximum of the reward functional,

Φ(Mmax,0,T) = Eg(τ∗, S(τ∗)) = sup

τ ∈Mmax,0,T

Eg(τ, S(τ )), (2.1) where g(t, x) : R+×R+�→ R+is a measurable real valued function. Note that

the problem in (2.1) consists of two tasks: 1) evaluate the reward functional as explicit as possible, 2) to find an optimal stopping time τ∗ at which the

supremum is attained.

The books by Cont and Tankov [20], Rolski, Schmidli, Schmidt and Teugels [77], Schoutens [79] and Shiryaev [84] gives a comprehensive overview of models of stochastic processes used in finance and insurance.

A typical example for financial applications is where the underlying price process is modeled by the geometric Brownian motion,

S(t) = S(0)eµt+σW (t), t≥ 0 (2.2) where µ ∈ R, σ > 0 are real numbers, W (t) is a standard Brownian motion and S(0) is a positive constant.

We define the reward function w(t, x) as w(t, x) = sup

τ ∈Mmax,t,T

(12)

2

Γ Γc

T τ∗

Figure 1: An illustration of a stopping domain Γ, an optimal stopping time τ∗ and

a continuation domainΓc.

where Mmax,t,T is the set of all Markov stopping times t ≤ τ ≤ T (i.e.

such that the event {τ ≤ s} depends only on trajectory S(u), t ≤ u ≤ s for t ≤ s ≤ T ). The reward function in (2.3) is connected to the reward functional as

Φ(Mmax,0,T) = w(0, S(0)). (2.4)

As was shown in Chow, Robbins and Siegmund [18], and Shiryaev [83], in the case of a Markov process S(t), the optimal stopping time exists and has, under some minor conditions, the following form,

τ∗= inf{t ≥ 0 : S(t) ∈ Γt} ∧ T, (2.5)

where A ∧ B = min{A, B} and Γtis the optimal stopping domain at time t

defined as

Γt= {x : g(t, x) = w(t, x)}, 0 ≤ t ≤ T.

Thus, the optimal stopping time is the first moment when the process S(t) enters the stopping domain

Γ = {Γt,0 ≤ t ≤ T }.

This result is illustrated in Figure 1.

A similar formulation of the optimal stopping problem can be given for discrete time processes. Then, we assume that the underlying price process

(13)

2 OPTIMAL STOPPING 3 0 10 20 30 40 50 0 5 10 15 20 25 30 35 Underlying price Payoff 0 10 20 30 40 50 0 10 20 30 40 50 60 70 Underlying price Payoff

Figure 2: Left: Example of a standard put payoff function with strike price K = 35. Right: Example of a piecewise linear payoff function with strike prices K1= 35 and K2= 15.

follows a discrete time Markov process S(n), n = 0, 1, . . . ,. See for example the books by F¨ollmer and Schied [34] and Pliska [74] for discrete time models used in finance. The reward functional for the optimal stopping problem then has the following structure:

Φ(MN n) = Eg(τ ∗, S(τ)) = sup τ ∈MN n Eg(τ, S(τ )), (2.6) where MN

n is the set of all discrete Markov stopping times n ≤ τ ≤ N (i.e.

such that event {τ ≤ k} depends only on trajectory S(r), n ≤ r ≤ k for n≤ k ≤ N).

For financial applications in discrete time, a typical example is when the underlying price process is modeled by a geometric random walk,

S(n) = S(n − 1)eµ+σWn, n= 1, . . . , N, (2.7)

where µ ∈ R, σ > 0 are real numbers, Wn is a sequence of i.i.d. standard

normally distributed random variables and S(0) is a positive constant. The set of all asset prices at moment n such that the optimal strategy is to stop form an optimal stopping set.

Let us define the reward function wn(x) as

wn(x) = sup

τ ∈MN

n

E[g(n, S(n + 1))|S(n) = x], (2.8) The results by Chow, Robbins and Siegmund [18] and Shiryaev [83] tell that the optimal expected reward can be found by using the recursion (2.13),

(14)

4

i.e.,

Φ(MN

0) = Ew0(S(0)). (2.9)

The optimal stopping time maximizing (2.6) is the first time the price process hitting the stopping domain Γ

τ∗= min{k ≥ 0 : S(k) ∈ Γk} ∧ N. (2.10)

where

Γn= {x : g(n, x) = wn(x)}, n = 0, 1, . . . , N. (2.11)

The optimal stopping domain is a collection of optimal stopping sets and has the following structure:

Γ = {Γ0,Γ1, . . . ,ΓN}, (2.12)

In discrete time setting, a backward recursion can be applied, where the reward received if the process is stopped at moment n, is compared with the expected reward if the choice is to continue the process, i.e. the process is not stopped at moment n. This procedure starts at maturity time N and moves backward in time. The backward recursion has the following structure,

   wN(x) = g(N, x), wn(x) = max{g(n, x), E[wn+1(S(n + 1))|S(n) = x]}, n= N − 1, . . . , 0. (2.13)

In discrete time, Snell [90] formulated a general optimal stopping prob-lem for discrete time stochastic processes. He characterized the solution as the smallest supermartingale (also called the Snell envelope) dominating g(tn, S(n)), n = 0, . . . , N. We also mention the results in discrete time

of Darling, Liggett and Taylor [24], Dynkin [29] and Shiryaev [83]. For continuous time models, there exist few explicit solutions for the optimal stopping problem (2.1). They are usually attained by solving a parabolic free-boundary problem; see for example Dochviri [27], McKean [64], Salmi-nen [78] and Van Moerbecke [92]. The book by Peskir and Shiryaev [73] treats optimal stopping problems in both discrete and continuous time.

2.1

Optimal stopping in mathematical finance

A typical example of optimal stopping problem in mathematical finance are the American type options. Such an option is a contract that depends upon an underlying asset; it can be a stock, a currency, spot energy price, etc. The dynamics of the underlying asset is modeled by a stochastic price process.

(15)

2 OPTIMAL STOPPING 5 5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 55 Days Asset price 5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 55 Days Asset price

Figure 3: Left: An example of a stopping domain for a standard American put option having payoff function g(t, x) = e−rt

(35 − x)+. Right: An example of a

stopping domain for an American option having a non-standard payoff function given by g(t, x) = e−rtmax(20 + 3.5 · (15 − x), max(35 − x, 0)),

where r is daily interest rate corresponding to a yearly interest rate of 1%. Note in the picture to the right, at day 28 the stopping domain is not monotone increasing, this is an example of a classification error. Classification errors are studied in paperD.

The payoff function g(t, x) determines the reward received from the op-tion. A standard payoff function has the form g(t, x) = e−rt

(x − K)+ or

g(t, x) = e−rt

(K − x)+ for so-called ”call” or ”put” options.

We present stopping domains in Figure 3 that are connected to the payoff functions illustrated in Figure 2.

The special feature of an American type option is that the holder can choose to exercise the option at any moment of time up to time T and to receive the reward determined by the payoff function at that moment of time. For general mathematical finance and option theory, the books by Bj¨ork [6], Musiela and Rutkowski [66] and Shiryaev [84] may be used as reference publications.

As in general optimal stopping problems, explicit or numerical methods can be applied. McKean [64] showed that the optimal stopping problem for finding the price of an American option could be transformed into a partial differential equation with a free boundary. In Bensoussan [5] the existence of solutions for such problems were proved.

Kim [49] showed that the American option price can due to the possibility of early exercise, be decomposed into the sum of the corresponding European

(16)

6

price plus an early exercise premium, i.e.

Cam(S(t), t) = Ceu(S(t), t) + Ψ(S(t), t, b(t)). (2.14)

Where Cam, Ceudenote the price of American and European option

respec-tively, Ψ(s, t, b) denote the early exercise premium function and b(t) is the optimal stopping boundary. For American call and put options, written on stocks that pays continuous dividend Kim [49] uses the decomposition (2.14) to give explicit solutions for the perpetual call and put option and for the case when the interest and the dividend rate are equal to zero.

In the paper by Jacka [41], the American put option is studied. It is showed that the stopping domains are of one threshold structure and that the optimal stopping boundary b(t) is continuous. However, for non-standard payoff functions the stopping domain can have a multi-threshold structure. This case is studied in the papers by J¨onsson [44], J¨onsson, Kukush, and Silvestrov [45, 46] and Kukush and Silvestrov [51].

Carr, Jarrow and Myneni [16] continued the work on the early exer-cise premium but used an estimate of the stopping boundary instead of the boundary itself. That is to make the early exercise premium in (2.14) easier to approximate numerically.

Other works we like to mention that treat the American put option ana-lytically are Ekstr¨om and Tysk [32], Geske and Johnson [36] and Lamberton and Villeneuve [54].

But in most cases, the only solution available to the problem is to find a ε-optimal solution based on an appropriate discrete time approximation of the model.

In discrete time setting, the reward function in (2.13) can be evaluated by using three different approaches.

The first one is based on numerical solving the corresponding partial differential equation using finite difference methods. This method can solve the American option pricing problem for some constrains and is efficient when the dimension of the problem is low. This method were firstly used by Brennan and Schwartz [12] by using a finite difference method to price American put options with finite time horizon. The authors consider both the case with a dividend paying stock and the case without.

The second method is based on the use of the Monte Carlo method. An advantage of the Monte Carlo method is that it allows one to investigate the structure of the optimal stopping domains. One of the first to apply this method in mathematical finance was Boyle [8].

Unfortunately this method does require very comprehensive calculations and thus shows a slow performance. However, the Monte Carlo method be-comes more competitive when the dimension of the price process is high. As

(17)

2 OPTIMAL STOPPING 7

examples when Monte Carlo method are used for pricing of American op-tions we mention Broadie and Glasserman [14] where the authors propose a stochastic tree method where the tree is firstly constructed using simulation and then backward in the tree the option price is calculated. The draw-back with this method is that the number of exercise moments should be comparably small, otherwise this method show a slow performance. A more recent contribution is Longstaff and Schwartz [56] that proposes a technique to first generate several trajectories and then, backward, for each time step, determine the expected continuation profit by regression and then compare it with the stopping profit. A survey regarding valuation of American options using Monte Carlo simulation can be found in Glasserman [37]. Other works where Monte Carlo simulation is used are J¨onsson [44], J¨onsson, Kukush, and Silvestrov [45, 46], Kukush and Silvestrov [51], Schoutens and Symens [80], Silvestrov, Galochkin and Malyarenko [87] and Silvestrov, Galochkin and Sibirtsev [88]. We apply a Monte Carlo method in paper D in order to study the stopping domain for American options with knock out domains.

The third approach is based on direct recursion calculations in (2.13), and can be applied in the cases, where the underlying price process is a simple one, for example, a binomial or trinomial sum-process. This method is effective if one can build the corresponding tree model satisfying so-called recombining conditions that imply not more than a polynomial rate of growth of a number of nodes in the corresponding tree as a function of the number of steps. Also convergence of tree approximation rewards to the corresponding optimal rewards for continuous time models should be proven.

Parkinson [71] is one of the first to apply this method where a tree method was used to approximate the American put price. However, the correspond-ing convergence problem is not treated in this paper.

Several authors have used tree models to price American type option, we would like to mention Boyle [9], Breen [11] and Cox, Ross and Rubin-stein [22]. Works that are related with convergence and error analysis are Jiang and Dai [42], Lamberton [52] and Leisen [55]. The latter work proves convergence of the binomial tree method for American options. They also prove existence and convergence of the optimal exercise boundary in the bi-nomial tree approximation. A numerical comparison of different tree models is performed in Joshi [43].

We apply this method in papers A, B, C and E to the analysis of models for the reselling of European options, American type options for multivariate price processes, mean-reverse price processes used to model the stochastic dynamics of energy prices and discrete time American options with knock out domains.

(18)

8

will always be less or equal to the American option price. Chaudhury [17] and Chung and Chang [19] give upper bounds for the American option price and can thus give an interval where the true American option price should lie.

More recent work about American options has been done by Carr [15] and Pascucci [72].

The results of study of American options with non-standard payoff func-tions can be found in the papers by J¨onsson [44], J¨onsson, Kukush, and Silvestrov [45, 46] and Kukush and Silvestrov [51]. It should be noted that the optimal stopping domains for such options can possess a complex multi-threshold structure as shown in Figure 3 and Figure 4.

The stopping domains for American options with multiple underlying price processes are studied in Broadie and Detemple [13]. Just like the study with non-standard payoff functions, this can also give a complex multi-threshold structure.

A survey concerning American contingent claims can be found in the book by Detemple [26].

Some other examples, where optimal stopping is used in option theory relate to the Russian option introduced by Shepp and Shiryaev [82]. Kifer [48] introduced the game option, also known as the Israeli option, where both the holder and the issuer is allowed to exercise the option. It has also been investigated by Ekstr¨om [30].

2.2

American knock out options

The difference between knock out options and ordinary American options is that the possibility that the option will be worthless if the underlying process S(t) enters some special domain called a knock out domain.

American barrier options of knock out type are a special case of the American knock out option.

In Lundgren [57, 58] and papers D and E we study general American knock out options in discrete time. In this case the knock out domain H has the following structure H = {Hn: n ≥ 0}, where Hn are some measurable

subsets of R+ = (0, ∞).

A barrier option of knock out type is a special case of a knock out option when H1= H2= · · · = HN.

(19)

2 OPTIMAL STOPPING 9 0 10 20 30 40 50 0 10 20 30 40 50 60 70 Underlying price Payoff 0 5 10 15 20 25 30 35 40 45 50 0 5 10 15 20 25 30 35 40 45 50 Days Asset price 0 5 10 15 20 25 30 35 40 45 50 0 5 10 15 20 25 30 35 40 45 50 Days Asset price 0 10 20 30 40 50 0 5 10 15 20 25 30 35 40 45 50 Days Asset price

Figure 4: Up left: The considered piecewise linear payoff function g(t, x) = e−rtmax(20 + 3.5 · (15 − x), max(35 − x, 0)). Up right: The consid-ered knock out domain given by Hn= {x : 0 ≤ x ≤ 0.03n2− 1.5n + 27},

n = 1, . . . , 50. Down left: The stopping domain for the case without knock out domain. Down right: The optimal stopping domain for the combination of payoff function and knock out domain.

time setting takes the following form: ΦH(MT0) = Eg(τ ∗, S(τ))χ(τ∗ < τH) = sup τ ∈MN 0 Eg(τ, S(τ ))χ(τ < τH), (2.15) where τH= min{n ≥ 0 : S(n) ∈ Hn} ∧ N. (2.16)

To solve this discrete optimization problem, an analogue of the backward recursion of type (2.13) can be applied. A description of the optimal stopping domain and optimal stopping moments similar to (2.9) and (2.10) can be given.

It is worth noting that the actual form of the optimal stopping domains for American options of the knock out type may be much more complex

(20)

10

than for ordinary American type options, as can be seen in Figure 4. Earlier works by Dayanik and Karatzas [25], Gau, Huang and Subrahmanyam [35] and Karatzas and Wang [47] consider American options with constant bar-rier. Novikov, Frishling and Kordzakhia [69] and Roberts and Shortland [76] consider time dependent barrier options.

The work in papers D and E differs from previous works, firstly because the knock out domain is allowed to be any subset of the positive half line. Sec-ondly, we study optimal stopping domains for different combinations of payoff functions and knock out domains, where we give examples when the optimal stopping domains possess complex multi-threshold structure. Thirdly, we show convergence of option rewards for discrete time American type options with knock out domains.

2.3

Reselling of European options

In the reselling problem, it is assumed that an investor has bought an option of European type. The difference with American options is that European ones can only be exercised at maturity T . There, exists, however, the possi-bility to resell the option on the open market. The reselling problem is the problem of finding the optimal time τ∗ that maximizes the expected reward

for the investor and of finding the corresponding optimal expected reward. It is well known that the market price of European option at moment t deviates from the theoretical price. In fact, the market price randomly oscillates around the theoretical price. In this case, the implied volatility σ(t) corresponding to the market price is used. It respectively randomly oscillates around σ.

The use of the market price of the European option C(t, S(t), σ(t)) is an approach used in practice. In this case, C(t, S, σ) may be interpreted as some kind of utility reward function commonly recognized and used by market agents for evaluation of market option prices.

Such models with stochastic volatility have for example been investigated by Hull and White [40] and in the recent paper by Ekstr¨om and Tysk [33].

The use of market option prices actualizes the problem of reselling for European options. In this case, it is assumed that an owner of an option can resell the option at some stopping time from the class MT which includes

all stopping times 0 ≤ τ ≤ T that are Markov moments with respect to the filtration Ft = σ((S(s), σ(s)), s ≤ t) generated by the vector process

(S(t), σ(t)). It is worth to note that the process σ(t) is indirectly observable as an implied volatility corresponding to the observable market price of an option.

(21)

2 OPTIMAL STOPPING 11

that the option is already bought. In this case, the owner is only interested in finding the optimal expected reward for reselling the option. This implies that the investors estimation of the actual trend of the price process should be involved is the analysis. This analysis is realized in the paper.

Similar problem as the reselling problem for European call options have been investigated by Ekstr¨om, Lindberg, Tysk and Wanntorp [31], where the authors find the optimal time to sell a call spread. Henderson and Hobson [38] are considering the problem of selling an asset instead of an option.

In the paper by Kukush, Mishura and Shevchenko [50], the authors con-sider a model of the reselling of European options, where they assume that the price process and the stochastic volatility follow a bivariate geometrical Brownian motion with correlated noise terms. They show that there is no arbitrage in this model and consider methods for estimating parameters of the model. The drawback of modeling volatility by a geometrical Brownian motion is that the implied volatility and consequently the market price can infinitely deviate, respectively, from the corresponding historical volatility and the theoretical Black-Scholes price.

We assume that a price process for an underlying asset is represented by a geometric Brownian motion S(t), while the stochastic implied volatility is represented, as mentioned above, by a mean-reverse exponential Ornstein-Uhlenbeck process σ(t). It is also assumed that the noise terms for process S(t) and σ(t) are correlated.

The problem of optimal reselling of the European option is imbedded in the problem of optimal execution of the American type option for the two-dimensional process (S(t), σ(t)) having the payoff function given by g(t, s, σ) = e−rtC(t, s, σ), where C(t, s, σ) is the price of the European option given by

the Black-Scholes formula.

The corresponding optimization problem has the following form,

Φ(Mmax,0,T) = sup

τ ∈Mmax,0,T

Ee−rτC(τ, S(τ ), σ(τ )), (2.17) where Mmax,0,T is the set of all Markov stopping times τ ≤ T , (T is the

maturity of the option) for the two-dimensional process (S(t), σ(t)). It is worth noting that the stochastic implied volatility σ(t) is an indirectly ob-servable process since it is determined via the Black-Scholes formula for im-plied volatility upon the observable underlying asset price and option market price.

It should be noted that the structure of the two-dimensional stopping domain Γ is very complicated. The problem of a detailed description of the optimal stopping domain remains open.

(22)

12

The study of the reselling problem is relatively new, the other works in this area is Kukush, Mishura and Shevchenko [50].

Usually the market price oscillates around the theoretical price. That can be described by a mean-reverse model for stochastic volatility. Such a model was proposed in paper B and also studied in paper C, where an exponential mean-reverse Ornstein-Uhlenbeck process [70] is used to model the stochastic implied volatility. This stochastic volatility model is similar to the model studied in Heston [39], where the stochastic volatility is assumed to follow a Cox-Ingersoll-Ross process.

In the papers B and C we firstly introduce the described reselling model and show how the problem can be transformed to an optimal stopping prob-lem for American options. Secondly, we show how to build approximate algorithms for finding optimal reselling reward. Thirdly, we show conver-gence of the algorithm to the considered model, and finally we illustrate the model with numerical examples.

3

Convergence of Option Rewards

We study the convergence of optimal expected rewards for American type op-tions. In the seminal paper by Cox, Ross and Rubinstein [22] the convergence of European option values for the binomial tree model to the Black-Scholes values, developed in Black and Scholes [7], have been shown. It should, however, be mentioned that the use of binomial or trinomial trees has been treated before, see for example the paper by Parkinson [71] or the book by Sharpe [81]. Other work using discrete models of tree type, related to our interest, that we would like to mention: Boyle and Lau [10] discuss the prob-lem of pricing barrier options using binomial trees and show how to fit the binomial tree to overcome this problem. Nelson and Ramaswamy [67] pro-pose a binomial tree for a mean reverting model, in their model the jump probabilities depends upon distance to the long term mean. This should be compared with the discrete trinomial tree proposed in paper A for a mean reverting process, where the jump probabilities depend upon time. Nelson and Ramaswamy also show convergence of their model for European option prices, partly by applying results from the book by Stroock and Varadhan [91]. For American option prices they however state: ”Unfortunately, we have not been able to extend this argument to the American case rigorously.”

The main problem of convergence of American option prices is that, in general, the limit of stopping times is not a stopping time. Aldous [1, 3] give compactness condition for weak convergence of stopping times, where he looks on the difference in increments of the process. This can be

(23)

com-3 CONVERGENCE OF OPTION REWARDS 13

pared with the compactness condition used in paper A where the logarithmic increments are considered instead. Aldous [2] proved convergence of reward functions for a quasi-left continuous c`adl`ag process by introducing a pre-diction process he can define weak extended convergence of processes and filtration. Similar studies have been done by Lamberton and Pag`es [53].

The first paper showing convergence of American option rewards and stopping times were done by Amin and Khanna [4]. The difference to our work is that they assumed the sequence of underlying asset values converges weakly to a diffusion. They also show results when the limiting process is a diffusion with discrete jumps at fixed dates. We only require the price processes to be a Markov process. Other minor differences are that they have a fixed starting point and only consider discrete time processes with uniform time step. Dupuis and Wang [28] also consider diffusion type price processes and show that the rate of convergence in such setting is ε for the reward function and √εfor the stopping times, where ε is the grid spacing. Coquet and Toldo [21] show convergence with one dimensional c`adl`ag process as price process. The book by Prigent [75] cover results on weak convergence in financial markets by using martingale based methods; the methods are applied to both American and European type options.

In the paper by Silvestrov, J¨onsson and Stenberg [89] the convergence of American option rewards for Markov type price processes modulated by a semi-Markov index is considered. Similar convergence problems are consid-ered in paper A for multivariate Markov price processes. In our convergence studies, we use general methods developed in Silvestrov [85, 86] for limit theorems for randomly stopped stochastic processes.

In conclusion we also would like to mention the works by Cutland, Kopp, Willinger and Wyman [23], Mulinacci and Pratelli [65] and Nieuwenhuis and Vellekoop [68] that are related to convergence of option rewards.

Compared with earlier works, we show convergence of option rewards in a more general setting in paper A. First of all, we consider an underlying k-dimensional price process of Markov type. Secondly we show convergence for a general payoff function with the only requirement that it has power type upper bounds. Thirdly, we require weak convergence of the transition probabilities.

3.1

Convergence of option rewards for multivariate

price processes

We consider, for every ε ≥ 0, a general k-dimensional Markov price process �

S(ε)(t) = (S(ε)

1 (t), . . . , S (ε)

(24)

inter-14

ested in American type options with payoff functions g(ε)(t, �s) which have

not more than a polynomial rate of growth. We also assume that the transi-tion characteristics of the price process and the payoff functransi-tion depend on a small perturbation parameter ε ≥ 0 and converge to the corresponding limit characteristics as ε → 0.

Let M(ε)max,0,T be a set of all Markov stopping times τ(ε) ≤ T for an

American option with maturity T . The object of our interest is the following reward functional

Φ(M(ε)max,0,T) = sup

τ(ε)∈M(ε) max,0,T

Eg(ε)(τ(ε), �S(ε)(ε))). (3.1)

In paper A, we give general conditions for convergence of optimal ex-pected rewards Φ(M(ε)max,0,T) for perturbed multivariate Markov price

pro-cesses.

The conditions of convergence include the following assumptions: a) Not more than polynomial rate of growth for the payoff function. b) Point wise convergence of payoff functions g(ε)(t, �s) as ε → 0.

c) Weak convergence of local transition probabilities for price process �

S(ε)

(t) as ε → 0.

d) Moment compactness condition for the price processes based on some modulus ∆(ε)(·) in the form of relation lim

ε→0∆(ε)(c) → 0 as c → 0.

We approximate the optimal expected reward functional Φ(M(ε)max,0,T) by

the optimal expected reward functional Φ(M(ε)Π,T) for a so-called Bermudian

option that gives the holder the right to exercise the option only at discrete time points on some partition Π = �t0 = 0 < t1 < · · · < tn = T � on the

interval [0, T ]. The class of all Markov stopping times connected with the Bermudian option M(ε)Π,T ⊆ M

(ε)

max,0,T.

Under these conditions we first get so-called skeleton approximations,

Φ(M(ε)max,0,T) − Φ(M

(ε)

Π,T) ≤ φ(∆

(ε)(d(Π))), (3.2)

where φ(x) ≥ 0 is some explicitly given continuous function such that φ(x) → 0 as x → 0, and d(Π) = max1≤i≤nti− ti−1 is the maximal time step in

partition Π.

Second, we prove convergence of the optimal expected reward for the Bermudian options,

Φ(M(ε)Π,T) → Φ(M (0)

(25)

3 CONVERGENCE OF OPTION REWARDS 15

Finally, by choosing a sequence of partitions Π = ΠN such that d(ΠN) →

0 as N → ∞ and combining asymptotic relations (3.2) and (3.3) we get under the conditions listed above, the desirable convergence result,

Φ(M(ε)max,0,T) → Φ(M

(0)

max,0,T) as ε → 0. (3.4)

The method used to prove convergence of American option rewards is illustrated in Figure 5. Φ(M(ε)max,0,T) ⇒ Φ(M (ε) ΠN,T) ⇓ ⇓ Φ(M(0)ΠN,T) ⇐ Φ(M(0)max,0,T)

Figure 5: Scheme of convergence for multivariate option rewards. First to get skele-ton approximations as in the upper row, second to get convergence in the Bermudian option rewards, combining these two results the convergence of American option rewards can be achieved.

3.2

Convergence in the reselling and related models

Our general convergence results for American type options are illustrated in papers A, B, C and E by applications to approximation tree algorithms for evaluation of optimal expected rewards for reselling and some other related models with multivariate exponential diffusion price processes.

The optimal reward functional for the reselling model given in (2.17) is approximated in this method by the optimal rewards for the corresponding approximation binomial-trinomial tree model.

It is a non-trivial problem. One of the difficulties in constructing approx-imate tree models for exponential multivariate diffusion process is connected with a so-called recombining condition. This condition means that the pro-cess which has first a jump up and then a jump down should end up in the same point as if it has first a jump down and then a jump up, for any node in the tree. Figure 6 illustrates this condition.

A straight forward fitting of expectations and correlation matrices for in-crements of the multivariate log-price process and vector summands for the corresponding multivariate tree model yields jumps depending upon the po-sition of the corresponding node in the tree. This can cause the recombining

(26)

16

condition to be violated. Consequently, the corresponding tree has an expo-nential rate of growth for the number of nodes as a function of number of jumps, and, therefore, the corresponding approximation algorithm becomes ineffective.

In the reselling model, a mean-reverse diffusion process is used to model dynamics of stochastic volatilities and, thus, the difficulty described above takes place. In Figure 7 the two type of trees used to approximate the reselling problem is illustrated.

The method proposed in papers A, B and C, in order to avoid this problem, is based on a possibility of representation of the corresponding multivariate log-price diffusion process as a non-random transformation of a more simple multivariate Gaussian process with independent increments, inhomogeneous in time. This makes it possible to include the transformation function in the corresponding payoff function and to fit a tree model to the corresponding multivariate process with independent increments. This fitting is a simpler problem and the jumps in the tree can be chosen with values independent of positions of the process in such way that the recombining condition holds.

Examples of tree approximations considered in papers A, B, C and E are:

a) Reselling model, where the price of underlying asset and the implied volatility are modeled by correlated price processes given respectively, by a geometric Brownian motion and a mean-reverse Ornstein-Uhlenbeck process. The approximation trees are of binomial-trinomial type with values of jumps independent of time and given by σ√ε and −σ√ε for the binomial model and u√ε,0, −u√εfor the trinomial model ap-proximating, respectively, for the Brownian motion and the

Ornstein-Figure 6: Left: an example of an element of a binomial tree with violated recombining condition. Right: an example of an element of a recombining binomial tree.

(27)

3 CONVERGENCE OF OPTION REWARDS 17

Figure 7: Left: an example of a recombing binomial tree. Right: an example of a recombing trinomial tree.

Uhlenbeck process. The jump probabilities depend upon time. The payoff function is given by the Black-Scholes formula.

b) A price process is modeled by a bivariate geometric Brownian mo-tion with a general payoff funcmo-tion. The approximamo-tion trees are of binomial-binomial type with values of jumps independent of time and given by σ1√εand −σ1√εfor the first component and σ2√εand −σ2√ε

for the second component. The jump probabilities do not depend upon time. General payoff functions are considered. As an example, we would like to mention the payoff function in the model of exchange of assets, where the payoff function has the form g(t, �s) = e−rt(s

1− s2).

c) A price process is modeled by the mean-reverse Ornstein-Uhlenbeck process used to model energy prices. The approximation tree has a trinomial form. The values of jumps are u√ε,0, −u√ε. The jump probabilities depend upon time. The payoff function may be of the standard form e−rt

[s − K]+ or of a general form g(t, s).

d) A price process is approximated by a binomial tree, where it is assumed that a knock out domain are connected to the option. The problem is solved by an reformulation of payoff function and transition proba-bilities so that the knock out domain is included in the optimization problem and the conditions of convergence still hold.

In all cases, we build the backward recursion algorithms for finding the optimal expected rewards for the corresponding tree models and prove con-vergence of these rewards of optimal expected rewards for corresponding continuous time models.

(28)

18

4

Summary of the papers

The present thesis is based on the following papers: A: Lundgren, and Sil-vestrov [61] ,B: Lundgren, SilSil-vestrov, and Kukush [63], C: Lundgren and Silvestrov [62], D Lundgren [59] and E Lundgren [60].

In paper A, general convergence results are given for American option re-wards for multivariate Markov price processes and options with general payoff functions. These results are used in paper B to prove convergence of approx-imation algorithms for the reselling model. In paper B, the optimal stopping problem related to the model of reselling European options is studied, ap-proximation tree algorithm is built, and convergence of the apap-proximation is proved. In paper C, the reselling problem are further investigated. Here we essentially improve results of paper B and give the complete solution of the approximation problem that is to build up an effective approximation algo-rithm for evaluation of optimal reward functional for the reselling model. In paper D, the American option with knock out domains is considered. The structure of the stopping domains is investigated by Monte Carlo simulation. In paper E, we show convergence of the American option rewards for the knock out option by applying the results of paper A.

4.1

Paper A

In paper A, general multi-dimensional Markov price processes are consid-ered. American type options with payoff functions having power type upper bounds are studied for such price processes. We assume that the transi-tion characteristics of the price process and the payoff functransi-tion depend on a small perturbation parameter ε ≥ 0 and converge to the corresponding limit characteristics as ε → 0. Convergence results for American type options for perturbed price processes are presented for models with continuous and dis-crete time. Also, asymptotical uniform skeleton approximations connecting the continuous and discrete time reward functional are obtained.

Two examples of applications for the methods developed in the paper are considered. In the first example, we build an approximation tree algorithm for an American type option for the price process represented by a bivariate Brownian motion. In particular, the model of exchange of assets is consid-ered. In the second example, we build an approximation tree algorithm for exponential mean-reverse Ornstein-Uhlenbeck processes. Such processes are used to model stochastic dynamics of energy prices.

(29)

4 SUMMARY OF THE PAPERS 19

4.2

Paper B

In paper B we study the reselling problem for European options. A buyer pays a price C0 when buying an European option. Since the option is of the

European type we can only exercise it at maturity and receive CT. There

exists, however, the possibility to resell the option at any moment of time τ ∈ [0, T ] and to receive a market price Cm

τ . It is well known that the

market price Cm

τ fluctuates from the theoretical price Cτ. We model this by

introducing a stochastic implied volatility that is modeled by an exponential mean-reverse Ornstein-Uhlenbeck process. We are interested in finding the value of the optimal expected reward. This problem can be transformed to the problem of finding the optimal expected reward for an American option having payoff function depending upon both the underlying price and the implied volatility and given, under fixed values of price and volatility, by the Black-Scholes formula.

To solve this problem, we construct an effective approximation bivari-ate binomial-trinomial tree algorithm. In order to construct a re-combining binomial-trinomial tree having quadratic rate of growth of number of nodes as a function of number of steps, we move some time dependencies from the processes to the payoff function, making the payoff function more com-plex but simplifying the processes to the model of Gaussian processes with independent increments.

The results of paper A on the convergence of optimal expected rewards for multivariate Markov price processes are applied to show convergence of the optimal expected rewards for approximation binomial-trinomial tree to the optimal expected rewards for the corresponding continuous bivariate process.

4.3

Paper C

In paper C we continue the study of the reselling problem for European options. Here we consider the same setting as in paper B, but we essentially improve results and give the complete solution of the approximation problem. That is, we show how to build up an effective approximation algorithm for evaluation of optimal reward functional for the reselling model.

In order to be able to guarantee that the jump probabilities satisfies the conditions pn,+, pn,·, pn,−≥ 0 and pn,++ pn,·+ pn,−= 1 an upper bound can

be needed to formulate on the correlation coefficient between the underlying process S(t) and the implied volatility process σ(t). By dividing the time interval into k parts, we can construct a tree that has k different values of jumps. Then the correlation coefficient needs to satisfy the condition |ρ| < e−αT

(30)

20

Numerical example is given where we illustrate how the reselling reward depends upon changes in the parameters of the model.

4.4

Paper D

In paper D, price processes S(t) with log-normal multiplicative increments are studied. The structure of the stopping domains for different types of pay-off functions and knock out domains is studied using Monte Carlo simulation. The payoff functions are usually assumed to be convex and non-decreasing. In particular, the models with the standard payoff function, piece-wise linear payoff functions, and quadratic payoff functions are investigated. knock out domains under consideration are represented by unions of different rectan-gular forms.

The structure studies of the stopping domains are based on Monte Carlo simulation. The idea is to impose a grid structure of asset prices for n = 1, . . . , N. For each asset price sn,jon the grid, we use Monte Carlo simulation

to determine if the point sn,j belongs to the stopping set Γnor to the

contin-uation set Γc

n. The values of the future reward are estimated by averages of

rewards accumulated for a large number of simulated trajectories originating from the point sn,j. The value of the estimated expected future reward are

compared with the reward received if exercising the option at moment n. The algorithm described above determines each point in the grid if it be-longs to the stopping domain or the continuation domain. By its estimation nature, it can classify a point to belong to the stopping domain, while it should belong to the continuation domain and vice versa. Probabilities of misclassification are also studied. We study the case with and without knock out domain and for moments n = N − 1 and n = 1.

4.5

Paper E

In paper E the American option with knock out domains are considered. In order to show convergence of the reward functional the problem are re-formulated in such a way that the convergence results in paper A can be applied.

The method used in paper E connected with construction of a modulated price process and payoff function such that the optimization problem for knock out American options reduce to the optimization problem for standard American options (without knock out domains) for a modified price process. This makes it possible to apply the general convergence results for Amer-ican type options in paper A.

(31)

5 SAMMANFATTNING P˚A SVENSKA 21

This is illustrated by two examples where we firstly consider the case when price process is assumed to follow a discrete time geometric random walk. To this process we assume an American put option with constant barrier is connected. In the second example, we consider the same type of option as in the first example but the price process is assumed to follow a discrete time mean reverting process instead.

5

Sammanfattning p˚

a svenska

Avhandlingen behandlar optimala stopproblem inom finansmarknaden. Det typiska exemplet p˚a ett optimalt stopproblem inom finansv¨arden ¨ar den Amerikanska optionen. Den ger ¨agaren r¨atten att n¨ar som helst k¨opa/s¨alja underliggande till ett f¨orutbest¨amt pris.

Att best¨amma priset f¨or en Amerikansk option kan dock vara komplicerat. Detta beror p˚a vilken process man antar att den underliggande varan f¨oljer, vilken typ av avkastningsfunktion optionen har etc.

Ett s¨att att v¨ardera Amerikanska optioner ¨ar att skatta prisfunktionen och den underliggande processen med en enklare diskret process. F¨or att veta att skattningen ¨ar korrekt s˚a beh¨over vi konvergensresultat. I artikel A s˚a ger vi generella villkor f¨or avkastningsfunktionen, den underliggande processen och dess ¨overg˚angssannolikheter f¨or den diskreta skattningen f¨or att kunna garantera konvergens.

I artiklarna B och C s˚a studerar vi n¨ar det ¨ar optimalt att s¨alja en Europeisk option. Vi visar att det problemet ¨overensst¨ammer med problemet n¨ar man ska l¨osa in en Amerikansk option som beror p˚a hur en underliggande aktie och optionens implicita volatilitet. Vi visar hur en diskret algoritm kan konstrueras f¨or att skatta v¨ardet av en s˚adan strategi och visar att den diskreta modellen konvergerar genom att anv¨anda resultaten fr˚an artikel A. I artikel D s˚a tittar vi p˚a Amerikanska optioner som har n˚agot f¨orb-judet omr˚ade f¨or underliggande. Det vill s¨aga att om underliggande tr¨affar det f¨orbjudna omr˚adet s˚a blir kontraktet v¨ardel¨ost. Ett specialfall av detta ¨ar den Amerikanska barri¨ar optionen av knock out typ. N¨ar man studerar optimala stopproblem s˚a s¨ager man att det finns tv˚a omr˚aden som underlig-gande kan befinna sig i. Dels ett omr˚ade d¨ar det ¨ar optimalt att forts¨atta att ¨aga optionen, ett s˚a kallat forts¨attningsomr˚ade. Dels ett omr˚ade d¨ar det ¨ar optimalt att l¨osa in optionen och ta emot avkastningen, det s˚a kallade stop-pomr˚adet. Vi studerar stoppomr˚adet som h¨anger samman med knock out optionen f¨or olika typer av avkastningsfunktioner och f¨orbjudna omr˚aden. F¨or att hitta stoppomr˚adet s˚a anv¨ander vi oss utav Monte Carlo simulering. Artikel E behandlar hur modellen f¨or den Amerikanska knock out

(32)

op-22

tionen ska formuleras f¨or att konvergensresultaten fr˚an artikel A ska kunna anv¨andas ¨aven i detta fall f¨or att kunna garantera konvergens.

References

[1] Aldous, D. (1978), Stopping times and tightness, Annals of Probability, 6, no.2, 335–340.

[2] Aldous, D. (1981), Weak convergence of stochastic processes for processes viewed in the strasbourg manner, Unpublished Manuscript, Statistical Laboratory, Cambridge University.

[3] Aldous, D. (1989), Stopping times and tightness II, Annals of Probability, 17, no.2, 586–595.

[4] Amin, K., Khanna, A. (1994), Convergence of American option values from discrete- to continuous-time financial models, Mathematical Finance, 4, no. 4, 289–304.

[5] Bensoussan, A. (1984), On the theory of option pricing, Acta Applied Mathematics, 2, 139–158.

[6] Bj¨ork, T. (2004), Arbitrage Theory in Continuous Time, 2nd edition, Oxford University Press.

[7] Black, F., Scholes, M. (1973), The pricing of options and copotrate liabilities, The Journal of Political Economy, 81, no. 3, 637–654.

[8] Boyle, P.P. (1977), Options: a Monte Carlo approach, Journal of Financial Economics, 4, 323–338.

[9] Boyle, P.P. (1986), Option valuation using a three-jump process, Interna-tional Options Journal, 3, 7–12.

[10] Boyle, P.P., Lau, S.H. (1994), Bumping up against the barrier with the binomial method, The Journal of Derivatives, 1, 6–14.

[11] Breen, R. (1991), The accelerated binomial option pricing model, The Journal of Financial and Quantitative Analysis, 26, no. 2, 153–164.

[12] Brennan, M.J., Schwartz, E.S. (1977), The valuation of American put options, The Journal of Finance, 32, no. 2, 449–462.

[13] Broadie, M., Detemple, J. (1997), The valuation of American options on multiple assets, Mathematical Finance, 7, no. 3, 241–286.

[14] Broadie, M., Glasserman, P. (1997), Pricing American-style securities by simulation, Journal of Economics Dynamics and Control, 21, 1323–1352.

(33)

REFERENCES 23

[15] Carr, P. (1998), Randomization and the American put, The Review of Financial Studies11, no. 3, 597–626.

[16] Carr, P., Jarrow, R., Myneni, R. (1992), Alternative characterizations of American put Options, Mathematical Finance, 2, no. 2, 87–100.

[17] Chaudhury, M. (2007), Upper bounds for American options, Research in Finance, 23, 161–191.

[18] Chow, Y.S., Robbins, H., Siegmund, D. (1971), The Theory of Optimal Stopping, Houghton Mifflin, Boston.

[19] Chung, S-L., Chang, H-C. (2007), Generalized analytical upper bounds for American option prices, Journal of Financial and Quantitative Analysis, 42, no. 1, 209–228.

[20] Cont, R., Tankov, P. (2004), Financial Modelling With Jump Processes, Chapman & Hall/CRC.

[21] Coquet, F., Toldo, S. (2007), Convergence of values in optimal stopping and convergence in optimal stopping times, Electronic Journal of Probability, 12, no. 8, 207–228.

[22] Cox, J., Ross, S., Rubinstein, M. (1979), Option price: A simplified approach, Journal of Financial Economics, 7, 229–263.

[23] Cutland, N.J., Kopp, P.E., Willinger, W., Wyman, M.C. (1997), Convergence of Snell envelopes and critical prices in the American put, In: Dempster, A., Pliska, S. (Eds.), Mathematics of Derivative Securities, Cambridge University Press, 61–87.

[24] Darling, D. A., Liggett, T., Taylor, H. M. (1972), Optimal stopping for partial sums, The Annals of Mathematical Statistics, 43, no. 4, 1363–1368. [25] Dayanik, S., Karatzas, I. (2003), On the optimal stopping problem for

one-dimensional diffusions, Stochastic Processes and their Applications, 107, 173–212.

[26] Detemple, J. (2005), American-Style Derivatives: Valuation and Computa-tion, Chapman & Hall/CRC.

[27] Dochviri, V.M. (1988), On optimal stopping with incomplete data, Lecture Notes in Mathematics, 1299, 64–68, Springer.

[28] Dupuis, P., Wang, H. (2005), On the convergence from discrete to continuous time in an optimal stopping problem, The Annals of Applied Probability, 15, no. 2, 1339–1366.

(34)

24

[29] Dynkin, E. B. (1965), The optimum choice of the instant for stopping a Markov process, Soviet Mathematics Doklady, 4, 627–629.

[30] Ekstr¨om, E. (2006), Properties of game options, Mathematical Methods of Operations Research, 63, 221–238.

[31] Ekstr¨om, E., Lindberg, C., Tysk, J., Wanntorp, H. (2009), Optimal liquida-tion of a call spread, Preprint Uppsala University, Sweden.

[32] Ekstr¨om, E., Tysk, J. (2006), The American put is log-concave in the log-price, Journal of Mathematical Analysis and Applications, 314, no.2, 710–723.

[33] Ekstr¨om, E., Tysk, J. (2010), The Black-Scholes equation in stochastic volatility models, Journal of Mathematical Analysis and Applications, 368, 498–507.

[34] F¨ollmer, H., Schied, A. (2004), Stochastic Finance - An Introduction in Discrete Time, 2nd edition. de Gruyter Studies in Mathematics 27, Walter de Gruyter.

[35] Gau, B., Huang, J.Z., Subrahmanyam, M. (2000), The valuation of American barrier options using the decomposition technique, Journal of Economic, Dynamics and Control, 24, 1783–1827.

[36] Geske, R., Johnson, H. E. (1984), The American put option valued analyti-cally, Journal of Finance, 39, no. 5, 1511–1524.

[37] Glasserman, P. (2003), Monte Carlo Methods in Financial Engineering Springer Verlag, Berlin.

[38] Henderson, V., Hobson, D. (2008), An explicit solution for an optimal stopping/optimal control problem which models an asset sale, Annals of Applied Probabability, 18, 1681–1705.

[39] Heston, S.L. (1993), A closed-form solution for options with stochastic volatility with applications to bond and currency options, Review of Finan-cial Studies, 6, 327–343.

[40] Hull, J., White, A. (1987), The pricing of options on assets with stochastic volatilities Journal of Finance, 42, no. 2, 281–300.

[41] Jacka, S.D. (1991), Optimal stopping and the American put, Mathematical Finance, 1, 1–14.

[42] Jiang, L., Dai, M. (1999), Convergence of binomial tree method for American options, In: Hua, C. and Rodino, L. (Eds.), Partial Differential Equations and their Applications, World Scientific, 106–118.

(35)

REFERENCES 25

[43] Joshi, M.S. (2009), The convergence of binomial trees for pricing the American put, Journal of Risk, 11, no. 4, 87–108.

[44] J¨onsson, H. (2001), Monte Carlo studies of American type call options with discrete time, Theory of Stochastic processes, 7(23), no. 1-2, 163–189. [45] J¨onsson, H., Kukush, A.G., Silvestrov, D.S. (2004), Threshold structure of

optimal stopping strategies for American type option. I, Teoriya ˘Imovirnoste˘ı ta Matemmatychna Statystyka, 71, 82–92. (English translation in Theory of Probability and Mathematical Statistics, 71, 93–103).

[46] J¨onsson, H., Kukush, A.G., Silvestrov, D.S. (2005), Threshold structure of optimal stopping strategies for American type option. II, Teoriya

˘

Imovirnoste˘ı ta Matemmatychna Statystyka, 72, 42–53. (English translation in Theory of Probability and Mathematical Statistics, 72, 47–58).

[47] Karatzas, I., Wang, H. (2000), A barrier option of American type, Applied Mathematics and Optimization, 42, no. 3, 259–279.

[48] Kifer, Y. (2000), Game options, Finance and Stochastics, 4, no. 4, 443–463. [49] Kim, I. J. (1990), The analytic valuation of American options, Rewiev of

Financial Studies, 3, 547–572.

[50] Kukush, A.G., Mishura, Y.S., Shevchenko, G.M. (2006), On reselling of european option, Theory of Stochastic Processes, 12(28), no. 3-4, 75–87. [51] Kukush, A.G., Silvestrov, D.S. (2004), Optimal pricing for American type

options with discrete time, Theory of Stochastic Processes, 10(26), no. 1-2, 72–96.

[52] Lamberton, D. (1998), Error estimates for the binomial approximation of American put option, Annals of Applied Probability, 8, no. 1, 206–233. [53] Lamberton, D., Pag`es, G. (1990), Sur l’approximation des r`eduites, Annales

de l’Institut Henri Poincar´e, 26, no. 2, 331–355.

[54] Lamberton, D., Villeneuve, S. (2003), Critical price near maturity for an American option on a dividend-paying stock, Annals of Applied Probability, 13, no. 2, 800–815.

[55] Leisen, D.P. (1998), Pricing the American put option, a detailed convergence analysis for binomial models, Journal of Economic Dynamics and Control, 22, 1419–1444.

[56] Longstaff, F.A., Schwartz, E.S. (2001), Valuing American options by simulation: A simple least squares approach, Review of Financial Studies, 14, 113–147.

(36)

26

[57] Lundgren, R. (2007), Optimal stopping domains for discrete time knock out American options, In: Skiadas, C.H. (Eds.), Recent Advances in Stochastic Modelling and Data Analysis, World Scientific, 613–620.

[58] Lundgren, R. (2007), Structure of optimal stopping domains for American options with knock out domains, Theory of Stochastic Processes, 13(29), no. 4, 98–129.

[59] Lundgren R. (2010), Convergence of American knock out options in discrete time, Research Report 2010-1, School of Education, Culture, and Communi-cation, Division of Applied Mathematics, M¨alardalen University, Sweden, 18 pages.

[60] Lundgren R. (2011), Simulation studies of stopping domains for American knock out options, Journal of Statistical Planning and Inference, (to appear). [61] Lundgren, R., Silvestrov, D. (2009), Convergence of option rewards for multivariate price processes, Research Report 2009:10 Department of Mathematics, Stockholm University, Sweden, 53 pages, ISSN: 1650-0377. [62] Lundgren, R., Silvestrov, D. (2011), Optimal stopping and reselling of

European options, In: Rykov, V., Balakrishan, N., Nikulin, M. (Eds.) Mathematical and Statistical Models and Methods in Reliability, Birkh¨auser, Boston, 378-394.

[63] Lundgren, R., Silvestrov, D.S., Kukush, A.G. (2008), Reselling of options and convergence of option rewards, Journal of Numerical and Applied Mathematics, 96(1), 149–172.

[64] MacKean, H.P. (1965), Appendix: A free boundary problem for the heat equation arising from a problem in mathematical economics, Industrial Management Review, 6, 32–39.

[65] Mulinacci, S., Pratelli, M. (1998), Functional convergence of Snell envelopes: applications to American options approximations, Finance and Stochastics, 2, 311-327.

[66] Musiela,M., Rutkowski,M. (2005), Martingale Methods in Financial Mod-elling, 2nd edition. Springer-Verlag Berlin Heidelberg.

[67] Nelson, D., Ramaswamy, K. (1990), Simple binomial processes as diffusion approximations in financial models, Review of Financial studies, 3, 393–430. [68] Nieuwenhuis, J.W., Vellekoop, M.H. (2004), Weak convergence of tree methods, to price options on defaultable assets, Decisions in Economics and Finance, 27, no. 2, 87–107.

(37)

REFERENCES 27

[69] Novikov, A., Frishling, V., Kordzakhia, N. (2003), Time-dependent barrier options and boundary crossing probabilities, Georgian Mathematical Journal, 10 no. 2, 325–334.

[70] Ornstein, L., Uhlenbeck, G. (1930), On the theory of Brownian motion, Physical Review, 36, 823–841.

[71] Parkinson, M. (1977), Option pricing: The American put, Journal of Business, 50, no. 1, 21–36.

[72] Pascucci, A. (2008), Free boundary and optimal stopping problems for American asian options, Finance and Stochastics, 12, 21–41.

[73] Peskir, G. and Shiryaev, A.N. (2006), Optimal Stopping and Free Boundary Value Problems, Birkh¨auser.

[74] Pliska, S.R. (1997), Introduction to Mathematical Finance, Blackwell. [75] Prigent, J.L. (2003), Weak Convergence of Financial Markets, Springer

Verlag, New York.

[76] Roberts, G.O., Shortland, C.F. (1997), Pricing barrier options with time-dependent coefficients, Mathematical Finance, 7 no. 1, 83–93.

[77] Rolski, T., Schmidli, H., Schmidt, V., Teugels, J. (1999), Stochastic Processes for Insurance and Finance, Wiley.

[78] Salminen, P (1985), Optimal stopping of one-dimensional diffusions, Mathe-matische Nachrichten, 124 85–101.

[79] Schoutens, W. (2003), L´evy Processes in Finance, Wiley.

[80] Schoutens, W., Symens, S. (2003), The Pricing of Exotic Options by Monte-Carlo Simulations in a Levy Market with Stochastic Volatility, International Journal for Theoretical and Applied Finance6, no. 8, 839–864.

[81] Sharpe, W.F. (1978), Investments, 2nd edition. Prentice-Hall.

[82] Shepp, L.A., Shiryaev, A.N. (1993), The Russian option: reduced regret, The Annals of Applied Probability, 3 no. 3, 631–640.

[83] Shiryaev, A.N. (1976), Sequential Analysis. Optimal Stopping Rules, Nauka, Moscow, (English edition: Optimal Stopping Rules, Springer, New York, 1978.)

[84] Shiryaev, A.N. (1999), Essentials of Stochastic Finance - Facts, Models, Theory, World Scientific, Singapore.

[85] Silvestrov, D. (1974), Limit Theorems for Composite Random Functions, Vysshaya Shkola and Izdatel’stvo Kievskogo Universiteta, Kiev.

(38)

28

[86] Silvestrov, D. (2004), Limit Theorems for Randomly Stopped Stochastic Processes, Springer-Verlag, London.

[87] Silvestrov, D.S., Galochkin V.G., Malyarenko A.A. (2001), OPTAN - a pilot program system for analysis of options, Theory of Stochastic Processes, 7(23), no. 1-2, 291–300.

[88] Silvestrov, D.S., Galochkin V.G., Sibirtsev V.G. (1999), Algorithms and programs for optimal Monte Carlo pricing of American type options, Theory of Stochastic Processes5(21), no. 1-2, 175–187.

[89] Silvestrov, D., J¨onsson, H., Stenberg, F. (2008, 2009), Convergence of option rewards for Markov type price processes modulated by stochastic indices I, II, Teoriya ˘Imovirnoste˘ı ta Matemmatychna Statystyka, 79, 149–165, 80, 138–155, (English translation in Theory of Probability and Mathematical Statistics, 79, 153–170, 80, 153–172).

[90] Snell, J.L. (1952), Applications of Martingale system theorems, Transactions of the American Mathematical Society, 2, no. 73, 293–312.

[91] Stroock, D., Varadhan, S. (1979), Multidimensional Diffusion Processes, Springer Verlag, New York.

[92] Van Moerbeke, P. (1976), On optimal stopping and free boundary problems, Archive for Rational Mechanics and Analysis, 60, 101–148.

Figure

Figure 1: An illustration of a stopping domain Γ, an optimal stopping time τ ∗ and a continuation domain Γ c .
Figure 2: Left: Example of a standard put payoff function with strike price K = 35.
Figure 3: Left: An example of a stopping domain for a standard American put option having payoff function g(t, x) = e −rt (35 − x) +
Figure 4: Up left: The considered piecewise linear payoff function g(t, x) = e −rt max(20 + 3.5 · (15 − x), max(35 − x, 0))
+4

References

Related documents

46 Konkreta exempel skulle kunna vara främjandeinsatser för affärsänglar/affärsängelnätverk, skapa arenor där aktörer från utbuds- och efterfrågesidan kan mötas eller

Från den teoretiska modellen vet vi att när det finns två budgivare på marknaden, och marknadsandelen för månadens vara ökar, så leder detta till lägre

Av tabellen framgår att det behövs utförlig information om de projekt som genomförs vid instituten. Då Tillväxtanalys ska föreslå en metod som kan visa hur institutens verksamhet

In this work an adaptive technique [8] is used to solve the Black-Scholes equation along with using a principal component analysis (PCA) and asymp- totic expansion to reduce the

The difference is that a European option only gives the owner the right to exercise at a specific time in the future, whereas an American option can be exercised at any time up

Industrial Emissions Directive, supplemented by horizontal legislation (e.g., Framework Directives on Waste and Water, Emissions Trading System, etc) and guidance on operating

Here we will find the general analytical solution of American-style Asian option under jump diffusion process, for the case of floating strike and we will end by studying the

From Figure 2 one can see that the three-piece exponential (EP3) bound- ary of Ju (1998) ’tracks’ the ’true’ early exercise boundary rather well, whereas the four-piece