arxiv: v2 [q-fin.tr] 11 Apr 2013
|
|
- Hubert Harry Dorsey
- 6 years ago
- Views:
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 - Sylvain DELARE Université Paris Diderot, LPMA) - Christian ROBER Université Lyon 1, Laboratoire SAF) - Mathieu ROSENBAUM
More informationConvergence 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 informationStatistical 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 informationStochastic 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 informationFunctional 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 informationBrownian 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 informationNonparametric 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 informationPoint 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 informationOutline. 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 informationA 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 informationSample 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 informationIntroduction 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 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 informationLecture 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 informationRegular 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 informationPROBABILITY: 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 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 informationDuration-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 informationMASSACHUSETTS 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 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 informationAdditive 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 informationAn 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 informationLiquidation 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 informationStochastic 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 informationSTAT 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 informationMonte-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 informationPROBABILITY: 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 informationOn 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 informationExercises 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 informationStochastic 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 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 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 informationLimit 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 informationMultivariate 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 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 informationSelected 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 informationLecture 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 informationStatistical 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 informationNon-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 informationFunctional 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 informationOptimal 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 informationLecture 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 information1. 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 informationWeak 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 informationTopics 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 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 informationProofs 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 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 information1. 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 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 informationGaussian 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 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 informationWEAK 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.
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 informationStrong 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 informationP (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 informationThe 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 informationMarkov 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 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 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 informationOn 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 informationA 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 informationLIMITS 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 informationA 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 informationSubsampling 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 informationSTAT 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 informationMinimization 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 informationFiltrations, 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 informationarxiv: 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 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 informationis 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 informationSUPPLEMENT 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 informationMATH 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 informationOrder 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 informationNotes 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 informationA 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 informationGENERALIZED 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 informationCross-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 informationLeast 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 informationMonte 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 information1 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 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 informationEstimation 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 informationLimit 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 informationExponential 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 informationContagious 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 informationModeling 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 informationMinimal 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 informationBivariate 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 informationSample 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 informationRecent 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 informationOther 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 informationON 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 informationUNCERTAINTY 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 information9 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 informationSome 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 informationParameter 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 informationBranching 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 informationA 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 informationSPATIAL 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