arxiv: v2 [q-fin.tr] 11 Apr 2013

Size: px
Start display at page:

Download "arxiv: v2 [q-fin.tr] 11 Apr 2013"

Transcription

1 Estimating the efficient price from the order flow: a Brownian Cox process approach arxiv: v2 [q-fin.r] 11 Apr 213 Sylvain Delattre LPMA, Université Paris Diderot Paris 7) delattre@math.univ-paris-diderot.fr Christian Y. Robert SAF, Université Lyon 1 christian.robert@univ-lyon1.fr Mathieu Rosenbaum LPMA, Université Pierre et Marie Curie Paris 6) mathieu.rosenbaum@upmc.fr April 15, 213 Abstract At the ultra high frequency level, the notion of price of an asset is very ambiguous. Indeed, many different prices can be defined last traded price, best bid price, mid price,... ). hus, in practice, market participants face the problem of choosing a price when implementing their strategies. In this work, we propose a notion of efficient price which seems relevant in practice. Furthermore, we provide a statistical methodology enabling to estimate this price form the order flow. Key words: Efficient price, order flow, response function, market microstructure, Cox processes, fractional part of Brownian motion, non parametric estimation, functional limit theorems. 1 Introduction 1.1 What is the high frequency price? he classical approach of mathematical finance is to consider that the prices of basic products future, stock,... ) are observed on the market. In particular, their values are used in order to price complex derivatives. Since options traders typically rebalance their portfolio once or a few times a day, such derivatives pricing problems typically occur at the daily scale. When working at the ultra high frequency scale, even pricing a basic product, that is assigning a price to it, becomes a challenging issue. Indeed, one has access to trades 1

2 and quotes in the order book so that at a given time, many different notions of price can be defined for the same asset: last traded price, best bid price, best ask price, mid price, volume weighted average price,... his multiplicity of prices is problematic for many market participants. For example, market making strategies or brokers optimal execution algorithms often require single prices of plain assets as inputs. Choosing one definition or another for the price can sometimes lead to very significantly different outcomes for the strategies. his is for example the case when the tick value the minimum price increment allowed on the market) is rather large. Indeed, this implies that the prices mentioned above differ in a non negligible way. In practice, high frequency market participants are not looking for the fair economic value of the asset. What they need is rather a price whose value at some given time summarizes in a suitable way the opinions of market participants at this time. his price is called efficient price. Hence, this paper aims at providing a statistical procedure in order to estimate this efficient price. 1.2 Ideas for the estimation strategy In this paper, we focus on the case of large tick assets. We define them as assets for which the bid-ask spread is almost always equal to one tick. Our goal is then to infer an efficient price for this type of asset. Naturally, it is reasonable to assume that the efficient price essentially lies inside the bid-ask spread but we wish to say more. In order to retrieve the efficient price, the classical approach is to consider the imbalance of the order book, that is the difference between the available volumes at the best bid and best ask levels, see for example [4]. Indeed, it is often said by market participants that the price is where the volume is not. Here we consider a dynamic version of this idea through the information available in the order flow. More precisely, we assume that the intensity of arrival of the limit order flow at the best bid or the best ask level depends on the distance between the efficient price and the considered level: if this distance is large, the intensity should be high and conversely. hus, we assume the intensity can be written as an increasing deterministic function of this distance. his function is called the order flow response function. In our approach, a crucial step is to estimate the response function in a non parametric way. hen, this functional estimator is used in order to retrieve the efficient price. Note that it is also possible to use the buy or sell market order flow. In that case, the intensity of the flow should be high when the distance is small. Indeed, in this situation, market takers are not loosing too much money with respect to the efficient price) when crossing the spread. 1.3 Organization of the paper he paper is organized as follows. he model and the assumptions are described in Section 2. Particular properties of the efficient price are given in Section 3 and the main statistical procedure is explained in Section 4. he theorems about the response function can be found in Section 5 and the limiting behavior of the estimator of the efficient price is given 2

3 in Section 6. One numerical illustration can be found in Section 7 and a conclusion is given in Section 8. Finally the proofs are relegated to Section 9. 2 he model 2.1 Description of the model We assume the tick size is equal to one, meaning that the asset can only take integer values. Moreover, the efficient price is given by P t = P + σw t, t [, ], where W t ) t is a Brownian motion on some filtered probability space Ω, F, F t ) t, P) and P is a F -measurable random variable, independent of W t ) t and uniformly distributed on [p, p + 1), with p N. Note that such a simple dynamics for the efficient price is probably still reasonable at our high frequency scale. Let Y t be the fractional part of P t, denoted by {P t }, that is: Y t = {P t } = P t P t. o fix ideas and without loss of generality, we focus in the rest of the paper on the limit order flow at the best bid level. We assume that when a limit order is posted at time t at the best bid level, its price B t is given by P t. herefore at time t, the efficient price is P t = B t + Y t. We denote by N t the total number of limit orders posted over [, t]. In order to translate the fact that the intensity of N t should be an increasing function of Y t, we assume that N t ) t is a Cox process also known as doubly stochastic Poisson process) with arrival intensity at time t given by µhy t ), where µ is a positive constant and where the response function h : [, 1) R + satisfies the following assumption: Assumption H1: Response function he function h is C 1 on [, 1), with derivative h for which there exist some positive constants c 1 and c 2 such that for any x in [, 1), c 1 h x) c 2. Moreover, we have the following identifiability condition: 1 hx)dx = 1. Remark that H1 implies in particular that h is bounded on [, 1). he response function expresses the influence of the efficient price on the order flow. In particular, the limiting case 3

4 where h is constant corresponds to orders arriving according to a standard Poisson process. his simple arrival mechanism is often assumed in practice for technical convenience, see for example [4], [5]. However, it is not very realistic. Nevertheless, it is clear that our model is still probably too simple, Y t being the only source of fluctuation of the intensity. For example, a more reasonable version should probably allows for seasonalities in the flow and the possible influence of exogenous variables. We assume that we observe the point process N t ) on [, ]. o get asymptotic properties, we let tend to infinity. It will be also necessary to assume that µ = µ depends on. More precisely, we have the following assumption: Assumption H2: Asymptotic setting For some ε >, as tends to infinity, 5/2+ε /µ. Note that up to scaling modifications, we could also consider a setting where is fixed and the tick size is not constant equal to 1 but tends to zero. 3 Properties of the process Y t he intensity of the order flow is µ hy t ). his process Y t ) t has several nice properties. 3.1 Markov property First, recall that if U is uniformly distributed on [, 1] and X is a real-valued random variable, which is independent of U then {U + X} is also uniformly distributed on [, 1]. hus, since P is uniformly distributed on [, 1], we obtain that Y t ) is a stationary Markov process. hen, from the properties of the sample paths of the Brownian motion, it is clear that Y t ) is Harris recurrent and therefore, for f a bounded positive measurable function, almost surely, 1 1 lim fy s )ds = fs)ds, + see for example heorem 3.12 in [7]. 3.2 Regenerative property Beyond the Markov property, the process Y t ) also enjoys a regenerative property. We define the sequence of stopping time ν n ) n N the following way: ν =, ν 1 = inf {t > : P t N} and for n 2: ν n = inf{t > ν n 1 : P t = P νn 1 ± 1} = inf{t > ν n 1 : W t = W νn 1 ± 1/σ}. 1) he cycles Y t+νn ) t<νn+1 ν n are independent and identically distributed for n 1. hus Y t ) is a regenerative process. Note that ν 2 ν 1 has the same law as τ 1 = inf{t, W t = 1/σ}. 4

5 he Laplace transform of τ 1 is given by E[e γτ 1 ] = 1/ cosh 2γ)/σ ), γ. From this expression, we get E[τ 1 ] = 1/σ 2 and E[τ1 2] = 5/3σ4 ). hus for f a bounded positive measurable function, by heorem VI.3.1 in [1], we get that almost surely 1 1 ν2 lim fy s )ds = + E[ν 2 ν 1 ] E[ fy t )dt ] = σ 2 E τ1 f{σw t })dt ). ν 1 In particular, this implies that σ 2 E [ τ 1 f{σw t })dt ] = 1 fs)ds. Furthermore, by heorem VI.3.2, we get the following convergence in distribution as : 1 with Z f a centered random variable defined by Z f = fy t )dt σ 2 E [ τ 1 f{σw t })dt ]) d N, σ 2 Var[Z f ]), 2) τ1 f{σw t })dt τ 1 σ 2 E [ τ 1 f{σw t })dt ]. We end this section with two particular cases for the function f which will be useful in the following. First, if f = h, using the identifiability condition in Assumption H1, we get that almost surely: herefore, lim + 1 hy t )dt = σ 2 E [ τ 1 h{σw t })dt ] = 1. Z h = = τ1 1/σ 1/σ I{Wt<}h1 + σw t ) 1) + I {Wt>}hσW t ) 1) ) dt I{u<} h1 + σu) 1) + I {u>} hσu) 1))L 1/σ,1/σ u)du, where L 1/σ,1/σ is the local time of W t ) t stopped at τ 1. In the same way, for fx) = I {hx) t}, t [, h1 )), we have So we define Z e t) = Z f by 1 I {hs) t} ds = h 1 t). Z e t) = = τ1 1/σ 1/σ I{Ws<,h1+σW s) t} + I {Ws>,hσW s) t} h 1 t) ) ds I{u<,h1+σu) t} + I {u>,hσu) t} h 1 t) ) L 1/σ,1/σ u)du. 5

6 In fact, it will be useful to see Z e t)) t [,h1 )) as a process indexed by t. hus, we introduce the covariance function ρt 1, t 2 ) = Cov[Z e t 1 ), Z e t 2 )] = E[Z e t 1 )Z e t 2 )]. 3) Note that this function and Var[Z h ] can be computed as multiple integrals on [ 1/σ, 1/σ] 2. Indeed, for u, v) [ 1/σ, 1/σ] 2 an explicit expression for the bivariate Laplace transform of L 1/σ,1/σ u), L 1/σ,1/σ v) ) is available, see Formula I ) in [3]. Using differentiation, from this expression, one can easily derive 1 E[L 1/σ,1/σ u)l 1/σ,1/σ v)] which enables to deduce ρt 1, t 2 ) and Var[Z h ]. 4 Estimating the flow response function and its inverse We want to estimate the flow response function. Before estimating h, we need to estimate µ. Since E [ N ] [ 1 = E hy t )dt ] = 1 E[hY t )]dt = 1, µ it is natural to propose ˆµ = N 4) as an estimator for µ. We have the following proposition, whose proof is given in Section 9.1. Proposition 1 If µ as, we have ˆµ µ 1 ) d N, σ 2 Var[Z h ]). We now explain how the estimator of the flow response function is built. Let k be a known deterministic sequence of positive integers. hen define for j = 1,..., k ˆθ j = k N j/k N j 1)/k ˆµ = k N N j/k N j 1)/k ). Now remark that ˆθ j is approximately equal to 1 µ /k µ /k i=1 Nj 1)/k +i/µ N j 1)/k +i 1)/µ ). Conditional on the path of Y t ), the variables in the sum are independent and if /k is small enough, they approximately follow a Poisson law with parameter µ hy j 1)/k )/k. herefore, if moreover µ /k is sufficiently large, one can expect that ˆθ j is close to hy j 1)/k ). 1 We skip the explicit writing of E[L 1/σ,1/σ u)l 1/σ,1/σ v)], the expression being quite heavy. 6

7 In fact, throughout the paper, we assume that k is chosen so that for some p >, as tends to infinity, p+1/2 k 1/2,. µ k p/2 Note that since 5/2+ε /µ such a sequence k exists. he ˆθ j introduced above are k estimators of quantities of the form hu j ), u j [, 1]. However, we do not have access to the values of the u i. Nevertheless, we know that they are uniformly distributed on [, 1]. We therefore rank the ˆθ j : ˆθ 1) ˆθ 2)... ˆθ k ). For u [, 1), we define the estimator of hu) the following way: ĥu) = ˆθ uk +1). Note that the identifiability condition holds since 1 ĥu)du = 1 k k hen, the estimator of h 1 is naturally defined by: ĥ 1 t) = 1 k k ˆθ j = 1. Indeed, ĥ is the right continuous generalized inverse of ĥ 1. I {ˆθj t}. 5) 5 Limit theorems for the response function We now give the theorems associated to ĥ 1 and ĥ. For b >, we write D[, b) resp. D[, b]), for the space of càdlàg functions from [, b) resp. [, b]) into R. We define the convergence in law in D[, b) as the convergence in law of the restrictions of the stochastic processes to any compact [, a], < a < b, for the Skorohod topology J 1. We have the following result for ĥ 1 : heorem 1 Under H1 and H2, as tends to infinity, we have ĥ 1 ) h 1 ) ) d σg ) ) h h 1 ) ) h1 ) σgv)dv, in D[, h1 )), where G ) is a continuous centered Gaussian process with covariance function ρ defined by 3). Note that although there are k terms in the sum defining h 1 in 5), the rate of convergence in heorem 1 is which is slower than k. his is due to the strong dependence between the ˆθ j within each cycle and the fact that the number of cycles is of order. he same type of result holds also for ĥ: 7

8 heorem 2 Under H1 and H2, as tends to infinity, we have in D[, 1). ) d ĥ ) h ) h )σg h ) ) h1 ) + h ) σgv)dv, Note that both in heorem 1 and heorem 2, the covariance functions of the limiting processes can be computed explicitly, see Section Estimation of the efficient price Let t, ]. A nice corollary of heorem 1 is the possibility of estimating Y t, which is the fractional part of the efficient price at time t, and therefore the efficient price itself provided t is an order arrival time). Indeed, we have an estimator ĥ 1 of h 1 and we know how to estimate hy t ), namely we take 2 ĥy t ) = k N t N t /k ˆµ herefore, it is natural to define an estimator Y t ) of Y t by Ŷ t = ĥ 1 ĥy t ) ). he following corollary shows that this estimator is indeed a good estimator of the efficient price Y t. Corollary 1 For t >, under H1 and H2, as tends to infinity, we have Ŷt Y t ) d σg hyt ) ),. with G independent of Y t. hus, thanks to this approach it is possible to retrieve the efficient price from the order flow, with the usual parametric rate. Remark that the law of σg hy t ) ) can be explicitly computed. Indeed, conditional on Y t we have a centered Gaussian random variable with variance σ 2 ρ hy t ), hy t ) ). hen we can use the fact that Y t is independent of G and uniformly distributed. 7 One numerical illustration In order to give an idea of the way our procedure behaves, we give in this section one short numerical illustration. We take σ = 1, = 5, µ = 1 and k = 15. We consider 2 Of course we can take any other time interval of length /k containing t. Also, one could probably consider a smoothed version of ĥyt) using a kernel rather than a rectangle window. 8

9 two different cases 3 for the function h: hu) = 2u and hu) = 4u 3. We give in Figure 1 examples of estimator of the order flow response function h, over one sample path in each case. Averages of estimators together with 95% confidence intervals obtained from 1 Monte Carlo simulations are drawn in Figure u u Figure 1: One sample path for hu) = 2u left) and hu) = 4u 3 right). In red the function h and in black its estimation u u Figure 2: 1 Monte Carlo simulations for hu) = 2u left) and hu) = 4u 3 right). In red the function h, in black its estimation, and in green the 95 % confidence intervals. We see that in each case, the function h is indeed fairly well estimated. here is of course a positive bias around zero, especially in the first case. his is in particular due to the fact that our estimators cannot be negative, whereas in both cases h) =. 3 Note that strictly speaking, in the second case, the function h does not fully satisfy our assumptions since the derivative in zero is equal to zero. 9

10 8 Conclusion In this paper, we take advantage of the dynamic information from the order flow in order to estimate an efficient price, whereas the classical approach is to use the static information from the imbalance of the order book. his is done through a simple model which follows the principle that the price is where the volume is not. Indeed we assume that the intensity of the flow at the considered level of the order book best bid or best ask) is a deterministic function order flow response function) of the distance between this level and the efficient price. his function is estimated in a non parametric way and the estimation of the efficient price is derived from it. Remark that in our approach, the arrival of the orders at the best bid and best ask levels are connected through the influence of the efficient price. However, in the estimating procedures, the bid and ask sides are treated separately. An interesting extension of this work, left for further research, would be to find a suitable way to combine these estimators, or even to consider a more intricate dependence structure between the order flows at the best bid and best ask levels. 9 Proofs All the proofs are done under H1 and H2. In the following, c denotes a constant that may vary from line to line, and even in the same line. 9.1 Proof of Proposition 1 We give here an elementary proof of Proposition 1. However, note that this result can also be obtained as a by product of Corollary 2 shown in Section 9.3. he characteristic function of ˆµ /µ 1) satisfies E [ exp it ˆµ /µ 1) )] = exp it [ [ N )] ] )E E Y exp it, µ where E Y denotes the expectation conditional of the path of Y t ). Since conditional on the path of Y t ), N follows a Poisson random variable with parameter the characteristic function is equal to Now µ hy t )dt, exp it [ )E exp e itµ )] ) 1 1)µ hy t )dt. e itµ ) 1 1 itµ ) 1 t 2 µ ) 2. 1

11 Using the boundedness of h and the fact that µ goes to infinity, we get that the limit as tends to infinity of the characteristic function is given by [ lim E exp it ) 1 hy t )dt 1 ))]. + From Equation 2), this is equal to 9.2 An auxiliary result exp σ2 t 2 2 Var[Zh ] ). We now give an auxiliary result from which heorem 1 and heorem 2 will be easily deduced. For j = 1,..., k, we set θ j = and define the function ĥe by for u [, 1). Note that We also define which satisfies k µ N j/k N j 1)/k ) = ˆµ µ ˆθj, ĥ e u) = θ uk +1) = ĥu) ˆµ µ, ˆµ µ = N µ = 1 k We have the following proposition: k ĥ 1 e t) = 1 k θ j = 1 k I {θj t}, ĥ 1 ˆµ ) t) = ĥ 1 e t. µ ĥ e u)du. Proposition 2 As tends to infinity, we have ĥ 1 e ) h 1 ) ) d σg ), in D[, h1 )). By heorem 15.1 in [2], it is enough to prove that the finite dimensional distributions converge and that a tightness criterion holds. he proof is splitted into two Lemmas: Lemma 1 for the convergence of the finite dimensional distributions and Lemma 2 for the C-tightness. Lemma 1 For any n N, for any t 1,..., t n ) [, h1 )) n, as, ĥ 1 e t i ) h 1 t i ) ) d N, σ 2 ν) i=1,...,n where ν i,j = ρ t i, t j ), with ρ the covariance function defined in 3). 11

12 Proof. We give the proof for n = 1. he other cases are obtained the same way by linear combinations of the components of ĥ 1 e t i ) h 1 t i ) ) i=1,...,n. We write ĥ 1 e t) h 1 t) = , with 1 = 1 k 2 = 1 3 = 1 k k I k k { I {θj t} 1 k k I k { j /k j 1) /k hy u)du t} 1 I {hyu) t}du h 1 t). We study each component separately. We have 1 = 1 k j /k j 1) /k hy u)du t}, k j/k k v j,k, j 1)/k I {hyu) t}du, with v j,k, = or - v j,k, = 1 if θ j t < k j/k j 1)/k hy t )dt, - v j,k, = 1 if k j/k j 1)/k hy t )dt t < θ j. For ε >, we have E[ v j,k, ] A 1 + A 2, with A 1 = P [ k j/k A 2 = E [ I {vj,k, =±1}I { t k j 1)/k hy u )du t ε ], j /k j 1) /k hy t)dt >ε } ]. We have k j/k hy t )dt hy j 1)/k ) c k j 1)/k j/k j 1)/k Y t Y j 1)/k dt. Now remark that from Markov inequality and the fact that Y s is uniformly distributed, we easily get that for any p >, δ >, λ > : P [ sup Y t Y s λ ] c p δ p/2 λ p + δ 1/2 ). herefore, for p >, t [s,s+δ] P [ k j/k j 1)/k hy t )dt hy j 1)/k ) ε ] c /k ) p/2 12 ε p + c/k ) 1/2.

13 hus A 1 P [ hy j 1)/k ) t ] /k ) p/2 2ε + c + c/k ) 1/2. Using the stationary distribution of Y j 1)/k we obtain that for ε small enough, A 1 h 1 t+2ε ) h 1 t 2ε )+c /k ) p/2 +c/k ) 1/2 cε +c /k ) p/2 +c/k ) 1/2. We now turn to A 2. Remarking that v j,k, = 1 implies ε p ε p ε p we get θ j k j/k hy t )dt t k j 1)/k A 2 P [ θj k j/k j/k j 1)/k hy t )dt, j 1)/k hy t )dt > ε ]. d Recall that conditional on the path of Y t ), θ j = Zk /µ ), where Z is a Poisson random variable with parameter j/k µ hy t )dt. j 1)/k hus we obtain Using that we finally get k A 2 c µ ε 2. E[ 1 ] 1 k k E[ v j,k, ], E[ p+1)/2 k 1 ] cε + c + c + c. µ ε 2 k ε p kp/2 We take ε = ζ /, with ζ tending to zero. hen E[ 1 ] cζ + c p+1/2 ζ p kp/2 + c k µ ζ 2 + c k. hanks to the assumptions on k, we can find a sequence ζ tending to such that E[ 1 ] goes to. We now turn to 2. We have 2 = 1 = 1 k j/k I{ k j 1)/k k j/k j 1)/k w j,k,s, ds I ) j /k hy {hy j 1) /k u)du t} s) t} ds 13

14 with w j,k,s, = or - w j,k,s, = 1 if k j/k j 1)/k hy u )du t < hy s ) - w j,k,t, = 1 if hy s ) t < k j/k j 1)/k hy u )du. Now remark that w j,k,s, = 1 implies hys ) k j/k hy u )du t k j 1)/k j/k j 1)/k hy u )du. Using the same kind of computations as for 1, for ε positive and small enough and p >, this leads to P [ w j,k,s, = 1 ] cε + c /k ) p/2 ε + c/k ) 1/2. )p hen, in the same way as for 1, we easily obtain that E[ 2 ] goes to. Eventually, from 2), we get that 3 converges in law towards a centered Gaussian random variable with variance σ 2 Var[Z e t)]. We now give the second Lemma needed to prove Proposition 2. Lemma 2 Let α t) = 1 he sequence α t) ) t<h1 ) is tight in D[, h1 )). I {hys) t}ds h 1 t) ). Proof. Recall first the following classical C-tightness criterion, see heorem 15.6 in [2]: If for p >, p 1 > 1 and all t 1, t 2 < h1 ) E [ α t 1 ) α t 2 ) p] c t 1 t 2 p ) 1, then α t) ) t h1) is tight in D[, h1 )). We now use the sequence of stopping time ν n ) n N defined by 1). Let n = inf{i : ν i }. From heorem II.5.1 and heorem II.5.2 in [6], we have E[ν n ] c 6) 14

15 and and We now define Y i t 1, t 2 ) = 1 ν i E[ n σ 2 )2 ] c. 7) ν i 1 I {t1 <hy t) t 2 }dt 1 σ 2 h 1 t 2 ) h 1 t 1 ) )) Ỹ n +1t 1, t 2 ) = 1 νn I {t1 <hy t) t 2 }dt. Remark that for i 2, the Y i t 1, t 2 ) are centered see Section 3.2) and iid. Moreover, using the occupation formula together with a aylor expansion, we get with E [ Y 2 t 1, t 2 ) ) [ 2] νi 1 E I {t1 <hy t) t 2 }dt ) ] 2 c 1 t 2 t 1 2 E[L ) 2 ], ν i 1 L = sup L 1/σ,1/σ u) ). u [ 1/σ,1/σ] From the Ray-Knight version of Burkhölder-Davis-Gundy inequality, see [7], we know that all polynomial moments of L are finite and thus We show in the same way that E [ Y 2 t 1, t 2 ) ) 2] c 1 t 2 t 1 2. E [ Y 1 t 1, t 2 ) ) 2] + E [Ỹn +1t 1, t 2 ) ) 2] c 1 t 2 t 1 2. We now use the preceding inequalities in order to show that the tightness criterion holds. We have α t 1 ) α t 2 ) = B 1 + B 2 B 3 + B 4, with hus By heorem I.5.1 in [6], B 1 = Y 1 t 1, t 2 ), B 2 = n i=2 Y i t 1, t 2 ), B 3 = Ỹn +1t 1, t 2 ), B 4 = n 1 σ 2 1) h 1 t 2 ) h 1 t 1 ) ). E [ E [ α t 1 ) α t 2 ) 2] c n i=2 4 E[ B i 2 ]. i=1 Y i t 1, t 2 ) 2 ] ce[n ]E [ Y 2 t 1, t 2 ) ) 2]. 15

16 Moreover, using 6), we obtain E[n ]E [ Y 2 t 1, t 2 ) 2] c t 1 t 2 2. Finally, we easily get from 7) that E[B 4 ) 2 ] c t 1 t 2 2. and therefore the tightness criterion is satisfied. 9.3 Proofs of heorem 1 and heorem 2 In order to prove heorem 1 and heorem 2, we start with the following corollary of Proposition 2: Corollary 2 As, ĥ 1 e ) h 1 ), ĥe ) h ), converges in law towards σg ), h )σg h ) ), 1 h1 ) in D[, h1 )) D[, 1) R for the product topology). ) ĥ e u)du 1 ) σgv)dv, Proof. he result follows directly from Proposition 2 together with heorem in [8] and the continuous mapping theorem. with We now give the proof of heorem 1. We write ĥ 1 ) h 1 ) ) = V 1 + V 2, V 1 = ĥ 1 e ˆµ /µ ) h 1 ˆµ /µ ) ), V 2 = h 1 ˆµ /µ ) h 1 ) ). From Proposition 2 together with the continuity of the composition map, see heorem in [8], we have d V1 σg ), in D[, h1 )). Now recall that 1 ĥ e u)du = ˆµ /µ. hus using Corollary 2 together with heorem in [8], we get V2 d ) h h 1 ) ) h1 ) σgv)dv, in D[, h1 )). he convergences of V 1 and V 2 taking place jointly, this ends the proof of heorem 1. From heorem 1, heorem 2 follows from heorem in [8]. 16

17 9.4 Proof of Corollary 1 Let ĥey t ) be the oracle estimate of hy t ) defined the same way as ĥy t) but with µ instead of ˆµ. Now we use that ĥ 1 µ ) = ĥ 1 e ) ˆµ and the fact that the cycles on which ĥey t ) is defined are negligible in the estimation of h 1 ) and µ. hen, from Proposition 2, we have ĥ 1 µ ) h 1 ) ) ) d, ˆµ ĥey t ) σg ), hyt ) ) in D[, h1 )) R +, with G independent of Y t. Using heorem in [8] and the fact that ĥy t ) = µ ˆµ ĥ e Y t ), we derive ĥ 1 ĥy t ) ) h 1 ĥ e Y t ) )) d σg hyt ) ). Now recall that h 1 ĥ e Y t ) ) h 1 hy t ) ) c ĥ e Y t ) hy t ). and that conditional on the path of Y t ), ĥey t ) d = Zk /µ ), where Z is a Poisson random variable with parameter µ t t /k hy u )du. herefore, using a aylor expansion and the fact that P [ sup Y u Y t 1 ] c/k ) 1/2, u [t /k,t] we get that [ ĥe E Y t ) hy t ) ] is smaller than ce [ sup Y u Y t ] + c k / µ ) ) 1/2 c/k ) 1/2 + c k / µ ) ) 1/2. u [t /k,t] Eventually, since / k and k /µ tend to zero, we get E [ ĥe Y t ) hy t ) ] c/ k + ck /µ ) 1/2, which concludes the proof. Acknowledgements We are grateful to Renaud Drappier and Sebouh akvorian from BNP Paribas for their very interesting comments. We also thank the referee for his helpful remarks. 17

18 References [1] Asmussen, S. 23). Applied probability and queues. 2nd ed. Springer. [2] Billingsley, P. 1968). Convergence of probability measures. Wiley, New York. [3] Borodin, A. and Salminen, P. 23). Handbook of Brownian motion: facts and formulae. Birkhauser. [4] Cont, R. and de Larrard, A. 213). Price dynamics in a Markovian limit order book market. SIAM J. Financial Math. 4, [5] Farmer, J.D., Patelli, P. and Zovko, I.I. 25). he predictive power of zero intelligence in financial markets. PNAS 12, [6] Gut, A. 1988). Stopped random walks. Limit theorems and applications. Springer- Verlag, New York. [7] Revuz, D. and Yor, M. 25). Continuous martingale and Brownian motion. Springer, third edition. [8] Whitt, W. 22). Stochastic-Process Limits. Springer, New York. 18

Estimating the efficient price from the order flow : a Brownian Cox process approach

Estimating the efficient price from the order flow : a Brownian Cox process approach Estimating the efficient price from the order flow : a Brownian Cox process approach - Sylvain DELARE Université Paris Diderot, LPMA) - Christian ROBER Université Lyon 1, Laboratoire SAF) - Mathieu ROSENBAUM

More information

Convergence at first and second order of some approximations of stochastic integrals

Convergence at first and second order of some approximations of stochastic integrals Convergence at first and second order of some approximations of stochastic integrals Bérard Bergery Blandine, Vallois Pierre IECN, Nancy-Université, CNRS, INRIA, Boulevard des Aiguillettes B.P. 239 F-5456

More information

Statistical Inference for Ergodic Point Processes and Applications to Limit Order Book

Statistical Inference for Ergodic Point Processes and Applications to Limit Order Book Statistical Inference for Ergodic Point Processes and Applications to Limit Order Book Simon Clinet PhD student under the supervision of Pr. Nakahiro Yoshida Graduate School of Mathematical Sciences, University

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

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico

More information

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory

More information

Nonparametric regression with martingale increment errors

Nonparametric regression with martingale increment errors S. Gaïffas (LSTA - Paris 6) joint work with S. Delattre (LPMA - Paris 7) work in progress Motivations Some facts: Theoretical study of statistical algorithms requires stationary and ergodicity. Concentration

More information

Point Process Control

Point Process Control Point Process Control The following note is based on Chapters I, II and VII in Brémaud s book Point Processes and Queues (1981). 1 Basic Definitions Consider some probability space (Ω, F, P). A real-valued

More information

Outline. A Central Limit Theorem for Truncating Stochastic Algorithms

Outline. A Central Limit Theorem for Truncating Stochastic Algorithms Outline A Central Limit Theorem for Truncating Stochastic Algorithms Jérôme Lelong http://cermics.enpc.fr/ lelong Tuesday September 5, 6 1 3 4 Jérôme Lelong (CERMICS) Tuesday September 5, 6 1 / 3 Jérôme

More information

A Barrier Version of the Russian Option

A Barrier Version of the Russian Option A Barrier Version of the Russian Option L. A. Shepp, A. N. Shiryaev, A. Sulem Rutgers University; shepp@stat.rutgers.edu Steklov Mathematical Institute; shiryaev@mi.ras.ru INRIA- Rocquencourt; agnes.sulem@inria.fr

More information

Sample path large deviations of a Gaussian process with stationary increments and regularily varying variance

Sample path large deviations of a Gaussian process with stationary increments and regularily varying variance Sample path large deviations of a Gaussian process with stationary increments and regularily varying variance Tommi Sottinen Department of Mathematics P. O. Box 4 FIN-0004 University of Helsinki Finland

More information

Introduction General Framework Toy models Discrete Markov model Data Analysis Conclusion. The Micro-Price. Sasha Stoikov. Cornell University

Introduction General Framework Toy models Discrete Markov model Data Analysis Conclusion. The Micro-Price. Sasha Stoikov. Cornell University The Micro-Price Sasha Stoikov Cornell University Jim Gatheral @ NYU High frequency traders (HFT) HFTs are good: Optimal order splitting Pairs trading / statistical arbitrage Market making / liquidity provision

More information

On the quantiles of the Brownian motion and their hitting times.

On 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 information

Lecture 17 Brownian motion as a Markov process

Lecture 17 Brownian motion as a Markov process Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is

More information

Regular Variation and Extreme Events for Stochastic Processes

Regular Variation and Extreme Events for Stochastic Processes 1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for

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

A MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT

A 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 information

Duration-Based Volatility Estimation

Duration-Based Volatility Estimation A Dual Approach to RV Torben G. Andersen, Northwestern University Dobrislav Dobrev, Federal Reserve Board of Governors Ernst Schaumburg, Northwestern Univeristy CHICAGO-ARGONNE INSTITUTE ON COMPUTATIONAL

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 7 9/25/2013

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 7 9/25/2013 MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 7 9/5/013 The Reflection Principle. The Distribution of the Maximum. Brownian motion with drift Content. 1. Quick intro to stopping times.

More information

Practical conditions on Markov chains for weak convergence of tail empirical processes

Practical 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 information

Additive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535

Additive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535 Additive functionals of infinite-variance moving averages Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535 Departments of Statistics The University of Chicago Chicago, Illinois 60637 June

More information

An essay on the general theory of stochastic processes

An essay on the general theory of stochastic processes Probability Surveys Vol. 3 (26) 345 412 ISSN: 1549-5787 DOI: 1.1214/1549578614 An essay on the general theory of stochastic processes Ashkan Nikeghbali ETHZ Departement Mathematik, Rämistrasse 11, HG G16

More information

Liquidation in Limit Order Books. LOBs with Controlled Intensity

Liquidation in Limit Order Books. LOBs with Controlled Intensity Limit Order Book Model Power-Law Intensity Law Exponential Decay Order Books Extensions Liquidation in Limit Order Books with Controlled Intensity Erhan and Mike Ludkovski University of Michigan and UCSB

More information

Stochastic integral. Introduction. Ito integral. References. Appendices Stochastic Calculus I. Geneviève Gauthier.

Stochastic integral. Introduction. Ito integral. References. Appendices Stochastic Calculus I. Geneviève Gauthier. Ito 8-646-8 Calculus I Geneviève Gauthier HEC Montréal Riemann Ito The Ito The theories of stochastic and stochastic di erential equations have initially been developed by Kiyosi Ito around 194 (one of

More information

STAT 331. Martingale Central Limit Theorem and Related Results

STAT 331. Martingale Central Limit Theorem and Related Results STAT 331 Martingale Central Limit Theorem and Related Results In this unit we discuss a version of the martingale central limit theorem, which states that under certain conditions, a sum of orthogonal

More information

Monte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan

Monte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan Monte-Carlo MMD-MA, Université Paris-Dauphine Xiaolu Tan tan@ceremade.dauphine.fr Septembre 2015 Contents 1 Introduction 1 1.1 The principle.................................. 1 1.2 The error analysis

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

On Lead-Lag Estimation

On Lead-Lag Estimation On Lead-Lag Estimation Mathieu Rosenbaum CMAP-École Polytechnique Joint works with Marc Hoffmann, Christian Y. Robert and Nakahiro Yoshida 12 January 2011 Mathieu Rosenbaum On Lead-Lag Estimation 1 Outline

More information

Exercises in stochastic analysis

Exercises in stochastic analysis Exercises in stochastic analysis Franco Flandoli, Mario Maurelli, Dario Trevisan The exercises with a P are those which have been done totally or partially) in the previous lectures; the exercises with

More information

Stochastic Calculus for Finance II - some Solutions to Chapter VII

Stochastic Calculus for Finance II - some Solutions to Chapter VII Stochastic Calculus for Finance II - some Solutions to Chapter VII Matthias hul Last Update: June 9, 25 Exercise 7 Black-Scholes-Merton Equation for the up-and-out Call) i) We have ii) We first compute

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

On the submartingale / supermartingale property of diffusions in natural scale

On 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 information

Limit theorems for multipower variation in the presence of jumps

Limit theorems for multipower variation in the presence of jumps Limit theorems for multipower variation in the presence of jumps Ole E. Barndorff-Nielsen Department of Mathematical Sciences, University of Aarhus, Ny Munkegade, DK-8 Aarhus C, Denmark oebn@imf.au.dk

More information

Multivariate Normal-Laplace Distribution and Processes

Multivariate Normal-Laplace Distribution and Processes CHAPTER 4 Multivariate Normal-Laplace Distribution and Processes The normal-laplace distribution, which results from the convolution of independent normal and Laplace random variables is introduced by

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

Selected Exercises on Expectations and Some Probability Inequalities

Selected Exercises on Expectations and Some Probability Inequalities Selected Exercises on Expectations and Some Probability Inequalities # If E(X 2 ) = and E X a > 0, then P( X λa) ( λ) 2 a 2 for 0 < λ

More information

Lecture 12. F o s, (1.1) F t := s>t

Lecture 12. F o s, (1.1) F t := s>t Lecture 12 1 Brownian motion: the Markov property Let C := C(0, ), R) be the space of continuous functions mapping from 0, ) to R, in which a Brownian motion (B t ) t 0 almost surely takes its value. Let

More information

Statistical inference for ergodic point processes and application to Limit Order Book

Statistical inference for ergodic point processes and application to Limit Order Book Available online at www.sciencedirect.com ScienceDirect Stochastic Processes and their Applications 127 (217) 18 1839 www.elsevier.com/locate/spa Statistical inference for ergodic point processes and application

More information

Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings. Krzysztof Burdzy University of Washington

Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings. Krzysztof Burdzy University of Washington Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings Krzysztof Burdzy University of Washington 1 Review See B and Kendall (2000) for more details. See also the unpublished

More information

Functional central limit theorem for super α-stable processes

Functional central limit theorem for super α-stable processes 874 Science in China Ser. A Mathematics 24 Vol. 47 No. 6 874 881 Functional central it theorem for super α-stable processes HONG Wenming Department of Mathematics, Beijing Normal University, Beijing 1875,

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

Lecture 21 Representations of Martingales

Lecture 21 Representations of Martingales Lecture 21: Representations of Martingales 1 of 11 Course: Theory of Probability II Term: Spring 215 Instructor: Gordan Zitkovic Lecture 21 Representations of Martingales Right-continuous inverses Let

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

Weak quenched limiting distributions of a one-dimensional random walk in a random environment

Weak quenched limiting distributions of a one-dimensional random walk in a random environment Weak quenched limiting distributions of a one-dimensional random walk in a random environment Jonathon Peterson Cornell University Department of Mathematics Joint work with Gennady Samorodnitsky September

More information

Topics in fractional Brownian motion

Topics in fractional Brownian motion Topics in fractional Brownian motion Esko Valkeila Spring School, Jena 25.3. 2011 We plan to discuss the following items during these lectures: Fractional Brownian motion and its properties. Topics in

More information

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER

ON 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 information

Proofs of the martingale FCLT

Proofs of the martingale FCLT Probability Surveys Vol. 4 (2007) 268 302 ISSN: 1549-5787 DOI: 10.1214/07-PS122 Proofs of the martingale FCLT Ward Whitt Department of Industrial Engineering and Operations Research Columbia University,

More information

arxiv: v1 [math.pr] 24 Sep 2018

arxiv: 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 information

1. Aufgabenblatt zur Vorlesung Probability Theory

1. Aufgabenblatt zur Vorlesung Probability Theory 24.10.17 1. Aufgabenblatt zur Vorlesung By (Ω, A, P ) we always enote the unerlying probability space, unless state otherwise. 1. Let r > 0, an efine f(x) = 1 [0, [ (x) exp( r x), x R. a) Show that p f

More information

Consistency of the maximum likelihood estimator for general hidden Markov models

Consistency 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 information

Gaussian Processes. 1. Basic Notions

Gaussian Processes. 1. Basic Notions Gaussian Processes 1. Basic Notions Let T be a set, and X : {X } T a stochastic process, defined on a suitable probability space (Ω P), that is indexed by T. Definition 1.1. We say that X is a Gaussian

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

WEAK VERSIONS OF STOCHASTIC ADAMS-BASHFORTH AND SEMI-IMPLICIT LEAPFROG SCHEMES FOR SDES. 1. Introduction

WEAK VERSIONS OF STOCHASTIC ADAMS-BASHFORTH AND SEMI-IMPLICIT LEAPFROG SCHEMES FOR SDES. 1. Introduction WEAK VERSIONS OF STOCHASTIC ADAMS-BASHFORTH AND SEMI-IMPLICIT LEAPFROG SCHEMES FOR SDES BRIAN D. EWALD 1 Abstract. We consider the weak analogues of certain strong stochastic numerical schemes considered

More information

(6, 4) Is there arbitrage in this market? If so, find all arbitrages. If not, find all pricing kernels.

(6, 4) Is there arbitrage in this market? If so, find all arbitrages. If not, find all pricing kernels. Advanced Financial Models Example sheet - Michaelmas 208 Michael Tehranchi Problem. Consider a two-asset model with prices given by (P, P 2 ) (3, 9) /4 (4, 6) (6, 8) /4 /2 (6, 4) Is there arbitrage in

More information

Strong approximation for additive functionals of geometrically ergodic Markov chains

Strong approximation for additive functionals of geometrically ergodic Markov chains Strong approximation for additive functionals of geometrically ergodic Markov chains Florence Merlevède Joint work with E. Rio Université Paris-Est-Marne-La-Vallée (UPEM) Cincinnati Symposium on Probability

More information

P (A G) dp G P (A G)

P (A G) dp G P (A G) First homework assignment. Due at 12:15 on 22 September 2016. Homework 1. We roll two dices. X is the result of one of them and Z the sum of the results. Find E [X Z. Homework 2. Let X be a r.v.. Assume

More information

The Asymptotic Theory of Transaction Costs

The Asymptotic Theory of Transaction Costs The Asymptotic Theory of Transaction Costs Lecture Notes by Walter Schachermayer Nachdiplom-Vorlesung, ETH Zürich, WS 15/16 1 Models on Finite Probability Spaces In this section we consider a stock price

More information

Markov Chains. Arnoldo Frigessi Bernd Heidergott November 4, 2015

Markov Chains. Arnoldo Frigessi Bernd Heidergott November 4, 2015 Markov Chains Arnoldo Frigessi Bernd Heidergott November 4, 2015 1 Introduction Markov chains are stochastic models which play an important role in many applications in areas as diverse as biology, finance,

More information

The Azéma-Yor Embedding in Non-Singular Diffusions

The 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 information

Lecture 2. We now introduce some fundamental tools in martingale theory, which are useful in controlling the fluctuation of martingales.

Lecture 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 information

On pathwise stochastic integration

On pathwise stochastic integration On pathwise stochastic integration Rafa l Marcin Lochowski Afican Institute for Mathematical Sciences, Warsaw School of Economics UWC seminar Rafa l Marcin Lochowski (AIMS, WSE) On pathwise stochastic

More information

A Note on the Central Limit Theorem for a Class of Linear Systems 1

A Note on the Central Limit Theorem for a Class of Linear Systems 1 A Note on the Central Limit Theorem for a Class of Linear Systems 1 Contents Yukio Nagahata Department of Mathematics, Graduate School of Engineering Science Osaka University, Toyonaka 560-8531, Japan.

More information

LIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974

LIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974 LIMITS FOR QUEUES AS THE WAITING ROOM GROWS by Daniel P. Heyman Ward Whitt Bell Communications Research AT&T Bell Laboratories Red Bank, NJ 07701 Murray Hill, NJ 07974 May 11, 1988 ABSTRACT We study the

More information

A Central Limit Theorem for Fleming-Viot Particle Systems Application to the Adaptive Multilevel Splitting Algorithm

A Central Limit Theorem for Fleming-Viot Particle Systems Application to the Adaptive Multilevel Splitting Algorithm A Central Limit Theorem for Fleming-Viot Particle Systems Application to the Algorithm F. Cérou 1,2 B. Delyon 2 A. Guyader 3 M. Rousset 1,2 1 Inria Rennes Bretagne Atlantique 2 IRMAR, Université de Rennes

More information

Subsampling Cumulative Covariance Estimator

Subsampling Cumulative Covariance Estimator Subsampling Cumulative Covariance Estimator Taro Kanatani Faculty of Economics, Shiga University 1-1-1 Bamba, Hikone, Shiga 5228522, Japan February 2009 Abstract In this paper subsampling Cumulative Covariance

More information

STAT Sample Problem: General Asymptotic Results

STAT Sample Problem: General Asymptotic Results STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative

More information

Minimization of ruin probabilities by investment under transaction costs

Minimization of ruin probabilities by investment under transaction costs Minimization of ruin probabilities by investment under transaction costs Stefan Thonhauser DSA, HEC, Université de Lausanne 13 th Scientific Day, Bonn, 3.4.214 Outline Introduction Risk models and controls

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

arxiv: v1 [math.pr] 11 Jan 2013

arxiv: v1 [math.pr] 11 Jan 2013 Last-Hitting Times and Williams Decomposition of the Bessel Process of Dimension 3 at its Ultimate Minimum arxiv:131.2527v1 [math.pr] 11 Jan 213 F. Thomas Bruss and Marc Yor Université Libre de Bruxelles

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

is a Borel subset of S Θ for each c R (Bertsekas and Shreve, 1978, Proposition 7.36) This always holds in practical applications.

is a Borel subset of S Θ for each c R (Bertsekas and Shreve, 1978, Proposition 7.36) This always holds in practical applications. Stat 811 Lecture Notes The Wald Consistency Theorem Charles J. Geyer April 9, 01 1 Analyticity Assumptions Let { f θ : θ Θ } be a family of subprobability densities 1 with respect to a measure µ on a measurable

More information

SUPPLEMENT TO PAPER CONVERGENCE OF ADAPTIVE AND INTERACTING MARKOV CHAIN MONTE CARLO ALGORITHMS

SUPPLEMENT TO PAPER CONVERGENCE OF ADAPTIVE AND INTERACTING MARKOV CHAIN MONTE CARLO ALGORITHMS Submitted to the Annals of Statistics SUPPLEMENT TO PAPER CONERGENCE OF ADAPTIE AND INTERACTING MARKO CHAIN MONTE CARLO ALGORITHMS By G Fort,, E Moulines and P Priouret LTCI, CNRS - TELECOM ParisTech,

More information

MATH 6605: SUMMARY LECTURE NOTES

MATH 6605: SUMMARY LECTURE NOTES MATH 6605: SUMMARY LECTURE NOTES These notes summarize the lectures on weak convergence of stochastic processes. If you see any typos, please let me know. 1. Construction of Stochastic rocesses A stochastic

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

Notes on Random Variables, Expectations, Probability Densities, and Martingales

Notes on Random Variables, Expectations, Probability Densities, and Martingales Eco 315.2 Spring 2006 C.Sims Notes on Random Variables, Expectations, Probability Densities, and Martingales Includes Exercise Due Tuesday, April 4. For many or most of you, parts of these notes will be

More information

A slow transient diusion in a drifted stable potential

A slow transient diusion in a drifted stable potential A slow transient diusion in a drifted stable potential Arvind Singh Université Paris VI Abstract We consider a diusion process X in a random potential V of the form V x = S x δx, where δ is a positive

More information

GENERALIZED RAY KNIGHT THEORY AND LIMIT THEOREMS FOR SELF-INTERACTING RANDOM WALKS ON Z 1. Hungarian Academy of Sciences

GENERALIZED RAY KNIGHT THEORY AND LIMIT THEOREMS FOR SELF-INTERACTING RANDOM WALKS ON Z 1. Hungarian Academy of Sciences The Annals of Probability 1996, Vol. 24, No. 3, 1324 1367 GENERALIZED RAY KNIGHT THEORY AND LIMIT THEOREMS FOR SELF-INTERACTING RANDOM WALKS ON Z 1 By Bálint Tóth Hungarian Academy of Sciences We consider

More information

Cross-Selling in a Call Center with a Heterogeneous Customer Population. Technical Appendix

Cross-Selling in a Call Center with a Heterogeneous Customer Population. Technical Appendix Cross-Selling in a Call Center with a Heterogeneous Customer opulation Technical Appendix Itay Gurvich Mor Armony Costis Maglaras This technical appendix is dedicated to the completion of the proof of

More information

Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises

Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,

More information

Monte Carlo Integration I [RC] Chapter 3

Monte Carlo Integration I [RC] Chapter 3 Aula 3. Monte Carlo Integration I 0 Monte Carlo Integration I [RC] Chapter 3 Anatoli Iambartsev IME-USP Aula 3. Monte Carlo Integration I 1 There is no exact definition of the Monte Carlo methods. In the

More information

1 Continuous-time chains, finite state space

1 Continuous-time chains, finite state space Université Paris Diderot 208 Markov chains Exercises 3 Continuous-time chains, finite state space Exercise Consider a continuous-time taking values in {, 2, 3}, with generator 2 2. 2 2 0. Draw the diagramm

More information

Stochastic Differential Equations.

Stochastic 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 information

Estimation of the Bivariate and Marginal Distributions with Censored Data

Estimation of the Bivariate and Marginal Distributions with Censored Data Estimation of the Bivariate and Marginal Distributions with Censored Data Michael Akritas and Ingrid Van Keilegom Penn State University and Eindhoven University of Technology May 22, 2 Abstract Two new

More information

Limit theorems for Hawkes processes

Limit theorems for Hawkes processes Limit theorems for nearly unstable and Mathieu Rosenbaum Séminaire des doctorants du CMAP Vendredi 18 Octobre 213 Limit theorems for 1 Denition and basic properties Cluster representation of Correlation

More information

Exponential martingales: uniform integrability results and applications to point processes

Exponential martingales: uniform integrability results and applications to point processes Exponential martingales: uniform integrability results and applications to point processes Alexander Sokol Department of Mathematical Sciences, University of Copenhagen 26 September, 2012 1 / 39 Agenda

More information

Contagious default: application of methods of Statistical Mechanics in Finance

Contagious default: application of methods of Statistical Mechanics in Finance Contagious default: application of methods of Statistical Mechanics in Finance Wolfgang J. Runggaldier University of Padova, Italy www.math.unipd.it/runggaldier based on joint work with : Paolo Dai Pra,

More information

Modeling Ultra-High-Frequency Multivariate Financial Data by Monte Carlo Simulation Methods

Modeling Ultra-High-Frequency Multivariate Financial Data by Monte Carlo Simulation Methods Outline Modeling Ultra-High-Frequency Multivariate Financial Data by Monte Carlo Simulation Methods Ph.D. Student: Supervisor: Marco Minozzo Dipartimento di Scienze Economiche Università degli Studi di

More information

Minimal Sufficient Conditions for a Primal Optimizer in Nonsmooth Utility Maximization

Minimal Sufficient Conditions for a Primal Optimizer in Nonsmooth Utility Maximization Finance and Stochastics manuscript No. (will be inserted by the editor) Minimal Sufficient Conditions for a Primal Optimizer in Nonsmooth Utility Maximization Nicholas Westray Harry Zheng. Received: date

More information

Bivariate Uniqueness in the Logistic Recursive Distributional Equation

Bivariate Uniqueness in the Logistic Recursive Distributional Equation Bivariate Uniqueness in the Logistic Recursive Distributional Equation Antar Bandyopadhyay Technical Report # 629 University of California Department of Statistics 367 Evans Hall # 3860 Berkeley CA 94720-3860

More information

Sample of Ph.D. Advisory Exam For MathFinance

Sample of Ph.D. Advisory Exam For MathFinance Sample of Ph.D. Advisory Exam For MathFinance Students who wish to enter the Ph.D. program of Mathematics of Finance are required to take the advisory exam. This exam consists of three major parts. The

More information

Recent results in game theoretic mathematical finance

Recent results in game theoretic mathematical finance Recent results in game theoretic mathematical finance Nicolas Perkowski Humboldt Universität zu Berlin May 31st, 2017 Thera Stochastics In Honor of Ioannis Karatzas s 65th Birthday Based on joint work

More information

Other properties of M M 1

Other properties of M M 1 Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected

More information

ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS

ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)

More information

UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE

UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.

More information

9 Brownian Motion: Construction

9 Brownian Motion: Construction 9 Brownian Motion: Construction 9.1 Definition and Heuristics The central limit theorem states that the standard Gaussian distribution arises as the weak limit of the rescaled partial sums S n / p n of

More information

Some SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen

Some SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen Title Author(s) Some SDEs with distributional drift Part I : General calculus Flandoli, Franco; Russo, Francesco; Wolf, Jochen Citation Osaka Journal of Mathematics. 4() P.493-P.54 Issue Date 3-6 Text

More information

Parameter Estimation

Parameter Estimation Parameter Estimation Consider a sample of observations on a random variable Y. his generates random variables: (y 1, y 2,, y ). A random sample is a sample (y 1, y 2,, y ) where the random variables y

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

A one-level limit order book model with memory and variable spread

A one-level limit order book model with memory and variable spread A one-level limit order book model with memory and variable spread Jonathan A. Chávez-Casillas José E. Figueroa-López September 8, 216 Abstract Motivated by Cont and Larrard (213) s seminal Limit Order

More information

SPATIAL STRUCTURE IN LOW DIMENSIONS FOR DIFFUSION LIMITED TWO-PARTICLE REACTIONS

SPATIAL STRUCTURE IN LOW DIMENSIONS FOR DIFFUSION LIMITED TWO-PARTICLE REACTIONS The Annals of Applied Probability 2001, Vol. 11, No. 1, 121 181 SPATIAL STRUCTURE IN LOW DIMENSIONS FOR DIFFUSION LIMITED TWO-PARTICLE REACTIONS By Maury Bramson 1 and Joel L. Lebowitz 2 University of

More information