Hybrid Atlas Model. of financial equity market. Tomoyuki Ichiba 1 Ioannis Karatzas 2,3 Adrian Banner 3 Vassilios Papathanakos 3 Robert Fernholz 3
|
|
- Randolph Chandler
- 5 years ago
- Views:
Transcription
1 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, New York 3 INTECH, Princeton November 2009
2 2 Outline Introduction Hybrid Atlas model Martingale Problem Stability Effective dimension Rankings Long-term growth relations Portfolio analysis Stochastic Portfolio Theory Target portfolio Universal portfolio Conclusion
3 3 Flow of Capital Figure: Capital Distribution Curves (Percentage) for the S&P 500 Index of 1997 (Solid Line) and 1999 (Broken Line).
4 Log-Log Capital Distribution Curves Figure: Capital distribution curves for 1929 (shortest curve) (longest curve), every ten years. Source Fernholz( 02). What kind of models can describe this long-term stability?
5 A Model of Rankings [Hybrid Atlas model] Capital process X := {(X 1 (t),..., X n (t)), 0 t < }. Order Statistics: X (1) (t) X (n) (t) ; 0 t <. Log capital Y := log X : Y (1) (t) Y (n) (t) ; 0 t <. Dynamics of log capital: d Y (k) (t) = (γ + γ i + g k ) d t + σ k d W i (t) if Y (k) (t) = Y i (t) ; for 1 i, k n, 0 t <, where W ( ) is n dim. B. M. company name i kth ranked company Drift ( mean ) γ i g k Diffusion ( variance ) σ k > 0 Banner, Fernholz & Karatzas ( 05), Chatterjee & Pal ( 07, 09), Pal & Pitman ( 08). 5
6 A Model of Rankings [Hybrid Atlas model] Capital process X := {(X 1 (t),..., X n (t)), 0 t < }. Order Statistics: X (1) (t) X (n) (t) ; 0 t <. Log capital Y := log X : Y (1) (t) Y (n) (t) ; 0 t <. Dynamics of log capital: d Y (k) (t) = (γ + γ i + g k ) d t + σ k d W i (t) if Y (k) (t) = Y i (t) ; for 1 i, k n, 0 t <, where W ( ) is n dim. B. M. company name i kth ranked company Drift ( mean ) γ i g k Diffusion ( variance ) σ k > 0 Banner, Fernholz & Karatzas ( 05), Chatterjee & Pal ( 07, 09), Pal & Pitman ( 08). 5
7 A Model of Rankings [Hybrid Atlas model] Capital process X := {(X 1 (t),..., X n (t)), 0 t < }. Order Statistics: X (1) (t) X (n) (t) ; 0 t <. Log capital Y := log X : Y (1) (t) Y (n) (t) ; 0 t <. Dynamics of log capital: d Y (k) (t) = (γ + γ i + g k ) d t + σ k d W i (t) if Y (k) (t) = Y i (t) ; for 1 i, k n, 0 t <, where W ( ) is n dim. B. M. company name i kth ranked company Drift ( mean ) γ i g k Diffusion ( variance ) σ k > 0 Banner, Fernholz & Karatzas ( 05), Chatterjee & Pal ( 07, 09), Pal & Pitman ( 08). 5
8 6 Illustration (n = 3) of interactions through rank Time X 3 (t)>x 1 (t)>x 2 (t) X 1 (t)>x 3 (t)>x 2 (t) x 1 > x 3 > x 2 x 3 > x 1 > x 2 x 1 > x 2 > x 3 X 1 (t)>x 2 (t)>x 3 (t) x 3 > x 2 > x 1 x 2 > x 1 > x 3 X 3 X 2 X 1 x 2 > x 3 > x 1 Paths in R + Time. A path in different wedges of R n. Symmetric group Σ n of permutations of {1,..., n}. For n = 3, {(1, 2, 3), (1, 3, 2), (2, 1, 3), (2, 3, 1), (3, 1, 2), (3, 2, 1)}.
9 6 Illustration (n = 3) of interactions through rank Time X 3 (t)>x 1 (t)>x 2 (t) X 1 (t)>x 3 (t)>x 2 (t) x 1 > x 3 > x 2 x 3 > x 1 > x 2 x 1 > x 2 > x 3 X 1 (t)>x 2 (t)>x 3 (t) x 3 > x 2 > x 1 x 2 > x 1 > x 3 X 3 X 2 X 1 x 2 > x 3 > x 1 Paths in R + Time. A path in different wedges of R n. Symmetric group Σ n of permutations of {1,..., n}. For n = 3, {(1, 2, 3), (1, 3, 2), (2, 1, 3), (2, 3, 1), (3, 1, 2), (3, 2, 1)}.
10 Vector Representation d Y (t) = G(Y (t))d t + S(Y (t))dw (t) ; 0 t < Σ n : symmetric group of permutations of {1, 2,..., n}. For p Σ n define wedges (chambers) R p := {x R n : x p(1) x p(2) x p(n) }, R n = p Σn R p, (the inner points of R p and R p are disjoint for p p Σ n ), Q (i) k : = {x R n : x i is ranked kth among (x 1,..., x n )} = {p : p(k)=i} R p ; 1 i, k n, n j=1 Q(j) k = R n = n l=1 Q(i) l and R p = n k=1 Q(p(k)) k. G(y) = p Σ n (g p 1 (1) + γ 1 + γ,..., g p 1 (n) + γ n + γ) 1 Rp (y), S(y) = p Σ n diag(σ p 1 (1),..., σ p 1 (n) ) }{{} s p 1 Rp (y) ; y Rn. 7
11 Vector Representation d Y (t) = G(Y (t))d t + S(Y (t))dw (t) ; 0 t < Σ n : symmetric group of permutations of {1, 2,..., n}. For p Σ n define wedges (chambers) R p := {x R n : x p(1) x p(2) x p(n) }, R n = p Σn R p, (the inner points of R p and R p are disjoint for p p Σ n ), Q (i) k : = {x R n : x i is ranked kth among (x 1,..., x n )} = {p : p(k)=i} R p ; 1 i, k n, n j=1 Q(j) k = R n = n l=1 Q(i) l and R p = n k=1 Q(p(k)) k. G(y) = p Σ n (g p 1 (1) + γ 1 + γ,..., g p 1 (n) + γ n + γ) 1 Rp (y), S(y) = p Σ n diag(σ p 1 (1),..., σ p 1 (n) ) }{{} s p 1 Rp (y) ; y Rn. 7
12 8 Martingale Problem Theorem [Krylov( 71), Stroock & Varadhan( 79), Bass & Pardoux( 87)] Suppose that the coefficients G( ) and a( ) := SS ( ) are bounded and measurable, and that a( ) is uniformly positive-definite and piecewise constant in each wedge. For each y 0 R n there is a unique one probability measure P on C([0, ), R n ) such that P(Y 0 = y 0 ) = 1 and t f (Y t ) f (Y 0 ) L f (Y s ) d s ; 0 0 t < is a P local martingale for every f C 2 (R 2 ) where L f (x) = 1 n 2 i,j=1 a ij(x)d ij f (x) + n i=1 G i(x)d i f (x) ; x R n. This implies that the hybrid Atlas model is well-defined.
13 9 Model Market capitalization X follows Hybrid Atlas model: the log capitalization Y i = log X i of company i has drift γ + g k + γ i and volatility σ k, when company i is k th ranked, i.e., Y Q (i) k for 1 k, i n. d Y i (t) = ( γ + + n k=1 n k=1 g k 1 (i) Q (Y (t)) + γ i )d t k σ k 1 (i) Q (Y (t))d W i (t) ; 0 t <. k
14 10 Model assumptions Market capitalization X follows Hybrid Atlas model: the log capitalization Y i = log X i of company i has drift γ + g k + γ i and volatility σ k, when company i is k th ranked, i.e., Y Q (i) k for 1 k, i n. Assume σ k > 0, (g k, 1 k n), (γ i, 1 i n) and γ are real constants with stability conditions n n g k + γ i = 0, k=1 i=1 k (g l +γ p(l) ) < 0, k = 1,..., n 1, p Σ n. l=1 γ i = 0, 1 i n, g 1 = = g n 1 = g < 0, g n = (n 1)g > 0. γ i = 1 (2i)/(n + 1), 1 i n, g k = 1, k = 1,..., n 1, g n = n 1.
15 10 Model assumptions Market capitalization X follows Hybrid Atlas model: the log capitalization Y i = log X i of company i has drift γ + g k + γ i and volatility σ k, when company i is k th ranked, i.e., Y Q (i) k for 1 k, i n. Assume σ k > 0, (g k, 1 k n), (γ i, 1 i n) and γ are real constants with stability conditions n n g k + γ i = 0, k=1 i=1 k (g l +γ p(l) ) < 0, k = 1,..., n 1, p Σ n. l=1 γ i = 0, 1 i n, g 1 = = g n 1 = g < 0, g n = (n 1)g > 0. γ i = 1 (2i)/(n + 1), 1 i n, g k = 1, k = 1,..., n 1, g n = n 1.
16 10 Model assumptions Market capitalization X follows Hybrid Atlas model: the log capitalization Y i = log X i of company i has drift γ + g k + γ i and volatility σ k, when company i is k th ranked, i.e., Y Q (i) k for 1 k, i n. Assume σ k > 0, (g k, 1 k n), (γ i, 1 i n) and γ are real constants with stability conditions n n g k + γ i = 0, k=1 i=1 k (g l +γ p(l) ) < 0, k = 1,..., n 1, p Σ n. l=1 γ i = 0, 1 i n, g 1 = = g n 1 = g < 0, g n = (n 1)g > 0. γ i = 1 (2i)/(n + 1), 1 i n, g k = 1, k = 1,..., n 1, g n = n 1.
17 Model Summary The log-capitalization Y = log X follows d Y i (t) = ( γ + + n k=1 n k=1 g k 1 (i) Q (Y (t)) + γ i )d t k σ k 1 (i) Q (Y (t))d W i (t) ; 0 t < k where σ k > 0, (g k, 1 k n), (γ i, 1 i n) and γ are real constants with stability conditions n n g k + γ i = 0, k=1 i=1 k (g l +γ p(l) ) < 0, k = 1,..., n 1, p Σ n. l=1
18 Stochastic stability The average Y ( ) := n i=1 Y i( ) / n of log-capitalization: dy (t) = γ d t + 1 n n n k=1 σ k 1 (i) Q (Y (t))d W i (t) k i=1 }{{} d B k (t) is a Brownian motion with variance rate n k=1 σ2 k /n2 drift γ by the Dambis-Dubins-Schwartz Theorem. Proposition Under the assumptions the deviations Ỹ ( ) := (Y 1 ( ) Y ( ),..., Y n ( ) Y ( )) from the average are stable in distribution, i.e., there is a unique invariant probability measure µ( ) such that for every bounded, measurable function f we have the Strong Law of Large Numbers 1 T lim f (Ỹ T T (t)) d t = f (y)µ(d y), 0 Π where Π := {y R n : y y n = 0}. a.s.
19 Stochastic stability The average Y ( ) := n i=1 Y i( ) / n of log-capitalization: dy (t) = γ d t + 1 n n n k=1 σ k 1 (i) Q (Y (t))d W i (t) k i=1 }{{} d B k (t) is a Brownian motion with variance rate n k=1 σ2 k /n2 drift γ by the Dambis-Dubins-Schwartz Theorem. Proposition Under the assumptions the deviations Ỹ ( ) := (Y 1 ( ) Y ( ),..., Y n ( ) Y ( )) from the average are stable in distribution, i.e., there is a unique invariant probability measure µ( ) such that for every bounded, measurable function f we have the Strong Law of Large Numbers 1 T lim f (Ỹ T T (t)) d t = f (y)µ(d y), 0 Π where Π := {y R n : y y n = 0}. a.s.
20 13 Average occupation times Especially taking f ( ) = 1 Rp ( ) or 1 (i) Q ( ), we define from µ k the average occupation time of X in R p or Q (i) k : 1 T θ p := µ(r p ) = lim 1 Rp (X(t))d t T T θ k,i := µ(q (i) k ) = lim T 1 T 0 T 0 1 (i) Q (X(t))d t, 1 k, i n, k since 1 Rp (Ỹ ( )) = 1 R p (X( )) and 1 (i) Q (X( )) = 1 (i) k Q (Ỹ ( )). By k definition 0 θ k,i = {p Σ n:p(k)=i} θ p 1 for 1 k, i n, n l=1 θ l,i = n j=1 θ k,j = 1 for 1 k, i n. What is the invariant distribution µ?
21 Attainability One-dimensional Brownian motion attains the origin infinitely often. Two-dimensional Brownian motion does not attain the origin. Does the process X( ) attain the origin? t t X(t) = X(0) + b(x(s)) d s + σ(x(s))d W (s) 0 0 where b and σ are bounded measurable functions. Friedman( 74), Bass & Pardoux( 87).
22 15 Effective Dimension Let us define effective dimension ED( ) by ED(x) = trace(a(x)) x 2 x A(x)x where A( ) = σ( )σ( ). ; x R n \ {0}, Proposition Suppose X(0) 0. If inf x R n \{0} ED(x) 2, then X( ) does not attain the origin. If sup x R n \{0} ED(x) < 2 and if there is no drift, i.e., b( ) 0, then X( ) attains the origin. Exterior Dirichlet Problem by Meyers and Serrin( 60). Removal of drift by Girsanov s theorem. If there is drift, take [trace(a(x))+x b(x)] x 2 x A(x)x.
23 16 Triple collision Now consider triple collision: { Xi (t) = X j (t) = X k (t) for some t > 0, 1 i < j < k n }. What is the probability of triple collision? Fix i = 1, j = 2, k = 3. Let us define the sum of squared distances: s 2 (x) := (x 1 x 2 ) 2 +(x 2 x 3 ) 2 +(x 3 x 1 ) 2 = x DD x ; x R n, where (n 3) matrix D is defined by D := (d 1, d 2, d 3 ) with d 1 := (1, 1, 0,..., 0), d 2 := (0, 1, 1, 0,..., 0), d 3 := ( 1, 0, 1,..., 0). Z := {x R n : s(x) = 0 }.
24 Define the local effective dimension: R(x) := trace(d A(x)D) x DD x x DD A(x)DD x ; x R n \ Z. Proposition Suppose s(x(0)) 0. If inf x R n \Z R(x) 2, then P(X 1 (t) = X 2 (t) = X 3 (t) for some t 0) = 0. If sup x R n \Z R(x) < 2 and if there is no drift, i.e., b( ) 0, then P(X 1 (t) = X 2 (t) = X 3 (t) for some t 0) = 1. R( ) 2 for n dim. BM, i.e., A( ) I. If there is drift, take [trace(d A(x)D)+x DD b(x)] x DD x x DD A(x)DD x. Idea of Proof: a comparison with Bessel process with dimension two. 17
25 Define the local effective dimension: R(x) := trace(d A(x)D) x DD x x DD A(x)DD x ; x R n \ Z. Proposition Suppose s(x(0)) 0. If inf x R n \Z R(x) 2, then P(X 1 (t) = X 2 (t) = X 3 (t) for some t 0) = 0. If sup x R n \Z R(x) < 2 and if there is no drift, i.e., b( ) 0, then P(X 1 (t) = X 2 (t) = X 3 (t) for some t 0) = 1. R( ) 2 for n dim. BM, i.e., A( ) I. If there is drift, take [trace(d A(x)D)+x DD b(x)] x DD x x DD A(x)DD x. Idea of Proof: a comparison with Bessel process with dimension two. 17
26 Define the local effective dimension: R(x) := trace(d A(x)D) x DD x x DD A(x)DD x ; x R n \ Z. Proposition Suppose s(x(0)) 0. If inf x R n \Z R(x) 2, then P(X 1 (t) = X 2 (t) = X 3 (t) for some t 0) = 0. If sup x R n \Z R(x) < 2 and if there is no drift, i.e., b( ) 0, then P(X 1 (t) = X 2 (t) = X 3 (t) for some t 0) = 1. R( ) 2 for n dim. BM, i.e., A( ) I. If there is drift, take [trace(d A(x)D)+x DD b(x)] x DD x x DD A(x)DD x. Idea of Proof: a comparison with Bessel process with dimension two. 17
27 18 Rankings Recall Y (1) ( ) Y (2) ( ) Y (n) ( ). Let us denote by Λ k,j (t) the local time accumulated at the origin by the nonnegative semimartingale Y (k) ( ) Y (j) ( ) up to time t for 1 k < j n. Theorem[Banner & Ghomrasni ( 07)] For a general class of semimartingale Y ( ), the rankings satisfy d Y (k) (t) = n i=1 1 (i) Q (Y (t))d Y i (t) k + (N k (t)) 1[ n l=k+1 k 1 ] d Λ k,l (t) d Λ l,k (t) l=1 where N k (t) is the cardinality {i : Y i (t) = Y (k) (t)}.
28 19 Rankings Recall Y (1) ( ) Y (2) ( ) Y (n) ( ). Let us denote by Λ k,j (t) the local time accumulated at the origin by the nonnegative semimartingale Y (k) ( ) Y (j) ( ) up to time t for 1 k < j n. Lemma Under the non-degeneracy condition σ k > 0 for k = 1,..., n, dy (k) (t) = ( γ + g k n i=1 for k = 1,..., n, 0 t T. ) γ i 1 (i) Q (Y (t)) d t + σ k d B k (t) k ] [ d Λ k,k+1 (t) d Λ k 1,k (t) Idea of Proof: a comparison with a Bessel process with dimension one to show Λ k,l ( ) 0, k l 2..
29 19 Rankings Recall Y (1) ( ) Y (2) ( ) Y (n) ( ). Let us denote by Λ k,j (t) the local time accumulated at the origin by the nonnegative semimartingale Y (k) ( ) Y (j) ( ) up to time t for 1 k < j n. Lemma Under the non-degeneracy condition σ k > 0 for k = 1,..., n, dy (k) (t) = ( γ + g k n i=1 for k = 1,..., n, 0 t T. ) γ i 1 (i) Q (Y (t)) d t + σ k d B k (t) k ] [ d Λ k,k+1 (t) d Λ k 1,k (t) Idea of Proof: a comparison with a Bessel process with dimension one to show Λ k,l ( ) 0, k l 2..
30 20 Long-term growth relations Proposition Under the assumptions we obtain the following long-term growth relations: Y i (T ) log X i (T ) lim = lim = γ T T T T Thus the model is coherent: = lim T log n i=1 X i(t ) T 1 lim T T log µ i(t ) = 0 a.s.; i = 1,..., n a.s. where µ i ( ) = X i ( )/(X 1 ( ) + + X n ( )). Moreover, n g k θ k,i + γ i = 0 ; i = 1,..., n. k=1
31 n g k θ k,i + γ i = 0 ; i = 1,..., n. k=1 The log-capitalization Y grows with rate γ and follows d Y i (t) = ( γ + + n k=1 n k=1 g k 1 (i) Q (Y (t)) + γ i )d t k σ k 1 (i) Q (Y (t))d W i (t) ; 0 t < k for i = 1,..., n. The ranking (Y (1) ( ),..., Y (n) ( )) follows dy (k) (t) = ( γ + g k n i=1 for k = 1,..., n, 0 t <. ) γ i 1 (i) Q (Y (t)) d t + σ k d B k (t) k ] [ d Λ k,k+1 (t) d Λ k 1,k (t).
32 Semimartingale reflected Brownian motions The adjacent differences (gaps) Ξ( ) := (Ξ 1 ( ),..., Ξ n ( )) where Ξ k ( ) := Y (k) ( ) Y (k+1) ( ) for k = 1,..., n 1 can be seen as a semimartingale reflected Brownian motion (SRBM): Ξ(t) = Ξ(0) + ζ(t) + (I n Q)Λ(t) where ζ( ) := (ζ 1 ( ),..., ζ n ( )), Λ( ) := (Λ 1,2 ( ),..., Λ n 1,n ( )), ζ k ( ) := n i=1 0 1 (i) Q (Y (s)) d Y (s) k n i=1 0 1 (i) Q (Y (s)) d Y (s) k+1 for k = 1,..., n 1, and Q is an (n 1) (n 1) matrix with elements
33 23 Q := 0 1/2 1/2 0 1/2 1/ /2 1/2 0 Thus the gaps Ξ k := Y (k) ( ) Y (k+1) ( ) follow. Ξ(t) = Ξ(0) + ζ(t) }{{} + (I n Q)Λ(t) }{{} semimartingale reflection part In order to study the invariant measure µ, we apply the theory of semimartingale reflected Brownian motions developed by M. Harrison, M. Reiman, R. Williams and others. In addition to the model assumptions, we assume linearly growing variances: σ 2 2 σ2 1 = σ2 3 σ2 2 = = σ2 n σ 2 n 1.
34 Invariant distribution of gaps and index Let us define the indicator map R n x p x Σ n such that x p x (1) x p x (2) x p(n), and the index process P t := p Y (t). Proposition Under the stability and the linearly growing variance conditions the invariant distribution ν( ) of (Ξ( ), P ) is ( ν(a B) = n 1 q Σ n k=1 ) 1 λ 1 q,k exp( λ p, z )d z p Σ A n for every measurable set A B where λ p := (λ p,1,..., λ p,n 1 ) is the vector of components λ p,k := 4( k l=1 g l + γ p(l) ) σk 2 + > 0 ; p Σ n, 1 k n 1. σ2 k+1 Proof: an extension from M. Harrison and R. Williams ( 87).
35 25 Average occupation time Corollary The average occupation times are ( θ p = n 1 q Σ n k=1 λ 1 q,k ) 1 n 1 j=1 λ 1 p,j and θ k,i = {p Σ n:p(k)=i} θ p for p Σ n and 1 k, i n. (1,10) (10,10) 0.15 If all γi = 0 and σ 2 1 = σ2 n, then θ k,i = n 1 for 1 k, i n. Heat map of θk,i when n = 10, σ 2 k = 1 + k, g k = 1 for k = 1,..., 9, g 10 = 9, and γ i = 1 (2i)/(n + 1) for i = 1,..., n (k,i) (1,1) (10,1)
36 26 Market weights come from Pareto type Corollary The joint invariant distribution of market shares µ (i) ( ) := X (i) ( )/(X 1 ( ) + + X n ( )) ; i = 1,..., n has the density (m 1,..., m n 1 ) = λ p,1 λ p,n 1 θ p p Π n m λ p,1+1 1 m λ p,2 λ p, m λ p,n 1 λ p,n 2 +1 n 1 0 < m n m n 1... m 1 < 1, m n = 1 m 1 m n 1. This is a distribution of ratios of Pareto type disribution. m λ p,n 1+1 n,
37 Expected capital distribution curves From the expected slopes E ν [ log µ (k) log µ (k 1) log(k+1) log k ] = Eν (Ξ k ) log(1+k 1 ) obtain expected capital distribution curves. e 1 e 1 we iv e 5 e 5 v Weight iii Weight e 10 i ii e 10 Rank Rank n = 5000, gn = c (2n 1), g k = 0, 1 k n 1, γ 1 = c, γ i = 2c, 2 i n, σ 2 k = k 10 5, 1 k n. (i) c = 0.02, (ii) c = 0.03, (iii) c = (iv) c = 0.02, g 1 = 0.016, g k = 0, 2 k n 1, g n = (0.02)(2n 1) , (v) g1 = = g 50 = 0.016, g k = 0, 51 k n 1, g n = (0.02)(2n 1)
38 28 Empirical data Variance Rate Historical capital distribution curves. Data Source: Fernholz ( 02) Rank Growing variances
39 29 Capital Stocks and Portfolio Rules Market X = ((X 1 (t),..., X n (t)), t 0) of n companies log X i(t ) T T n X i (0) = G i (t)dt + S i,ν (t)dw ν (t), 0 0 ν=1 with initial capital X i (0) = x i > 0, i = 1,..., n, 0 T <. Define a ij ( ) = n ν=1 S iν( )S jν ( ) and A ij ( ) = 0 a ij(t)d t. Long only Portfolio rule π and its wealth V π. Choose π n + := {x R n : x i = 1, x i 0 } invest π i V π of money to company i for i = 1,..., n, i.e., π i V π /X i share of company i : n dv π (t) = i=1 V π (0) = w. π i (t)v π (t) dx i (t), 0 t <, X i (t)
40 29 Capital Stocks and Portfolio Rules Market X = ((X 1 (t),..., X n (t)), t 0) of n companies log X i(t ) T T n X i (0) = G i (t)dt + S i,ν (t)dw ν (t), 0 0 ν=1 with initial capital X i (0) = x i > 0, i = 1,..., n, 0 T <. Define a ij ( ) = n ν=1 S iν( )S jν ( ) and A ij ( ) = 0 a ij(t)d t. Long only Portfolio rule π and its wealth V π. Choose π n + := {x R n : x i = 1, x i 0 } invest π i V π of money to company i for i = 1,..., n, i.e., π i V π /X i share of company i : n dv π (t) = i=1 V π (0) = w. π i (t)v π (t) dx i (t), 0 t <, X i (t)
41 29 Capital Stocks and Portfolio Rules Market X = ((X 1 (t),..., X n (t)), t 0) of n companies log X i(t ) T T n X i (0) = G i (t)dt + S i,ν (t)dw ν (t), 0 0 ν=1 with initial capital X i (0) = x i > 0, i = 1,..., n, 0 T <. Define a ij ( ) = n ν=1 S iν( )S jν ( ) and A ij ( ) = 0 a ij(t)d t. Long only Portfolio rule π and its wealth V π. Choose π n + := {x R n : x i = 1, x i 0 } invest π i V π of money to company i for i = 1,..., n, i.e., π i V π /X i share of company i : n dv π (t) = i=1 V π (0) = w. π i (t)v π (t) dx i (t), 0 t <, X i (t)
42 30 Portfolios and Relative Arbitrage Market portfolio: Take π(t) = m(t) = (m 1 (t),, m n (t)) n + where m i (t) = X i (t) X 1 + X n (t), i = 1,..., n, 0 t <. Diversity weighted portfolio: Given p [0, 1], take π i (t) = (m i (t)) P p n j=1 (m for i = 1,..., n, 0 t <. j (t)) p Functionally generated portfolio (Fernholz ( 02) & Karatzas ( 08)). A portfolio π represents an arbitrage opportunity relative to another portfolio ρ on [0, T ], if P(V π (T ) V ρ (T )) = 1, P(V π (T ) > V ρ (T )) > 0. Can we find an arbitrage opportunity π relative to m?
43 30 Portfolios and Relative Arbitrage Market portfolio: Take π(t) = m(t) = (m 1 (t),, m n (t)) n + where m i (t) = X i (t) X 1 + X n (t), i = 1,..., n, 0 t <. Diversity weighted portfolio: Given p [0, 1], take π i (t) = (m i (t)) P p n j=1 (m for i = 1,..., n, 0 t <. j (t)) p Functionally generated portfolio (Fernholz ( 02) & Karatzas ( 08)). A portfolio π represents an arbitrage opportunity relative to another portfolio ρ on [0, T ], if P(V π (T ) V ρ (T )) = 1, P(V π (T ) > V ρ (T )) > 0. Can we find an arbitrage opportunity π relative to m?
44 31 Constant-portfolio For a constant-proportion π( ) π, [ n { V π Aii (t) ( (t) = w exp π i 2 +log Xi (t) )} 1 X i (0) 2 for 0 t <. i=1 n ] π i A ij (t)π j i,j=1 Here A ij ( ) = 0 a ij(t)d t and a ij ( ) = n ν=1 S iν( )S jν ( ), d Y (t) = G(Y (t))d t + S(Y (t))dw (t) ; 0 t <, G(y) = p Σ n (g p 1 (1) + γ 1 + γ,..., g p 1 (n) + γ n + γ) 1 Rp (y), S(y) = p Σ n diag(σ p 1 (1),..., σ p 1 (n) ) }{{} s p 1 Rp (y) ; y Rn.
45 32 Target Portfolio(Cover( 91) & Jamshidian( 92)) [ n { V π Aii (t) ( ( ) = w exp π i 2 +log Xi ( ) )} 1 X i (0) 2 i=1 n ] π i A ij ( )π j i,j=1 Target Portfolio Π (t) maximizes the wealth V π (t) for t 0: V (t) := max V π (t), π n + Π (t) := arg max V π (t), π n + where by Lagrange method we obtain Π i (t) = (2A ii (t) n j=1 1 ) 1 [ 2 n 2 A jj (t) n j= ( A ii (t) log Xi (t) ) ; 0 t <. X i (0) 1 ( A jj (t) log Xj (t) )] X j (0)
46 Asymptotic Target Portfolio Under the hybrid Atlas model with the assumptions 1 v(π) := lim T T log V π (T ) = γ + 1 ( n n π i a ii π i a ) ii π j 2 i=1 i=1 }{{} γπ where (a ij ) 1 i n is the (i,i) element of a 1 := lim T T T 0 (a ij (t)) 1 i,j n d t = p Σ n θ p s p s p. Asymptotic target portfolio maximizes the excess growth γπ : ( n n π := arg max π i a π n ii π i a ) ii π j. + We obtain π i = 1 [ 1 n 2 ( n 2 a ii 1 a j=1 jj i=1 ) 1 ] 33 i=1 = lim Π i (t) ; i = 1,..., n. t
47 Asymptotic Target Portfolio Under the hybrid Atlas model with the assumptions 1 v(π) := lim T T log V π (T ) = γ + 1 ( n n π i a ii π i a ) ii π j 2 i=1 i=1 }{{} γπ where (a ij ) 1 i n is the (i,i) element of a 1 := lim T T T 0 (a ij (t)) 1 i,j n d t = p Σ n θ p s p s p. Asymptotic target portfolio maximizes the excess growth γπ : ( n n π := arg max π i a π n ii π i a ) ii π j. + We obtain π i = 1 [ 1 n 2 ( n 2 a ii 1 a j=1 jj i=1 ) 1 ] 33 i=1 = lim Π i (t) ; i = 1,..., n. t
48 34 Universal Portfolio(Cover( 91) & Jamshidian( 92)) Universal portfolio is defined as π Π i ( ) := n i V π ( )dπ + V π ( )dπ, 1 i n, V bπ ( ) = n + V π ( )dπ n +. dπ n + Proposition Under the hybrid Atlas model with the model assumptions, 1 lim T T log V bπ (T ) V π (T ) = lim 1 T T log V bπ (T ) V (T ) = 0 P a.s.
49 34 Universal Portfolio(Cover( 91) & Jamshidian( 92)) Universal portfolio is defined as π Π i ( ) := n i V π ( )dπ + V π ( )dπ, 1 i n, V bπ ( ) = n + V π ( )dπ n +. dπ n + Proposition Under the hybrid Atlas model with the model assumptions, 1 lim T T log V bπ (T ) V π (T ) = lim 1 T T log V bπ (T ) V (T ) = 0 P a.s.
50 35 Conclusion Ergodic properties of Hybrid Atlas model Diversity weighted portfolio, Target portfolio, Universal portfolio. Further topics: short term arbitrage, generalized portfolio generating function, large market (n ), numéraire portfolio, data implementation. References: 1. arxiv: arxiv: (to appear in Annals of Applied Probability) Tomoyuki Ichiba (UCSB)
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 informationEQUITY MARKET STABILITY
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,
More informationSome 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 informationLOCAL 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 informationA 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 informationIOANNIS 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 informationConvex 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 informationTRANSPORT 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 informationWalsh 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 informationCompeting 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 informationEXPLOSIONS 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 informationALEKSANDAR 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 informationGeometry 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 informationOn 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 informationGeometry 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 informationStochastic 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 informationExact 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 informationVOLATILITY 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 informationRESEARCH 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 informationA 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 informationThe 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 informationOptimal Arbitrage under Constraints
Optimal Arbitrage under Constraints Subhankar Sadhukhan Submitted in partial fulfillment of the Requirements for the degree of Doctor of Philosophy in the Graduate School of Arts and Sciences COLUMBIA
More informationDIVERSITY 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 informationMulti-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form
Multi-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form Yipeng Yang * Under the supervision of Dr. Michael Taksar Department of Mathematics University of Missouri-Columbia Oct
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationThomas 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 informationSensitivity analysis of the expected utility maximization problem with respect to model perturbations
Sensitivity analysis of the expected utility maximization problem with respect to model perturbations Mihai Sîrbu, The University of Texas at Austin based on joint work with Oleksii Mostovyi University
More informationOPTIONAL 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 informationON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER
ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose
More informationGeneralized Hypothesis Testing and Maximizing the Success Probability in Financial Markets
Generalized Hypothesis Testing and Maximizing the Success Probability in Financial Markets Tim Leung 1, Qingshuo Song 2, and Jie Yang 3 1 Columbia University, New York, USA; leung@ieor.columbia.edu 2 City
More informationP i [B k ] = lim. n=1 p(n) ii <. n=1. V i :=
2.7. Recurrence and transience Consider a Markov chain {X n : n N 0 } on state space E with transition matrix P. Definition 2.7.1. A state i E is called recurrent if P i [X n = i for infinitely many n]
More informationSEPARABLE 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 informationLarge 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 informationModel-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 informationStrong Markov property of determinantal processes associated with extended kernels
Strong Markov property of determinantal processes associated with extended kernels Hideki Tanemura Chiba university (Chiba, Japan) (November 22, 2013) Hideki Tanemura (Chiba univ.) () Markov process (November
More informationAsymptotic Coupling of an SPDE, with Applications to Many-Server Queues
Asymptotic Coupling of an SPDE, with Applications to Many-Server Queues Mohammadreza Aghajani joint work with Kavita Ramanan Brown University March 2014 Mohammadreza Aghajanijoint work Asymptotic with
More informationHJB 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 informationOn the dual problem of utility maximization
On the dual problem of utility maximization Yiqing LIN Joint work with L. GU and J. YANG University of Vienna Sept. 2nd 2015 Workshop Advanced methods in financial mathematics Angers 1 Introduction Basic
More informationStrong Markov property of determinatal processes
Strong Markov property of determinatal processes Hideki Tanemura Chiba university (Chiba, Japan) (August 2, 2013) Hideki Tanemura (Chiba univ.) () Markov process (August 2, 2013) 1 / 27 Introduction The
More informationOptional decomposition for continuous semimartingales under arbitrary filtrations *
Electron. Commun. Probab. 2 (215), no. 59, 1 1. DOI: 1.1214/ECP.v2-49 ISSN: 183-589X ELECTRONIC COMMUNICATIONS in PROBABILITY Optional decomposition for continuous semimartingales under arbitrary filtrations
More informationOn 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 informationResearch Article Existence and Uniqueness Theorem for Stochastic Differential Equations with Self-Exciting Switching
Discrete Dynamics in Nature and Society Volume 211, Article ID 549651, 12 pages doi:1.1155/211/549651 Research Article Existence and Uniqueness Theorem for Stochastic Differential Equations with Self-Exciting
More informationReversible Coalescing-Fragmentating Wasserstein Dynamics on the Real Line
Reversible Coalescing-Fragmentating Wasserstein Dynamics on the Real Line Vitalii Konarovskyi Leipzig university Berlin-Leipzig workshop in analysis and stochastics (017) Joint work with Max von Renesse
More informationModel-free portfolio theory and its functional master formula
Model-free portfolio theory and its functional master formula Alexander Schied Leo Speiser Iryna Voloshchenko arxiv:166.3325v3 [q-fin.pm] 8 Dec 216 First version: June 1, 216 This version: December 8,
More informationBranching Brownian motion seen from the tip
Branching Brownian motion seen from the tip J. Berestycki 1 1 Laboratoire de Probabilité et Modèles Aléatoires, UPMC, Paris 09/02/2011 Joint work with Elie Aidekon, Eric Brunet and Zhan Shi J Berestycki
More informationOptimal Stopping Problems and American Options
Optimal Stopping Problems and American Options Nadia Uys A dissertation submitted to the Faculty of Science, University of the Witwatersrand, in fulfilment of the requirements for the degree of Master
More informationOn the quantiles of the Brownian motion and their hitting times.
On the quantiles of the Brownian motion and their hitting times. Angelos Dassios London School of Economics May 23 Abstract The distribution of the α-quantile of a Brownian motion on an interval [, t]
More informationBERNOULLI ACTIONS AND INFINITE ENTROPY
BERNOULLI ACTIONS AND INFINITE ENTROPY DAVID KERR AND HANFENG LI Abstract. We show that, for countable sofic groups, a Bernoulli action with infinite entropy base has infinite entropy with respect to every
More informationSmall-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 informationarxiv: v1 [math.pr] 24 Sep 2018
A short note on Anticipative portfolio optimization B. D Auria a,b,1,, J.-A. Salmerón a,1 a Dpto. Estadística, Universidad Carlos III de Madrid. Avda. de la Universidad 3, 8911, Leganés (Madrid Spain b
More informationLecture 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 informationSupplement A: Construction and properties of a reected diusion
Supplement A: Construction and properties of a reected diusion Jakub Chorowski 1 Construction Assumption 1. For given constants < d < D let the pair (σ, b) Θ, where Θ : Θ(d, D) {(σ, b) C 1 (, 1]) C 1 (,
More informationOptimal 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 informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationON THE DIFFUSION OPERATOR IN POPULATION GENETICS
J. Appl. Math. & Informatics Vol. 30(2012), No. 3-4, pp. 677-683 Website: http://www.kcam.biz ON THE DIFFUSION OPERATOR IN POPULATION GENETICS WON CHOI Abstract. W.Choi([1]) obtains a complete description
More informationGradient Estimation for Attractor Networks
Gradient Estimation for Attractor Networks Thomas Flynn Department of Computer Science Graduate Center of CUNY July 2017 1 Outline Motivations Deterministic attractor networks Stochastic attractor networks
More informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationSupermodular 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 informationBessel-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 informationRandom Process Lecture 1. Fundamentals of Probability
Random Process Lecture 1. Fundamentals of Probability Husheng Li Min Kao Department of Electrical Engineering and Computer Science University of Tennessee, Knoxville Spring, 2016 1/43 Outline 2/43 1 Syllabus
More informationTHE ASYMPTOTIC ELASTICITY OF UTILITY FUNCTIONS AND OPTIMAL INVESTMENT IN INCOMPLETE MARKETS 1
The Annals of Applied Probability 1999, Vol. 9, No. 3, 94 95 THE ASYMPTOTIC ELASTICITY OF UTILITY FUNCTIONS AND OPTIMAL INVESTMENT IN INCOMPLETE MARKETS 1 By D. Kramkov 2 and W. Schachermayer Steklov Mathematical
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
More informationn! (k 1)!(n k)! = F (X) U(0, 1). (x, y) = n(n 1) ( F (y) F (x) ) n 2
Order statistics Ex. 4. (*. Let independent variables X,..., X n have U(0, distribution. Show that for every x (0,, we have P ( X ( < x and P ( X (n > x as n. Ex. 4.2 (**. By using induction or otherwise,
More informationPractical conditions on Markov chains for weak convergence of tail empirical processes
Practical conditions on Markov chains for weak convergence of tail empirical processes Olivier Wintenberger University of Copenhagen and Paris VI Joint work with Rafa l Kulik and Philippe Soulier Toronto,
More informationUnderstanding Regressions with Observations Collected at High Frequency over Long Span
Understanding Regressions with Observations Collected at High Frequency over Long Span Yoosoon Chang Department of Economics, Indiana University Joon Y. Park Department of Economics, Indiana University
More informationarxiv: 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 informationOptimal transportation and optimal control in a finite horizon framework
Optimal transportation and optimal control in a finite horizon framework Guillaume Carlier and Aimé Lachapelle Université Paris-Dauphine, CEREMADE July 2008 1 MOTIVATIONS - A commitment problem (1) Micro
More informationOn the submartingale / supermartingale property of diffusions in natural scale
On the submartingale / supermartingale property of diffusions in natural scale Alexander Gushchin Mikhail Urusov Mihail Zervos November 13, 214 Abstract Kotani 5 has characterised the martingale property
More informationUnspanned Stochastic Volatility in the Multi-Factor CIR Model
Unspanned Stochastic Volatility in the Multi-Factor CIR Model Damir Filipović Martin Larsson Francesco Statti April 13, 2018 Abstract Empirical evidence suggests that fixed income markets exhibit unspanned
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationMaximum pressure policies for stochastic processing networks
Maximum pressure policies for stochastic processing networks Jim Dai Joint work with Wuqin Lin at Northwestern Univ. The 2011 Lunteren Conference Jim Dai (Georgia Tech) MPPs January 18, 2011 1 / 55 Outline
More informationPoisson 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 informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationFE 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 informationOn 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 informationShort-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 informationTwo Brownian Particles with Rank-Based Characteristics and Skew-Elastic Collisions
Two Brownian Particles with ank-based Characteristics and Skew-Elastic Collisions arxiv:16.435v [math.p] 3 Aug 1 E. OBET FENHOLZ TOMOYUKI ICHIBA IOANNIS KAATZAS August 19, 1 Abstract We construct a two-dimensional
More informationIntroduction Optimality and Asset Pricing
Introduction Optimality and Asset Pricing Andrea Buraschi Imperial College Business School October 2010 The Euler Equation Take an economy where price is given with respect to the numéraire, which is our
More informationBrownian 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 informationModelling of Dependent Credit Rating Transitions
ling of (Joint work with Uwe Schmock) Financial and Actuarial Mathematics Vienna University of Technology Wien, 15.07.2010 Introduction Motivation: Volcano on Iceland erupted and caused that most of the
More informationQUANTITATIVE FINANCE RESEARCH CENTRE
QUANTITATIVE FINANCE RESEARCH CENTRE QUANTITATIVE FINANCE RESEARCH CENTRE Research Paper 241 December 2008 Viability of Markets with an Infinite Number of Assets Constantinos Kardaras ISSN 1441-8010 VIABILITY
More informationIntroduction to self-similar growth-fragmentations
Introduction to self-similar growth-fragmentations Quan Shi CIMAT, 11-15 December, 2017 Quan Shi Growth-Fragmentations CIMAT, 11-15 December, 2017 1 / 34 Literature Jean Bertoin, Compensated fragmentation
More informationConvex Geometry. Carsten Schütt
Convex Geometry Carsten Schütt November 25, 2006 2 Contents 0.1 Convex sets... 4 0.2 Separation.... 9 0.3 Extreme points..... 15 0.4 Blaschke selection principle... 18 0.5 Polytopes and polyhedra.... 23
More informationTrading 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 informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationConsistency of the maximum likelihood estimator for general hidden Markov models
Consistency of the maximum likelihood estimator for general hidden Markov models Jimmy Olsson Centre for Mathematical Sciences Lund University Nordstat 2012 Umeå, Sweden Collaborators Hidden Markov models
More informationA variational formula for risk-sensitive control
A variational formula for risk-sensitive control Vivek S. Borkar, IIT Bombay May 9, 2018, IMA, Minneapolis PART I: Risk-sensitive reward for Markov chains (joint work with V. Anantharam, Uni. of California,
More informationPath Decomposition of Markov Processes. Götz Kersting. University of Frankfurt/Main
Path Decomposition of Markov Processes Götz Kersting University of Frankfurt/Main joint work with Kaya Memisoglu, Jim Pitman 1 A Brownian path with positive drift 50 40 30 20 10 0 0 200 400 600 800 1000-10
More informationA relative entropy characterization of the growth rate of reward in risk-sensitive control
1 / 47 A relative entropy characterization of the growth rate of reward in risk-sensitive control Venkat Anantharam EECS Department, University of California, Berkeley (joint work with Vivek Borkar, IIT
More informationThe Azéma-Yor Embedding in Non-Singular Diffusions
Stochastic Process. Appl. Vol. 96, No. 2, 2001, 305-312 Research Report No. 406, 1999, Dept. Theoret. Statist. Aarhus The Azéma-Yor Embedding in Non-Singular Diffusions J. L. Pedersen and G. Peskir Let
More informationDirichlet uniformly well-approximated numbers
Dirichlet uniformly well-approximated numbers Lingmin LIAO (joint work with Dong Han Kim) Université Paris-Est Créteil (University Paris 12) Hyperbolicity and Dimension CIRM, Luminy December 3rd 2013 Lingmin
More informationSpectral properties of Markov operators in Markov chain Monte Carlo
Spectral properties of Markov operators in Markov chain Monte Carlo Qian Qin Advisor: James P. Hobert October 2017 1 Introduction Markov chain Monte Carlo (MCMC) is an indispensable tool in Bayesian statistics.
More informationHeavy Tailed Time Series with Extremal Independence
Heavy Tailed Time Series with Extremal Independence Rafa l Kulik and Philippe Soulier Conference in honour of Prof. Herold Dehling Bochum January 16, 2015 Rafa l Kulik and Philippe Soulier Regular variation
More informationLecture 2. We now introduce some fundamental tools in martingale theory, which are useful in controlling the fluctuation of martingales.
Lecture 2 1 Martingales We now introduce some fundamental tools in martingale theory, which are useful in controlling the fluctuation of martingales. 1.1 Doob s inequality We have the following maximal
More informationSMSTC (2007/08) Probability.
SMSTC (27/8) Probability www.smstc.ac.uk Contents 12 Markov chains in continuous time 12 1 12.1 Markov property and the Kolmogorov equations.................... 12 2 12.1.1 Finite state space.................................
More informationNumerical 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 informationRuin probabilities in multivariate risk models with periodic common shock
Ruin probabilities in multivariate risk models with periodic common shock January 31, 2014 3 rd Workshop on Insurance Mathematics Universite Laval, Quebec, January 31, 2014 Ruin probabilities in multivariate
More informationMarch 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 informationMean-Variance Hedging for Continuous Processes: New Proofs and Examples
Mean-Variance Hedging for Continuous Processes: New Proofs and Examples Huyên Pham Équipe d Analyse et de Mathématiques Appliquées Université de Marne-la-Vallée 2, rue de la Butte Verte F 9366 Noisy-le-Grand
More information