EQUITY MARKET STABILITY

Size: px
Start display at page:

Download "EQUITY MARKET STABILITY"

Transcription

1 EQUITY MARKET STABILITY Adrian Banner INTECH Investment Technologies LLC, Princeton (Joint work with E. Robert Fernholz, Ioannis Karatzas, Vassilios Papathanakos and Phillip Whitman.) Talk at WCMF6 conference, Santa Barbara September 27, 2014 For Information Purposes Only 1 / 26

2 Model for an equity market Equity market framework with capitalizations modeled by continuous semimartingales; e.g., Bachelier, Samuelson... dx i (t) X i (t) = α i(t) dt + N ξ iν (t) dw ν (t), i = 1,..., N, ν=1 where W ( ) = (W 1 ( ),, W n ( )) is B.M. in dimension N N. Mean Return rates α( ) = (α 1 ( ),..., α N ( )) Volatility rates ξ( ) = (ξ iν ( )) 1 i N,1 ν N Covariance rates σ( ) = ξ( )ξ ( ), (N N) matrix given by σ ij (t) := N ξ iν (t)ξ jν (t) = ν=1 1 X i (t)x j (t) d dt X i, X j (t). Both α( ) and σ( ) are assumed to be locally integrable. 2 / 26

3 Logarithmic representation Define log capitalizations Y ( ) by so that the equation becomes Y i ( ) := log X i ( ), dx i (t) X i (t) = α i(t) dt + dy i (t) = γ i (t) dt + N ξ iν (t) dw ν (t) ν=1 N ξ iν (t) dw ν (t), ν=1 where the growth rate of the i th asset is γ i (t) := α i (t) 1 2 σ ii(t). 3 / 26

4 Market weights Define the total market capitalization X (t) := X 1 (t) + + X N (t). Consider now the relative weights of the various assets, in terms of capitalization: µ i (t) := X i(t) X (t) = X i (t), i = 1,..., N. X 1 (t) + + X N (t) These are strictly positive numbers, add up to 1, and are determined on the basis of the capitalizations of the market s various assets at any given time. They are perfectly observable, at all times. 4 / 26

5 Ranked log capitalizations Descending order statistics: max Y i( ) =: Y (1) ( ) Y (N) ( ) := min Y i( ). 1 i N 1 i N Ties are resolved lexicographically, by resorting to the lowest index (name) i. Similar ordering for the capitalizations X (k) ( ) = exp(y (k) ( )) and for the ordered market weights µ (k) ( ) = X (k) ( ) X 1 (t) X N (t), k = 1,, N. 5 / 26

6 Capital distribution, U.S. equity market % 1% 10 bps Market Weight 1 bp 0.1 bp 0.01 bp bp bp bp ,000 10,000 Stocks Ranked by Capitalization 6 / 26

7 Implications of market stability Many rules-based strategies ( smart beta ) are seen to outperform cap-weighted indexes over the long term. Equity market stability is critical to this outperformance. Example: a continuously rebalanced equal-weight portfolio E. (N.B. over $10B invested in equal-weight portfolios.) Its relative value process versus the market M is (Fernholz 2002): d log(v E (t)/v M (t)) = ds(t) + θ(t) dt where θ( ) is a.s. positive and S(t) := 1 N N log µ i (t). A stable capital distribution implies that S( ) is mean-reverting, hence the relative performance is positive over sufficiently long time periods. i=1 7 / 26

8 Collision local times The ordered processes Y (k) are themselves semimartingales (as we shall see), and for 1 k < l N, the gap ( ) ( ) G (k,l) X(k) ( ) µ(k)( ) ( ) := Y (k) ( ) Y (l) ( ) = log = log 0 X (l) ( ) µ (l) ( ) is a continuous, nonnegative semimartingale. We shall denote by Λ (k,l) ( ) := L G (k,l) ( ) its local time at the origin (collision local time of order l k + 1). For a nonnegative continuous semimartingale G( ), recall its local time process L G ( ) = {G(t)=0} dg(t) = lim ε 0 2ε 0 1 {0 G(t)<ε} d G (t). 8 / 26

9 Reasonable conditions for stability A: No single stock crashes relative to the entire market (coherence): lim T 1 T log µ i(t ) = 0, a.s., i = 1,..., N. B: Convergence in distribution: ( µ(1) (T ),..., µ (N) (T ) ) (M 1,, M N ), as T. C: Strong laws of large numbers for local times: lim T 1 2T Λ(k,k+1) (T ) =: λ k > 0, a.s., k = 1,, N. 9 / 26

10 A result of Fernholz Suppose that the eigenvalues of σ( ) are bounded from above and below a.s., and that γ i (t) = γ(t) for all i = 1,..., n (common growth rate). Then lim T 1 T µ T (1) (t) dt = 1. 0 This is a degenerate form of stability in which a single stock occupies virtually the entire market capitalization. This is not consistent with empirical observations. We seek a simple model with stability but not concentration into a subset of the stocks. 10 / 26

11 First-order models With given real numbers γ, γ 1,..., γ N and σ 1 > 0,..., σ N > 0, and independent Brownian motions W 1 ( ),..., W N ( ), we consider the system of SDEs with rank-based interactions dy i (t) = ( γ + ) N γ k 1 {Yi (t)=y (k) (t)} dt k=1 + N σ k 1 {Yi (t)=y (k) (t)}dw i (t). k=1 Bass & Pardoux (1987) show that this system has a weak solution, unique in distribution (well-posed martingale problem). 11 / 26

12 Pathwise uniqueness Ichiba, Karatzas & Shkolnikov (2013) show that there is pathwise uniqueness, thus also strength, for the above system up until the first time 3 particles collide simultaneously. For this never to happen, it is necessary and almost sufficient (IKS op. cit.) that the mapping {1,, N} k σ 2 k (0, ) be concave. Sarantsev has recently proved that concavity is indeed sufficient. However, the weaker property Λ (k,l) ( ) 0 for l k 2 is always satisfied. 12 / 26

13 In such a system, the ranked log-capitalizations are rather simple to describe: Y (k) (t) = Y (k) (0) + (γ + γ k )t + σ k B k (t) + 1 ( ) Λ (k,k+1) (t) Λ (k 1,k) (t). 2 In other words, they behave as independent Brownian motions B k ( ) := N i=1 0 1 {Yi (t)=y (k) (t)}dw i (t) (by the P. Lévy theorem) with drift, bouncing off each other with the help of the two-by-two collision local times. Alternative Picture: Brownian motion with drift, reflecting off the walls of the Weyl chamber W = {(y 1,..., y N ) R N : y 1 y 2 y N }. 13 / 26

14 Stochastic stability This system is stable if γ 1 < 0,..., γ γ N 1 < 0, γ γ N = 0. Most basic example ( Atlas Configuration): γ = g > 0, γ 1 = = γ N 1 = g, γ N = (N 1)g. Under this condition, we have the SLLN for local times lim T 1 2T Λ(k,k+1) (T ) = λ k := (γ γ k ) > 0, a.s. for k = 1,..., N / 26

15 Under this stability condition, the vector process of gaps ( ) G (k,k+1) µ(k) ( ) ( ) = Y (k) ( ) Y (k+1) ( ) = log, k = 1,..., N 1 µ (k+1) ( ) converges in distribution to a random vector (G 1,..., G N 1 ) with exponential marginals P (G k > u) = e ur k, u > 0 where r k := 4λ k σ 2 k + σ2 k+1 = 1 E(G k ). The components G 1,..., G N 1 of this vector are independent, if the variances are either all the same or grow linearly by rank: σ 2 2 σ 2 1 = = σ 2 N σ2 N / 26

16 Rank-based variances, U.S., (smoothed) Variance Rank 16 / 26

17 Estimated γ k, U.S. Equity market, g k Rank 17 / 26

18 Size-based approaches The study of rank-based models such as first-order models leads to local times. Can they be avoided? Consider a simplistic size-based model is of the form dy i (t) = G(Y i (t)) dt + S(Y i (t)) dw i (t) for suitable G, S. This is not realistic, since there may be no stability in (Y 1 ( ),..., Y N ( )). Refine by choosing G and S so that (Y 1 ( ),..., Y N ( )) has an invariant distribution, then set X i (t) = M(t) exp Y i (t) for some market growth factor M(t), independent of Y i (t). Alternatively, use the weight-based variant dy i (t) = G(µ i (t)) dt + S(µ i (t)) dw i (t). 18 / 26

19 Example: volatility-stabilized model With α 0, the weight-based model (Fernholz & Karatzas 2004) dy i (t) = α 2µ i (t) dt + 1 µi (t) dw i(t) is known as a volatility-stabilized model. The coordinates can be expressed as time-changed square Bessel processes of dimension 2(1 + α). Coherence holds for α > 0 but fails for α = 0. The invariant distribution of market weights is Dirichlet (Pal 2011). Many other properties are known (e.g. Goia 2009, Shkolnikov 2011). 19 / 26

20 Example: logarithmic potential Stabilization from the previous model follows from unbounded volatility (and growth rate, for α > 0). Consider instead the size-based model dy i (t) = ( g + σ2 2Y i (t) for g > 0, σ > 0, α 0. Coordinates are now independent. ) dt + σ exp( αy i (t)) dw i (t), With λ = 2g/σ 2, the stable density G( ) on (0, ) is { C exp ( 2αx λ G(x) = 2α e2αx + Ei(2αx) ), if α > 0 λ 2 xe λx, if α = 0. Very good agreement with observed data. 20 / 26

21 A jump model A potential barrier in the form of an unbounded growth rate may seem unrealistic (similar to Atlas stock). Consider instead a jump model: dy i (t) = γ dt + σ dw i (t) if Y i (t) > 0; when Y i (t) = 0, it is kicked back into (0, ) at random location η chosen (independently) according to the distribution F (x) = P(η x), 0 < x <. Interpretation: If the relative stock capitalization is too small, the stock exits the market (bankruptcy). It is immediately replaced by a new stock (IPO). This gives a constant number of stocks for convenience; a reinterpretation as a Poisson point process would allow a variable number of stocks. 21 / 26

22 Stability in the jump model As described in Gihman & Skorokhod (1972), the system has a stable distribution G( ) if ( x ) e 2γx e 2γy (1 F (y)) dy dx <. 0 Question: how can we ensure that F G? 0 Assume the kicks are absolutely continuous w.r.t. Lebesgue measure. Then 1 F (x) = for some pdf f ( ) on (0, ). x f (u) du 22 / 26

23 It develops that where f (x) = N f (x) 0 N f (u) du, x ( N f (x) := e 2γx e 2γy If f ( ) is C 2, it must then solve the ODE 0 y ) f (u) du dy. f (x) + 2γf (x) + f (0 + )f (x) = 0, x (0, ). A solution (not known to be unique!) is given by f (x) = ( γ σ 2 ) 2 xe γx/σ 2, leading to a Gamma distribution with parameter 2, bearing more than a passing resemblance to the analogous distribution for the logarithmic potential with constant variance (α = 0). 23 / 26

24 Summary Macroscopic models should be consistent with observations of stability and growth rates. All the models considered here have a stabilization mechanism at the bottom (Atlas stock, logarithmic potential, bankruptcy/ipo). Stabilization may also occur at the top, e.g. antitrust. Many open questions and generalizations, e.g.: Understand ranks and local time properties for non-rank-based models Study properties of canonical portfolios, e.g. equal-weight, universal, in these models Better linkage of common features between these models Thank you! 24 / 26

25 A few references A. Banner, E.R. Fernholz & I. Karatzas (2005) Atlas models of equity markets. Annals of Applied Probability 15, A. Banner & R. Ghomrasni (2008) Local times of ranked continuous semimartingales. Stochastic Processes & Applications 118, R. Bass & E. Pardoux (1987) Uniqueness for diffusions with piecewise constant coëfficients. Probability Theory & Related Fields 76, S. Chatterjee & S. Pal (2010) A phase transition behavior for Brownian motions interacting through their ranks. Probability Theory & Related Fields 147, / 26

26 E.R. Fernholz (2002) Stochastic Portfolio Theory. Springer-Verlag, New York. T. Ichiba, I. Karatzas & M. Shkolnikov (2013) Strong solutions of stochastic equations with rank-based coëfficients. Probability Theory & Related Fields 56, S. Pal & M. Shkolnikov (2013) Concentration of measure for Brownian particle systems interacting through their ranks. Probability Theory & Related Fields, to appear. M. Shkolnikov (2012) Large systems of diffusions interacting through their ranks. Stochastic Processes & Their Applications 122, / 26

Hybrid Atlas Model. of financial equity market. Tomoyuki Ichiba 1 Ioannis Karatzas 2,3 Adrian Banner 3 Vassilios Papathanakos 3 Robert Fernholz 3

Hybrid Atlas Model. of financial equity market. Tomoyuki Ichiba 1 Ioannis Karatzas 2,3 Adrian Banner 3 Vassilios Papathanakos 3 Robert Fernholz 3 Hybrid Atlas Model of financial equity market Tomoyuki Ichiba 1 Ioannis Karatzas 2,3 Adrian Banner 3 Vassilios Papathanakos 3 Robert Fernholz 3 1 University of California, Santa Barbara 2 Columbia University,

More information

Stable Diffusions Interacting through Their Ranks, as Models of Large Equity Markets

Stable Diffusions Interacting through Their Ranks, as Models of Large Equity Markets 1 Stable Diffusions Interacting through Their Ranks, as Models of Large Equity Markets IOANNIS KARATZAS INTECH Investment Management, Princeton Department of Mathematics, Columbia University Joint work

More information

Some Aspects of Universal Portfolio

Some Aspects of Universal Portfolio 1 Some Aspects of Universal Portfolio Tomoyuki Ichiba (UC Santa Barbara) joint work with Marcel Brod (ETH Zurich) Conference on Stochastic Asymptotics & Applications Sixth Western Conference on Mathematical

More information

LOCAL TIMES OF RANKED CONTINUOUS SEMIMARTINGALES

LOCAL TIMES OF RANKED CONTINUOUS SEMIMARTINGALES LOCAL TIMES OF RANKED CONTINUOUS SEMIMARTINGALES ADRIAN D. BANNER INTECH One Palmer Square Princeton, NJ 8542, USA adrian@enhanced.com RAOUF GHOMRASNI Fakultät II, Institut für Mathematik Sekr. MA 7-5,

More information

A Second-Order Model of the Stock Market

A Second-Order Model of the Stock Market A Second-Order Model of the Stock Market Robert Fernholz INTECH Conference in honor of Ioannis Karatzas Columbia University, New York June 4 8, 2012 This is joint research with T. Ichiba, I. Karatzas.

More information

IOANNIS KARATZAS Mathematics and Statistics Departments Columbia University

IOANNIS KARATZAS Mathematics and Statistics Departments Columbia University STOCHASTIC PORTFOLIO THEORY IOANNIS KARATZAS Mathematics and Statistics Departments Columbia University ik@math.columbia.edu Joint work with Dr. E. Robert FERNHOLZ, C.I.O. of INTECH Enhanced Investment

More information

Competing Brownian Particles

Competing Brownian Particles Competing Brownian Particles Andrey Sarantsev University of California, Santa Barbara February 9, 2018 Andrey Sarantsev California State University, Los Angeles 1 / 16 Brownian Motion A Brownian motion

More information

EXPLOSIONS AND ARBITRAGE IOANNIS KARATZAS. Joint work with Daniel FERNHOLZ and Johannes RUF. Talk at the University of Michigan

EXPLOSIONS AND ARBITRAGE IOANNIS KARATZAS. Joint work with Daniel FERNHOLZ and Johannes RUF. Talk at the University of Michigan EXPLOSIONS AND ARBITRAGE IOANNIS KARATZAS Department of Mathematics, Columbia University, New York INTECH Investment Management LLC, Princeton Joint work with Daniel FERNHOLZ and Johannes RUF Talk at the

More information

TRANSPORT INEQUALITIES FOR STOCHASTIC

TRANSPORT INEQUALITIES FOR STOCHASTIC TRANSPORT INEQUALITIES FOR STOCHASTIC PROCESSES Soumik Pal University of Washington, Seattle Jun 6, 2012 INTRODUCTION Three problems about rank-based processes THE ATLAS MODEL Define ranks: x (1) x (2)...

More information

Convex polytopes, interacting particles, spin glasses, and finance

Convex polytopes, interacting particles, spin glasses, and finance Convex polytopes, interacting particles, spin glasses, and finance Sourav Chatterjee (UCB) joint work with Soumik Pal (Cornell) Interacting particles Systems of particles governed by joint stochastic differential

More information

Numerical Methods with Lévy Processes

Numerical Methods with Lévy Processes Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:

More information

Geometry of functionally generated portfolios

Geometry of functionally generated portfolios Geometry of functionally generated portfolios Soumik Pal University of Washington Rutgers MF-PDE May 18, 2017 Multiplicative Cyclical Monotonicity Portfolio as a function on the unit simplex -unitsimplexindimensionn

More information

On a class of stochastic differential equations in a financial network model

On a class of stochastic differential equations in a financial network model 1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University

More information

Walsh Diffusions. Andrey Sarantsev. March 27, University of California, Santa Barbara. Andrey Sarantsev University of Washington, Seattle 1 / 1

Walsh Diffusions. Andrey Sarantsev. March 27, University of California, Santa Barbara. Andrey Sarantsev University of Washington, Seattle 1 / 1 Walsh Diffusions Andrey Sarantsev University of California, Santa Barbara March 27, 2017 Andrey Sarantsev University of Washington, Seattle 1 / 1 Walsh Brownian Motion on R d Spinning measure µ: probability

More information

The multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications

The multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications The multidimensional Ito Integral and the multidimensional Ito Formula Eric Mu ller June 1, 215 Seminar on Stochastic Geometry and its applications page 2 Seminar on Stochastic Geometry and its applications

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

University Of Calgary Department of Mathematics and Statistics

University Of Calgary Department of Mathematics and Statistics University Of Calgary Department of Mathematics and Statistics Hawkes Seminar May 23, 2018 1 / 46 Some Problems in Insurance and Reinsurance Mohamed Badaoui Department of Electrical Engineering National

More information

Branching Processes II: Convergence of critical branching to Feller s CSB

Branching Processes II: Convergence of critical branching to Feller s CSB Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied

More information

arxiv: v1 [q-fin.mf] 26 Sep 2018

arxiv: v1 [q-fin.mf] 26 Sep 2018 Trading Strategies Generated by Path-dependent Functionals of Market Weights IOANNIS KARATZAS DONGHAN KIM arxiv:189.1123v1 [q-fin.mf] 26 Sep 218 September 27, 218 Abstract Almost twenty years ago, E.R.

More information

Optimal portfolio strategies under partial information with expert opinions

Optimal portfolio strategies under partial information with expert opinions 1 / 35 Optimal portfolio strategies under partial information with expert opinions Ralf Wunderlich Brandenburg University of Technology Cottbus, Germany Joint work with Rüdiger Frey Research Seminar WU

More information

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,

More information

ALEKSANDAR MIJATOVIĆ AND MIKHAIL URUSOV

ALEKSANDAR MIJATOVIĆ AND MIKHAIL URUSOV DETERMINISTIC CRITERIA FOR THE ABSENCE OF ARBITRAGE IN DIFFUSION MODELS ALEKSANDAR MIJATOVIĆ AND MIKHAIL URUSOV Abstract. We obtain a deterministic characterisation of the no free lunch with vanishing

More information

Geometry and Optimization of Relative Arbitrage

Geometry and Optimization of Relative Arbitrage Geometry and Optimization of Relative Arbitrage Ting-Kam Leonard Wong joint work with Soumik Pal Department of Mathematics, University of Washington Financial/Actuarial Mathematic Seminar, University of

More information

Tools of stochastic calculus

Tools of stochastic calculus slides for the course Interest rate theory, University of Ljubljana, 212-13/I, part III József Gáll University of Debrecen Nov. 212 Jan. 213, Ljubljana Itô integral, summary of main facts Notations, basic

More information

Affine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009

Affine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009 Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine

More information

SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM. 1. Introduction This paper discusses arbitrage-free separable term structure (STS) models

SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM. 1. Introduction This paper discusses arbitrage-free separable term structure (STS) models SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM DAMIR FILIPOVIĆ Abstract. This paper discusses separable term structure diffusion models in an arbitrage-free environment. Using general consistency

More information

Stochastic Calculus February 11, / 33

Stochastic Calculus February 11, / 33 Martingale Transform M n martingale with respect to F n, n =, 1, 2,... σ n F n (σ M) n = n 1 i= σ i(m i+1 M i ) is a Martingale E[(σ M) n F n 1 ] n 1 = E[ σ i (M i+1 M i ) F n 1 ] i= n 2 = σ i (M i+1 M

More information

Stochastic Volatility and Correction to the Heat Equation

Stochastic Volatility and Correction to the Heat Equation Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century

More information

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009 A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance

More information

RESEARCH STATEMENT ANDREY SARANTSEV

RESEARCH STATEMENT ANDREY SARANTSEV RESEARCH STATEMENT ANDREY SARANTSEV My research in Stochastic Finance and Stochastic Analysis is focused on the following topics: 1. Interacting Brownian particles with rank-dependent dynamics. This model

More information

Translation Invariant Experiments with Independent Increments

Translation Invariant Experiments with Independent Increments Translation Invariant Statistical Experiments with Independent Increments (joint work with Nino Kordzakhia and Alex Novikov Steklov Mathematical Institute St.Petersburg, June 10, 2013 Outline 1 Introduction

More information

March 16, Abstract. We study the problem of portfolio optimization under the \drawdown constraint" that the

March 16, Abstract. We study the problem of portfolio optimization under the \drawdown constraint that the ON PORTFOLIO OPTIMIZATION UNDER \DRAWDOWN" CONSTRAINTS JAKSA CVITANIC IOANNIS KARATZAS y March 6, 994 Abstract We study the problem of portfolio optimization under the \drawdown constraint" that the wealth

More information

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3 Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................

More information

Supermodular ordering of Poisson arrays

Supermodular ordering of Poisson arrays Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore

More information

Order book resilience, price manipulation, and the positive portfolio problem

Order book resilience, price manipulation, and the positive portfolio problem Order book resilience, price manipulation, and the positive portfolio problem Alexander Schied Mannheim University Workshop on New Directions in Financial Mathematics Institute for Pure and Applied Mathematics,

More information

On countably skewed Brownian motion

On countably skewed Brownian motion On countably skewed Brownian motion Gerald Trutnau (Seoul National University) Joint work with Y. Ouknine (Cadi Ayyad) and F. Russo (ENSTA ParisTech) Electron. J. Probab. 20 (2015), no. 82, 1-27 [ORT 2015]

More information

Exact Simulation of Multivariate Itô Diffusions

Exact Simulation of Multivariate Itô Diffusions Exact Simulation of Multivariate Itô Diffusions Jose Blanchet Joint work with Fan Zhang Columbia and Stanford July 7, 2017 Jose Blanchet (Columbia/Stanford) Exact Simulation of Diffusions July 7, 2017

More information

Generalized Gaussian Bridges of Prediction-Invertible Processes

Generalized Gaussian Bridges of Prediction-Invertible Processes Generalized Gaussian Bridges of Prediction-Invertible Processes Tommi Sottinen 1 and Adil Yazigi University of Vaasa, Finland Modern Stochastics: Theory and Applications III September 1, 212, Kyiv, Ukraine

More information

Some Tools From Stochastic Analysis

Some Tools From Stochastic Analysis W H I T E Some Tools From Stochastic Analysis J. Potthoff Lehrstuhl für Mathematik V Universität Mannheim email: potthoff@math.uni-mannheim.de url: http://ls5.math.uni-mannheim.de To close the file, click

More information

Ernesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010

Ernesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010 Optimal stopping for Hunt and Lévy processes Ernesto Mordecki 1 Lecture III. PASI - Guanajuato - June 2010 1Joint work with Paavo Salminen (Åbo, Finland) 1 Plan of the talk 1. Motivation: from Finance

More information

Stochastic Portfolio Theory: An Overview

Stochastic Portfolio Theory: An Overview Stochastic Portfolio Theory: An Overview ROBERT FERNHOLZ INTECH One Palmer Square Princeton, NJ 8542, USA bob@enhanced.com IOANNIS KARATZAS Departments of Mathematics and Statistics Columbia University

More information

DIVERSITY AND RELATIVE ARBITRAGE IN EQUITY MARKETS

DIVERSITY AND RELATIVE ARBITRAGE IN EQUITY MARKETS DIVERSITY AND RELATIVE ARBITRAGE IN EQUITY MARKETS ROBERT FERNHOLZ INTECH, One Palmer Square Princeton, NJ 8542 bob@enhanced.com IOANNIS KARATZAS Departments of Mathematics and Statistics Columbia University,

More information

1. Stochastic Processes and filtrations

1. Stochastic Processes and filtrations 1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S

More information

An Analytic Method for Solving Uncertain Differential Equations

An Analytic Method for Solving Uncertain Differential Equations Journal of Uncertain Systems Vol.6, No.4, pp.244-249, 212 Online at: www.jus.org.uk An Analytic Method for Solving Uncertain Differential Equations Yuhan Liu Department of Industrial Engineering, Tsinghua

More information

OPTIONAL DECOMPOSITION FOR CONTINUOUS SEMIMARTINGALES UNDER ARBITRARY FILTRATIONS. Introduction

OPTIONAL DECOMPOSITION FOR CONTINUOUS SEMIMARTINGALES UNDER ARBITRARY FILTRATIONS. Introduction OPTIONAL DECOMPOSITION FOR CONTINUOUS SEMIMARTINGALES UNDER ARBITRARY FILTRATIONS IOANNIS KARATZAS AND CONSTANTINOS KARDARAS Abstract. We present an elementary treatment of the Optional Decomposition Theorem

More information

An Uncertain Control Model with Application to. Production-Inventory System

An Uncertain Control Model with Application to. Production-Inventory System An Uncertain Control Model with Application to Production-Inventory System Kai Yao 1, Zhongfeng Qin 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 2 School of Economics

More information

Trading strategies generated by Lyapunov functions

Trading strategies generated by Lyapunov functions Finance Stoch 27) 2:753 787 DOI.7/s78-7-332-8 Trading strategies generated by Lyapunov functions Ioannis Karatzas,2 Johannes Ruf 3 Received: 2 July 26 / Accepted: 3 January 27 / Published online: 25 May

More information

{σ x >t}p x. (σ x >t)=e at.

{σ x >t}p x. (σ x >t)=e at. 3.11. EXERCISES 121 3.11 Exercises Exercise 3.1 Consider the Ornstein Uhlenbeck process in example 3.1.7(B). Show that the defined process is a Markov process which converges in distribution to an N(0,σ

More information

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias José E. Figueroa-López 1 1 Department of Statistics Purdue University Seoul National University & Ajou University

More information

On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem

On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family

More information

Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process

Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Mátyás Barczy University of Debrecen, Hungary The 3 rd Workshop on Branching Processes and Related

More information

Poisson random measure: motivation

Poisson random measure: motivation : motivation The Lévy measure provides the expected number of jumps by time unit, i.e. in a time interval of the form: [t, t + 1], and of a certain size Example: ν([1, )) is the expected number of jumps

More information

Bessel-like SPDEs. Lorenzo Zambotti, Sorbonne Université (joint work with Henri Elad-Altman) 15th May 2018, Luminy

Bessel-like SPDEs. Lorenzo Zambotti, Sorbonne Université (joint work with Henri Elad-Altman) 15th May 2018, Luminy Bessel-like SPDEs, Sorbonne Université (joint work with Henri Elad-Altman) Squared Bessel processes Let δ, y, and (B t ) t a BM. By Yamada-Watanabe s Theorem, there exists a unique (strong) solution (Y

More information

Model-free portfolio theory and its functional master formula

Model-free portfolio theory and its functional master formula Model-free portfolio theory and its functional master formula arxiv:66.335v [q-fin.pm] 7 Jun 6 Alexander Schied, Leo Speiser, and Iryna Voloshchenko Department of Mathematics University of Mannheim 683

More information

Optimal exit strategies for investment projects. 7th AMaMeF and Swissquote Conference

Optimal exit strategies for investment projects. 7th AMaMeF and Swissquote Conference Optimal exit strategies for investment projects Simone Scotti Université Paris Diderot Laboratoire de Probabilité et Modèles Aléatories Joint work with : Etienne Chevalier, Université d Evry Vathana Ly

More information

Economics 672 Fall 2017 Tauchen. Jump Regression

Economics 672 Fall 2017 Tauchen. Jump Regression Economics 672 Fall 2017 Tauchen 1 Main Model In the jump regression setting we have Jump Regression X = ( Z Y where Z is the log of the market index and Y is the log of an asset price. The dynamics are

More information

STA 256: Statistics and Probability I

STA 256: Statistics and Probability I Al Nosedal. University of Toronto. Fall 2017 My momma always said: Life was like a box of chocolates. You never know what you re gonna get. Forrest Gump. There are situations where one might be interested

More information

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Robert J. Elliott 1 Samuel N. Cohen 2 1 Department of Commerce, University of South Australia 2 Mathematical Insitute, University

More information

I forgot to mention last time: in the Ito formula for two standard processes, putting

I forgot to mention last time: in the Ito formula for two standard processes, putting I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy

More information

1 Brownian Local Time

1 Brownian Local Time 1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =

More information

VOLATILITY AND ARBITRAGE

VOLATILITY AND ARBITRAGE The Annals of Applied Probability 218, Vol. 28, No. 1, 378 417 https://doi.org/1.1214/17-aap138 Institute of Mathematical Statistics, 218 VOLATILITY AND ARBITRAGE BY E. ROBERT FERNHOLZ,IOANNIS KARATZAS,,1

More information

Stat410 Probability and Statistics II (F16)

Stat410 Probability and Statistics II (F16) Stat4 Probability and Statistics II (F6 Exponential, Poisson and Gamma Suppose on average every /λ hours, a Stochastic train arrives at the Random station. Further we assume the waiting time between two

More information

Lecture 4: Introduction to stochastic processes and stochastic calculus

Lecture 4: Introduction to stochastic processes and stochastic calculus Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London

More information

2012 NCTS Workshop on Dynamical Systems

2012 NCTS Workshop on Dynamical Systems Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,

More information

Conditional Distributions

Conditional Distributions Conditional Distributions The goal is to provide a general definition of the conditional distribution of Y given X, when (X, Y ) are jointly distributed. Let F be a distribution function on R. Let G(,

More information

Mean-field dual of cooperative reproduction

Mean-field dual of cooperative reproduction The mean-field dual of systems with cooperative reproduction joint with Tibor Mach (Prague) A. Sturm (Göttingen) Friday, July 6th, 2018 Poisson construction of Markov processes Let (X t ) t 0 be a continuous-time

More information

Generalized Gaussian Bridges

Generalized Gaussian Bridges TOMMI SOTTINEN ADIL YAZIGI Generalized Gaussian Bridges PROCEEDINGS OF THE UNIVERSITY OF VAASA WORKING PAPERS 4 MATHEMATICS 2 VAASA 212 Vaasan yliopisto University of Vaasa PL 7 P.O. Box 7 (Wolffintie

More information

FE 5204 Stochastic Differential Equations

FE 5204 Stochastic Differential Equations Instructor: Jim Zhu e-mail:zhu@wmich.edu http://homepages.wmich.edu/ zhu/ January 20, 2009 Preliminaries for dealing with continuous random processes. Brownian motions. Our main reference for this lecture

More information

Multi-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes

Multi-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes Multi-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes Anatoliy Swishchuk Department of Mathematics and Statistics University of Calgary Calgary, Alberta, Canada Lunch at the Lab Talk

More information

Simulation methods for stochastic models in chemistry

Simulation methods for stochastic models in chemistry Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation

More information

Local vs. Nonlocal Diffusions A Tale of Two Laplacians

Local vs. Nonlocal Diffusions A Tale of Two Laplacians Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2

More information

Exercises. T 2T. e ita φ(t)dt.

Exercises. T 2T. e ita φ(t)dt. Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.

More information

HJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011

HJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011 Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance

More information

Representing Gaussian Processes with Martingales

Representing Gaussian Processes with Martingales Representing Gaussian Processes with Martingales with Application to MLE of Langevin Equation Tommi Sottinen University of Vaasa Based on ongoing joint work with Lauri Viitasaari, University of Saarland

More information

A class of globally solvable systems of BSDE and applications

A class of globally solvable systems of BSDE and applications A class of globally solvable systems of BSDE and applications Gordan Žitković Department of Mathematics The University of Texas at Austin Thera Stochastics - Wednesday, May 31st, 2017 joint work with Hao

More information

Polynomial Preserving Jump-Diffusions on the Unit Interval

Polynomial Preserving Jump-Diffusions on the Unit Interval Polynomial Preserving Jump-Diffusions on the Unit Interval Martin Larsson Department of Mathematics, ETH Zürich joint work with Christa Cuchiero and Sara Svaluto-Ferro 7th General AMaMeF and Swissquote

More information

Stability and Sensitivity of the Capacity in Continuous Channels. Malcolm Egan

Stability and Sensitivity of the Capacity in Continuous Channels. Malcolm Egan Stability and Sensitivity of the Capacity in Continuous Channels Malcolm Egan Univ. Lyon, INSA Lyon, INRIA 2019 European School of Information Theory April 18, 2019 1 / 40 Capacity of Additive Noise Models

More information

Filtrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition

Filtrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition Filtrations, Markov Processes and Martingales Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition David pplebaum Probability and Statistics Department,

More information

Jump-type Levy Processes

Jump-type Levy Processes Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic

More information

Optimal investment with high-watermark fee in a multi-dimensional jump diffusion model

Optimal investment with high-watermark fee in a multi-dimensional jump diffusion model Optimal investment with high-watermark fee in a multi-dimensional jump diffusion model Karel Janeček Zheng Li Mihai Sîrbu August 2, 218 Abstract This paper studies the problem of optimal investment and

More information

The Uses of Differential Geometry in Finance

The Uses of Differential Geometry in Finance The Uses of Differential Geometry in Finance Andrew Lesniewski Bloomberg, November 21 2005 The Uses of Differential Geometry in Finance p. 1 Overview Joint with P. Hagan and D. Woodward Motivation: Varadhan

More information

Properties of an infinite dimensional EDS system : the Muller s ratchet

Properties of an infinite dimensional EDS system : the Muller s ratchet Properties of an infinite dimensional EDS system : the Muller s ratchet LATP June 5, 2011 A ratchet source : wikipedia Plan 1 Introduction : The model of Haigh 2 3 Hypothesis (Biological) : The population

More information

Thomas Knispel Leibniz Universität Hannover

Thomas Knispel Leibniz Universität Hannover Optimal long term investment under model ambiguity Optimal long term investment under model ambiguity homas Knispel Leibniz Universität Hannover knispel@stochastik.uni-hannover.de AnStAp0 Vienna, July

More information

Economics 883 Spring 2016 Tauchen. Jump Regression

Economics 883 Spring 2016 Tauchen. Jump Regression Economics 883 Spring 2016 Tauchen Jump Regression 1 Main Model In the jump regression setting we have X = ( Z Y where Z is the log of the market index and Y is the log of an asset price. The dynamics are

More information

Exact Simulation of Diffusions and Jump Diffusions

Exact Simulation of Diffusions and Jump Diffusions Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction

More information

Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations.

Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics

More information

A. Bovier () Branching Brownian motion: extremal process and ergodic theorems

A. Bovier () Branching Brownian motion: extremal process and ergodic theorems Branching Brownian motion: extremal process and ergodic theorems Anton Bovier with Louis-Pierre Arguin and Nicola Kistler RCS&SM, Venezia, 06.05.2013 Plan 1 BBM 2 Maximum of BBM 3 The Lalley-Sellke conjecture

More information

Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model

Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model Xiaowei Chen International Business School, Nankai University, Tianjin 371, China School of Finance, Nankai

More information

Notes, March 4, 2013, R. Dudley Maximum likelihood estimation: actual or supposed

Notes, March 4, 2013, R. Dudley Maximum likelihood estimation: actual or supposed 18.466 Notes, March 4, 2013, R. Dudley Maximum likelihood estimation: actual or supposed 1. MLEs in exponential families Let f(x,θ) for x X and θ Θ be a likelihood function, that is, for present purposes,

More information

Multilevel Monte Carlo for Lévy Driven SDEs

Multilevel Monte Carlo for Lévy Driven SDEs Multilevel Monte Carlo for Lévy Driven SDEs Felix Heidenreich TU Kaiserslautern AG Computational Stochastics August 2011 joint work with Steffen Dereich Philipps-Universität Marburg supported within DFG-SPP

More information

Vast Volatility Matrix Estimation for High Frequency Data

Vast Volatility Matrix Estimation for High Frequency Data Vast Volatility Matrix Estimation for High Frequency Data Yazhen Wang National Science Foundation Yale Workshop, May 14-17, 2009 Disclaimer: My opinion, not the views of NSF Y. Wang (at NSF) 1 / 36 Outline

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

MFM Practitioner Module: Quantitative Risk Management. John Dodson. September 23, 2015

MFM Practitioner Module: Quantitative Risk Management. John Dodson. September 23, 2015 MFM Practitioner Module: Quantitative Risk Management September 23, 2015 Mixtures Mixtures Mixtures Definitions For our purposes, A random variable is a quantity whose value is not known to us right now

More information

Prof. Erhan Bayraktar (University of Michigan)

Prof. Erhan Bayraktar (University of Michigan) September 17, 2012 KAP 414 2:15 PM- 3:15 PM Prof. (University of Michigan) Abstract: We consider a zero-sum stochastic differential controller-and-stopper game in which the state process is a controlled

More information

On Minimal Entropy Martingale Measures

On Minimal Entropy Martingale Measures On Minimal Entropy Martingale Measures Andrii Andrusiv Friedrich Schiller University of Jena, Germany Marie Curie ITN Outline 1. One Period Model Definitions and Notations Minimal Entropy Martingale Measure

More information

Stochastic differential equation models in biology Susanne Ditlevsen

Stochastic differential equation models in biology Susanne Ditlevsen Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential

More information

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral

More information

OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS

OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS APPLICATIONES MATHEMATICAE 29,4 (22), pp. 387 398 Mariusz Michta (Zielona Góra) OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS Abstract. A martingale problem approach is used first to analyze

More information

SDE Coefficients. March 4, 2008

SDE Coefficients. March 4, 2008 SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)

More information

THE MALLIAVIN CALCULUS FOR SDE WITH JUMPS AND THE PARTIALLY HYPOELLIPTIC PROBLEM

THE MALLIAVIN CALCULUS FOR SDE WITH JUMPS AND THE PARTIALLY HYPOELLIPTIC PROBLEM Takeuchi, A. Osaka J. Math. 39, 53 559 THE MALLIAVIN CALCULUS FOR SDE WITH JUMPS AND THE PARTIALLY HYPOELLIPTIC PROBLEM ATSUSHI TAKEUCHI Received October 11, 1. Introduction It has been studied by many

More information