arxiv: v3 [q-fin.cp] 16 Nov 2015

Transcription

1 arxiv: v3 [q-fin.cp] 16 Nov 2015 The leas squares mehod for opion pricing revisied Maciej Klimek and Marcin Piera November 18, 2015 Absrac I is shown ha he he popular leas squares mehod of opion pricing converges even under very general assumpions. This subsanially increases he freedom of creaing differen implemenaions of he mehod, wih varying levels of compuaional complexiy and flexible approach o regression. I is also argued ha in many pracical applicaions even modes non-linear exensions of sandard regression may produce saisfacory resuls. This claim is illusraed wih examples. Key words: leas squares opion pricing, Snell envelopes, opimal sopping, approximaion of condiional expecaion, American opions, baske opions, Mone Carlo simulaion, LIBOR marke model, Heson-Nandi model. MSC 2010: Primary - 91G20, 91G60, 93E24; Secondary - 60G40, 62J02. 1 Inroducion For over a decade several varians of he so called leas squares mehod of American opion pricing have been widely used by financial praciioners and a he same ime sudied by researchers. The origins of he mehod can be found in he work of Carriere [5], Tsisiklis, Van Roy [23] (see also [22]), Longsaff, Schwarz [17] and Clémen, Lamberon, Proer [7]. Basically, he mehod seeks a way of approximaing condiional expecaions needed in he valuaion process eiher direcly as in [17] and [7], or indirecly hrough he value funcion as in [23]. A modificaion of he algorihm from [17] was sudied in [7] from he poin of view of he convergence of he mehod. Subsequenly, several papers on his subjec have been published we will menion jus a few of hem relaed o he presen aricle. Glasserman and Yu [11] invesigaed in 2004 he convergence of he leas squares like mehods, where basically he necessary condiional expecaions are approximaed by finie linear combinaions of approximaing funcions. More specifically, hey looked ino he problem of accuracy of esimaions when he number of approximaing funcions and he number Deparmen of Mahemaics, P.O.Box 480, Uppsala Universiy, Uppsala, Sweden; Maciej.Klimek@mah.uu.se Insiue of Mahemaics, Jagiellonian Universiy, ul. Łojasiewicza 6, Kraków. Poland; Marcin.Piera@im.uj.edu.pl 1

2 of simulaed rajecories increase. They assumed ha he underlying is a mulidimensional Markov process. The raher pessimisic oucome, from he pracical poin of view, is ha for polynomials as he approximaing funcions and for he convenional (resp. geomeric) Brownian moion as he underlying, he number of required pahs may grow exponenially in he degree (resp. he square of he degree) of he polynomials. Glasserman and Yu remarked ha similar propery may hold also for more general approximaing funcions (wih he number of approximaing funcions replacing he maximal degree). Also in 2004, Senof [21] analyzed and exended he convergence resuls presened in [7]. In paricular, he has considered he problem of choosing he opimal number of regressors in relaion o he number of simulaed rajecories. In 2005, Egloff [9] proposed an exension o he original Longsaff-Schwarz [17] as well as Tsisiklis-Van Roy ([22], [23]) algorihms by reaing he opimal sopping problem for mulidimensional discree ime Markov processes as a generalized saisical learning problem. His resuls also improve hose from [7]. Egloff commens ha despie very good performance of leas squares algorihms in some pracical calculaions, precise esimaes of he saisical quaniies involved in hese procedures may be difficul, leading o some less impressive performance in oher cases. Zanger [24] proposed in 2009 anoher exension o he leas squares mehod by considering fairly arbirary subses of informaion spaces as he approximaing ses. He has also produced some new and ineresing convergence resuls showing in paricular ha someimes he exponenial dependence on he number of ime seps can be avoided. I should be menioned ha he leas squares approach can be also seen as par of he sochasic mesh framework proposed by Broadie and Glasserman ([3], [4]; see also [16] and [10]). I should also be observed ha wo feaures seem o be common o he aricles menioned above. Firsly, he underlying is assumed o be Markovian. Secondly, he convergence raes of he mehod, in all is incarnaions, are no encouraging from he compuaional poin of view. In he presen paper, we exend he Clémen, Lamberon, Proer approach [7] o a fairly general seing for he regression approximaing condiional expecaions. Also, in a naural way, he underlying does no have o be Markovian and he pay-offs are allowed o be pah-dependen. While lack of Markov propery can be easily circumnavigaed in oher ways, his always implies addiional compuaional cos. Obviously, by aiming a beer approximaion of condiional expecaion, he poenial compuaional complexiy increases considerably. However, he main advanage of relaxaion of he assumpions is he increase in freedom o cusomize he mehod. Moreover, we would like o argue ha he leas squares mehods should be seen as a general framework leading o a variey of specific implemenaions. The main reason is essenially he fac ha he informaion space for condiional expecaion, or in oher words is range, is in many ineresing cases infinie dimensional. Ineviably, in hese cases any approximaion of condiional expecaions, or value funcions depending on condiional expecaions, has o involve significanly resricive exrinsic assumpions o make pracical compuaions possible. While general convergence resuls are necessary o moivae he overall approach and some compuaional complexiy may be addressed along he lines of [19], i is mos likely ha he fuure developmens will evolve in he direcion of simplified ime-series models. I is quie conceivable ha an alernaive source of realism and numerical efficiency could exploi he advances in boh ime-series analysis and frame heory (see e.g. [14]). The empirical basis for such speculaions comes from he fac ha in many real problems even by aking only a few 2

3 non-linear regressors, and someimes ignoring lack of he Markov propery, one migh arrive a saisfacory resuls from he pracical poin of view. There seem o be much anecdoal evidence coming from he financial indusry supporing he las saemen and in his paper we provide furher corroboraing evidence in he form of hree empirical examples. The maerial is organized as follows. The inroducion is followed by a shor review of consequences of he classic Dobrushin-Minlos heorem, which can lead o viable numerical approximaions of condiional expecaions. Afer recalling briefly how Snell envelopes are used in pricing of American-syle opions, we show ha he mehods proposed by Clémen, Lamberon and Proer [7] can be exended o cover he case of American syle opions wih a very general approach o regression. The seing includes pah dependen pay-offs and a non-markovian mulidimensional underlying. This is followed by hree compuaional examples illusraing he viabiliy of he mehod under raher resricive assumpions. Firs, we presen pricing of a one year Eurodollar American pu and call opions wih differen srike prices. Then, we use he leas squares approach o price a 1.5 monh American pu opion, whose payoff funcion depends on wo marke indices, namely DAX and EUROSTOXX50. Finally, we use he leas squares algorihm o price wo 1.5 monh American pu opions, whose payoff funcion is based on a single marke index under he assumpion ha he underlyings can be described by he Heson-Nandi GARCH(1,1) model [12]. Again, we will use EUROSTOXX50 and DAX indices as he respecive underlying insrumens. 2 Approximaion of condiional expecaion In his secion we will inroduce some basic noaion and recall a classic resul of Dobrushin and Minlos [8], which provides moivaion, as well as a choice of pracical recipes, for approximaion of condiional expecaions via he so called admissible projecion sysems. The Dodrushin- Minlos heorem shows a specific example of such approximaion, bu of course here exis infiniely many non-polynomial consrucions ha would have he same propery. Le (Ω,F,P) be a probabiliy space. Since we will be dealing only wih random variables of finie variance, we can rely on he Hilber space geomery in addressing he issues of ineres (see [20]). A closed subspace S L 2 (Ω,F,P) is said o be probabilisic if i conains consans and is closed wih respec o aking he maximum of wo of is elemens, i.e. if X,Y S, hen X Y S. For any non-empy se X L 2 (Ω,F,P), is laice envelope La(X) is defined as he smalles probabilisic subspace of L 2 (Ω,F,P) conaining X. Moreover, if X = {X 1,...,X n } and B n denoes he σ-algebra of Borel ses in R n, hen i is no difficul o prove ha La(X) = L 2 (Ω,σ(X),P) = L 2 (Ω,(X 1,...,X n ) 1 (B n ),P). The laer is someimes referred o as he informaion space generaed by X 1,...,X n. Even if X consiss of jus one scalar random variable, La(X) is ypically infinie-dimensional. Since i is also he range of he orhogonal projecion E[ X 1,...,X n ], i would be desirable from he numerical sandpoin o be able o approximae such projecions wih projecions ono smaller finie-dimensional vecor spaces using available leas squares algorihms. However, approximaing an orhogonal projecion wih infinie-dimensional range by projecions ono finie dimensional subspaces 3

4 makes mos error esimaes useless, unless he naure of he projeced objecs is somehow known beforehand. In order o consruc finie-dimensional approximaion of condiional expecaion one could use he following heorem, which is a sligh reformulaion of a resul of Dobrushin and Minlos [8]. Theorem 2.1. Le (Ω,F,P) be a probabiliy space and le α > 0. Le P n denoe he space of all polynomials of n real variables. If X 1,...,X n are random variables such ha e X j L α (Ω,F,P) for j = 1,...,n, hen: (a) P(X 1,...,X n ) L p (Ω,F,P) for any polynomial P P n and p [1, ); (b) he vecor space {P(X 1,...,X n ) : P P n } is dense in L p (Ω,σ(X 1,...,X n ),P) for every p [1, ). I should be noed ha he converse o par (a) is false as shown in he following example. Example 2.2. Le n = 1 and le X 1 be a discree valued random variable wih probabiliy mass funcion Since for any q 1 and α > 0 P[X 1 = m] = m=1 1 m lnm m=1 m q m lnm < and 1, m N. m lnm m=1 e αm =, mlnm he propery (a) from Theorem 2.1 is saisfied bu e X 1 L α (Ω,F,P). If he probabiliy measure P has a bounded suppor, in R n, hen he assumpion of he Dobrushin-Minlos heorem is rivially saisfied. In fac, in his special case he conclusion of he heorem follows direcly from he Sone-Weiersrass Theorem. I is also easy o see ha if X is Gaussian, hen e X L 1. However, if X is lognormal, hen is momen generaing funcion does no exis in he inerval (0, ) and hence e α X L α for all α > 0. In concree applicaions, he condiion e X L α can someimes be achieved by changing he probabiliy disribuion of very large values of X. For insance, his can be accomplished by runcaion of probabiliy disribuion or some direc aenuaion of he random variable X. Anoher possibiliy is he use of suiable weigh funcions. In his conex, he Dubrushin- Minlos heorem can be used o jusify he densiy par in he consrucion of several classic polynomial bases in spaces of square inegrable funcions associaed wih he names of Jacobi, Gagenbauer, Legendre, Chebyshev, Laguerre and Hermie (see e.g. [6]). Le V be an informaion space generaed by random variables X 1,...,X n. Suppose ha one can furnish a sequence of Borel funcions q m : R n R, wih m N, such ha he se {q m (X 1,...,X n ) : m N} is linearly dense in V (e.g. wih he help of he Dobrushin-Minlos heorem). Then he condiional expecaion operaor E[ X 1,...,X n ] is he poinwise limi of he sequence of projecions ono linear spaces V m = {q k (X 1,...,X n ) : 1 k m} as m ր. This observaion leads o an auxiliary concep of admissible projecion sysems. 4

5 Given a discree ime filraion {,Ω} = F 0 F 1... F T F in he probabiliy space (Ω,F,P), we define an admissible projecion sysem as a family of orhogonal projecions P m : L 2 (Ω,F,P) L 2 (Ω,F,P), where = 1,...,T and m N, wih ranges V 1 V 2 V 3..., whose union is dense in L 2 (Ω,F,P) for each value of. Noe ha for any such sysem and for any fixed, we ge poinwise convergence of he projecions P m o E[ F ]. However, his is no a norm convergence unless he underlying sequence of subspaces becomes consan afer finiely many seps. I is well known ha Snell envelopes are useful in valuaion of American pu opions in discree ime models (see e.g. [18], p.127). They also furnish he main heoreical ingredien of he leas squares opion pricing algorihm which is he main opic of his paper. The sandard use of Snell envelopes can be easily exended o provide pricing algorihms for more general American syle opions, ha is opions ha allow execuion a any ime prior o mauriy, bu wih a wide variey of pay-off paerns. For a given probabiliy space (Ω,F,P), le (F ) T =0 be a filraion, where F 0 = {,Ω} and F T = F. Assume ha an adaped sochasic process (Z ) T =0 is inegrable. One could look a (Z ) as he inrinsic value process, ha is he (discouned) value of execuing some Americansyle opion a ime. The Snell envelope of (Z ) is defined as he adaped process U such ha U T = Z T and U = max(z,e[u +1 F ]) for {0,...,T 1}. The value U 0 corresponds o he he price of he opion associaed wih pay-offs given by (Z ). Indeed (U ) can be seen as he applicaion of dynamic programming principle o he opimal sopping problem sup{ez ν : ν C0 T}, where CT 0 denoe he se of all sopping imes wih values in he se {0,1,...,T} (cf. [15] and references herein, for basic properies of Snell envelopes and applicaions o pricing American-syle opions). The dynamic programming principle could be also rewrien in erms of he series of sopping imes (τ ), defined recursively by puing τ T = T and τ = 1 {Z E[Z τ+1 F ]} +τ +1 1 {Z<E[Z τ+1 F ]}, = 1,...,T 1. In paricular, we ge U = E[Z τ F ] and consequenly, τ 0 is opimal for (Z ). The key elemen in any numerical implemenaion of Snell envelopes is he abiliy o approximae he condiional expecaion operaor. Excep for he finie case, one has o deal wih infinie-dimensional spaces of random variables. Some elucidaion seems o be in order here. Given an admissible projecion sysem (P m ), for a fixed m N we define he sopping imes τ [m] by recursion, puing τ [m] T τ [m] = T and = 1 {Z P m(z τ [m] )} +τ [m] {Z <P m(z τ [m] +1 )}, = 1,...,T 1. Then he following heorem generalizes a resul due o Clémen, Lamberon and Proer (see Theorem 3.1 in [7]): [ ] Theorem 2.3. If (P m ) is an admissible projecion sysem, hen lim m E Z [m] τ F = E[Z τ F ] for = 1,...,T, where he convergence is in L 2. 5

6 Proof: Despie a much more general seing we have adoped here, we can use sandard properies of projecions in Hilber spaces and proceed as in [7]. Obviously, he above consideraions remain valid for vecor valued sochasic processes. 3 The leas squares mehod of opion pricing Assuming ha he filraion is generaed by a discree ime mulivariae sochasic process, we will show how o use Mone Carlo mehods o approximae numerically he value for he opimal sopping for a given adaped process (Z ), i.e. how o approximae he Snell envelope (U ) of ha process. To do so, given an admissible projecion sysem, we basically need o approximae [ numerically E Z τ [m] ] for m N, due o Theorem 2.3 and he fac ha U 0 = max(z 0,E[Z τ1 ]). In wha follows we will denoe he se of all real (m n)-marices by R m n wih he convenion ha R m = R 1 m. Throughou he secion we will use noaion and mehods similar o hose inroduced in [7] bu adaped o our less resricive assumpions. Suppose ha (X ) T =0 is a discree ime d-dimensional sochasic process on he probabiliy space (Ω,F,P), wih X 0 being a consan. This process is mean o represen he prices of he underlying asses for an American syle opion we wish o valuae. Le X = (X 1,...,X T ) : Ω R d T and le F = σ(x 0,...,X ) = σ(x 1,...,X ) for = 1,...,T. Given a family of Borel funcions f : R d (+1) R +, where = 0,...,T, we define Z = f (X 0,...,X ) for = 0,...,T. This sequence represens suiably discouned inrinsic prices of he opion we wan o consider. Such a general choice of funcions f expands he poenial applicabiliy well beyond American pu opions. Nex, we need o chose an admissible projecion sysem for he filraion associaed wih X. This is equivalen o choosing for each {1,...,T} a suiable sequence of Borel funcions q k : R d T R, where k N, which depend only on he firs column variables, and are such ha he sequence {q k(x)} k N is linearly dense and linearly independen in he space L 2 (Ω,σ(X 1,...,X ),P). Then, we can selec an increasing sequence of inegers (k m ) m N, such ha he spaces V m = Lin{q k(x) : k = 1,...,k m} and he orhogonal projecions P m : L 2 (Ω,σ(X),P) V m have all he righ properies. The symbol Lin denoes he linear envelope of he given se of vecors. If he sopping imes τ [m] are defined as in he previous secion, hen for some α m R km 1 we have ( ) P m Z [m] τ = e m (X)α m, +1 where he mapping e m is given by he formula e m = (q 1,...,qkm ) : R d T R km. In view of our assumpions, he Gram marix of he componens [ of e m (X) (wih ]] respec o he inner produc (Y 1,Y 2 ) E[Y 1 Y 2 ]), ha is he marix A m = E [q(x)q i j (X) R km km, 1 i,j k m 6

7 is inverible and hence α m = (A m ) 1 [ ] E Z [m] τ q 1(X) +1 [. ] E Z [m] τ q km (X) +1. Given a number N, ( he nex sep i) o use Mone-Carlo simulaion o generae independen rajecories X (n) = X (n) 1,...,X (n) T R d T of he process X, for n = 1,2,...,N. Each simulaion has he fixed saring poin X (n) 0 = X 0 R d 1. ( ) [ ] Define Z (n) := f X (n) 0,...,X (n) and le Ẑ = Z (1),...,Z (N) R N 1. This column vecor consiss simply of he values a ime of all simulaed rajecories of he process Z. Define also q(x k (1) ) V (m,n) = Lin. : k = 1,...,k m RN 1 q(x k (N) ) and P (m,n) = Proj V (m,n) : R N 1 R N 1 wih respec o he inner produc x,y N, where x,y denoes he sandard scalar produc. Noe ha e m V (m,n) (X(1) ) = Lin he columns of he marix. R N km RN 1. e m (X (N) ) If we define he sopping imes τ [m] by leing τ [m] T = T and he formula τ [m] = 1 {Z P m(z τ [m] )} +τ [m] +1 hen for some α m R km 1 we have P m +1 1 {Z <P m(z τ [m] +1 ( ) Z [m] τ = e m (X)α m. +1 )}, = 1,...,T 1, Similarly, if we define he approximaive sopping imes τ n,m,n and by puing by requiring ha τ n,m,n T = T τ n,m,n = 1 { Z (n) π n [ P (m,n) for = 1,...,T 1, ]} +τ n,m,n +1 1 { [ (Ẑτ n,m,n ) Z (n) <π n +1 P (m,n) ]}, (Ẑτ n,m,n ) +1 7

8 where π n : R N 1 R is he projecion on he n-h coordinae, hen for some α (m,n) R km 1 we have P (m,n) Z (1) τ 1,m,N +1. Z (N) τ N,m,N +1 = e m (X (1) ). e m (X(N) ) α (m,n). Le A (m,n) denoe he (k m k m )-Gram marix associaed wih he columns of he marix e m (X(1) ). e m (X(N) ), (wih respec o he inner produc x,y N ). Then his is simply he Gram marix esimaor for he given sample. Then α (m,n) is a soluion of he equaion A (m,n) α (m,n) = 1 N By he Law of Large Numbers A (m,n) he marix A (m,n) e m (X(1) ). e m (X (N) ) Z (1) τ 1,m,N +1. Z (N) τ N,m,N +1. a.s. A m as N, and hence for sufficienly large N is inverible (almos surely). In his case α (m,n) = 1 N ( A (m,n) ) 1 e m (X (1) ). e m (X (N) ) For convenience we will wrie α m = ( α m 1,...,αm T 1) and α (m,n) = objecs are k m (T 1)-marices. Z (1) τ 1,m,N +1. Z (N) τ N,m,N +1 The nex heorem is a direc exension of Theorem 3.2 and Lemma 3.2 from [7]. Theorem 3.1. Wih he above noaion, as N, 1 N N n=1 Z (n) τ n,m,n [ a.s. E Z τ [m] ], = 1,...,T.. ( ) α (m,n) 1,...,α (m,n) T 1. Boh 8

9 Proof: Define B = {(a m,z,x) : z < e m (x)a m } R km (T 1) R T R d T for = 1,...,T 1, where a m = (a m 1,...,am T 1 ), z = (z 1,...,z T ), and x = (x 1,...,x T ). By B c we will denoe he complemen of B. We define an auxiliary funcion F : R km (T 1) R T R d T R, by recursion by puing F T (a m,z,x) = z T and F (a m,z,x) = z 1 B c + F +1 (a m,z,x)1 B for = 1,...,T 1. I is easy o see ha F (a m,z,x) = z 1 B c + T 1 s=+1 z s 1 B... B s 1 B c s +z T1 B... B T 1 for = 1,...,T 1. Moreover, F (a m,z,x) is independen ofa m 1,...,am 1, F (α m,z,x) = Z [m] τ and F (α (m,n),z (n),x (n) ) = Z (n). For = 2,...,T define also wo oher auxiliary funcions τ n,m,n G (a m,z,x) = F (a m,z,x)e m 1 (x) and ψ (a m ) = E[G (a m,z,x)]. Using his noaion, one can see ha for = 1,...,T 1: α m = (A m ) 1 ψ +1 (α m ); (1) = (A (m,n) ) 1 1 N G +1 (α (m,n),z (n),x (n) ). N (2) α (m,n) n=1 The following esimae is a higher-dimensional counerpar of Lemma 3.1 in [7] and can be derived along he same lines as ha lemma: [ T T 1 ] F (a,z,x) F (ã,z,x) z s 1 { zs e m s (x)ãs em s (x) ãs as }, (3) s= where 1 T 1, a = (a 1,...,a T 1 ) R km (T 1), ã = (ã 1,...,ã T 1 ) R km (T 1), z R T and x R d T. Using (1, 2, 3), and under he echnical assumpion ha P(e m (X)αm = Z ) = 0, he reasoning from [7] can be easily modified o work wihin our more general seup. In general, his addiional echnical requiremen can be fulfilled by using approximaion of he conrac funcions f by funcions wih probabilisically negligible fibers and by inroducion of small amoun of random noise perurbing he probabiliy disribuion of X. Theorems 2.3 and 3.1 provide a recipe for approximaion of E[Z τ1 ] and hence also U 0 = max(z 0,E[Z τ1 ]), as required. s= 4 Examples In his secion we show hree examples of applicaions of he above leas squares algorihm. The firs example covers American call and pu opions wrien on Eurodollar fuures, which are assumed o conform o he Brace-Gaarek-Musiela model [2]. Nex we price baske and dual-srike American pu opions for EUROSTOXX50 and DAX indices, under he sandard 9

10 bivariae Brownian dynamics. Finally, we show how o price univariae American pu opions, boh for EUROSTOXX50 and DAX indices, assuming ha he dynamics of he underlyings could be expressed using he Heson-Nandi GARCH(1,1) model [12]. We have decided no o include convergence speed analysis, as i would make he presened examples much more complicaed (e.g. he proper variance reducion echnique is a crucial sep for any marke implemenaion), wihou adding much o he conclusions drawn in his paper. We refer o [7, 1], and references herein, for deailed analysis abou convergence speed in univariae Markovian case. For ransparency, we use only he sandard models for parameer esimaion and Mone Carlo simulaion. In paricular, only he prices of he underlyings are used for calibraion purposes and no Mone Carlo variance reducion echnique is implemened. Neverheless, we presen he (smoohed) densiy funcion of simulaed prices for every example (see Figures 2, 4 and 5) o give some insigh ino he accuracy of our implemenaion. I should be noed ha while our examples are raher sraighforward, he accuracy seems o be saisfacory. This allow us o be opimisic abou he leas squares algorihm approach for opion pricing, even when he dynamics of he underlyings is complicaed and no heoreical price is known. Our implemenaion of he leas squares algorihm is based on he in-he-money realizaions o speed up he convergence and reduce he number of polynomials needed o achieve sufficien level of accuracy. I is worh menioning ha in real-world models, o improve he convergence rae and he speed of he algorihm one migh use addiional informaion available from he marke (e.g. prices of various derivaives, based on he same underlying insrumens) as well as various modificaions of he sandard Mone Carlo algorihm (for more advanced models cf. [1] or [10] and references herein). All compuaions were done using R (64-bi). In paricular we have used he libraries fopions (for Heson-Nandi parameer calibraion, CRR prices and Mone Carlo simulaion), orhopolynom (for differen base funcions in L-S algorihm), imeseries (for marke daa handling) and Rsge (for parallel compuaions). 4.1 Eurodollar opions In his subsecion we use he leas squares algorihm o price one year Eurodollar American pu and call opions wih differen srike prices, given he real-marke daily prices of he Eurodollar fuures. I should be noed ha he sandard Black-Scholes model canno be used when he opion price is based on more han one LIBOR rae (e.g. when he opion s lifeime is longer han 3 monhs). This is due o he fac ha forward raes over consecuive ime inervals are relaed o each oher and canno all be log-normal under he same spo risk-neural measure. Consequenly, models of such insrumens in he sandard risk-neural seing are based on non-markovian dynamics. A. Brace, D. Gaarek and M. Musiela [2] proposed a model which can overcome his inconvenience (BGM Model) by uilizing a forward arbirage-free risk-neural measure. In he lieraure, i is also referred o as he LIBOR Marke Model (LMM). I is worh menioning ha he dynamics of ineres raes described in BGM model is very closely relaed 10

11 o he Heah-Jarrow-Moron (HJM) Model. Nex we will presen a brief overview of he BGM model, followed by some basic informaion concerning he seup of he leas squares algorihm The Brace-Gaarek-Musiela Model. The Brace-Gaarek-Musiela model is a sochasic model of ime-evoluion of ineres raes. I will be used here o simulae he (Mone Carlo) pahs of LIBOR fuures. We will now ouline a simplified version of he model ha suis our framework, and we will make some commens on he esimaion procedure. Le T 0 = 0 and T i = T i for i = 1,2,3,4. In realiy he daes of expiraion for he consecuive Eurodollar fuures differ slighly from 90 days. This migh poenially have an impac on he resuls, especially when we consider shor erm opions. Neverheless, we will use he heoreical values for simpliciy. Le L 0 be a spo LIBOR rae and le L i : [0,T i ] Ω R be he i-h forward LIBOR rae. Assuming d sources of randomness, he dynamics of he i-h LIBOR rae can be described by he equaion i dlogl i () = δ j L j () 1+δ j L j () σ j() σ i() σ i ()d+σ i ()dw Q Spo (), 2 j=i() where [0,T i ], δ i = T i+1 T i = 3/12 is he lengh of he accrual period of he i-h LIBOR forward rae, σ i () : [0,T i ] Ω R d is he insananeous volailiy of he i-h LIBOR forward rae, i() denoes he index of he bond (corresponding o he appropriae Eurodollar fuure) which is firs o expire a ime, and finally, W Q Spo() is a sandard (d-dimensional) Brownian moion under he spo LIBOR measure Q Spo (see [13] for more deails). We are assuming here ha he sources of randomness are independen of each oher and ha he proper dependency srucure is modelled wih σ i. For he Mone Carlo simulaion we will use a sandard Euler discreizaion of he above SDE, wih he ime sep = 1 360, i.e. i logl i () = δ j L j () 1+δ j L j () σ j() σ i() σ i () +σ i ()ǫ, (4) 2 j=i() where ǫ N(0,I) is a d-dimensional sandard normally disribued random vecor. In our implemenaion we will use d = 3. To calibrae he model we need o define he funcions σ i (), for i = 1,2,3,4. We will assume ha σ i () (for i = 1,2,3,4) is ime homogeneous, i.e., ha here exiss a funcion λ = (λ 1,λ 2,λ 3 ) : [0,T] R 3 such ha σ i () = λ(t i ) for [0,T i ] and i = 1,2,3,4. We will provide values of λ(t i ) for i = 1,2,3,4 and assume ha λ() = λ(t i ) for [T i 1,T i ]. We will apply he Principal Componens Analysis (PCA) o Eurodollar fuures daa o approximae he values of he 4 3 marix Λ = [λ j (T i )]. In oher words, we base our esimaion process on he correlaion beween he Eurodollar fuures. The difficuly wih calibraion of PCA is ha Eurodollar fuures have fixed mauriy daes, and so for a given T we can monior a conrac wih volailiy λ(t) only once per hree monhs. To overcome his, we will use linear inerpolaion of he quoed prices of Eurodollar fuures (which is in fac a common marke pracice). I should be noed ha we need he L 5 prices o perform such inerpolaion. Using 11

12 ha approach o Eurodollar fuures prices, we obain he values of conracs wih volailiy λ(t i ) (i=1,2,3,4) for every rading day. We also use linear inerpolaion of forward LIBOR raes for he days when he marke is no operaing (i.e. we inerpolae he conrac prices using known quoes from he las rading day before and he nex rading day afer he dae in quesion). Because of ha assumpion, in order o conduc he PCA and o esimae σ i (for i = 1,2,3,4), we will need (for each day) he prices of he five Eurodollar fuures closes o delivery. Le us now commen on he PCA esimaion process. We assume ha λ j (T i ) = Θ i s j α i,j d, k=1 s2 k α2 i,k for i = 1,2,3,4 and j = 1,2,3. Here, s 2 j denoes he variance of he j-h facor compued by PCA (wih s 2 1 s2 2 s2 3 ), α i,j measures he influence of he j-h facor when he ime o mauriy is in he period [T i 1,T i ] and Θ i := 3 j=1 s jα i,j is he oal volailiy in he i-h period. We also assume ha he facors are uncorrelaed and ha he relaive influence of every facor is 1 (i.e. for j 1,j 2 {1,2,3} we have 4 i α i,j 1 α i,j2 = 0 if j 1 j 2 and 4 i α i,j 1 α i,j2 = 1 if j 1 = j 2 ). Combining (4) and he parameers from he PCA we will be able simulae Eurodollar fuures pahs Seup, daa deails and he leas squares mehod parameers. We wish o price he quarerly Eurodollar American call and pu opions EDZ2 (GEZ2 in Globex noaion; i means ha he underlying insrumen is he December 2012 Eurodollar fuure). The firs rade day for EDZ2 is December 13, 2010, and he expiraion dae is December 17, We will esimae he value of several such pu and call opions during he period from December 20, 2011 o January 20, 2012, wih differen srike prices - ranging from o While he values of he American call opions could be compued wihou he use of he leas squares algorihm, because hey coincide wih he European calls, we will calculae hem anyway o provide more insigh ino how he parameers are fied o he marke daa. In oher words, we wish o check empirically if he differences beween he marke prices and he compued prices are he resul of badly fied model parameers or are due o a problem wih accuracy of he leas squares algorihm. For he calibraion purposes we will use he daily closing prices of Eurodollar fuures and he spo LIBOR rae. Given a dae, we will use a period of he same lengh as he ime o mauriy of he opion (i.e. if he opion lifeime is 300 days, hen we ake las 300 days daa before ime o calibrae our model). The leas squares algorihm needs several inpus. As he funcions generaing he informaion abou he pas we use he sandard exponenially weighed Laguerre polynomials of degree no greaer han 3. Our implemenaion is based on he Mone Carlo simulaion of he L 4 values obained using (4). The algorihm needs also formulas for he ineres rae (for he purpose of discouning) in wo insances. Firsly, o discoun he values of opions from one period o anoher (in he recursive sep-by-sep par). Secondly, o compue he final price of he opion (i.e. o discoun he opimal prices from every simulaion o ime T 0 = 0). While he 12

13 second ineres rae could be associaed wih sandard spo LIBOR rae, he firs one mus be based on he evoluion of asses (i.e. for every pah in he Mone Carlo run, one mus esimae separaely he spo rae a ime using he prices of Eurodollar conracs) Esimaion and numerical resuls. In his subsecion we presen deailed esimaion resuls for he dae December 20, A similar procedure has been also conduced for all remaining days under consideraion. Assuming he Brace-Gaarek-Musiela dynamics and aking ino accoun he Eurodollar fuures closing prices during he period from December 26, 2010 o December 20, 2011 we have conduced PCA and obained as he esimae of Λ. To price several pu and call opions wih differen srike prices and he closing dae falling on December 20, 2011, we have generaed 1000 Mone Carlo simulaions of size 10,000 and using he leas squares algorihm we have obained esimaed opion prices for differen srike prices. The resuls are presened in Table 1. The Mone Carlo disribuions of he prices of he Eurodollar pu and call opions wih he srike price can be seen in Figure 2. Figure 1 shows examples of 100 Mone Carlo pahs, ogeher wih he acual realizaion of he process. Similar analysis has been performed for all days from December 21, 2011 o January 20, During ha period, EDZ2 was he fourh closes o delivery Eurodollar Fuure. In Fig. 3 we can see he dynamics of he original pu and call opion prices, he sample means from 1000 simulaions (each of size 10,000) and he lower and upper 5% quaniles for he pu and call opions wih srike price The values of he mean and sandard deviaion of he simulaed prices of he opions, as well as he corresponding marke prices of he opions, can be seen in Table 2. We have chosen he srike price because he mean volume of ransacions was highes in he considered period. I is ineresing o noe ha he esimaed prices corresponding o his srike price say consisenly higher han he marke price (see Table 2), which migh be he resul of he fac ha he opion was paricularly acively raded (see Table 1, where he marke price is in mos cases lower han he esimaed price). I should also be noed ha he value of σ in Table 1 is highes for he srike price equal o 99.50, which may explain he ineres in he opion wih his paricular srike price. 4.2 Baske and dual-srike opions In his subsecion we will use he leas squares algorihm o price 1.5 monh baske and dualsrike American pu opions whose payoff funcions are based on wo marke indices, namely DAX and EUROSTOXX50. The laer will be denoed by he symbol EUR for breviy. We will assume ha he underlying insrumens follow he sandard bivariae Brownian dynamics. Unforunaely, bivariae opions are usually over-he-couner (OTC) insrumens, so i is difficul o find marke daa for such opions. Neverheless, we could do a parial comparison wih 13

14 price price day Figure 1: Examples of 100 Mone Carlo pahs for he L 4 conac (for December 20, 2011) and he realized pah (red) during he firs 100 and 300 days. day densiy Mean MC Marke price densiy Mean MC Marke price opion price Figure 2: The smoohed densiies of he simulaed prices of he pu (lef) and he call (righ) opion on December 20, 2011, wih srike price The disribuion is based on 1000 Mone Carlo runs, each of size 10,000. opion price he relevan one-dimensional sandard American pu opions based on DAX and EUR. As was he case wih he previous example, we sar wih some background informaion DAX and EUROSTOXX50 Indices. The univariae sandard American pu opions based on DAX and EUR are raded on he Eurex Exchange. In fac, he underlyings are no indices bu exchange-raded funds (ETF), which are acively raded on he German sock marke (Deusche Börse Group). The DAX and EUR indices are highly correlaed, chiefly due o he inclusion of some common socks in heir baskes. The esimaed value of Pearson s linear correlaion coefficien for he period from Ocober 23, 2012 o January 08, 2013 is equal o Some conagion beween hese indices migh poenially occur, bu in such a shor period of ime his aspec is negligible. In general, he issue of conagion could be addressed by adoping models wih differen dynamics (e.g. of he mulivariae GARCH variey). Such 14

15 Table 1: The esimaed prices of he Eurodollar opions on December 20, 2011, based on 1000 simulaions (each of size 10,000). Here µ denoes he sample mean of he 1000 prices obained wih MC simulaion, while σ denoes he sample sandard deviaion. Dae: Dec. 20, 2011 Pu Call Srike price Marke price µ σ Marke price µ σ pu opion price mean MC price Marke price MC quaniles call opion price mean MC price Marke price MC quaniles Dec 19 Dec 26 Jan 02 Jan 09 Jan 16 Dec 19 Dec 26 Jan 02 Jan 09 Jan 16 dae Figure 3: The esimaed prices (wih confidence level 90%) and hisorical prices of he pu (lef) and call (righ) opion prices wih srike price 99.50, during he period from December 20, 2011 o January 20, dae approach would be very closely relaed o he Heson and Nandi opion pricing model [12], which is he mehodology we will adop in he las example. 15

16 Table 2: Simulaed and hisorical prices of he Eurodollar opions wih srike price based on 1000 Mone Carlo simulaions (each of size 10,000). As before, µ denoes he sample mean, whereas σ denoes he sample sandard deviaion of he prices calculaed form he simulaions. Dae Pu Call Marke price µ σ Marke price µ σ Dec.20, Dec.21, Dec.22, Dec.23, Dec.27, Dec.28, Dec.29, Dec.30, Jan.03, Jan.04, Jan.05, Jan.06, Jan.09, Jan.10, Jan.11, Jan.12, Jan.13, Jan.17, Jan.18, Jan.19, Jan.20, Baske and dual-srike opions. As has been already saed, baske and dualsrike opions are mainly OTC derivaives. In his example we will consider a bivariae American pu opion. The payoff funcions a ime, for a bivariae baske American pu opion (1) and dual-srike American pu opion (2) is given by p (1) () = max ( K 1 S 1 (),K 2 S 2 (),0 ) (, p (2) K1 +K 2 () = max S ) 1()+S 2 (),0, 2 2 where S 1 () and S 2 () are he prices of he firs and he second underlying a ime, respecively, and K 1, K 2 are he srike prices Model seup, daa deails and implemenaion parameers. We will be assuming ha he price process(s 1 (),S 2 ()) is modeled by a 2-dimensional geomeric Brownian moion, wih he insananeous correlaion coefficien and insananeous sandard deviaions for he processes logs 1 and logs 2 denoed by ρ,σ 1 and σ 2, respecively. We will consruc 16

17 a bivariae baske and dual-srike American pu opions based on 1 DAX ETF share and 2.5 EUR ETF shares (o have similar srike prices in boh cases). We will price Baske and Dualsrike pu opions on January 08, 2013 wih he expiraion dae March 16, 2013 (o make i comparable o exising univariae opions). The opion lifeime will be 49 business days. The srike prices will range from 65 o 75 and from 66 o 76, for he firs and he second srike price, respecively. To esimae ρ,σ 1 and σ 2 we will use he las 50 observaions of he price of ETF (DE) DAX and ETF (DE) EUROSTOXX50. Choosing a relaively shor ime inerval for calibraion purposes is quie common in pracice (e.g. his is he case wih he esimaion of he VIX volailiy index). As in he previous case, we will need wo inpus for he leas squares algorihm: an ineres rae (for discouning) and appropriae basis funcions. Because he opion lifeime is shor, we will assume ha he ineres rae is consan and equal o r = 1.50% (he ECB ineres rae on January 08, 2013). Moreover, we will use he following exponenially weighed polynomials of wo variables o perform he regression: e 1 2, e x 2 x, e y 2 y, e (x+y) 4 xy, e (x+y) 4 xy 2, e (x+y) 4 x 2 y, e (x+y) 4 x 2 y Esimaion and numerical resuls. The esimaed (annualized) covariance marix gives us he values ρ = 0.920, σ 1 = and σ 2 = Using hese numbers, we run 100,000 Mone Carlo simulaions (each of size 49). Nex, using he leas squares algorihm we compue he prices of he baske and dual-srike American pu opions for differen srike prices. We also compue he leas squares prices for he univariae American pu opions based on 1 DAX ETF share and 2.5 EUR ETF share. Apar from he marke daa, we also presen he heoreical price according o he Cox-Ross-Rubinsein model (CRR) as i is used by he Eurex Exchange o quoe opion prices when no rading akes place. I should be noed ha he volume of ransacion of American pu opions is very low, so unforunaely he marke price is jus for comparison purposes. Also, he leas squares price should be compared wih he CRR price raher han he marke price (as i is compued under he compaible assumpions abou he asse dynamics). The prices (obained using single 100,000 Mone Carlo run) can be seen in Table 3. The columns wih names DAX and EUR denoe he sandard univariae pu opions (i.e. wih 1 DAX ETF and 2.5 EUR ETF share as he underlying, respecively). We have also performed muliple Mone Carlo runs (1000), each of size 10,000 for he baske and dual-srike opions wih he srike prices K 1 = K 2 = 70. The corresponding Mone Carlo densiy funcion could be seen in Figure 4 (his could provide some informaion abou he model and/or he Mone Carlo bias). 4.3 The Heson-Nandi model In he las example we will use he leas squares algorihm o price wo 1.5 monh American pu opions whose payoff is based on a single marke index. We will use daa from he previous example, i.e., we will price opions wrien on DAX and EUR indices. We will assume ha he 17

18 Table 3: Prices of he opions according o hisorical sock marke daa, he CRR model and he leas squares algorihm. Here S 0 = (68.05,69.72), r = 1.50%, T = 49/252, σ 1 = 0.133, σ 2 = 0.119, ρ = Srike price Marke price CRR price leas squares price EUR DAX EUR DAX EUR DAX EUR DAX Baske Dual-Srike dynamics of he underlying insrumens could be described wih he Heson-Nandi GARCH model [12]. Le S denoe he price of he underlying. Using he Heson-Nandi GARCH dynamics, we assume ha he log-reurns of he random process S could be described by formula logs = r daily +λσ 2 +σ ǫ, wih σ 2 = ω +βσ 2 1 +α(ǫ 1 γσ 1 ) 2, where denoes he daily backward difference, he parameer r daily denoes daily risk-free ineres rae, (λ,ω,β,α,γ) are model parameers and ǫ is he sandard Gaussian whie noise. In addiion, we will assume ha here is no asymmery in he model, i.e. γ = 0. 18

19 densiy ,000 MC price densiy ,000 MC price opion price opion price Figure 4: The smoohed densiies of he leas squares prices of he baske (lef) and he dualsrike (righ) American pu opions for he srike prices K 1 = K 2 = 70. The verical lines depic he sample mean of he leas squares prices. If we use he sandard Heson-Nandi dynamics (wih he objecive probabiliy measure) hen he discouning par of he leas squares algorihm will be pah dependen. In order o avoid his complicaion we will swich o he risk-neural measure and use he risk-neural dynamic of he underlying reurn. The risk neural process is obained simply by replacing (previously esimaed) parameers λ and γ wih ( 0.5) and (γ+λ+0.5), respecively (see [12] for deails). Moreover, we will use he long run expeced sandard deviaion from Heson Nandi model for comparison purposes (see [12]): ω +α σ HN = 1 β αγ2. (5) The EUR and DAX daa will be used again as he underlyings. As before, he opions expiraion dae will be March 16, 2013 and we will price hem on January 08, 2013 (hus he opion lifeime T will be 49 business days). The weighed Laguerre polynomials of degree no greaer han 3 will serve as he base polynomials for he regression procedure. As before, we will assume ha (annualised) risk free rae is equal o r = 1.50% and pu r daily = r/252 (as here are 252 rading days each year). Using he las 50 prices of 2.5 EUR ETF and 1 DAX ETF shares, we have obained wo ses of parameers: λ ω α β EUR DAX The iniial values of he underlying are and 69.72, respecively, in he EUR and DAX case. The (annualized) volailiies obained from (5) are equal o and 0.137, respecively. The mean sample prices of American pu opions obained from en simulaions (each consising of 100,000 Mone Carlo pahs) can be seen in Table 4. We also presen he heoreical European pu opion prices according o Heson-Nandi model [12], as well as he American pu 19

20 opions prices and early exercise premiums (i.e. he differences beween he prices of American and European pu opions) according o Cox-Ross-Rubinsein model (CRR), wih volailiies obained from (5). Boh models are presened for comparison purposes. Moreover, we perform muliple Mone Carlo runs (1000), each of size 10,000, o calculae prices of he American pu opions wih he srike price 70 (boh for EUR and DAX). Smoohed simulaed probabiliy densiy funcions are ploed in Figure 5. Table 4: Prices of he EUR and DAX American pu opions according o he leas squares algorihm (L-S), compared wih he acual marke prices, CRR model prices and he Heson- Nandi European pu opion prices. EA denoes he early exercise premium. EUR American pu opions Srike price Marke price CRR price CRR EA H-N price L-S price DAX American pu opions Srike price Marke price CRR price CRR EA H-N price L-S price densiy ,000 MC price densiy ,000 MC price opion price opion price Figure 5: The smoohed disribuions of he leas squares prices of he EUR (lef) and DAX (righ) American pu opions. The verical lines correspond o he sample mean of he leas squares prices. 20

21 5 Concluding remarks We have shown ha he widely used leas squares approach o Mone Carlo based pricing of American opions remains valid under very general and flexible choice of assumpions. In paricular, convergence o he heoreical price obained via Snell envelopes remains rue wih a highly adapable seup for approximaion of condiional expecaions. Of course one should be aware ha he compuaional cos of liberalizaion of he assumpions may be poenially very high. However, a growing body of empirical evidence indicaes ha in many pracical applicaions even relaively limied non-linear exensions of sandard regression may produce saisfacory resuls, as illusraed also by our hree examples. The relaxaion of he assumpions of he mehod should be seen primarily as increase in freedom of choice of seings for a specific implemenaion of he algorihm, which wih careful choices may neverheless reain compuaional viabiliy. Acknowledgmens: The second auhor acknowledges he suppor by Projec operaed wihin he Foundaion for Polish Science IPP Programme "Geomery and Topology in Physical Models" co-financed by he EU European Regional Developmen Fund, Operaional Program Innovaive Economy References [1] C.J. Beveridge and M.S. Joshi. Pracical policy ieraion: generic mehods for obaining rapid and igh bounds for Bermudan exoic derivaives using Mone Carlo simulaion. Universiy of Melbourne, Faculy of Business & Economics, Cenre for Acuarial Sudies, [2] A. Brace, D. Gąarek, and M. Musiela. The marke model of ineres rae dynamics. Mahemaical Finance, 7(2): , [3] M. Broadie and P. Glasserman. Pricing American-syle securiies using simulaion. Journal of Economic Dynamics and Conrol, 21: , [4] M. Broadie and P. Glasserman. A sochasic mesh mehod for pricing high-dimensional American opions. Technical repor, Columbia Business School, New York, A revised version of a working paper from [5] J.F. Carriere. Valuaion of early-exercise price for opions using simulaions and nonparameric regression. Insurance: Mahemaics and Economics, 19:19 30, [6] T.S. Chihara. An inroducion o orhogonal polynomials. Gordon and Breach, New York, [7] E. Clémen, D. Lamberon, and P. Proer. An analysis of a leas squares regression mehod for American opion pricing. Finance and Sochasics, 6: ,

22 [8] R.L. Dobrushin and R.A. Minlos. Polynomials in linear random funcions. Russian Mah. Surveys, 32:2:71 127, [9] D. Egloff. Mone Carlo algorihms for opimal sopping and saisical learning. Ann. Appl. Probab., 15(2): , [10] P. Glasserman. Mone Carlo mehods in financial engineering. Springer, New York, [11] P. Glasserman and B. Yu. Number of pahs versus number of basis funcions in American opion pricing. Ann. Appl. Probab., 14(4): , [12] S.L. Heson and S Nandi. A closed-form garch opion valuaion model. Review of Financial Sudies, 13(3): , [13] F. Jamshidian. Libor and swap marke models and measures. Finance and Sochasics, 1 (4): , [14] M. Klimek, M. Masuura, and Y. Okabe. Sochasic flows and finie block frames. Journal of Mahemaical Analysis and Applicaions, 342(2): , [15] D. Lamberon and B. Lapeyre. Inroducion o sochasic calculus applied o finance. CRC Press, [16] G. Liu and L.J. Hong. Revisi of sochasic mesh mehod for pricing American opions. Operaions Research Leers, 37(6): , [17] F.A. Longsaff and E.S. Schwarz. Valuing American opions by simulaion: a simple leas squares approach. Review of Financial Sudies, 14: , [18] S. Pliska. Inroducion o Mahemaical Finance - Discree Time Models. Blackwell Publishers, Oxford, [19] J. Rus. Using randomizaion o break he curse of dimensionaliy. Economerica, 65(3): , [20] C.G. Small and D.L. McLeish. Hilber space mehods in probabiliy and saisical inference. Wiley, New York, [21] L. Senof. Convergence of he leas squares Mone Carlo approach o American opion valuaion. Managemen Sci., 50(9): , [22] J. N. Tsisiklis and B. Van Roy. Opimal sopping of Markov processes: Hilber space heory, approximaion algorihms, and an applicaion o pricing high-dimensional financial derivaives. IEEE Trans. Auoma. Conrol, 44(10): , [23] J. N. Tsisiklis and B. Van Roy. Regression mehods for pricing complex American-syle opions. IEEE Trans. on Neural Neworks, 12(4): , [24] D. Z. Zanger. Convergence of a leas-squares Mone Carlo algorihm for bounded approximaing ses. Appl. Mah. Finance, 16(1-2): ,