How to Model a Financial Bubble Mathematically Lecture 1 Philip Protter, Columbia University Statistical Methods and Applications in Insurance and Finance CIMPA School April 12, 2013 A Little History • Stock markets date back to at least 1531 in Antwerp, Belgium • There are over 150 stock market exchanges world wide, of which the most significant count 103: • There are 34 in Europe, including 5 in the U.K. and 4 in France • There are 20 in North America (7 in Canada, 12 in the U.S., and 1 in Mexico) • There are 3 principal stock exchanges in Africa: Johannesburg, Lagos, and Casablanca Why do Stock Markets Exist? • In the U.S., for example, the railroads needed vast amounts of capital to build their tracks, and created the need for a stock exchange • The Dow Jones Industrial Average officially began in 1896 • In 1884, 12 years earlier, its predecessor began: Customer’s Afternoon Letter which contained 11 stocks, 9 of which were railroads • In 1885, there were 12 railroads and 2 industrials in the Dow Jones letter • In 1886, there were 10 railroads and 2 industrials in the Dow Jones letter The Dow Jones Industrial Average in January, 1896 *American Sugar Chicago, Milwaukee & St. Paul Chicago, Rock Island & Pacific Delaware, Lackawanna & Western Missouri Pacific Union Pacific Chicago, Burlington & Quincy Chicago & North Western Delaware & Hudson Canal Louisville & Nashville Northern Pacific preferred *Western Union *Indicates an industrial (not a railroad) The initial Dow Jones Industrial Average without Railroads (May 26, 1896) American Cotton Oil American Tobacco Distilling & Cattle Feeding Laclede Gas North American U.S. Leather preferred American Sugar Chicago Gas General Electric National Lead Tennessee Coal & Iron U.S. Rubber North American was replaced by US. Cordage Preferred, and Distilling & Cattle Feeding became American Spirits, in August, 1896 Basic Mathematical Models for Asset Pricing Finance • Let S = (St )0≤t≤T represent the (nonnegative) price process of a risky asset (e.g., the price of a stock, a commodity such as “pork bellies,” a currency exchange rate, etc.) • The present is often thought of as time t = 0. One is interested in the unknown price at some future time T , and thus ST constitutes a “risk.” • Example: An American company contracts at time t = 0 to deliver machine parts to Germany at time T . Then the unknown price of Euros at time T (in dollars) constitutes a risk for that company. • In order to reduce this risk, one may use “derivatives”: one can purchase — at time t = 0 — the right to buy Euros at time T at a price that is fixed at time 0, and which is called the “strike price.” • If the price of Euros is higher at time T , then one exercises this right to buy the Euros, and the risk is removed. This is one example of a derivative, called a call option. • More generally, a derivative is any financial security whose value is derived from the price of another asset, financial security, or commodity. Call and Put Options • A call option with strike price K , and payoff at time T can be represented mathematically as C = (ST − K )+ where x + = max(x, 0). • Analogously, the payoff to a put option with strike price K at time T is P = (K − ST )+ and this corresponds to the right to sell the security at price K at time T . • Calls and Puts are related, and we have ST − K = (ST − K )+ − (K − ST )+ . This relation is known as put – call parity. More complicated simple options • We can calls and puts as building blocks for more complicated derivatives. • For example, if V = max(K , ST ) then V = ST + (K − ST )+ = K + (ST − K )+ . • More generally, if f : R+ → R+ is convex then f (x) = f (0) + f+0 (0)x + Z ∞ (x − y )+ µ(dy ) (1) 0 where f+0 (x) is the right continuous version of the (mathematical) derivative of f , and µ is a positive measure on R with µ = f 00 , where the mathematical derivative is in the generalized function sense. Thus if f is convex, and if V = f (ST ) is our financial derivative, then V is a portfolio consisting of a continuum of European call options: Z ∞ V = f (0) + f+0 (0)ST + (ST − K )+ µ(dK ). 0 Other kinds of derivatives • We can also have path dependent derivatives. V = F (S)T = F (St ; 0 ≤ t ≤ T ) which are functionals of the paths of S. • For example if S has c` adl`ag paths (c`adl`ag is a French acronym for “right continuous with left limits”) then F : D → R+ , where D is the space of functions f : [0, T ] → R+ which are right continuous with left limits. The time value of money • Inflation makes money worth less as time goes on • Deflation makes it worth more • Evaluating a claim that pays off $D at time T , when current time is zero, can be done in time T dollars, or in time 0 dollars; if we use time T dollars for the payoff, but time 0 dollars for the evaluation, we must discount the payoff by the rate of inflation (deflation) • Suppose we have $D at time 0, and invest it in a bank which pays interest rate r for one time unit (eg, one year). After one year, we have $(D + rD). • If we are paid interest every 3 months , or 1/4 year, and leave the interest in the bank, we have $D + Dr /4 after the first quarter, $D(1 + r /4)2 after the second, and $D(1 + r /4)4 after one year. • If we compound n times in one year and leave the money in the bank, we have $D(1 + r /n)n • Taking limits limn→∞ $D(1 + r /n)n = $De r ; for t time units analogously the limit = $De rt , which solves the ODE dRt = r; Rt R0 = D • In general if r is a stochastic process (rt )t≥0 , then Z Rt = D + t rs Rs ds 0 ⇒ Rt = De Rt 0 rs ds . A simple Portfolio • A simple portfolio has a varying quantity of shares of a stock, plus a varying amount of money in a liquid, risk-free money account. • The value of a portfolio, V , depends on the trading strategy a for stocks, and b for the money account • A trading strategy is a vector of stochastic processes (a, b) • Following a strategy (a, b) gives a dynamic portfolio value process: Vt (a, b) = at St + bt Rt . • A trading strategy (a, b) is called self-financing if Z at St + bt Rt = a0 S0 + b0 R0 + t Z as dSs + 0 t bs dRs 0 Comments on self-financing Z at St + bt Rt = a0 S0 + b0 R0 + t Z as dSs + 0 t bs dRs (2) 0 • Intuitively Self-financing means that we do not consume money for other purposes, or add new money; we will soon give a heuristic justification of equation (2) • S is taken, by assumption, to have sample paths which are right continuous and have left limits (c`adl`ag), and R is continuous; hence the right side of (2) is at least c`adl`ag • This creates implicit restrictions on the illusory arbitrariness of the choice of a (predictable) and b (right continuous) • If r ≡ 0 then (Rt )t≥0 ≡ 1, hence dRt = 0 and (2) becomes Z t at St + bt = a0 S0 + b0 + as dSs (3) 0 • This means once we have chosen strategy a, then b is determined Heuristic justification of self-financing • Suppose a, b, S are all three semimartingales, and that R ≡ 1. Then we have: (at+dt − at )St+dt = −(bt+dt − bt ) (4) which says that the change in stock holdings creates a corresponding change in the money account. • Equation (4) becomes (at+dt − at )(St+dt − St ) + (at+dt − at )St = −(bt+dt − bt ) ≈ d[a, S]t + St− dat = −dbt (5) • By Integration by parts, (5) becomes Z at St = a0 S0 + b0 + t t Z as dSs + 0 Ss− das + [a, S]t 0 ⇒ d(at St ) − at− dSt = −dbt Z t ≡ at St + bt = a0 S0 + b0 + as dSs . 0 (6) What is Arbitrage? • In language: Arbitrage is the chance, no matter how small, to make a profit without taking any risk • Definition A model is arbitrage free on [0, 1] if there does not exist a self-financing strategy (a, b) such that V0 (a, b) = 0, VT (a, b) ≥ 0, P(VT (a, b) > 0) > 0. (7) • We want to convert this idea into useful mathematics • Folk Theorem: There is no arbitrage if and only if there exist a new probability Q, equivalent to P (ie, same sets of probability zero, written Q ∼ P), such that S is a martingale. • The above folk theorem is based on it being true in simple cases (eg, finite probability space Ω [J.M. Harrison & S.R. Pliska]) Martingales, local martingales, and sigma martingales • We assume given a complete, filtered probability space (Ω, F, F, P), where F = (Ft )t≥0 • A stochastic process M is a martingale if E (|Mt |) < ∞, and for s ≤ t, E (Mt |Fs ) = Ms a.s. • Martingales are insufficient; for example: • If X is a submartingale, we want a decomposition of X = M + A, where M is a martingale and A is an increasing, predictably measurable process. This is not true in general, instead we need the concept of local martingale. • If R tN is a martingale, we would like the stochastic integral Hs dNs to be a martingale, too. This is not true in general, 0 but instead (if X has continuous paths) it is a local martingale. • In R t general, if N is a martingale, then the stochastic integral Hs dNs is a σ martingale 0 Definitions • A stochastic process X with X0 = 0 is a local martingale if there exists a sequence of stopping times (Tn )n≥1 with limn→∞ Tn = ∞ a.s., such that Xt∧Tn is a martingale for every n ≥ 1 • A stochastic process X with X0 = 0 is a σ martingale if there exists Ra martingale M and a predictable process H such that t Xt = 0 Hs dMs for all t ≥ 0 • Note: Martingales ⊂ Local Martingales ⊂ σ Martingales • If X is a nonnegative (or just bounded from below) σ martingale, then it is a local martingale. So X ≥ 0 ⇒ Local Martingales = σ Martingales • Stochastic integration is closed for σ martingales • For continuous processes, stochastic integration is closed for local martingales One way local martingales can arise • Let X be the unique weak solution of the stochastic differential equation dXt = σ(Xt )dBt , X0 = 1 where B is a standard Brownian motion • (Blei-Engelbert) If there exists an α ∈ (0, 1) such that Z ∞ α 1 dy < ∞ σ(y )2 then X is a local martingale, and not a martingale. We call such a process a strict local martingale The Canonical Example of a Local Martingale • Let Bt = (Bt1 , Bt2 , Bt3 ) be standard 3 dimensional Brownian motion, with B0 = (1, 0, 0). • Let u : R3 \{(0, 0, 0)} → R+ be given by u(x) = 1 kx k • Xt = u(Bt ) is a positive real valued local martingale, with E (X0 ) = 1; X is called the inverse Bessel process • X is not a martingale, because one can show that lim E (Xt ) = 0 t→∞ and therefore it is not constant • The inverse Bessel process satisfies the SDE dXt = −(Xt )2 dBt , X0 = x0 > 0 The Canonical Example of a σ martingale • τ is an exponential r.v. with parameter λ = 1 • U is independent of τ and P(U = 1) = P(U = −1) = 12 • Xt = U1{t≥τ } ; then X is a martingale Rt • Let Hs = 1s for s > 0, and let Mt = 0 Hs dXs • Note that M has unbounded positive and negative jumps • E (|Mν |) = ∞ for every stopping time ν with P(ν > 0) > 0, so M is not a martingale, and not a local martingale, but M is in fact a σ martingale. Semimartingales and arbitrage • Suppose S has continuous paths and is a semimartingale with decomposition St = S0 + Mt + At , with M0 = A0 = 0, and Q ∼ P; take dQ |Ft ) Z t = EP ( dP which is a martingale • By Girsanov’s theorem the decomposition of S under Q is given by Z St = (Mt − 0 t 1 d[Z , M]s ) + (At + Zs− Z 0 t 1 d[Z , M]s ); Zs− • Therefore if Z can be chosen so that Z t At = − 0 1 d[Z , M]s , Zs− we have that S is a Q-local martingale. • By the Kunita-Watanabe inequality, from (8) we have d[Z , M]t d[Z , Z ]t d[M, M]t (8) • Recall (8): Z At = − 0 t 1 d[Z , M]s , Zs− therefore we must have that dAt d[M, M]t in order for M to be martingale or local martingale under Q. • This is not always the case; for example by Tanaka’s formula, if Z St = 1 + |Bt | = 1 + t sign(Bs )dBs + L0t 0 = 1 + βt + L0 (B)t , then d[β, β]t = dt, but dL0t 6 dt. • Therefore 6 ∃ Q ∼ P for (9) such that S is a Q (local) martingale (9) What if S is continuous and not a semimartingale? • If there exists a Q ∼ P such that S is a local martingale under dP |Ft ) Q, then let Yt = EQ ( dQ • By Girsanov, Z St = (St + 0 t 1 d[Y , S]s ) − Ys Z 0 t 1 d[Y , S]s Ys is a P decomposition of S. Therefore S is a P semimartingale, a contradiction • Thus a necessary condition for Q ∼ P is that S be a P semimartingale Why are (local) martingales so important? • Martingales model fair gambling games • A price process which is a model under the risk neutral measure should have constant expectation • Martingales have the property that t 7→ E (Mt ) is constant • Theorem A stochastic process X is a martingale if and only if E (Mτ ) = E (M0 ) for every bounded stopping time τ . • Thus, M has constant expectation not just for fixed times, but for stopping times as well. The First Fundamental Theorem of Asset Pricing • First Version: J. M. Harrison and S. R. Pliska, circa 1979 showed that a finite probability space (Ω, F, (Sn )n=0,1,2,... , P) has No Arbitrage if and only if there exists another probability measure Q ∼ P such that S is a martingale • Second Version: David Kreps, circa 1981 realized that No Arbitrage was not a strong enough condition to guarantee such a result in a more general case. He created a new condition and called it No Free Lunch • Ignoring admissibility conditions for now, Kreps said that S admits a Free Lunch on [0, T ] if there exists a function f ∈ L∞ + (Ω, F, P) such that P(f > 0) > 0, and a net RT (fα )α∈I = (gα − hα )α∈I , with hα ≥ 0 and gα = 0 Hsα dSs , for admissible H α . And also fα → f in the Mackey topology on L∞ induced by L1 • The Mackey topology is often written as σ(L∞ , L1 ), which means that for a sequence (Xn )n≥1 ∈ L∞ , then Xn → X , if for any Y ∈ L1 , E (Xn Y ) → E (XY ). Economic intuition of No Free Lunch RT • Often we think of f as being of the form f = 0 Hs dSs • Kreps saw that f could not in general be restricted to this form for an admissible process H. (If it were, one could follow this trading strategy H and replicated f , and have classical arbitrage [starting with 0 and ending with f ≥ 0]) • But suppose f can be approximated by (fα )α∈I in a suitable topology • Let (hα )α∈I be the “errors” in the approximation, representing “money thrown away.” • No Free Lunch does not allow arbitrage, but it does allow arbitrage to exist in the limit Kreps’ Theorem Theorem (Kreps, 1981) A bounded process S = (St )0≤t≤T admits NFL if and only if there exists Q ∼ P such that S is a martingale under Q. This creates three immediate questions: 1. Can we replace [0, T ] with [0, ∞)? 2. What if S is not bounded? 3. What does convergence in nets mean vis `a vis an economics interpretation? The Four Fundamental Papers that Clarified the Issues Surrounding the First Fundamental Theorem 1. Harrison, J.M, Kreps, D.M. (1979) Martingales and Arbitrage in Multiperiod Securities Markets, Journal of Economic Theory 20, 381-408 2. Harrison, J.M, Pliska, S.R. (1981) Martingales and Stochastic Integrals in the Theory of Continuous Trading, Stochastic Processes and their Applications 11, 215-260 3. Kreps, D.M. (1981) Arbitrage and Equilibrium in Economics with infinitely many Commodities, Journal of Mathematical Economics 8, 15-35 4. Harrison, J.M, Pliska, S.R. (1983) A stochastic calculus model of continuous trading: Complete markets, Stochastic Processes and their Application 11, 313-316 13 years later: Delbaen and Schachermayer • Delbaen and Schachermayer, 1994: Convergence with nets is replaced with convergence of sequences; S bounded is replaced with S locally bounded, and M a martingale is replaced with M a local martingale • Delbaen and Schachermayer, 1998: The general case is treated, where S can be c`adl`ag, and does not have to be locally bounded, and M is replaced with a σ martingale. • Before we discuss these results, we need the concept of an admissible trading strategy The Doubling Strategy • Bet $1 at even money • Stop betting if you win and collect $1 net winnings; otherwise bet again, waging $2 • Stop if you win; you have now lost $1 and won $2, for a profit of $1; otherwise bet again, waging $4 • In general: stop whenever you win, otherwise bet again, doubling your last bet; your net winnings will be $1 • The probability is 1 that you will eventually win $1, so this is an arbitrage strategy, known as the doubling strategy Problems with the Doubling Strategy • Need to make an unlimited number of bets (time constraints) • Need “no fees” to make such bets (transaction costs) • Need to have a counterparty (liquidity) • But the above are practical problems; a theoretical problem is the need for infinite resources • We can eliminate the doubling strategy with an admissibility condition Admissibility Definition: Let S be a semimartingale, α R> 0. A predictable t process H is α-admissible if H0 = 0, and 0 Hs dSs ≥ −α, for all t ≥ 0. H is admissible if there exists an α > 0 such that H is α-admissible. Note: • We are implicity assuming That if H is admissible it is predictably measurable and is in the space of S-integrable processes • This condition of admissibility is intrinsically asymmetric: H can increase without bound, but is strictly limited in how much it can be negative The Kreps-Delbaen-Schachermayer Theory • We work on the semi-infinite time interval [0, ∞], on a filtered complete probability space (Ω, F, F, P), where F = (Ft )t≥0 . • We further assume we have a risky asset price process S = (St )t≥0 and that the spot interest rate r = 0 • A Contingent Claim is simply an FT measurable random variable; examples are C = (ST − K )+ , which is a call at strike price K and maturity time T ; another example is P = (K − ST )+ which is a Rput We let H · S denote the t stochastic integral process ( 0 Hs dSs )t≥0 • Note that a call C is unbounded if S is unbounded, but a put P is always bounded, irrespective of the behavior of S Definitions • We let L0+ denote finite-valued, nonnegative random variables (a.s.). We define K = {(H · S)∞ |H is admissible} Kα = {H · S)∞ |H is α − admissible} • No Arbitrage (NA): K ∩ L0+ = {0} • Intuition: Starting with nothing, the only nonnegative result we can end up with is identically 0; i.e., nothing • Next we define A0 = K − L0+ = {X = H · S)∞ |H is admissible, f ≥ 0, finite} A = A0 ∩ L∞ = {|X | ≤ k, some k : X = (H · S)∞ − f } • No Free Lunch (NFL) [Kreps]: A¯M ∩ L∞ + = {0}, where the ¯ M denotes closure in the Mackey topology σ(L1 , L∞ ) (·) • No Free Lunch with Vanishing Risk (NFLVR) [Delbaen-Schachermayer]: A¯ ∩ L∞ + = {0}, where the ∞ closure of A is in L , that is, the a.s. sup norm, as opposed to the Mackey closure of Kreps and NFL • Theorem: NFLVR is invariant under a change to an equivalent probability measure • NFLVR has become the accepted definition of no arbitrage; it is considered to be the “gold standard.” • However, we will see when we consider bubbles, that NFLVR is just a bit too weak. • The idea of No Dominance was introduced by Robert Merton in 1973, but largely forgotten No Dominance • Let P(S) be all probabilities equivalent to the underlying probability P such that if Q ∈ P(S) then S is a Qσ martingale. Let J = {J ∈ FT |J is bounded from below and sup EQ (S) < ∞} Q∈P(S) Λ(J)t = {the market price at time t of the contingent claim J} • Definition: An element D of J Q-dominates another element C of J if there exists a time t < T such that C − Λ(C )t Q{C − Λ(C )t ≤ D − Λ(D)t , for all t ≥ 0, Q a.s., and < D − Λ(D)t } > 0 for some t ≥ 0 • We say that the model has No Dominance (ND) under P if for any contingent claim C ∈ P(S), there does not exist another claim D in P(S) which dominates C • Theorem: If No Dominance holds for one Q ∈ P(S), then it holds for Q ∈ P(S) • Theorem: If for any H ∈ A we have Λ((H · S)T )0 = 0, then No Dominance implies (NA). • Theorem: If for any H ∈ A we have Λ((H · S)T )0 = 0 and Λ is lower semicontinuous on L∞ with the k · k norm, then No Dominance implies (NFLVR) End of Lecture 1 How to Model a Financial Bubble Mathematically Lecture 2 Philip Protter, Columbia University Statistical Methods and Applications in Insurance and Finance CIMPA School April 12, 2013 • We have a stochastic process (St )t≥0 modeling the price of a risky asset; S exists on a complete, filtered probability space (Ω, F, P, F), where F = (Ft )t≥0 • We discussed that S needs to be a semimartingale: a semimartingale is a process S that has a decomposition St = S0 + Mt + At , where M0 = A0 = 0 a.s., and M is a local martingale and A is a c`adl` gprocess with paths of finite variation on compact time intervals; a local martingale M is a process such that there exists an increasing sequence of stopping times (Tn )n≥0 with T0 = 0 a.s. and limn→∞ Tn = ∞ a.s. • Folk Theorem: There are no arbitrage opportunities if and only if there is an equivalent probability measure Q such that under Q the risky asset price process S is a martingale. • We let L0+ denote finite-valued, nonnegative random variables (a.s.). We define K = {(H · S)∞ |H is admissible} Kα = {H · S)∞ |H is α − admissible} • No Arbitrage (NA): K ∩ L0+ = {0} • Intuition: Starting with nothing, the only nonnegative result we can end up with is identically 0; i.e., nothing • Next we define A0 = K − L0+ = {X = H · S)∞ |H is admissible, f ≥ 0, finite} A = A0 ∩ L∞ = {|X | ≤ k, some k : X = (H · S)∞ − f } • No Free Lunch (NFL) [Kreps]: A¯ ∩ L∞ + = {0}, where the ¯ denotes closure in the Mackey topology, that is σ(L∞ , L1 ) (·) • No Free Lunch with Vanishing Risk (NFLVR) [Delbaen-Schachermayer]: A¯ ∩ L∞ + = {0}, where the ∞ closure of A is in L , that is, the a.s. sup norm, as opposed to the Mackey closure of Kreps and NFL • Theorem: NFLVR is invariant under a change to an equivalent probability measure • NFLVR has become the accepted definition of no arbitrage; it is considered to be the “gold standard.” • However, we will see when we consider bubbles, that NFLVR is just a bit too weak. • The idea of No Dominance was introduced by Robert Merton in 1973, but largely forgotten No Dominance • Let P(S) be all probabilities equivalent to the underlying probability P such that if Q ∈ P(S) then S is a Qσ martingale. Let J = {J ∈ FT |J is bounded from below and sup EQ (S) < ∞} Q∈P(S) Λ(J)t = is the market price at time t of the contingent claim J • Definition: An element D of J Q-dominates another element C of J if there exists a time t < T such that C − Λ(C )t Q{C − Λ(C )t ≤ D − Λ(D)t , for all t ≥ 0, Q a.s., and < D − Λ(D)t } > 0 for somet ≥ 0 We say that the model has No Dominance (ND) under P if for any contingent claim C ∈ P(S), there does not exist another claim D in P(S) which dominates C • Theorem: If No Dominance holds for one Q ∈ P(S), then it holds for Q ∈ P(S) • Theorem: If for any H ∈ A we have Λ((H · S)T )0 = 0, then No Dominance implies (NA). • Theorem: If for any H ∈ A we have Λ((H · S)T )0 = 0 and Λ is lower semicontinuous on L∞ with the k · k norm, then No Dominance implies (NFLVR) A General Framework • Let S be a semimartingale modeling a risky asset price process (so S ≥ 0) • Assume that NFLVR holds • Let D = (Dt )(t ≥ 0 be its cumulative cash flow of dividends • Assume the spot interest rate r ≡ 0 • We assume that the risky asset has a (finite valued) maturity, or lifetime, of the risky asset • Let Xτ = the terminal payoff, or liquidation value at time τ • Assume that Xτ ≥ 0 and Dt ≥ 0 • ∆Dt = Dt − Dt− , and if ∆Dt0 > 0, for some t0 , then St0 denotes the price ex dividend • Ex dividend refers to the trading of shares when a declared dividend belongs to the seller, rather than to the buyer. The Wealth Process • The wealth of the investor at time t is given by Z τ ∧t Wt = St + dDu + Xτ 1{t≥τ } 0 • We assume there exists a probability measure Q ∼ P such that W is a Q-local martingale • Note that since W ≥ 0 always, we do not have need of σ martingales Trading Strategies • A trading strategy is a vector process (πt , ηt )t≥0 • (πt )t≥0 is the trading strategy for the risky asset • (ηt )t≥0 is the (risk-free) money market trading strategy • Wtπ = πt St + ηt is the Wealth process corresponding to the strategy (π, η) Self-Financing Trading Strategies • W0π = 0 • The Value Process V corresponding to the strategy (π, η) is given by Vtπ,η Z = t πu dWu 0 • Let α > 0. A self-financing strategy π is α-admissible if Vtπ,η ≥ −α. The strategy π is admissible if it is α-admissible for some α > 0 The Case of Complete Markets • A market is complete for a class of contingent claims H if every claim X ∈ H can be perfectly hedged • How to interpret this statement in mathematics? • First we must discount for the time value of money; usually we work in a finite horizon, ie, a time interval [0, T ] • We can take, for example, H to be all bounded X ∈ FT ; or all random variables X in FT such that X /RT is bounded, RT where RT = 1 + 0 Rs rs ds, and r is the spot interest rate process • The Second Fundamental Theorem of Finance: A market under H is complete if and only if every X ∈ H can be perfectly hedged. That is, for any X ∈ H there exists a hedging strategy π such that Z X =α+ T πs dWs 0 This is equivalent to there being only one equivalent probability measure Q such that W is a Q local martingale. • The unique equivalent measure Q that turns W into a local martingale is called the risk neutral measure • Recall that we are operating under the assumption of NFLVR, so we know that at least one such Q exists • The price of such a claim X is now intuitively clear: if RT 1 X ∈ LQ (FT ) and 0 πs dWs is a martingale, then the price should be EQ (X ) = α The Problem of Unique Prices • A suicide strategy is a strategy σ such that W0σ = 1 but WTσ = 0 • Suicide strategies can lead to the non-uniqueness of prices • For example let Y be a contingent claim that suppose θ is a strategy such that Y = c + WTθ • The fair price of Y should be EQ (Y ) = c. But if one adds the suicide strategy σ, one gets Y = c + 1 + WTθ+σ which by the same reasoning should have price c + 1. Eliminating suicide strategies • The problem arises because a local martingale Z need not have constant expectation, and can even have ZT = 0 at a finite time T • Harrison and Pliska eliminate this possibility be restricting themselves to martingales, at the cost of generality • One can also eliminate suicide strategies with Merton’s No Dominance assumption, and then allow local martingales (recall that No Dominance implies NFLVR) NFLVR does not imply ND • Consider two risky assets maturing at time τ with payoff Xτ and Yτ , respectively. Suppose Xτ ≥ Yτ a.s. Then: Xt? = EQ {Xτ |Ft }1{t<τ } ≥ EQ {Yτ |Ft }1{t<τ } = Yt? • Let β be a nonnegative local martingale such that βτ = 0 a.s., and βt0 > Xt?0 − Yt?0 for some t0 • Such a non-trivial β exists, and it is unbounded, of necessity (if it were bounded, it would be a nonnegative [true] martingale with terminal value 0, and hence identically 0) • Suppose next the two risky asset prices are Xt = Xt? Yt = βt + Yt? • No Dominance is violated, because Yt0 > Xt0 for some t0 , Xτ ≥ Yτ • But NFLVR is not violated. Why? • Suppose we have a strategy designed to take advantage of this mis-pricing: Sell Y and buy X , say at time t0 , and then hold it to maturity. This gives a gain of βt0 > 0 at time t0 , and since Xτ − Yτ ≥ 0, with no other cash flows, we are safe to keep our gain βt0 ; this creates an apparent arbitrage • However, to do this, if t0 < u ≤ τ , then the value of the trading strategy is −Yu + Xu = βu + (Xu? − Yu? ) • Since β is of necessity unbounded, the admissibility condition eliminates this strategy of apparent arbitrage; doubling all over again! What is the difference between S being a martingale under the risk neutral measure, and a local martingale? • It turns out that this nuance can be used to give an explanation of financial bubbles • But first, let us describe the popular conception of a financial bubble US Stock Prices 1929 (Donaldson & Kamstra [1996]) Important Recent Bubbles • Minor crashes in the 1960s and 1980s • Junk bond financing led to the major crash of 1987 • Japanese housing bubble circa 1970 to 1989 • The “dot com” crash, from March 11th, 2000 to October 9th, 2002. Led by speculation due to the promise of the internet; The Nasdaq Composite lost 78% of its value as it fell from 5046.86 to 1114.11. • The 2008 US housing bubble and subprime mortgages NASDAQ Index 1998-2000 (Brunnermeir & Nagel) Current US Housing Price Trend (Center for Responsible Lending) Oil Futures (WTRG Economics) Oil Futures (WTRG Economics) Oil Futures (WTRG Economics) Oil Futures (WTRG Economics) Oil Futures (WTRG Economics) Oil Futures (WTRG Economics) We can also see the bubble in the entire world’s stock markets bursting. The next slide gives the results from the Morgan Stanley Capital International All Country World Index. Global Stock Markets (Greed and Fear, Oct 9, 2008) Our Framework for Modeling Bubbles • We assume NFLVR • We have a risky asset price process S with S ≥ 0, and the spot rate rt ≡ 0, so Rt ≡ 1 • We have a dividend process D, a wealth process W , and a liquidation vale Xτ • We are in a complete market framework, with a unique risk neutral measure Q • a fair price to pay for the risky asset at time t is the conditional expectation of the future cash flows (taken under the risk neutral measure) The Fundamental Price In complete markets with a finite horizon T , we use the risk neutral measure Q, and for t < T the fundamental price of the risky asset is defined to be: St? Z = EQ { T dDu + XT |Ft } t Definition (Bubble) A bubble in a static market for an asset with price process S is defined to be: β = S − S? Static Markets Theorem (Three types of bubbles) 1. β is a local martingale (which could be a uniformly martingale) if P(τ = ∞) > 0; 2. β is a local martingale but not a uniformly integrable martingale, if it is unbounded, but with P(τ < ∞) = 1; 3. β is a strict Q local martingale, if τ is a bounded stopping time. • Type 1 is akin to fiat money • Type 2 is tested in the empirical literature • Type 3 is essentially “new.” Type 3 are the most interesting! • Fiat money (Type 1 bubbles) has no intrinsic value, and become almost worthless • German hyper inflation of the 1920s • Zimbabwean hyperinflation of the current day 100 million German marks, 1923 100 trillion Zimbabwe dollars, 2009 Hyperinflation in Zimbabwe In 1980, 1ZWR=US $1 Month ZWR per USD Sept 2008 1, 000 Oct 2008 90, 000 Nov 2008 1, 200, 000 Mid Dec 2008 60, 000, 000 End Dec 2008 2, 000, 000, 000 Mid Jan 2009 1, 000, 000, 000, 000 2 Feb 2009 300, 000, 000, 000, 000 Zimbabwean Currency Follow-up • The Zimbabwe dollar was abandoned officially on April 12, 2009 • Since then it has been legal to use other currencies • Most common are reported to be the South African Rand, the Botswana Pula, The British pound sterling, and the US dollar Theorem (Bubble Decomposition) The risky asset price admits a unique decomposition S = S ? + (β 1 + β 2 + β 3 ) where 1. β 1 is a c`adl`ag nonnegative uniformly integrable martingale with limt→∞ βt1 = X∞ a.s. 2. β 2 is a c`adl`ag nonnegative NON uniformly integrable martingale with limt→∞ βt2 = 0 a.s. 3. β 3 is a c`adl`ag non-negative supermartingale (and strict local martingale) such that limt→∞ E {βt3 } = 0 and limt→∞ βt3 = 0 a.s. Ben Bernanke and the Federal Reserve End of Lecture 2 How to Model a Financial Bubble Mathematically Lecture 3 Philip Protter, Columbia University Statistical Methods and Applications in Insurance and Finance CIMPA School April 13, 2013 Review from Lecture 2 • Let S be a semimartingale modeling a risky asset price process (so S ≥ 0) • Assume that NFLVR holds • Let D = (Dt )t≥0 be its cumulative cash flow of dividends • Assume the spot interest rate r ≡ 0 • Let Xτ = the terminal payoff, or liquidation value at time τ • The wealth of the investor at time t is given by Z τ ∧t Wt = St + dDu + Xτ 1{t≥τ } 0 • We assume there exists a probability measure Q ∼ P such that W is a Q-local martingale • A trading strategy is a vector process (πt , ηt )t≥0 • W0π = 0 • The Value Process V corresponding to the strategy (π, η) is given by Vtπ,η = t Z πu dWu 0 • Let α > 0. A strategy π is α-admissible if Vtπ,η ≥ −α. The strategy π is admissible if it is α-admissible for some α > 0 • The Second Fundamental Theorem of Finance: A market under H is complete if and only if for every X ∈ H there exists a hedging strategy π such that Z X =α+ T πs dWs 0 This is equivalent to there being only one equivalent probability measure Q such that W is a Q local martingale. The Fundamental Price in a Complete Market Setting • Since markets are assumed complete, let Q ∼ P be the unique risk neutral measure • We define the fundamental price of the risky asset S, denoted S ? , to be the future discounted future cash flow one expects to get, conditional on current information • In mathematics, S ? is given by St? Z = EQ { T dDs + Xτ 1{τ <∞} |Ft }1{t<τ } t • The payoff at time t = ∞ does not contribute to S ? Theorem: St? is well defined. Moreover, lim S ? t→∞ t = 0 a.s. • Observe that W is a nonnegative Q supermartingale, so St? ∈ L1 (dQ), and the result follows the supermartingale convergence theorem, and the facts that (Dt )t≥0 and Xτ are nonnegative. • Note that in contrast, we cannot assume that St? is in L1 (dP) • Corollary: Z t∧τ ? ? Wt = St + dDu + Xτ 1{τ ≤t} 0 is a uniformly integrable martingale under Q, and Z τ ? W∞ = dDu + Xτ 1{τ <∞} 0 Bubbles • A bubble is defined to be a process β = (βt )t≥0 given by βt = St − St? • Note that βt ≥ 0 for all t, a.s. • Theorem: If there exists a non trivial bubble (ie, βt 6= 0 for some t > 0) for the risky asset price process S, then 1. If P(τ = ∞) > 0 then β is a local martingale without restrictions (it can even be a uniformly integrable martingale) 2. If P(τ < ∞) = 1, and β is unbounded, then β is a local martingale, and it cannot be a uniformly integrable martingale 3. If τ is bounded, and then β must be a strict local martingale Theorem (Bubble Decomposition) The risky asset price admits a unique decomposition S = S ? + (β 1 + β 2 + β 3 ) where 1. β 1 is a c`adl`ag nonnegative uniformly integrable martingale with limt→∞ βt1 = X∞ a.s. 2. β 2 is a c`adl`ag nonnegative NON uniformly integrable martingale with limt→∞ βt2 = 0 a.s. 3. β 3 is a c`adl`ag non-negative supermartingale (and strict local martingale) such that limt→∞ E {βt3 } = 0 and limt→∞ βt3 = 0 a.s. Examples • We call the three types of bubbles in the decomposition bubbles of Type 1, Type 2 and Type 3 • Example of a Type 1 bubble: Let St = 1, all t, 0 < t < ∞, and no dividends. This is an example of fiat money • In this case τ = ∞ a.s., and X∞ = 1, and Dt ≡ 0 all t ≥ 0. • Therefore ? S∞ Z t∧τ dDu + Xτ 1{t≥τ } |Ft = EQ 0 • Hence βt = St − St? = 1 =0 A second example of a Type 1 bubble • Let Bt = (Bt1 , Bt2 , Bt3 ), the three dimensional standard Brownian motion • k Bt k is called the Bessel process; Xt = 1 k Bt k is known as the inverse Bessel process • Assume there are no dividends, only the asset price. One can show that lim Xt = 0 t→∞ and that X is a local martingale; indeed, X satisfies the SDE dXt = −Xt2 dBt ; X0 = 1 where B is a Brownian motion • Also E (X0 ) = 1 and limt→∞ E (Xt ) = 0 Example of a Type 2 bubble • Let τ be a stopping time with P(τ > t) > 0, for all t > 0, and P(τ < ∞) = 1 • Let St? = 1{t<τ } , payoff 1 at time τ 1 − 1{τ ≤t} βt = P(τ > t) St = St? + βt • One can show that β is a martingale which is not uniformly integrable, and β∞ = 0 • So β is a bubble which is not uniformly integrable Example of a Type 3 bubble [A. Cox and D. Hobson, 2005] • Let T be a fixed (non random) time, and define St? = 1{[0,T )} (t), XT = 1 • Let the bubble be given by Z βt = 0 t √ βu dBu T −u • Then β is a strict local martingale, with Bt = 0; define St = St? + βt Historical example of an option • Aristotle, in his treatise Politics; Book 1, Part XI, writes of Thales of Miletus, a pre-Socratic Greek philosopher and one of the Seven Sages of Greece • Thales wanted to justify his beliefs in astronomy, which allowed him to predict (correctly, as it turned out) that there would be a bumper olive crop harvest (Source: Walter Schachermayer) • According to Aristotle, Thales “gave deposits for the use of all the olive-presses in Chios and Miletus, which he hired at a low price because no one bid against him. When the harvest-time came, and many were wanted all at once and of a sudden, he let them out at any rate which he pleased, and made a quantity of money.” An ancient olive press used to make oil Put-Call Parity in the Presence of Bubbles • A Call Option has the payoff structure at the maturity time T of (ST − K )+ and a put (K − ST )+ and a forward contract at strike price K and maturity time T has a payoff at time T of ST − K • Recall that trivially (ST − K )+ − (K − ST )+ = ST − K • Let Ct (K ), Pt (K ), and Vt (K ) be the market prices at time t and strike price K with common maturity time T of a call, a put, and a forward • Let Ct (K )? , Pt (K )? , and Vt (K )? be the fundamental prices at time t and strike price K with common maturity time T of a call, a put, and a forward • The traditional approach for a complete market for put call parity is to define the time t price of (for example) a European call to be EQ {(ST − K )+ |Ft } and then put-call parity follows from the linearity of conditional expectation • The issue of whether or not market prices agree with the conditional expectation prices is assumed to be true • With bubbles, the market prices of calls, puts, and forwards need not satisfy put-call parity Example of Put-Call Parity Failing • Let B i , 1 ≤ i ≤ 5 be five iid standard Brownian motions • Define Mt1 = exp(Bt1 − t/2) Z t Mi √ s dBsi , Mti = 1 + T −s 0 2≤i ≤5 • Consider a market with finite time horizon T • It is complete, given M i , 1 ≤ i ≤ 5 • M 1 is a uniformly integrable martingale, and the rest are strict local martingales on [0, T ] • Let St? = sup Mt1 ; P(K )t s≤t ? St = St? + Mt2 ; = P (K )t + Mt4 ; C (K )t = C ? (K )t + Mt3 V (K )t = V ? (K )t + Mt5 • All the traded securities in the this example have bubbles • Let δtC , δtP , and δtF be the bubbles parts of the market prices for the Call, Put, and Forward. • Under special conditions only (the absence of bubbles) do we have market price put-call parity: Ct (K ) − Pt (K ) = Ft (K ) if and only if δtF = δtC − δtP Ct (K ) − Pt (K ) = St − K if and only if δtS = δtC − δtP Implications for Models in the Black-Scholes Paradigm • To take advantage of these bubbles based on the convergence at time T , one needs only to short sell at least one asset • Such a strategy, however, is not admissible due to possible unbounded losses • By the Black-Scholes paradigm we mean a continuous risky asset price process under the now standard NFLVR structure • The important consequence is that in the presence of bubbles, the Black-Scholes formula need not hold • This is because the time t market price of a call option, Ct (K ), can differ from the price EQ {(ST − K )+ |Ft } Consequence for Black-Scholes Paradigm Models • Implied volatility from the B-S formula need not equal historical volatility; indeed, if there is a bubble, implied volatility should exceed historical volatility • However, if one assumes No Dominance, then the usual understanding of the Black-Scholes model applies • Another issue is Merton’s No Early Exercise Theorem • This theorem states that while an American call option with strike price K and maturity time T has the a priori impression of presenting more flexibility in the exercise of the option, in reality the optimal strategy is to exercise it at maturity T . Therefore the fair prices of an American call option and that of a European call option are the same • The proof of Merton’s theorem uses Jensen’s inequality and assumes the risky asset risk neutral price process is a martingale Under NFLVR and continuous complete markets, Merton’s No Early Exercise Theorem need not hold • We give an example where No Early Exercise fails to hold • Let Bt = (Bt1 , Bt2 , Bt3 ) be a standard Brownian motion with B0 = (1, 0, 0) • Recall that the inverse Bessel process is Xt = 1 k Bt k which is a strict local martingale • If X models a risky asset price process, then the price process is a bubble • (Xt )t≥0 is a uniformly integrable collection, E (X0 ) = 1, and limt→∞ Xt = 0 a.s. and in L1 • If X is a risk neutral (Q) martingale, then by Jensen’s inequality, t 7→ EQ {(Xt − K )+ } is monotone increasing • For the inverse Bessel process, with Soumik Pal, we have shown that the prices of European calls decrease as a function of time to expiration • That is, for S the inverse Bessel process, the function T 7→ E {(ST − K )+ } is monotone decreasing if K ≤ 12 , and otherwise it is initially increasing and then strictly decreasing for 2K + 1 −1 T ≥ K log . 2K − 1 • A similar results holds for all continuous strict local martingales with asymptotic behavior similar to that of the inverse Bessel process • This result is intuitive in the presence of bubbles, since in a bubble, the best strategy is to get in and out early, and not to wait a long time to liquidate your positions Bubble Decomposition Theorem [Bubble Decomposition]: St = St? + βt = St? + βt1 + βt2 + βt3 , is a unique decomposition such that • β 1 ≥ 0 is a uniformly integrable martingale with limt→∞ βt1 = X∞ a.s. • β 2 ≥ 0, is not a uniformly integrable martingale, but of course is a local martingale and is possibly a martingale, and limt→∞ βt2 = 0 a.s. • β 3 ≥ 0 is a strict local martingale such that limt→∞ βt3 = 0 a.s. and in L1 Why Does Short Selling Not Correct for Bubbles? • Two reasons are proposed in the literature: • The first is structural limitations: This is the limited ability and/or expensive cost to borrow an asset for short sales (eg, Duffie, Gˆarleanu, and Pederson [2002]) • As regards the first, in markets where short selling does not exist (especially the third world), there do not seem to be more bubbles • The second is the risk the short seller takes that the price will continue to go up (the danger of trying to predict a bubble) • In mathematics this translates into admissibility violations Two Problems with Complete Markets and Bubbles • What is nice is that the risk neutral measure Q is unique, and we therefore have a unique fundamental price • An undesirable property is the impossibility of bubble birth: A nonnegative local martingale cannot spring up after being zero; once a nonnegative local martingale reaches zero, it sticks at zero forever after • The biggest problem is that while bubbles make sense in complete markets under NFLVR, bubbles do not exist under No Dominance. This is serious, because we will see later we need No Dominance to establish fundamental put-call parity • Theorem: Under No Dominance, Type 2 and Type 3 bubbles do not exist in a complete market (with NFLVR) Proof that Bubbles Do Not Exist in Complete Markets under ND • Theorem: Under No Dominance, Type 2 and Type 3 bubbles do not exist in a complete market (with NFLVR) • Proof: For Type 2 and Type 3 bubbles, β∞ = 0. Let W be the wealth process corresponding to the risky asset price process S • There exist hedging processes π 1 and π 2 such that Wt? βt Z t + πu1 dWu Z t0 = β0 + πu2 dWu = W0? 0 • Let η 1 and η 2 make π 1 and π 2 self-financing, so that both π 1 and π 2 are admissible • We have two ways to generate W • The first way is buy and hold • The second way is to follow π 1 , obtaining W ? • The cost of the first position is W0 ≥ W0? , with W0 > W0? if there is a non-trivial bubble • That means that π 1 dominates the buy and hold strategy, which violates No Dominance; so β cannot exist. • We conclude: bubbles exist in a complete market under NFLVR [Lowenstein and Willard, Cox and Hobson], but cannot be born after time t = 0, create a Black-Scholes paradox, and violate put-call parity. • Bubbles do not exist in a complete market under No Dominance, which is stronger than NFLVR • The non existence of bubbles under ND solves the Black-Scholes paradigm paradox, for example • What happens in incomplete markets? In incomplete markets under No Dominance the argument showing bubbles do not exist, no longer applies Do Bubbles Exist in Incomplete Markets? • To discuss bubbles in incomplete markets, we need to decide what we mean by a fundamental price, since there is an infinite choice of risk neutral measures • There are five basic methods to choose such a measure • The first is Utility Indifference Pricing: Risk Neutral prices span an interval on the real line, and choosing the right price depends on the utility function of preferences of the agent selling the contingent claim • The Egocentric Method: Simply choose one arbitrarily • The Convenience Method: Choose a risk neutral measure • • • • that gives the price process mathematically nice properties: for example, makes it a Markov process, or even a L´evy process The Canonical Method: Find a reasonable criterion (eg, minimal variance of the error, minimal distance to the historical measure in a distance one chooses, minimal entropy) and let it determine the risk neutral measure The Ostrich Method: Prove results under a risk neutral measure already chosen; that is, you do not specify it, but pretend someone else has done so already All of these methods assume that, once chosen, the risk neutral measure is fixed and and never changes The current methods share two properties: 1. In one sense or another they can be considered to be ad hoc 2. Once a risk neutral measure is chosen, it remains fixed in the model for all time (which is often a finite horizon model, and time is modeled by an interval [0, T ]) • Instead of these approaches, one can let the market choose the risk neutral measure • En passant, the market is often wrong (otherwise we would not have bubbles!), but the issue is not “truth,” as it is, for example, in physics • The usual approach is as follows: • Begin with the risky asset price process and a (possibly random) savings rate; • By a change of num´ eraire argument, assume interest rates are zero; • Find and choose a risk neutral measure Q; • If g (XT ) is a financial (European style) derivative at time T , the price process is declared to be EQ {g (Xt )|Ft } for 0 ≤ t ≤ T. The approach of Jacod and myself: • Begin with the risky asset price process and a (possibly random) savings rate; • Next assume there are market given price processes for a large number of (European style) derivatives of the form g (XT ), where T varies; • Find the collection Q of risk neutral measures that make both the price process and all of the derivative price processes local martingales; • If there are enough derivative prices, the cardinality of Q might be one; • Alternatively the cardinality of Q might be zero or ∞ • compatibility issues are serious here. Other Approaches; Some key ones follow: • H. Dengler and R. Jarrow (1997): A simple jump model where the price process is pure jump, with two distinct jump sizes; two call options are used to complete the market. • B. Dupire (1997): A simple stochastic volatility model, where the stock price varies but the expiration time T is fixed. A PDE is obtained, and when (and if) solved, it gives a unique martingale measure; • E. Derman and I. Kani (1997): the stochastic volatility case, where the strike price K varies and T is fixed. Smoothness in K is assumed, and by twice differentiating obtain a density for the call option, leading to a PDE, which when (and if) solved leads to a martingale measure choice. • M. Schweizer and J. Wissel (2006), to appear in Math Finance: improved on previous ideas by using different maturities, creating a term structure of volatilities, within a Brownian based stochastic volatility framework. They find conditions for equivalent martingale measures to exist, and also a condition for it to be unique. How Do We Let the Market Choose the Risk Neutral Pricing Measure? • Begin with the risky asset price process S and a (possibly random) savings rate; • Next assume there are market-given price processes for a large number of (European style) derivatives of the form g (ST ), where T varies; • Find the collection Q of risk neutral measures that make both the price process and all of the derivative price processes local martingales • If there are enough derivative prices, the cardinality of Q might be one; • Alternatively the cardinality of Q might be zero or ∞ (compatibility issues are serious here). Regime Change • Recall that we have an infinite number of possible risk neutral measures • Let us assume that market has chosen a unique one determined by (an infinite number of) option prices • A fortiori the option prices must be internally consistent as well; else we would have no risk neutral measure matching, due to arbitrage opportunities (Delbaen-Schachermayer) • Given the nature of the market over time, it is unreasonable to assume that it stays with that risk neutral measure for all time • When the market changes from one risk neutral choice to a different one, we call it a regime change • One can construct regime changes via stopping times • Very recently Biagini, F¨ ollmer, and Nedelcu have proposed a regime change via a continuous flow of measures • This gives an elegant way for a bubble to be born gradually • Finally, one can also simply change the coefficients of the SDE governing the price process, either at a random time, or via a continuous flow The Fundamental Price for Incomplete Markets • What follows is based on work with Robert Jarrow and Kazuhiro Shimbo • Recall that for the case of complete markets with a finite horizon T , with risk neutral measure Q, and for t < T the fundamental price of the risky asset is defined to be: Z St? = EQ { t T dDu + XT |Ft } • In incomplete markets, if one Q is chosen by the market for all time (a “static market”), the definition is analogous. • In an incomplete market with an infinite horizon, we assume there exists a countable sequence of stopping times 0 = T0 < T1 < T2 < . . . increasing to ∞ a.s. which represents change times from one risk neutral probability to another • The stochastic interval [Ti , Ti+1 ) consists of the i th regime • In an incomplete market with an infinite horizon and regime change, the fundamental price of the risky asset with end time τ for the asset, t < τ , and for regime i at time t, is defined to be: Z τ St? = EQ i { dDu + Xτ 1{τ <∞} |Ft } t where Q i is the risk neutral measure chosen by the market • Note that Xτ 1{τ =∞ } is not included The Evaluation Measure • We can piece all of these measures Q i together to get one measure Q ? • Q ? need not be a risk neutral measure; (if it were, then in effect we would not have regime change) • Some people find Q ? not being a risk neutral measure, although it is equivalent, to be troubling • We call Q ? the evaluation measure, and write it Q t? to denote that it changes with the time t • The Fundamental Price can be written compactly: St? = X Z 1[Ti ,Ti+1 ) (t)EQ i { τ dDu + Xτ 1{τ <∞} |Ft } t i or using Q ? = Q t? Z St? = EQ t? { t τ dDu + Xτ 1{τ <∞} |Ft } End of Lecture 3 Thank you for your attention How to Model a Financial Bubble Mathematically Lecture 4 Philip Protter, Columbia University Statistical Methods and Applications in Insurance and Finance CIMPA School April 13, 2013 • At the end of Lecture 3 we explained how the market chooses one of the infinite number of risk neutral measures • The observer can tell this by observing an infinite number of option prices • The Fundamental Price can be written compactly: St? = X Z 1[Ti ,Ti+1 ) (t)EQ i { τ dDu + Xτ 1{τ <∞} |Ft } t i or using Q ? = Q t? St? = E Q t? Z { t τ dDu + Xτ 1{τ <∞} |Ft } • In incomplete markets, if one Q is chosen by the market for all time (a “static market”), the definition is analogous. • In an incomplete market with an infinite horizon, we assume there exists a countable sequence of stopping times 0 = T0 < T1 < T2 < . . . increasing to ∞ a.s. which represents change times from one risk neutral probability to another • The stochastic interval [Ti , Ti+1 ) consists of the i th regime • In an incomplete market with an infinite horizon and regime change, the fundamental price of the risky asset with end time τ for the asset, t < τ , and for regime i at time t, is defined to be: Z τ St? = EQ i { dDu + Xτ 1{τ <∞} |Ft } t where Q i is the risk neutral measure chosen by the market • Note that Xτ 1{τ =∞ } is not included The Evaluation Measure • We can piece all of these measures Q i together to get one measure Q ? • Q ? need not be a risk neutral measure; (if it were, then in effect we would not have regime change) • Some people find Q ? not being a risk neutral measure, although it is equivalent, to be troubling • We call Q ? the evaluation measure, and write it Q t? to denote that it changes with the time t • The Fundamental Price can be written compactly: St? = X Z 1[Ti ,Ti+1 ) (t)EQ i { τ dDu + Xτ 1{τ <∞} |Ft } t i or using Q ? = Q t? Z St? = EQ t? { t τ dDu + Xτ 1{τ <∞} |Ft } Bubble Birth • Recall the definition of a bubble: Defintion[Bubble]: A bubble in a static market for an asset with price process S is defined to be: βt = St − St? ≥ 0, t≥0 • A bubble in a dynamic market for t < τ in regime i is: Z βt = St − EQ t? { τ dDu + Xτ 1{τ <∞} |Ft } ≥ 0 t • Since we are in regime i, we have in this case Q t? = Q i . • If there are no bubbles for regime i − 1, a change to a new risk neutral measure at change time Ti might in effect create a bubble; if it does, we call this bubble birth Derivatives • Recall we touched on European Calls and Puts in Lecture 3 • For a European call with payoff (ST − K )+ at time T , we denote its market price at time t < T as Ct (K ) • Analogously for a put with payoff (K − ST )+ at time T , its time t market price is Pt (K ) • A forward with payoff ST − K has market price at time t denoted Vt (K ) • Their fundamental prices are given by Ct? (K ) = EQ t? {(ST − K )+ } Pt? (K ) = EQ t? {(K − ST )+ } Vt? (K ) = EQ t? {(ST − K )} • Theorem [Put-Call Parity]: We have put-call parity for fundamental prices Ct? (K ) − Pt? (K ) = Vt? (K ) • The proof follows from the linearity of the conditional expectation with respect to Q t? • Theorem: Pt (K ) = Pt? (K ) • The proof follows from the fact that it is bounded, plus the hypothesis of No Dominance • In contrast, if there is a bubble in the risky asset price S, then it is captured by the call option market price process Ct (K ) • Theorem: for K > 0 we have Ct (K ) − Ct? (K ) = St − EQ t? {ST |Ft } • The proof of this theorem follows from the observation regarding bubbles of forwards: Vt? (K ) = EQ t? {ST − K |Ft } = EQ t? {ST |Ft } − K = St? − K ≤ St − K = Vtf (K ) which means that a forward has a type 3 bubble of size δt3 = Vtf (K ) − Vtf ? (K ) = (St − K ) − (St? − K ) = St − St? • Put-call parity and the fact that Pt (K ) = Pt? (K ) imply Ct (K ) − Ct? (K ) = St − EQ t? {ST |Ft } • Note that with calls we are dealing with Type 3 bubbles, since they exist on the compact time interval [0, T ] Testing for Bubbles • Recall that we have put-call parity even in the presence of bubbles • This is in contrast to the views expressed in the literature • For example, Battalio and Schultz find no violations of put call parity during the internet stock price bubble, and argue that this evidence is inconsistent with an internet stock price bubble • However a “correct” model for the price operator EQ ( ·| Ft ), when applied to both call and put options on the same spot commodity, would give differential results in the presence of price bubbles • Puts would be priced correctly, but calls would not. If type 3 bubbles exist, this difference would be observable • Also, the mispricings would be independent of the moneyness of the options, but dependent on the time to maturity Existence • dXt = σ(Xt )dWt + b(t, Xt , Yt )dt dYt = g (Yt )dBt + f (Yt )dt • W and B are iid Brownian motions • The process X is a strict local martingale under a risk neutral measure Q if and only if Z ∞ x dx < ∞, any > 0 2 σ (x) • This model is incomplete and satisfies No Dominance Quotes • Federal Reserve Chairman Ben Bernanke: “It is extraordinarily difficult in real time to know if an asset price is appropriate or not.” Congressional Testimony, December, 2009 • New York Fed President William Dudley: “What I am proposing is that we try to identify bubbles in real time, try to develop tools to address those bubbles, try to use those tools when appropriate to limit the size of those bubbles and, therefore, try to limit the damage when those bubbles burst.” Interview with Planet Money, April 9, 2010 Detecting Bubbles - Joint Hypo with {S} • Specify a stochastic process dSs = σ(Ss )dZs + µ(Ss , Ys )ds; S0 = x where Z a standard BM, and Y is a diffusion driven by an independent Brownian motion • Under a risk neutral measure Q, this simplifies to Z t St S0 = + σ(Ss )dZs . Bt B0 0 • The process S is a strict local martingale if and only if Z ∞ x dx < ∞. σ(x)2 • Empirically test this integral, which we can do under the original probability measure P • Can independently test for fit of asset price evolution. • Results: Asset price bubbles existed during the Tech bubble. Bubble Detection Procedure • We assume the price process follows an SDE of the form dXt = σ(Xt )dWt + b(t, Xt , Yt )dt (1) • First we estimate σ from available historical data, using estimators from the literature, slightly adapted for our purposes • If the function σ exhibits some regularity of behavior, we attempt to extrapolate an extension of its behavior beyond the range of values assumed by X • If we assume σ is a parameterized family of functions, then we need only estimate the parameter, and this is routine • Otherwise, if we can show it is contained by a parameterized family via comparison theorems, we can extend it as well A simple model • Let S satisfy a stochastic differential equation of the form: dSt = σ ? (St )dBt + µ(St , Yt )dt; St S0 > 0 where Y can be (for example) an independent diffusion. • Our goal is to determine the function σ ? (x) from frequent observations of St as time evolves Natural Limitations of the Problem • By the assumption of the structure of S, we have St > 0 for all t ≥ 0 • We assume that σ ? and µ are such that a unique solution exists, and without explosions; for example, xσ ? (x) and xµ(x, y ) are Lipshitz continuous • Due to the lack of explosions, on a finite time interval [0, T ] the process S is bounded, ω by ω, so we can know xσ ? (x) only on a compact interval of R+ • We consider only the one dimensional problem; multiple dimensions are much harder, as is the case of stochastic volatility • The definitive paper on this subject is doe to Jean Jacod: Non-parametric kernel estimation of the coefficient of a diffusion, Scandinavian J. Statist., 27 (2000), pp. 83-96. A Key Tool: The Space Time Formula • Let X be a semimartingale with a local time process Lat for the level a; then [X , X ]ct ∞ Z = −∞ Lat da • More generally, Z t g (Xs− )d[X , X ]cs Z ∞ = 0 −∞ Lat g (a)da for any bounded, Borel function g • In the case of our SDE, S = X and d[S, S]ct = St2 σ ? (St )2 dt • Let us simplify notation by assuming σ(x) = xσ ? (x); then d[S, S]ct = σ 2 (St )dt The Idea of Florens-Zmirou • In this Brownian paradigm, LxT 1 = lim →0 2 T Z 0 1{|Ss −x|<} ds • The actual local time of interest is given by `xT 1 = lim →0 2 Z 0 T 1{|Ss −x|<} d[S, S]s • Since [S, S]s = σ 2 (Ss )ds, we have `xT = σ 2 (x)LxT , whence `xT = σ 2 (x) LxT and we can create an estimator for σ 2 (x) The Estimator of Florens-Zmirou • We consider two sets of sums: LnT (x) = n T X 1{|Sti −x|<hn } 2nhn i=1 n T X n 1{|Sti −x|<hn } n(Sti+1 − Sti )2 `T (x) = 2nhn i=1 where hn is a sequence of positive real numbers converging to 0 and satisfying some constraints. • This allows us to construct an estimator of σ(x) given by: Pn Sn (x) = i=1 1{|Sti −x|<hn } n(Sti+1 − Pn i=1 1{|Sti −x|<hn } Sti )2 (2) Consistency of the Estimator of Florens-Zmirou • Florens-Zmriou proved the following SLLN type theorem: Theorem If σ is bounded above and below from zero, has three continuous and bounded derivatives, and if (hn )n≥1 satisfies nhn → ∞ and nhn4 → 0 then Sn (x) is a consistent estimator of σ 2 (x). • The proof of this theorem is based on the expansion of the transition density. The choice of a sequence hn converging to 0 and satisfying nhn → ∞ and nhn4 → 0 allows one to show that LnT (x) and `nT (x) converge in L2 (dQ) to LT (x) and σ 2 (x)LT (x), respectively. Hence Sn (x) is a consistent estimator of σ 2 (x), for any x that has been visited by the diffusion. Confidence Intervals of the Estimator of Florens-Zmirou • Florens-Zmirou also proved another CLT type theorem, useful to obtain confidence intervals for the estimator Sn (x) of σ(x). Theorem p If moreover nhn3 → 0 then Nxn ( σSn2 (x) − 1) converges in (x) √ distributionPto 2Z where Z is a standard normal random variable and Nxn = ni=1 1{|Sti −x|<hn } . • The requirement that nhn4 → 0 requires impractical amounts of data • By introducing a smoothing kernel, we can relax the condition to nhn2 → ∞ The Smooth Kernel Approach • The idea underlying the smooth kernel estimator is to replace the kernel K (x) = 12 1{|x|≤1} by a smooth kernel φ, which is a 6 C R positive function with compact support and such that R+ φ = 1 • We are interested in some kind of convergence of the following quantities: Vnx = n−1 1 X S ni − x )n(S i+1 − S i )2 φ( n n nhn hn (3) n−1 1 X S ni − x φ( ) nhn hn (4) i=0 Lxn = i=0 to σ 2 (x)Lx and Lx respectively, where hn satisfies nhn2 → ∞. The Theorem • After massive calculations and estimates, we get: Theorem x If nhn2 → ∞ then Snx = VLnx converges in probability to σ 2 (x) and n provides a consistent estimator of σ 2 (x). • In the work of Jacod of 2000, he is able to take hn = √1n and moreover obtains a rate of convergence and an associated Central Limit Theorem, as well as treating the multidimensional case. Relaxing the restriction that σ > 0 • Let F = (Ft )t≥0 be the underlying Brownian filtration. • Let τp be the first exit of S from the interval ( p1 , p]. We use the theory of the enlargement of filtrations • Let Gt = σ(Ss ; s ≤ t) ∨ σ(τp ) • Or more generally, Gt? = σ(Ss ; s ≤ t) ∨ U, where U = 1{τp >T } • If we work on the space (Ω, G ? , Q, G? ), where Q is the conditional probability that σ(S) stays in the interval ( p1 , p], then we can obtain our previous results for Q; all works because this special event is in the filtration. • We piece the convergence in probability together. This works if, in practice, S never attains the value 0 Going Beyond the Visible Range • Suppose we can estimate σ(x) with reasonable accuracy for all x in some compact interval [a, b] of R+ , as we now know that we can • What does σ(x) look like outside of [a, b]? • This is, of course, intrinsically unknowable • But we need to know it for (for example) a theory of bubble detection • In particular, we would like to know the asymptotic properties of σ(x) as x % ∞ • What to do? First Method: Parametric Estimation • This approach uses the classic paper of V. Genon-Catalot and J. Jacod, AIHP, 1993 • The idea is to consider a parameterized family of coefficients σ(α, x) and then once knows α, one knows the function σ(α, x) for all x ∈ R+ • It is too much to expect σ to be in such a family, but it might be “squeezed” in between two members of the parameterized family • A measure of asymptotic growth of interest in the theory of bubbles is whether or not Z ∞ x dx < ∞ 2 σ(x) • If the above integral is finite, then there is a bubble • Basically, σ(x) has to increase asymptotically faster than f (x) = x for the integral to be finite The Comparison Theorem • Theorem (Comparison Theorem) Assume that dSt = σ(t, St )dWt and that there exist two functions Σ and σ such that : for all t and x, σ(x) ≤ σ(t, x) ≤ Σ(x) and such that σ, Σ and σ are continuous, locally H¨ older continuous with exponent 12 , then: R∞ R∞ x 1. if for all c > 0, c Σ2x(x) dx = ∞ then σ(x) 2 dx = ∞ R∞ x 2. if there exists c > 0 such that c σ2 (x) dx < ∞ then R∞ x dx < ∞ σ(x)2 Bubble Birth and Regime Change • One way we can have bubble birth is to incorporate the ideas of regime change • At a sequence of stopping times, the function σ in the volatility of equation (1) changes to a different σ • Another possibility for regime change is for the risk neutral measure Q to change to another one • Biagini, F¨ ollmer, and Nedelcu have a new result that allows the risk neutral measures to change via a continuous flow of equivalent measures • This gives a continuous change to risk neutral measure, giving an elegant method for bubble birth A Comparison Theorem Comparison Theorem: Assume that dSt = σ(t, St )dWt and that there exist two functions Σ and σ such that : ∀t, ∀x, σ(x) ≤ σ(t, x) ≤ Σ(x) and such that σ, Σ and σ are continuous, locally Holder continuous with exponent 12 , then: R∞ 1. If ∀c > 0, c Σ2x(x) dx = ∞ then S is a martingale. R∞ 2. If ∃c > 0, c σ2x(x) dx < ∞ then S is a strict local martingale. Another approach is that of Reproducing Kernel Hilbert Spaces (RKHS) Estimation of σ Method 1: Parametric Estimation of σ • We choose a class of volatility functions large enough to include many of the forms used in practice • An example: power functions σ(x) = σ ¯ x α , where σ ¯ and α are the unknown parameters to be estimated • Next we use non-parametric estimators for the function σ • If the parametric estimators and the non-parametric estimators are similar, then we know the behavior for large x through the parametric continuation, and we can test for convergence or divergence • If they are not comparable, then this method fails Estimation of σ, Continued Method 2: Reproducing Kernel Hilbert Spaces [RKHS] • Step 1: We first interpolate an estimate of σ within the bounded interval where we have observations • Consequence: We lose the irregularities of non parametric estimators • We then choose an RKHS family corresponding to this interpolated σ so as to optimize how close we are to σ on the bounded interval. We give a criterion for such an optimization. • We do this by choosing a certain extrapolating RKHS, which once again determines the tail behavior of σ • More precisely: Let (Hm )m∈N denote our family of RKHS spaces; Given any choice of m, call it m0 , we can interpolate perfectly • If we stop here, this method is just as arbitrary as parametric estimation and continuation Estimation of σ, Continued • Instead we use the data a second time to choose m in a reasonable way • Thus this method is sort of a two step bootstrap method • Define the Sobolev space: H n (D) = u ∈ L2 (D)st.∀k ∈ [1, n], u (k) ∈ L2 (D) where u (k) is the weak derivative of u. • The norm that is usually chosen is R ||u||2 = Pn k=0 D (u (k) )2 (x)dx. • Due to Sobolev inequalities, an equivalent and more R R appropriate norm is ||u|| = D u 2 (x)dx + where τ keeps things away from 0 1 τ 2n D (u (n) )2 (x)dx, • We denote by K the kernel function of H n (]a, b[), where in this case D =]a, b[. • We then can establish the form of K , which is simple but messy to state Estimation of σ, Continued • We next must decide if an extrapolation is required: If the natural extension of the estimated volatility σ shows the integral in question does not diverge as x → ∞ and remains bounded on R+ , no extrapolation is required, the process is a martingale • If however it appears that the σ(x) → ∞ as x → ∞, then more work is required • Extrapolate σ using our RKHS spaces and choose the one that minimizes the distance, using a weighted norm, away from σ ˆ on the bounded interval where observations exist. • In this way we obtain a function that interpolates all the available data points, but also remains as close as possible to the interpolated function on the observation interval Does this work with data? We test this method with data from the Dot Com Bubble Era (2000-2002): Recall we need to estimate σ and determine the convergence or divergence of the integral Z ∞ x dx, any > 0 σ(x)2 LastMinute.com Historical open/close data LastMinute.com Estimates of σ using parametric methods There is a bubble Infospace.com Historical open/close data Infospace.com Estimates of σ Infospace.com Estimates of σ using RKHS There is a bubble TheGlobe.com Historical open/close data TheGlobe.com Estimates of σ There is a bubble Geocities Historical open/close data Geocities Estimates of σ There is NOT a bubble eToys.com Historical open/close data eToys.com Estimates of σ using parametric methods The Test is inconclusive The Recent Case of LinkedIn Stock • LinkedIn shares jumped 109.4 percent on their first day of trade • Chicago Fed President Charles Evans said he was withholding judgment over whether a new dot-com bubble was under way • “I have no way of knowing that those aren’t just exactly the right valuations,” Evans told reporters • We used data from the IPO on May 19 until May 24 to show there is a bubble Tick Data of LinkedIn RKHS best estimate of σ The Recent Case of Gold Prices • There was a lot of speculation recently that gold prices were in a bubble • An example was the Dealbook article of August 30, 2011, by S.M. Davidoff, “How to Deflate a Gold Bubble (That Might Not Even Exist).” Gold Spot Prices RKHS best estimate of σ There is Not a Bubble Implications for Risk Management • Bubbles distort economic decisions of financial institutions, regulators, investment firms. • The tools for risk management - financial derivatives - are not priced nor hedged in the standard text book fashion. • Organizations that use the standard textbook formulations will mismanage risk. • Regulatory perspective - the determination of economic capital will be too low, because in the presence of bubbles, the drift on the relevant asset price evolutions are too large (too optimistic). • Investment decisions of pension funds and institutions will also be distorted because risky asset prices will look more attractive than they should. Thank you for your Attention

© Copyright 2018