arxiv: v4 [math.pr] 23 Nov 2015

Size: px
Start display at page:

Download "arxiv: v4 [math.pr] 23 Nov 2015"

Transcription

1 arxiv: v4 math.pr] 3 Nov 15 Volatility estimatio uder oe-sided errors with applicatios to limit order books Markus Bibiger, Moritz Jirak, Markus Reiß Markus Bibiger, Departmet of Ecoomics, Maheim Uiversity, L7, 3-5, Maheim Moritz Jirak, Markus Reiß, Istitut für Mathematik Humboldt-Uiversität zu Berli Uter de Lide Berli, Germay Abstract: For a semi-martigale X t, which forms a stochastic boudary, a rate-optimal estimator for its quadratic variatio X, X t is costructed based o observatios i the viciity of X t. The problem is embedded i a Poisso poit process framework, which reveals a iterestig coectio to the theory of Browia excursio areas. We derive 1/3 as optimal covergece rate i a high-frequecy framework with observatios i mea. We discuss a potetial applicatio for the estimatio of the itegrated squared volatility of a efficiet price process X t from itra-day order book quotes. AMS subject classificatios: Primary 6H3; secodary 6G55. Keywords ad phrases: Browia excursio area, limit order book, itegrated volatility, Feyma Kac, high-frequecy data, Poisso poit process, oparametric miimax rate. 1. Itroductio Cosider observatios Y i above a stochastic boudary X t, t, 1], which is formed by the graph of a cotiuous semi-martigale. The objective is to optimally recover the drivig characteristic X, X t of the boudary X t, give the observatios Y i. A quatificatio of the iformatio cotet i these observatios is o-trivial ad leads to itriguig mathematical questios. Oe motivatio for cosiderig this stochastic boudary problem stems from fiacial applicatios i the cotext of limit order books. From a microecoomic poit of view ask prices will typically lie above the efficiet market price. Here the uderlyig latet efficiet log-price of a stock X t, t, 1], observed over Fiacial support from the Deutsche Forschugsgemeischaft via SFB 649 Ökoomisches Risiko ad FOR 1735 Structural Iferece i Statistics: Adaptatio ad Efficiecy is gratefully ackowledged. The authors are grateful for helpful commets by the referee. 1

2 M. Bibiger, M. Jirak & M. Reiß/ iid Fig 1. Left: Microstructure oise model Y i = X i/ + ε i, i =,..., = 1, with ε i Exp5. Right: Poisso poit process model with itesity λ t,y = 51 y Xt with X t a Itô process. a tradig period like a day, serves as the boudary, whereas ask prices form the observatios Y i. Bid prices ca be hadled symmetrically ad idepedetly, which ca be used to validate the model. Climate physics provides aother example where semi-martigales appear as stochastic boudaries. Cosiderable efforts are devoted to uderstadig the drivig stochastic term for SDEs modelig the log-term temperature evolutio, see for istace 19] ad 8]. Oe key source for historical temperature data is give by aual tree rigs dedrochroology ad dedroclimatology, see e.g. 14], whose relatioship with temperature i a ideal eviromet is kow. For idividual trees oly sizes up to this ideal boudary are observed due to growth obstructios like limited utritio, leadig to deviatios modeled as egative observatio errors. As a prototype model, we cosider the cotiuous Itô semi-martigale X t = X + t a s ds + t σ s dw s, t, 1], 1.1 with possibly stochastic drift ad volatility coefficiets a s ad σ s, defied o a filtered probability space Ω, F, F t, P with a stadard F t -Browia motio W. Its total quadratic variatio X, X 1 = σ s ds is commoly called itegrated squared volatility. Sectio 3 provides a geeralizatio to models with jumps. A atural cotiuous-time embeddig of the boudary problem is i terms of a Poisso poit process PPP. Coditioal o X t, t, 1] we observe a PPP o, 1] R with itesity measure ΛA = R 1 A t, yλ t,y dt dy, where λ t,y = λ1y X t. 1. We deote by T j, Y j the observatios of that poit process, which are homogeeously dispersed above the graph of X t, t, 1], cf. Figure 1. Theoretically

3 M. Bibiger, M. Jirak & M. Reiß/ 3 ad also ituitively, iformatio o the stochastic boudary ca oly be recovered from the lowest observatio poits ad a homogeeous itesity away from the boudary is assumed for coveiece oly. A associated discrete-time regressio-type model, which explais well the differece to classical oise models, is defied by Y i = X t i with observatio times t i + ε i, i =,...,, ε i, ε i iid F λ, 1.3 ad a error distributio fuctio F λ satisfyig F λ x = λx 1 + O1, as x. 1.4 Oe atural parametric specificatio is ε i Expλ, cf. Figure 1. The oise is assumed to be idepedet of the sigal part X. I microstructure oise models for trasactio prices it is usually assumed that Eε i ] = holds, while here X t i determies the boudary of the support measure for Y i. I fact, if the boudary fuctio was piecewise costat, the by stadard PPP properties we would obtai the regressio-type model 1.3 with expoetial oise from the PPP-model 1. by takig local miima o those pieces. Here we show that uder so called high-frequecy asymptotics, the fudametal quatities i both models exhibit the same asymptotic behaviour, see Propositio 3. below. Compare also 9] for the stroger Le Cam equivalece i the case of smoother boudaries. We shall first cocetrate o the more uiversal PPP model which also allows for simpler scalig ad geometric iterpretatio. Local miima m,k of Y j for T j i some small itervals kh, k + 1h, 1] will form the basic quatities to recover the boudary, which by PPP properties leads to the study of Pm,k > x = E ñ exp k+1h kh ô X t + x + dt, x R, where A + = maxa,, ad its associated momets. For the fudametal case X t = σw t, this opes a iterestig coectio to the theory of Browia excursio areas ad also reveals the difficulty of this problem. It is well documeted i the literature, see e.g. 3], that o explicit form of the expectatio i the expressio above is available. Essetially oly double Laplace trasforms ad related quatities are kow, cf. Propositio 3.3 below ad the attached discussio. This makes the recovery of X, X 1 a itricate probabilistic questio. Still, we are able to prove that our estimator attais the rate 1/3. What is more, by iformatio-theoretic argumets we are able to derive a lower boud showig that the 1/3 -rate is ideed miimax optimal. A more direct proof seems out of reach because the Poisso part from the oise itertwies with the Gaussia martigale part i a way which reders the likelihood ad respective Helliger distaces difficult to cotrol, eve asymptotically. So far, the growig fiace literature o limit order books focusses o modelig ad empirical studies. Empirical cotributios as 5], 9] ad 3] have ivestigated price ad volume distributio, iter-evet duratios as well as the

4 M. Bibiger, M. Jirak & M. Reiß/ 4 Fig. Order price levels for Facebook asset NASDAQ from 1: to 1:3 o Jue d 14. Colored areas highlight spreads betwee differet bid ad ask levels from level 1 up to level 5, bid-ask spread is colored i dark red. structure of the order-flow. Probabilistic models proposed for a limit order book iclude poit process models, see 1], 3] ad ], with mutually excitig processes. Other models come from queuig theory, for istace 15], 33] ad 13], or stochastic optimal cotrol theory as 1]. The mai objective of most modelig approaches is to explai how market prices arise from the book. For the fiacial applicatio, this papers adopts a differet course with the focus o estimatig the latet volatility based o observatios from a limit order book. Cotrarily, to the regular microstructure oise model which costitutes the stadard setup for developig volatility estimators from trasactio data, see e.g. 1], 36], 4] ad ], amog may others, our model assumes oe-sided oise. The optimal covergece rate for volatility estimatio i the model with Gaussia or regular cetered oise ad observatios o a equidistat grid is 1/4, see 18]. Recetly, as iformatio from order books become more ad more available, researchers ad practitioers have sparked the discussio to which kid of observed prices estimatio methods should be applied. 16] discuss this poit ad the possibilities of mid-quotes, executed traded prices or micro-prices which are volume-weighted combiatios of bid ad ask order levels. Noe of these observed time series, however, is free from market microstructure corruptios ad the idea of a uderlyig efficiet price remais utouched. Figure visualizes the iformatio about the evolutio of prices provided by a limit order book for oe specific data set. The colored areas highlight differeces betwee the five best bid ad five best ask levels, the dark area i the ceter markig the bid-ask spread betwee best bid ad best ask. The idea is that a efficiet price should lie, at least most of the time, below the best ask ad symmetrically above the best bid ad that its distace to this stochastic frotier is homogeeous. Sim- Data provided by LOBSTER academic data powered by NASDAQ OMX.

5 M. Bibiger, M. Jirak & M. Reiß/ 5 ilar reasoig served as the fudamet of the order book model by 7] as well as for the dyamic tradig model by ]. Sice modelig i sciece, ecoomics ad particularly fiace is always a compromise betwee catchig major features ad too complex descriptive models, robustess to model misspecificatio is a key issue. Therefore we propose a simple modificatio of our estimator such that occasioal violatios of the cotiuous semi-martigale model do ot chage the asymptotic properties of the estimator. We shall show that for geeral violatios, i particular evoked by jumps of the efficiet price ad the volatility our adjusted estimatio method is robust. The remaider of the paper is orgaized as follows. I Sectio we preset a estimatio approach based o local order statistics whose asymptotic properties are explored i Sectio 3, alog with the robustificatio agaist violatios. I Sectio 4 we prove the lower boud for the miimax estimatio rate. A empirical example is performed i Sectio 5 which cocludes with a discussio. Proofs are provided i the Appedix.. Volatility estimatio based o local miima We costruct the itegrated squared volatility estimator i both models 1. ad 1.3. We partitio the uit iterval ito h 1 N equi-spaced bis Tk = kh, k + 1h, k =,..., h 1 1, with bi-widths h. For simplicity suppose that h N. As the bi-width gets smaller h, whereas the umber of observed values o each bi gets large, h. If we thik of a costat sigal locally o a bi observed with oe-sided positive errors, classical parametric estimatio theory motivates to use the bi-wise miimum as a estimator of the local sigal it the forms a sufficiet statistic uder expoetial oise or equivaletly i the PPP model. I the regressio-type model 1.3 with equidistat observatio times t i = i/, we therefore set m,k = mi Y i, I i I k = {kh, kh + 1,..., k + 1h 1}..1 k Equally, i the PPP model 1. the local miima are give by m,k = mi T j T k Y j, T k = kh, k + 1h.. The same symbol m,k is used i both models because the followig costructio oly depeds o the m,k. All results ad proofs will refer to the cocrete model uder cosideratio. Sice Varm,k X t λh holds i both models, the variace is much smaller tha for a estimator based o a local mea. Nevertheless, we may cotiue i the spirit of the pre-averagig paradigm, cf. ], ad iterpret m,k as a proxy for X t o Tk, which i a secod step is iserted i the realized variace expressio h 1 i=1 X kh X k 1h without oise. The use of a locally costat sigal approximatio X t = X kh + O P h 1/ o T is oly admissible, k

6 M. Bibiger, M. Jirak & M. Reiß/ 6 Fig 3. The poits idicate the fuctio Ψσ with K = 31.6 for small left ad moderate right values of σ. The calculatio is based o accurate Mote Carlo simulatios. The lies show close liear fuctios for compariso. however, if h is chose so small that h 1/ = oλh 1, which would result i a sub-optimal procedure. Rate-optimality ca be attaied if we balace the magitude λh 1 of bi-wise miimal errors due to oise with the rage h 1/ of the motio of X o the bi. This gives the order h λ 3, h 1 3 λ 3..3 I the PPP model 1. this atural choice of the bi-width also follows icely by a scalig argumet: Wt = h 1/ W ht defies a stadard Browia motio for t, 1] based o the values of W o, h ]; the correspodigly scaled PPP observatios T j, Ȳj with T j = h 1 T j, Ȳ j = h 1/ Y j have a itesity with desity λ t,y = λh 3/ 1y W t, which becomes idepedet of exactly for h = λ /3. I this balaced setup the law of the statistics m,k depeds o the motio of X as well as the error distributio i a o-trivial way. Still, the atural statistics to assess the quadratic variatio of the boudary process X are the squared differeces m,k m,k 1 betwee cosecutive local miima. I the PPP model ad with the choice h = K 3 λ 3 for some costat K >.4 the law of h 1/ m,k is idepedet of, h ad λ ad for X t = X k 1h + σ t k 1h dw s o Tk 1 T k, we may itroduce Ψ σ = h 1 E ] m,k m,k 1, k = 1,..., h Below we shall derive theoretical properties of Ψ ad i particular we shall see that it is ivertible as soo as K > is chose sufficietly large. Numerically, the fuctio Ψ ca be determied by stadard Mote Carlo simulatios, see Figure 3, ad is thus available. This paves the way for a momet-estimator

7 M. Bibiger, M. Jirak & M. Reiß/ 7 approach. I fact, km,k m,k 1 approximates Ψσ t dt with correspodig summatio ad itegratio itervals. Uder regularity assumptios o t σ t ad by the smoothess of Ψ show below, we have Ψ 1 1 lr / k=l 1r 1 /+1 m,k m,k 1 h 1 r σ lr 1 h,.6 where r 1 h is a coarse grid size with r h 1 N, r 1 N. This gives rise to the followig estimator of itegrated squared volatility IV = σ t dt i the PPP model 1. with bi-width.4: r h 1 IV h,r = Ψ 1 1 lr / k=l 1r 1 /+1 m,k m,k 1 h 1 r h r 1..7 I the regressio-type model 1.3 the correspodig secod momets still deped o ad we write explicitly Ψ σ = h 1 E ] m,k m,k 1, k = 1,..., h We shall see below that Ψ Ψ holds, but a o-asymptotic form of the volatility estimator from regressio-type observatios is give by r h 1 IV h,r = Ψ 1 lr 1 / k=l 1r 1 /+1 m,k m,k 1 h 1 r h r 1..9 For a parametric estimatio of σ t = σ = cost., we employ the global momet-type estimator IV h,h. Here, iversio of the etire sum of squared differeces is coducted. I the oparametric case of varyig σ t istead a localized estimator IV h,r, with r, r 1 h, is applied. A balace betwee a secod order term o each coarse iterval of order r ad a approximatio error cotrolled by a semi-martigale assumptio o σ t of order r 1 to the choice r h 1/ λ 1/3. h will lead 3. The covergece rate of the estimator I order to cetralize the local miima, we write m,k m,k 1 = R,k L,k, k = 1,..., h 1 1, 3.1 where R,k = m,k X kh ad L,k = m,k 1 X kh measure the distaces betwee the miima o bi Tk ad T k 1, respectively, to the cetral true value X kh betwee both bis. I our high-frequecy framework the drift is asymptotically egligible ad a regular volatility fuctio will be approximated by a piecewise costat fuctio o blocks of the coarse grid. I this

8 M. Bibiger, M. Jirak & M. Reiß/ 8 Fig 4. Distributios of bi-wise miima of the sigal process, oise ad the covolutio. Based o 1 simulated bis with σ = 1, ε i Exp5, h = 1. settig, where X t = X kh + σw t W kh ad σ is determiistic, we may ivoke time-reversibility of Browia motio to see that X t X kh, t Tk 1, ad X t X kh, t Tk, form idepedet Browia motios of variace σ such that R,k, L,k, k = l 1r 1 + 1,..., lr 1, are all idetically distributed ad there is idepedece wheever differet bis are cosidered but R,k ad L,k+1 are depedet. From.5 ad 3.1 we ifer Ψσ kh h = ER,k] + EL,k] ER,k ]EL,k ] = VarR,k, ad similarly for Ψ. The histogram i Figure 4 shows the distributio of R,k equivaletly L,k i the regressio model joitly with the associated histograms for mi i I X t X kh k i ad mi i I ε i. I this situatio the law of k R,k is give as the covolutio betwee a expoetial distributio ad the law of the miimum of Browia motio o the discrete grid Ik. The latter coverges to the law of the miimum of W o, 1], but the simulatios cofirm the kow feature that the laws deviate rather strogly aroud zero for moderate discretisatios. Let us state ad prove a slightly more geeral result. Propositio 3.1. Choose h accordig to.4. Cosider t Tk for fixed k ad suppose that X t = X kh + t kh σ dw s, t Tk. The i the PPP model 1. for all x R P h 1/ R,k > xσ = E exp Kσ ] x + W t + dt. 3. Proof. By coditioig o the Browia motio we ifer from the PPP properties of T j, Y j : P h 1/ R,k > xσ xσh 1/ +X kh W = exp λ t,y dt dy T k

9 M. Bibiger, M. Jirak & M. Reiß/ 9 = exp λσ T k xh 1/ W t W kh + dt. Notig that W s = h 1/ W k+sh W kh, s, 1], is agai a Browia motio, the result follows by rescalig ad takig expectatios. For the regressio-type model the survival fuctio is asymptotically of the same form. Propositio 3.. Choose h accordig to.4. Suppose that X t = X kh + t kh σ dw s, t Tk, for a fixed bi umber k. The i the regressio-type model 1.3 for all x R lim P h 1/ R,k > xσ = E exp Kσ ] x + W t + dt. 3.3 The approximatio error due to o-costat σ ad drift is cosidered i detail i Appedix A.1 ad proved to be asymptotically egligible. This way, the asymptotic aalysis of our estimatio problem leads ito the theory of Browia excursio areas. Let R t be a real radom variable distributed as lim h 1/ R, th 1. The law of R t determies Ψσt via VarR t = 1 Ψσ t. 3.4 The Feyma Kac formula gives a coectio of the right-had side i Propositio 3.1 to a parabolic PDE based o the heat semigroup for Browia motio. We ca prove the followig explicit result o the Laplace trasform which determies the distributio of R t, t, 1]. Propositio 3.3. The Laplace trasform i t of E exp ϑ t ] x + W s + ds with ϑ R satisfies the followig idetity: ñ E exp st t ô ϑ x + W s + ds dt = ϑ 3 ζs x, ϑ, with ζ s x, ϑ = ζ s, x, ϑ1, x + ζ s,+ x, ϑ1, x defied by the fuctios ζ s,+ x, ϑ = π ϑ 1/3 Gi ϑ /3 s sgi ϑ /3 s + ϑ /3 s 1/ sai ϑ /3 s ϑ 1/3 Ai ϑ /3 s Ai ϑ 1/3 x + ϑ /3 s + πgi ϑ 1/3 x + ϑ /3 s, Ç ϑ /3 s 1/ Ai ϑ /3 s + ϑ 1/3 AI ϑ /3 s å ζ s, x, ϑ = sai ϑ /3 s ϑ 1/3 Ai ϑ /3 s s 1 ϑ /3 exp sx + s 1 ϑ /3,

10 M. Bibiger, M. Jirak & M. Reiß/ 1 where Ai is the Airy fuctio which is bouded o the positive half axis, Aix = π 1 cos t 3 /3 + xt dt, ad Gi is the Scorer fuctio bouded o the positive half axis ad we defie AIx = x Aiydy. Gix = π 1 si t 3 /3 + xt dt, This result geeralizes the Laplace trasform of the expoetial itegrated positive part of a Browia motio derived by 31]. Isertig x = ad settig ϑ = 1 reders the result by 31]. A iversio of the Laplace trasform i Propositio 3.3 i order to obtai a explicit form of the distributio fuctio ad the Ψ appears ufeasible as several experts vaily attempted to solve related problems, see 31] ad 3]. Exploitig the strog Markov property of Browia motio together with hittig times, we are able to circumvet this problem i our study of Ψσ, for details we refer to the Appedix. The observatio models 1. ad 1.3 as well as the semi-martigale model 1.1 for X might be idealized. I fiace, effects of surprise elemets ad iformatio processig might occasioally result i violatios of this model, for istace by price jumps. I such situatios, a regularizatio of Ψ 1 ca yield more robust estimatio results. We propose to trucate the estimator o the coarse grid by employig Ψ 1 τ = Ψ 1 τ 3.5 for some τ > istead of Ψ 1, givig the adjusted estimators IV h,r,τ IV h,r,τ i.7 ad.9. The trucatio level τ > is chose such that we ca guaratee sup t,1]\v σt τ almost surely where V deotes the uio of all violated blocks. I practice, ay over-estimated boud from idepedet historical data may work. Furthermore, observe that ay cotiuous process of fiite variatio Āt may corrupt the observatios via Y i = X t i +A t i +ε i without harmig our volatility estimator because it ca be icorporated as a drift ito the ew semi-martigale X t + A t. For order books, the corruptio A t i may accout for spreads due to market processig ad ivetory costs. We formulate ow the mai covergece results whose proofs are give i the Appedix. For that we impose some mild regularity o the drift ad diffusio coefficiet. Moreover, we eed that the fuctio Ψ is ivertible ad sufficietly regular, which by Propositio A.6 below is esured by a sufficietly large choice of K, but at least umerically seems to be the case for much smaller choices, cf. Figure 3. We work uder the geeral structural hypothesis that the volatility is a Itô semi-martigale with fiite activity jumps. This is a stadard assumptio i fiacial volatility estimatio, see e.g. 4] ad ], allowig for stochastic volatility with leverage. To remai cocise, we assume global coditios o the ad

11 M. Bibiger, M. Jirak & M. Reiß/ 11 characteristics, but extesios via localizatio techiques as i Sectio of 1] are clearly possible. Assumptio 3.4. We work i the stochastic volatility model with potetial jumps i X ad σ: X t = X + σ t = σ + t t a s ds + ã s ds + t t σ s dw s + σ s dw s + t t R x dµ X ds, dx, η s dw s + t R x dµ σ ds, dx. with fiite radom measures µ X, µ σ, i.e. µ X + µ σ, 1], R < almost surely. Assume that the volatility is uiformly bouded away from zero, i.e; if t 1 σ t σ almost surely for a determiistic costat σ >. The characteristics a s, ã s, σ s ad η s are progressively-measurable ad uiformly bouded. The costat K i the defiitio.4 of h is chose large eough that Propositio A.6 below applies. Theorem 3.5. Grat Assumptio 3.4 with cotiuous X from 1.1, choose h accordig to.4 ad r = κ 1/3 for some κ >. The the estimator.7 based o observatios from the PPP-model satisfies h IV,r σs ds = O P O the geeral Assumptio 3.4, utilizig the trucatio i sup t,1] σt τ almost surely, it also holds that 3.5 ad if IV h,r,τ σs ds = O P Based o the same strategy of proof we ca obtai a aalogous result for the regressio-type model. Corollary 3.6. Grat Assumptio 3.4, choose h accordig to.4 ad r = κ 1/3 for some κ >. The the estimator.9 based o observatios from the regressio-type model satisfies the same asymptotic properties as estimator.7 i Theorem 3.5. I the rate-optimal balaced setup there are three error cotributios of the same order: the implied observatioal oise o the bis kh, k+1h, the biwise approximatio of X, ad a secod order term o the coarse blocks arisig from the oliearity of Ψ. Their iterplay is o-trivial ad thus a geeral stable cetral limit theorem for the rescaled error does ot seem straight-forward. If we dropped the ambitio of rate-optimality, however, we could udersmooth or oversmooth by a differet choice of the block sizes h ad r such that oly oe or two error terms would prevail for which estimators with a simpler asymptotic distributio theory would be available. This is ot pursued here.

12 M. Bibiger, M. Jirak & M. Reiß/ 1 4. Lower boud for the rate of covergece Cosider our PPP-model 1.. We show that eve i the simpler parametric statistical experimet where X t = σw t, t, 1], ad σ > is ukow the optimal rate of covergece is 1/3 i a miimax sese. This lower boud for the parametric case the serves a fortiori as a lower boud for the geeral oparametric case. A lower boud for the discrete regressio-type model is obtaied i a similar way; i fact the proof is eve simpler, replacig the Poisso samplig T s j below by a determiistic desig of distace /3. Theorem 4.1. We have for ay sequece of estimators ˆσ of σ, from the parametric PPP-model for each σ >, the local miimax lower boud δ > : lim if if max ˆσ σ {σ,σ +δ 1/3 } P σ ˆσ σ δ 1/3 >, where the ifimum exteds over all estimators ˆσ based o the PPP-model 1. with λ = 1 ad X t = σw t. The law of the latter is deoted by P σ. The proof falls ito three mai parts. We first simplify the problem by cosiderig more iformative experimets. These reductios are give i the two steps below. The, i the third step we use bouds for the Helliger distace. The more techical step 3 is worked out i Appedix B. 1. A PPP with itesity Λ is obtaied as the sum of two idepedet PPPs with itesities Λ r ad Λ s, respectively, satisfyig Λ = Λ r + Λ s, see e.g. 5]. Hece, for b > the experimet of observig Ti r, Yr i i 1 from a PPP with regularised itesity desity y λ r t, y = Xt + /b 1 ad idepedetly Tj s, Ys j j 1 from a PPP with discotiuous itesity desity λ s = λ λ r is more iformative. We ow provide eve more iformatio by replacig Tj s, Ys j j 1 by Tj s, X Tj s j 1, the direct observatio of the martigale values at the radom times Tj s. A lower boud proved for observig Ti r, Yr i i 1 ad Tj s, X Tj s j 1 idepedetly thus also applies to the origial less iformative observatios.. Due to λ s t, ydt dy = /3b, we coclude that the times Tj s are give by a Poisso samplig of itesity /3b o, 1] ad there are a.s. oly fiitely may times Tj s j=1,...,j. Let us first work coditioally o Tj s ad put T s =, TJ+1 s = 1. All observatios of T i r, Yr i i 1 with Ti r Tj 1 s, T j s are trasformed via T r i, Y r i T r i T s j 1, Y r i Ti r X T j 1 s Tj s T s j 1 Tj s T j 1 s + X T i r T s j Tj s T j 1 s Notig that B t t/t B T, t, T ] defies a Browia bridge B,T o, T ], we thus obtai coditioally o Tj s for each j = 1,..., J + 1.

13 M. Bibiger, M. Jirak & M. Reiß/ 13 Fig 5. Estimated daily volatilities multiplied by 1 5 for Facebook, August 15. Days o which jumps have bee filtered are highlighted by -symbols. observatios of a PPP o, T s j T s j 1] with itesity desity λ j t, y = b 1 y σb,t s j T s j 1 t + 1. The trasformatio has redered the family of PPPs with itesity desities λ j j=1,...,j+1 idepedet by reducig the Browia motio to piecewise Browia bridges. Coditioally o Tj s we thus have idepedet observatios of Tj s, X Tj s j=1,...,j ad idepedet PPPs with itesity desities λ j j=1,...,j+1. By usig the latter more iformative experimet ad by choosig b 1/3 we show below that for a Poisso samplig Tj s j=1,...,j o, 1] of itesity /3b /3 of direct observatios X T s j as well as for idepedet observatios of PPPs, geerated by σ times a Browia bridge i-betwee the samplig poits Tj s j, we caot estimate at a better rate tha 1/3. This is accomplished by boudig the Helliger distace betwee the experimets for σ = σ ad σ = σ + δ 1/3. 5. Discussio For the applicatio to limit order books we model the relatioship betwee ask quotes ad a efficiet price process by a pure boudary model, ot takig ito accout the fie structure of order book dyamics. This agostic poit of

14 M. Bibiger, M. Jirak & M. Reiß/ 14 view seems attractive for statistical purposes because more complex models will usually require more data for the same estimatio accuracy ad are highly exposed to model misspecificatio. To check whether the semi-martigale boudary model leads to realistic results, we apply our estimator to limit order book data ad compare it to itegrated volatility estimators, which are commoly used for traded prices uder market microstructure oise. We cosider limit order book data of the Facebook asset FB traded at NASDAQ provided by LOBSTER academic data, recorded over the 1 tradig days i August 15. Empirical data aalysis with similar assets lead to comparable results. The August 15 time series has the advatage of startig with a relatively calm period before icorporatig a period of high tradig activity, which ca serve as a kid of stress test to the estimators. We estimate day-wise itegrated volatilities based o 1. Our estimator.9 with trucatio for the regressio-type model ad first level ask quotes log-prices ad symmetrically, but idepedetly bid quotes, called IV i the sequel. The average umber of ewly submitted best ask quotes per day,, i the cosidered period is about 1,. The average for bid quotes is similar, but the differece o certai days may be large. The maximal absolute differece is 3,537 i the cosidered period.. The local method of momet LMM estimator from 6], adjusted to possible jumps with the trucated versio from 8], ad log-prices from trades recostructed from the order book. The average umber of trades per day i the cosidered period is about 43,. The LMM is the asymptotically efficiet estimator for the stadard regular oise model with cetered oise ad we follow the implemetatio with a selectio of tuig parameter described i 7]. For the trucatio step we employ a global threshold τ = logh 1 pre h IV, with the pre-estimator of itegrated volatility obtaied from the first estimatio step of the two-stage adaptive LMM, for both approaches usig their differet bi-legths h ad sample sizes. h is chose i a data-drive way, for the LMM we arrive at about 1 bis per day ad for the IV at about 65 bis. We expect that comparative studies usig alterative estimators for the regular oise model, as e.g. realized kerels from 4] or preaveragig from ], would yield similar results. The results are preseted i Figure 5. It shows a rather close relatioship betwee the three sequeces of estimates. Estimates obtaied from bid ad ask quotes ca differ, but their differeces are very small. O August 4, 15, however, there is a large differece amog the estimators. O that day a flash crash maipulated the traded ad order quote prices challegig all market models. The huge differece amog the trucated estimates is due to the rougher time resolutio of bis for LMM which are equidistat i caledar time. The flash crash led to tremedous price movemets i very short time at the begiig of the tradig day alog with a huge tradig activity. I the regular oise model the LMM ascribes those movemets o its first two time bis to jumps ad tru-

15 M. Bibiger, M. Jirak & M. Reiß/ 15 cates, while our bi widths h are much smaller ad thus IV is still affected by this period because ot all bis are trucated. The same effect explais why the values of IV are sigificatly larger aroud that date. A priori, eve from a microecoomic perspective, it is ot clear whether the assumed efficiet price processes for the differet market mechaisms givig rise to bids, asks ad trades are the same or at least exhibit the same itegrated volatility. A proper statistical test for the latter hypothesis requires a simultaeous distributio theory for IV ad the LMM a CLT aloe is ot sufficiet, which is beyod the scope of the preset work ad a project i its ow right. For the LMM aloe, however, a feasible cetral limit theorem is available, see Theorem 4.4 of 6]. For the Facebook data set we have coducted a test o the hypothesis that the itegrated volatilities i the order book ad trasactio price models coicide, assumig idepedece of the estimators ad a Gaussia limit distributio where the variace of IV does ot exceed the oe of LMM. Applied to 1 tradig days ad at asymptotic level α = 5%, the test has accepted the ull o 14 days ad rejected o 7 days. This testig problem illustrates that more mathematical aalysis of the estimator s risk is highly desirable as well as a more profoud empirical study. Appedix A: Proofs of Sectio 3 Propositio 3. cosiders the simplified model where X t, t Tk, is approximated by X kh + t kh σ kh dw t. The resultig approximatio error is bouded withi Propositio A. for the PPP-model ad a aalogous proof carries over to the regressio-type model. I the sequel, we write A + = A1A, A = A 1A ad Z p = E Z p ] 1/p, p 1. Proof of Propositio 3.. By law ivariace of R,k with respect to k for X t = X + σw t, we ca simplify Ph 1/ R,k > xσ = P h 1/ mi X i/ X + ε i > xσ i=,...,h 1 = P mi W i/h + σ 1 h 1/ ε i > x, i=,...,h 1 where we used that h 1/ W t/h is aother Browia motio. We coditio o the drivig Browia motio W = W t, t, 1] ad obtai i terms of the distributio fuctio F λ of ε i : Ph 1/ R,k > xσ = E h 1 i= = E exp P ε i > σh 1/ x W i/h W ] h 1 i= log 1 F λ σh 1/ ] x W i/h.

16 M. Bibiger, M. Jirak & M. Reiß/ 16 The expasio 1.4 of F λ yields together with expadig the logarithm therefore Ph 1/ R,k > xσ = E exp h 1 σh 1/ λ i= ] x W i/h O1, where O1 is to be uderstood ω-wise ad holds uiformly over i ad wheever max t,1] x W t ω + is bouded. By the choice of h we have h 1/ λ = Kh 1 ad the itegrad is a Riema sum tedig almost surely to exp σk x W t + dt. Notig that a coditioal probability is always bouded by 1, the assertio follows by domiated covergece ad use of W = d W. Proof of Propositio 3.3. Throughout the proof, we drop the depedece o ϑ i ζ s x, ϑ, ζ s, x, ϑ ad ζ s,+ x, ϑ to lighte the otatio. We shall apply the Kac formula i the versio as i formulae 4.13 ad 4.14 of 4]. It coects the cosidered Laplace trasform with the solutio of a differetial equatio which becomes i our case: d ζ dx = sζ ϑ/3, x <, A.1a d ζ dx = ϑx + sζ ϑ /3, x >. A.1b Sice all assertios ecessary to apply the Kac formula are fulfilled, the Laplace trasform from above multiplied with a costat Lagragia ϑ /3 satisfies ñ E ϑ /3 exp st t ô ϑ x + W s + ds dt = ζ s x. The geeral solutio of A.1a is give by ζ s, x = A exp sx + ϑ /3 s 1, A.a with a costat A depedig o s but ot o x. Airy s fuctio Ai solves the homogeous differetial equatio of the type A.1b, whereas the Scorer fuctio Gi is a particular solutio of the ihomogeous equatio ζ xζ = π 1, both beig bouded o the positive real lie. Hece, a solutio asatz for A.1b is give by ζ s,+ x = BAi ϑ 1/3 x + ϑ /3 s + πgi ϑ 1/3 x + ϑ /3 s, A.b with a costat B. Cotiuity coditios o ζ ad dζ/dx at x = give rise to B = π ϑ 1/3 Gi ϑ /3 s sgi ϑ /3 s + ϑ /3 s 1/ sai ϑ /3 s ϑ 1/3 Ai ϑ /3 s.

17 M. Bibiger, M. Jirak & M. Reiß/ 17 I order to express A i a more cocise ad simple maer, we exploit the followig relatio for the Wroskia of Ai ad Gi: π Gi xaix Ai xgix = AIx = Aiy dy. A.3 A proof of the latter equality ca be foud i 35]. Thereby, we obtai Ç ϑ /3 s 1/ Ai ϑ /3 s + ϑ 1/3 AI ϑ /3 s å A = sai ϑ /3 s ϑ 1/3 Ai ϑ /3 s s 1 ϑ /3. This result cocludes the proof. x A.1. Asymptotic aalysis of the estimator Recall that due to Assumptio 3.4 we ca assume without loss of geerality that a, ã, σ, η C + a.s., ad that if t 1 σ t σ > a.s. Here, C + ad σ are absolute costats. From here o A B expresses shortly that A K B for two sequeces A, B ad some real costat K <. We use the otatio A a.s. B if this holds P-almost surely. Similarly, we write A = O a.s. B ad use A a.s. B for short otatio. We also write P k = P F kh ad aalogously for the coditioal expectatio. Moreover, we use k Z q = E k Z q ] 1/q, q 1. First, we establish Theorem 3.5 ad Corollary 3.6 o Assumptio 3.4 ad i absece of jumps i X ad σ for estimators.7 ad.9, respectively. Robustess of the trucated versios agaist violatios is proved at the ed of this sectio. As a first step, we aalyze the approximatio error assumig a locally costat volatility ad eglectig the drift. The we prove Theorem 3.5 exploitig properties of Ψ which are established i Appedix A.. We shall use the followig idetities for momets of real radom variables: E X ] = P X > x dx P X > x dx, A.4a E X ] = xp X > x dx + xp X > x dx. Lemma A.1. For ay fiite p > 1 ad x >, s t 1, we have A.4b P sup σr σ s xσ s F s a.s. t s /x x p/. s r t Proof of Lemma A.1. Sice ã C + a.s; we get that t s ã r d r C + t s almost surely. Hece usig Markov s ad Burkholder s iequality, for ay fiite p > 1 P sup σ r σ s z F s a.s. t s /z p/, z >, A.5 s r t

18 where we also used σ, η C + a.s. Sice we obtai that M. Bibiger, M. Jirak & M. Reiß/ 18 σ r σ s = σ r σ s + σ s σ r σ s, P sup σr σ s xσ s F s P sup σ r σ s s r t s r t»xσ s / F s + P sup σ r σ s x/4 F s. s r t Hece usig A.5 ad if s 1 σ s σ a.s., the claim follows. Propositio A.. Cosider h i.4 ad t Tk Å ã P mi Y j X kh > xσ kh h Fkh j T k Fkh ] Kσ kh = E exp x + W t + dt for fixed k. The + λ 1/3 Gx, a.s., where W t t 1 is a stadard Browia motio idepedet of F t t 1, Gx is determiistic ad x p Gx L 1 R for ay fiite p. If σ t is costat ad a t = for t Tk, the Gx =. Proof of Propositio A.. Throughout the proof, C, C 1,... deote positive, geeric costats that may vary from lie to lie. Propositio 3.1 already gives the last statemet i case of o drift ad bi-wise costat volatility. Let A z = Tk, z] ad z = xσ kh h. A.6 Let X t k = t kh σ s dw s ad A t k = t kh a s ds. The, usig basic properties of a PPP, it follows that Å ã ] ] P k mi Y j X kh > z = E k P Λ Az = X = E k exp ΛA z j T k ] = E k exp λ 1 { Xtk+ A tk y} dt dy A z ] = E k exp λ z X t k A t k + dt. Itroduce T k = λ V k = λ U k = λ T k k+1h T k z X t k A t k + dt, z σkh W t W kh + dt, kh t σ s σ kh dw s dt ad A k = λh max a t. kh t T k k+1h kh A.7

19 The we have the upper ad lower bouds M. Bibiger, M. Jirak & M. Reiß/ 19 V k U k A k T k V k + U k + A k. By scalig ad symmetry properties of Browia motio, we have that V k d = λσkh h 3/ W t + x + dt = Kσ kh W t + x + dt, with a stadard Browia motio W idepedet of F t t 1. I the sequel, we distiguish the two cases where x 1 ad x < 1. Case x 1: As a first objective, we derive a upper boud for E k exp yt k ], y >. To this ed, ote that by the Dambis-Dubis-Schwarz Theorem Thm. 4.6 i 4], o a possibly larger probability space extedig time ad processes from T = 1 to T =, there exists a Browia motio Wt idepedet of F kh such that X t k d = W Xk, Xk t. Lemma A.1 yields that for x R Ç t å P k sup σs σkh ds x σ kh a.s. h / x p / 1, p >. kh kh t k+1h Sice Xk, Xk t = t kh σs ds for t kh, we deduce P k sup X t k z / a.s. P k t T k sup t x σ kh Ç å a.s. P k sup Wt C σ kh x t 1 σ kh + + h / x p / 1 x σ kh a.s. P W 1 C» x + h / x p / 1, W t z / + h / x p / 1 A.8 with some C >. Next, observe that by the boudedess of a t, it follows that A k a.s. λh max t,1] a t a.s. λ 1/3. We thus obtai for y > the upper boud ] E k exp yt k a.s. P» W 1 C x + h / x p / 1 + exp C 1 x yh σ kh / + C y 1/3 + 1 x <, A.9 A.1

20 M. Bibiger, M. Jirak & M. Reiß/ where p arbitrarily large but fiite. Note that elemetary calculatios yield that A.1 also supplies a boud for E k exp yv k ]. Next, observe that t U k h λ sup σ s σ kh dw s def = U + k. kh kh t k+1h Applyig Burkholder s iequality, we get k+1h k U + k q a.s. λh σ s σ kh k ds kh hece aother applicatio of Burkholder s iequality yields that k U + k q a.s. λh 3/ h 1/ = λ 1/3 K 4/3, A.11 where { we also used that ã, σ, η < almost surely. Set U k, = U + k δ }. By the Markov iequality ad A.11, it follows for q 1 ad δ = λ 1/6 P k U c k, a.s. δ q λ q/3 = λ q/6. Usig the power series of expx ad A.9, we obtai by Cauchy-Schwarz ] E k exp T k 1 Uk, ] expu + E k exp V k + E k exp V k k + A k 1 ] 1 Uk, ] E k exp V k + k exp V k k expu + k + A k 1 1 Uk, ] Å k E k exp V k + C k 3 exp V k U + k + δ + λ ã, 1/3 1/ q/, A.1 with some costat C 3. From A.11, this is bouded by ] E k exp V k + C k 4 exp V k λ 1/3, C 4 >. A.13 O the other had, it follows from Cauchy-Schwarz that ] E k exp T k 1 U c» k, k exp T k P k Uk, c. A.14 Combiig the above, we thus coclude from A.1 y =, p large eough ad A.1 q large eough for some costat C 5 : ] ] E k exp T k a.s. E k exp V k +C k 5 exp V k λ 1/3 +» k exp T k P k U ck, a.s. E k exp ÅP» W 1 C 5 x + V k ] + C 5 λ 1/3 Å ã p /4 1 h + 1x +exp ã 1/ x h σ kh C 5. x

21 M. Bibiger, M. Jirak & M. Reiß/ 1 I the same maer oe obtais a lower boud ad hece the claim follows for x, sice if t 1 σ t σ > almost surely. Case x < 1: Let V,k z = { sup σ kh W t W kh < z /4 }, t T k X,k z = { sup X t k < z / }, t T k A.15 ad deote with V,k c z ad X,k c z their complemets. Observe that sice x 1 ad if s 1 σ s σ ad A t k a.s. h max a t a.s. h, t,1] we get = T k = V k for large eough o the set V,k z X,k z. Hece we obtai from A.8 that ] E k exp T k exp V k ] E k 1 V,k z X c,k z 1 V,k z a.s. P W 1 C 6» x + h / x p / 1, with some C 6 >. O the other had, usig very similar argumets as for the case x 1, oe derives that ] E k exp T k exp V k 1 V c,k z 1 a.s. ÅP» W λ 1/3 1 C 6 x Å ã p /4 1 h + +exp ã x h σ kh C 6 1/, x which completes the proof. Deote with R,k σ the versio of R,k where σt = σ is costat for t Tk T k 1, ad we use the same otatio L,k σ for L,k. It is apparet from the proof of Propositio A. that σ kh ca be replaced with ay σ k jh where j is fiite ad idepedet of. Lemma A.3. For p 1 we have h 1/ R p,k, h 1/ L p,k <. Proof of Lemma A.3. Usig the idetity E X p] = p x p 1 PX > x + 1 PX x dx, the claim follows from Propositio A. ad the tower property of coditioal expectatio. Lemma A.4. We have the equality E h 1 ] L,k R Fk 1h,k = Ψσ k 1h + O a.s. h 1/.

22 M. Bibiger, M. Jirak & M. Reiß/ Proof of Lemma A.4. Usig Propositio A. ad relatios A.4a, A.4b, oe readily computes that E ] ] L,k F k 1h = E L,k σ k 1h F k 1h + Oa.s. h, E L ] ],k F k 1h = E L,k σ k 1h F k 1h + Oa.s. h 3/, A.16 ad similarly E ] ] R,k F k 1h = E R,k σ k 1h F k 1h + Oa.s. h, E R ] ],k F k 1h = E R,k σ k 1h F k 1h + Oa.s. h 3/. A.17 Hece by A.16 ad A.17, we obtai E k 1 L,k R,k ] = E k 1 L,k σ k 1h ] + E k 1 R,k σ k 1h ] ] + O a.s. h 3/ E k 1 L,k R,k. A.18 It thus suffices to cosider the cross term i the last lie. Note that by the tower property of coditioal expectatios E k 1 L,k R,k ] = E k 1 L,k E k R,k ] ] = E k 1 L,k E k R,k ] E k 1 R,k σ k 1h ] + E k 1 R,k σ k 1h ] ] = E k 1 L,k E k R,k ] E k 1 R,k σ k 1h ] ] + E k 1 L,k E k 1 R,k σ k 1h ] ]. By Lemma A.3 ad A.17 E k 1 L,k E k R,k ] E k 1 R,k σ k 1h ] ] = Oa.s. h E k L,k ] The same argumets as above lead to E k 1 L,k E k 1 R,k σ k 1h ] ] = E k 1 L,k σ k 1h ] E k 1 R,k σ k 1h ] + O a.s. h 3/ = O a.s. h 3/. A.19. A. Sice E k 1 R,k σ k 1h ] = E k R,k σ k 1h ], we have by the tower property of coditioal expectatios E k 1 L,k σ k 1h ] E k 1 R,k σ k 1h ] = E k 1 E k L,k σ k 1h R,k σ k 1h ] ] = E k 1 L,k σ k 1h R,k σ k 1h ]. A.1 Combiig A.19, A. ad A.1, we obtai E k 1 L,k R,k ] = E k 1 L,k σ k 1h R,k σ k 1h ] + O a.s. h 3/ ad hece the claim follows.,

23 M. Bibiger, M. Jirak & M. Reiß/ 3 Lemma A.5. Let Ψx = Ψx, a,l = lr 1 /, s,l = lh r 1 = a,l h ad M k, = h 1 m,k m,k 1 r. We the have the followig upper bouds i ii iii h 1 r Ψ 1 a Ψσ,l s,l 1 r h 1 r a,l r k=a,l 1 +1 r h 1 σ s,l 1 h r 1 k=a,l 1 +1 Ψσ kh Ψσs,l 1 = O 1, r h 1 a,l k=a,l 1 +1 Ψσ kh Ψσ s,l 1 1 = O 1, σkh h r 1 1/3. 1 Proof of Lemma A.5. Note first that due to A.36 ad A.37 below, it follows that for ay p 1 sup Ψ σ t p, sup Ψ x <. A. t,1] x σ Case i: By Itô s formula ad the idepedece of W s ad Ws 3.4 Ψσ t Ψσ s = t Usig A., we thus obtai s Ψ σ r dσ r + 1 t s i Assumptio Ψ σ r σ r + η r dr. A.3 h 1 r Ψ 1 a Ψσ,l s,l 1 r Ψσkh k=a,l 1 +1 Ψσ 1 s,l 1 h 1 r Ψ 1 a Ψσ,l kh s,l 1 r Ψ σ r σ r dw r + η r dwr 1 + h 1 r r a,l k=a,l 1 +1 a,l kh k=a,l 1 +1 h 1 r Ψ 1 a Ψσ,l s,l 1 r k=a,l 1 +1 a,l h a,l k k+1h kh Ψ σ r σ r dw r + η r dw r 1 + O 1. Jese s iequality gives a boud with the orm for the squared 1 orm above ad Burkholder s iequality ad A. the yield the upper boud a,l k=a,l 1 +1 a,l k h = O 1. Combiig all bouds, we thus h 1 r r obtai i. Case ii ad iii ca be hadled i a very similar maer ad we omit the proofs.

24 M. Bibiger, M. Jirak & M. Reiß/ 4 Proof of Theorem 3.5 ad Corollary 3.6. Observe that Mk, = M k, E ] M k, F k 1h is a sequece of martigale differeces. Lemma A.3 yields that all momets of Mk, exist. Hece for ay idex set J {,..., h 1 1 }, Burkholder s iequality esures that for ay p 1 p M k, r» J, A.4 k J where J is the cardiality of the set J. Let ß M l = a,l k=a,l 1 +1 M k, + Ψ σ s,l 1 / >, ad M = h 1 1 M l. Propositio A.6 yields that Ψ x > for x >. The we obtai from the Markov iequality ad A.4 that P Åh 1 1 M c l l= ã h 1 1 l= h 1 1 l= P M c l p h 1 1 l= Ψ σ s,l 1 p a,l k=a,l 1 +1 M k, p r p/ = O 1, A.5 for p > 4. We are ow ready to proceed to the mai proof. From.7 it follows that IV h,r + r h 1 σ t dt = Ψ 1 r h 1 a,l k=a,l 1 +1 σ s,l 1 h r 1 σ t dt m,k m,k 1 h 1 r σs,l 1 h r 1. Cosider first the approximatio error i the quadratic variatio by settig the volatility locally costat o the blocks of the coarse grid. Due to Lemma A.5 iii, it suffices to cosider the error h 1 σlh h σt dt h 1 lh l 1h σ t σ lh dt, which by the triagle ad Burkholder s iequality is bouded by h 1 lh l 1h h + h ã dt 1/3. p

25 M. Bibiger, M. Jirak & M. Reiß/ 5 I order to boud the remaiig estimatio error r h 1 a,l Ψ 1 M k, σs,l 1 h r 1, k=a,l 1 +1 we use a Taylor expasio ad that the first two derivatives of Ψ 1 exist ad are bouded accordig to Propositio A.6 below. To this ed, deote with k,l, M, Ψ = E M k, F k 1h ] Ψσ kh + Ψσ kh Ψσ s,l 1. It the follows that Å a,l Ψ 1 k=a,l 1 +1 M k, Ψ σ s,l 1 + Ψ σ s,l 1 ã = σs,l 1 + Ψ 1 Ψσ s,l 1 Å a,l k=a,l 1 +1 Mk, + k,l, M, Ψ ã + 1 Ψ 1 Å a,l ξl Mk, + k,l, M, Ψ ã k=a,l 1 +1 def = σ s,l 1 + Ψ l,1 + Ψ l,, where ξ l Ψ σs,l 1 / > o the set Ml. We first deal with Ψ. To this l,1 ed, deote with Z l = a,l M k=a,l 1 +1 k,, which is a partial sum of martigale differeces. Hece, by Burkholder s iequality see also A.4, we obtai r h 1 Ψ 1 Ψσ s,l 1 Z l r h 1 Zl r h 1 = O 1. A.6 O the other had, we obtai from Lemma A.4 ad Lemma A.5 i that r h 1 Ψ 1 a Ψσ,l s,l 1 k=a,l 1 +1 Combiig A.6 ad A.7, we fid r h 1 k,l, M, Ψ = O P 1. A.7 Ψ 1 l,1 M = O P 1. A.8 I the same maer, but usig Lemma A.5 ii ad additioally Ψ 1 ξ l 1 Ml < by Propositio A.6 below, we obtai r h 1 Ψ 1 l, M 1 r h 1 Zl + a,l k=a,l 1 +1 k,l, M, Ψ = O1.

26 M. Bibiger, M. Jirak & M. Reiß/ 6 Sice P M c = O 1 by A.5, this suffices to guaratee that IV h,r ad usig aalogous bouds ad Propo- Based o a Taylor expasio for Ψ 1 sitio A.7, we obtai likewise IV h,r σ t dt = O P 1/3. σ t dt = O P 1/3. ad coclude Corollary 3.6. This completes the proof of Theorem 3.5 ad Corollary 3.6 i absece of jumps for estimators.7 ad.9, respectively. Fially, cosider the trucated estimators. Sice τ < ad Psup t,1]\v σ t τ = 1, the argumets above reveal that it suffices to show r h 1 τ Ψ 1 a,l k=a,l 1 +1 M k, σs lh,l 1 1 h r V r = O P 1/3, uiformly for N. However, sice we have that sup t 1 σt < almost surely ad r h 1 1 lh /r V is fiite almost surely, the left-had side above is bouded by r h 1 O P τ + sup σ t 1 lh /r V h r 1 = O P τ + sup σ t V 1/3 t 1 t 1 = O P 1/3. Therefore, it suffices to show that σ s1 s V ds = OP 1/3. From the claim follows. σs1 s V ds sup σs V r 1 h = O P 1/3, s 1 A.. Properties of Ψ It follows from Propositio A. ad A.4a, A.4b that for Ψσ = Ψσ h 1 E m,k m,k 1 ] = Ψσ + O λ 1/3. A.9 Havig uderstood the behaviour of Ψσ, aalogue properties of Ψσ readily follow. Let Hx = W t + x + dt. A.3

27 M. Bibiger, M. Jirak & M. Reiß/ 7 The by A.4a, A.4b, we derive Ψ σ Å = 4σ x E î e KσHxó + 1 E î e KσH xóã dx Å Å σ E î e KσHxó 1 + E î e KσH xóã dxã. A.31 Next, cosider the distributio o the egative half axis. With x <, we make the decompositio E ïe σk 1 ò Wt x dt = E ïe σk 1 Wt x dt 1 if W t x + 1 if W t x ò t 1 t 1 def = U 1 x + U x. Let T x be the first passage time of W to level x with desity f Tx t = x πt 3 e x /t, t, see 6.3 i Sectio.6 of 4]. From { T x 1 } = { if t 1 W t x } it follows from the strog Markov property of W that U 1 x = E ïe σk 1 ò Wt x dt s T x = s f Tx sds = Usig a time shift yields U 1 x = E ïe σk 1 s W t dt We the obtai that E ïe σk 1 ò Wt x dt = P if W t x + t 1 ò f Tx sds. E ïe σk1 s3/ 1 Wt dt òf Tx sds. E ïe σk1 s3/ 1 Wt dt òf Tx sds def = P if W t x + A x, for x <. A.3 t 1 Let IKσ, s = E ïe Kσ1 s3/ 1 Wt dt ò. The by A.3 = x 1 E î e KσH xó dx = x P if W t < x A x dx t 1 = 1 1 xp if W t < xdx t 1 IKσ, sds, IKσ, s xf Tx sdxds

28 M. Bibiger, M. Jirak & M. Reiß/ 8 sice xpif t 1 W t < xdx = 1. Likewise, it follows that î ó 1 E e KσH x dx = P if W t < x dx t 1 1 = π IKσ, sds. π We thus obtai Ψ σ Å = 4σ xe î e KσHxó dx Å σ E î e KσHxó dx π + π π ã IKσ, sds IKσ, sds ã IKσ, sds A.33 = σ Λ 1 σ Λ σ, A.34 with fuctioals Λ 1, Λ. I the sequel, we write k fx = k fx/ k x. The further aalysis of properties of Ψ is structured i several lemmas which combied imply the followig key propositio. Propositio A.6. Suppose that σ σ >, K Cσ for Cσ sufficietly large the exact value of Cσ follows from A.39. The we have uiformly for σ σ Ψ σ Å = 4σ 1 ã Ç σ 3 + O π Ψ σ Å = σ 1 ã + O π K 1 3 Ç σ 3 K 1 3 å å > ad >. A.35 Moreover, it holds that Usig the relatio sup Ψσ <. σ σ A.36 Ψ 1 ϱ = 1 Ψσ, Ψσ = ϱ, A.37 we get that the secod derivative is uiformly bouded for σ σ = Ψ 1 ϱ, i.e. sup Ψ 1 ϱ = sup Ψσ 3 <. A.38 ϱ ϱ Ψσ σ σ So far we have focused o results for Ψσ = Ψσ. Essetially the same results are valid for Ψ σ, which we state ow.

29 M. Bibiger, M. Jirak & M. Reiß/ 9 Propositio A.7. Itroduce Å ã B,1 = xp max i/h x dx i h 1 Å ã B, = P max i/h x i h 1 dx,» which satisfy B,1 1, B, π. The A.35 ad A.36 i Propositio A.6 remai valid if we replace Ψ σ with Ψ σ ad 1 π with B,1 B,. Likewise, A.38 also holds. Proof of Propositio A.6. We write shortly Λ 1 = Λ 1 σ, Λ = Λ σ. We have that ad Ψ σ = 4σ Λ 1 Λ + σ Λ 1 Λ Λ. Usig Lemmas A.8 ad A.11 from below, we obtai with A.34 Λ1 1 Å ã 6Kσ 3 Kσ def = R Λ1 logkσ Å ã Λ Kσ π π Kσ def = R Λ. π logkσ Moreover, applyig Lemmas A.9, A.1 ad A.11 yields Λ 1 6 Kσ Å ã Kσ 5 def = R Λ1, Λ σ logkσ π, ãå Λ Λ ÅR Λ + 4 Å 1 + π 1 K 1 3 σ π π σ def = R Λ. Kσ logkσ We thus obtai from the above that Ψ 4σ 1 ã 4σ ÅR Λ1 + R Λ R Λ + π π + 4σ R Λ1 + R Λ = O K 1 3 σ 3, Ψ σ 1 σ ÅR Λ1 + R Λ R Λ + π ã π = O K 1 3 σ 3 ã 5 ã A explicit sufficiet lower boud for K i terms of σ ca be computed from the two coditios 1 π > R Λ1 + R Λ R Λ + + σ R Λ1 + R Λ, π 1 π > R Λ 1 + R Λ R Λ +. A.39 π.

30 M. Bibiger, M. Jirak & M. Reiß/ 3 It remais to show the boudedess property for the first two derivatives of Ψ. By Lemma A.1, we have k J σ C 1 + σ k, k = 1,, A.4 where J σ = Λ 1 σ Λ σ ad C is a costat ot depedig o σ. Observe that hece the claim follows. Ψ = 4J + 6σ J + σ J, Proof of Propositio A.7. The proof ca be redirected to Propositio A.6 usig Propositio 3. ad a trucatio argumet for the itegrals over x. The correspodig computatios are very similar to those above ad the Lemmas give below. We therefore omit the details. Lemma A.8. For K >, p N, we obtai the followig decay behaviour of the momet itegrals: x p E î e KσHxó dx p+1 E Z p+1 ] + Γp + 1 Kσ p+1/3 p + 1 with Z N, 1. Proof of Lemma A.8. The followig useful relatio i terms of the N, 1- distributio fuctio Φ is derived from the law of the miimum of Browia motio: P T x l = 1 Φ x / l. A.41 The for < l < 1 x p E î e KσHxó dx By A.41 we have R 1 = + x p 1 Φx/ 4l dx = 4l p+1/ x p E 1 T x/ l ] dx x p E î e KσHx 1 T x/ > l ó dx def = R 1 + R. z p 1 Φz dz = 4l p+1/ E Z p+1 ]. p + 1 We further ote that T x/ > l implies Hx l x/ + x +dt = lx/, such that R x p e Kσlx/ dx = Å ã p+1 Γ p + 1. Kσl

31 Choosig l = Kσ /3, we obtai as asserted. x p E î e KσHxó dx Lemma A.9. Let K >. The M. Bibiger, M. Jirak & M. Reiß/ 31 p+1 E Z p+1 ] p p+1 Γ p + 1 Kσ p+1/3, E î KHxe KσHxó dx / π K 1/3 σ 4/3. Proof of Lemma A.9. We make the decompositio E î KσHxe KσHxó dx = + with some l >. Usig ye y 1 ad A.41, we obtai E î 1T x/ lkσhxe KσHxó dx E î 1T x/ > lkσhxe KσHxó dx, E î 1T x/ lkσhxe KσHxó dx 1 Φx/ 4l dx» = 8l/π. Now, usig ye y e y/ ad T x/ > l Hx lx/, we boud the other term by E î 1T x/ > lkσhxe KσHxó dx The choice l = Kσ /3 ad divisio by σ yield the claim. Lemma A.1. Let K >. The xe î KHxe KσHxó dx e Kσlx/4 dx = 4 Kσl. 6 K 1/ σ 3/. Proof of Lemma A.1. We proceed as i the proof of Lemma A.9 ad obtai for ay l > xe î KσHxe KσHxó dx x 1 Φx/ 4l + e Kσlx/4 dx The result follows with l = Kσ 1/. Lemma A.11. Let K σ 1. The Kσ IKσ, sds 3 logkσ = l + Kσl/4 1. /5, IKσ, sds σ 3 σ Kσ /5. logkσ

5. Best Unbiased Estimators

5. Best Unbiased Estimators Best Ubiased Estimators http://www.math.uah.edu/stat/poit/ubiased.xhtml 1 of 7 7/16/2009 6:13 AM Virtual Laboratories > 7. Poit Estimatio > 1 2 3 4 5 6 5. Best Ubiased Estimators Basic Theory Cosider agai

More information

1 Estimating sensitivities

1 Estimating sensitivities Copyright c 27 by Karl Sigma 1 Estimatig sesitivities Whe estimatig the Greeks, such as the, the geeral problem ivolves a radom variable Y = Y (α) (such as a discouted payoff) that depeds o a parameter

More information

Rafa l Kulik and Marc Raimondo. University of Ottawa and University of Sydney. Supplementary material

Rafa l Kulik and Marc Raimondo. University of Ottawa and University of Sydney. Supplementary material Statistica Siica 009: Supplemet 1 L p -WAVELET REGRESSION WITH CORRELATED ERRORS AND INVERSE PROBLEMS Rafa l Kulik ad Marc Raimodo Uiversity of Ottawa ad Uiversity of Sydey Supplemetary material This ote

More information

point estimator a random variable (like P or X) whose values are used to estimate a population parameter

point estimator a random variable (like P or X) whose values are used to estimate a population parameter Estimatio We have oted that the pollig problem which attempts to estimate the proportio p of Successes i some populatio ad the measuremet problem which attempts to estimate the mea value µ of some quatity

More information

Statistics for Economics & Business

Statistics for Economics & Business Statistics for Ecoomics & Busiess Cofidece Iterval Estimatio Learig Objectives I this chapter, you lear: To costruct ad iterpret cofidece iterval estimates for the mea ad the proportio How to determie

More information

Unbiased estimators Estimators

Unbiased estimators Estimators 19 Ubiased estimators I Chapter 17 we saw that a dataset ca be modeled as a realizatio of a radom sample from a probability distributio ad that quatities of iterest correspod to features of the model distributio.

More information

18.S096 Problem Set 5 Fall 2013 Volatility Modeling Due Date: 10/29/2013

18.S096 Problem Set 5 Fall 2013 Volatility Modeling Due Date: 10/29/2013 18.S096 Problem Set 5 Fall 2013 Volatility Modelig Due Date: 10/29/2013 1. Sample Estimators of Diffusio Process Volatility ad Drift Let {X t } be the price of a fiacial security that follows a geometric

More information

Lecture 9: The law of large numbers and central limit theorem

Lecture 9: The law of large numbers and central limit theorem Lecture 9: The law of large umbers ad cetral limit theorem Theorem.4 Let X,X 2,... be idepedet radom variables with fiite expectatios. (i) (The SLLN). If there is a costat p [,2] such that E X i p i i=

More information

Lecture 4: Parameter Estimation and Confidence Intervals. GENOME 560 Doug Fowler, GS

Lecture 4: Parameter Estimation and Confidence Intervals. GENOME 560 Doug Fowler, GS Lecture 4: Parameter Estimatio ad Cofidece Itervals GENOME 560 Doug Fowler, GS (dfowler@uw.edu) 1 Review: Probability Distributios Discrete: Biomial distributio Hypergeometric distributio Poisso distributio

More information

Combining imperfect data, and an introduction to data assimilation Ross Bannister, NCEO, September 2010

Combining imperfect data, and an introduction to data assimilation Ross Bannister, NCEO, September 2010 Combiig imperfect data, ad a itroductio to data assimilatio Ross Baister, NCEO, September 00 rbaister@readigacuk The probability desity fuctio (PDF prob that x lies betwee x ad x + dx p (x restrictio o

More information

14.30 Introduction to Statistical Methods in Economics Spring 2009

14.30 Introduction to Statistical Methods in Economics Spring 2009 MIT OpeCourseWare http://ocwmitedu 430 Itroductio to Statistical Methods i Ecoomics Sprig 009 For iformatio about citig these materials or our Terms of Use, visit: http://ocwmitedu/terms 430 Itroductio

More information

x satisfying all regularity conditions. Then

x satisfying all regularity conditions. Then AMS570.01 Practice Midterm Exam Sprig, 018 Name: ID: Sigature: Istructio: This is a close book exam. You are allowed oe-page 8x11 formula sheet (-sided). No cellphoe or calculator or computer is allowed.

More information

Sequences and Series

Sequences and Series Sequeces ad Series Matt Rosezweig Cotets Sequeces ad Series. Sequeces.................................................. Series....................................................3 Rudi Chapter 3 Exercises........................................

More information

5 Statistical Inference

5 Statistical Inference 5 Statistical Iferece 5.1 Trasitio from Probability Theory to Statistical Iferece 1. We have ow more or less fiished the probability sectio of the course - we ow tur attetio to statistical iferece. I statistical

More information

A random variable is a variable whose value is a numerical outcome of a random phenomenon.

A random variable is a variable whose value is a numerical outcome of a random phenomenon. The Practice of Statistics, d ed ates, Moore, ad Stares Itroductio We are ofte more iterested i the umber of times a give outcome ca occur tha i the possible outcomes themselves For example, if we toss

More information

FINM6900 Finance Theory How Is Asymmetric Information Reflected in Asset Prices?

FINM6900 Finance Theory How Is Asymmetric Information Reflected in Asset Prices? FINM6900 Fiace Theory How Is Asymmetric Iformatio Reflected i Asset Prices? February 3, 2012 Referece S. Grossma, O the Efficiecy of Competitive Stock Markets where Traders Have Diverse iformatio, Joural

More information

Introduction to Probability and Statistics Chapter 7

Introduction to Probability and Statistics Chapter 7 Itroductio to Probability ad Statistics Chapter 7 Ammar M. Sarha, asarha@mathstat.dal.ca Departmet of Mathematics ad Statistics, Dalhousie Uiversity Fall Semester 008 Chapter 7 Statistical Itervals Based

More information

Chapter 8. Confidence Interval Estimation. Copyright 2015, 2012, 2009 Pearson Education, Inc. Chapter 8, Slide 1

Chapter 8. Confidence Interval Estimation. Copyright 2015, 2012, 2009 Pearson Education, Inc. Chapter 8, Slide 1 Chapter 8 Cofidece Iterval Estimatio Copyright 2015, 2012, 2009 Pearso Educatio, Ic. Chapter 8, Slide 1 Learig Objectives I this chapter, you lear: To costruct ad iterpret cofidece iterval estimates for

More information

ECON 5350 Class Notes Maximum Likelihood Estimation

ECON 5350 Class Notes Maximum Likelihood Estimation ECON 5350 Class Notes Maximum Likelihood Estimatio 1 Maximum Likelihood Estimatio Example #1. Cosider the radom sample {X 1 = 0.5, X 2 = 2.0, X 3 = 10.0, X 4 = 1.5, X 5 = 7.0} geerated from a expoetial

More information

Asymptotics: Consistency and Delta Method

Asymptotics: Consistency and Delta Method ad Delta Method MIT 18.655 Dr. Kempthore Sprig 2016 1 MIT 18.655 ad Delta Method Outlie Asymptotics 1 Asymptotics 2 MIT 18.655 ad Delta Method Cosistecy Asymptotics Statistical Estimatio Problem X 1,...,

More information

Proceedings of the 5th WSEAS Int. Conf. on SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 17-19, 2005 (pp )

Proceedings of the 5th WSEAS Int. Conf. on SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 17-19, 2005 (pp ) Proceedigs of the 5th WSEAS It. Cof. o SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 7-9, 005 (pp488-49 Realized volatility estimatio: ew simulatio approach ad empirical study results JULIA

More information

0.1 Valuation Formula:

0.1 Valuation Formula: 0. Valuatio Formula: 0.. Case of Geeral Trees: q = er S S S 3 S q = er S S 4 S 5 S 4 q 3 = er S 3 S 6 S 7 S 6 Therefore, f (3) = e r [q 3 f (7) + ( q 3 ) f (6)] f () = e r [q f (5) + ( q ) f (4)] = f ()

More information

Sampling Distributions and Estimation

Sampling Distributions and Estimation Cotets 40 Samplig Distributios ad Estimatio 40.1 Samplig Distributios 40. Iterval Estimatio for the Variace 13 Learig outcomes You will lear about the distributios which are created whe a populatio is

More information

Online appendices from The xva Challenge by Jon Gregory. APPENDIX 10A: Exposure and swaption analogy.

Online appendices from The xva Challenge by Jon Gregory. APPENDIX 10A: Exposure and swaption analogy. APPENDIX 10A: Exposure ad swaptio aalogy. Sorese ad Bollier (1994), effectively calculate the CVA of a swap positio ad show this ca be writte as: CVA swap = LGD V swaptio (t; t i, T) PD(t i 1, t i ). i=1

More information

The Valuation of the Catastrophe Equity Puts with Jump Risks

The Valuation of the Catastrophe Equity Puts with Jump Risks The Valuatio of the Catastrophe Equity Puts with Jump Risks Shih-Kuei Li Natioal Uiversity of Kaohsiug Joit work with Chia-Chie Chag Outlie Catastrophe Isurace Products Literatures ad Motivatios Jump Risk

More information

SELECTING THE NUMBER OF CHANGE-POINTS IN SEGMENTED LINE REGRESSION

SELECTING THE NUMBER OF CHANGE-POINTS IN SEGMENTED LINE REGRESSION 1 SELECTING THE NUMBER OF CHANGE-POINTS IN SEGMENTED LINE REGRESSION Hyue-Ju Kim 1,, Bibig Yu 2, ad Eric J. Feuer 3 1 Syracuse Uiversity, 2 Natioal Istitute of Agig, ad 3 Natioal Cacer Istitute Supplemetary

More information

Lecture 4: Probability (continued)

Lecture 4: Probability (continued) Lecture 4: Probability (cotiued) Desity Curves We ve defied probabilities for discrete variables (such as coi tossig). Probabilities for cotiuous or measuremet variables also are evaluated usig relative

More information

Subject CT1 Financial Mathematics Core Technical Syllabus

Subject CT1 Financial Mathematics Core Technical Syllabus Subject CT1 Fiacial Mathematics Core Techical Syllabus for the 2018 exams 1 Jue 2017 Subject CT1 Fiacial Mathematics Core Techical Aim The aim of the Fiacial Mathematics subject is to provide a groudig

More information

AMS Portfolio Theory and Capital Markets

AMS Portfolio Theory and Capital Markets AMS 69.0 - Portfolio Theory ad Capital Markets I Class 6 - Asset yamics Robert J. Frey Research Professor Stoy Brook iversity, Applied Mathematics ad Statistics frey@ams.suysb.edu http://www.ams.suysb.edu/~frey/

More information

Standard Deviations for Normal Sampling Distributions are: For proportions For means _

Standard Deviations for Normal Sampling Distributions are: For proportions For means _ Sectio 9.2 Cofidece Itervals for Proportios We will lear to use a sample to say somethig about the world at large. This process (statistical iferece) is based o our uderstadig of samplig models, ad will

More information

Estimating Proportions with Confidence

Estimating Proportions with Confidence Aoucemets: Discussio today is review for midterm, o credit. You may atted more tha oe discussio sectio. Brig sheets of otes ad calculator to midterm. We will provide Scatro form. Homework: (Due Wed Chapter

More information

Monetary Economics: Problem Set #5 Solutions

Monetary Economics: Problem Set #5 Solutions Moetary Ecoomics oblem Set #5 Moetary Ecoomics: oblem Set #5 Solutios This problem set is marked out of 1 poits. The weight give to each part is idicated below. Please cotact me asap if you have ay questios.

More information

Maximum Empirical Likelihood Estimation (MELE)

Maximum Empirical Likelihood Estimation (MELE) Maximum Empirical Likelihood Estimatio (MELE Natha Smooha Abstract Estimatio of Stadard Liear Model - Maximum Empirical Likelihood Estimator: Combiatio of the idea of imum likelihood method of momets,

More information

. (The calculated sample mean is symbolized by x.)

. (The calculated sample mean is symbolized by x.) Stat 40, sectio 5.4 The Cetral Limit Theorem otes by Tim Pilachowski If you have t doe it yet, go to the Stat 40 page ad dowload the hadout 5.4 supplemet Cetral Limit Theorem. The homework (both practice

More information

EXERCISE - BINOMIAL THEOREM

EXERCISE - BINOMIAL THEOREM BINOMIAL THOEREM / EXERCISE - BINOMIAL THEOREM LEVEL I SUBJECTIVE QUESTIONS. Expad the followig expressios ad fid the umber of term i the expasio of the expressios. (a) (x + y) 99 (b) ( + a) 9 + ( a) 9

More information

An Empirical Study of the Behaviour of the Sample Kurtosis in Samples from Symmetric Stable Distributions

An Empirical Study of the Behaviour of the Sample Kurtosis in Samples from Symmetric Stable Distributions A Empirical Study of the Behaviour of the Sample Kurtosis i Samples from Symmetric Stable Distributios J. Marti va Zyl Departmet of Actuarial Sciece ad Mathematical Statistics, Uiversity of the Free State,

More information

Subject CT5 Contingencies Core Technical. Syllabus. for the 2011 Examinations. The Faculty of Actuaries and Institute of Actuaries.

Subject CT5 Contingencies Core Technical. Syllabus. for the 2011 Examinations. The Faculty of Actuaries and Institute of Actuaries. Subject CT5 Cotigecies Core Techical Syllabus for the 2011 Examiatios 1 Jue 2010 The Faculty of Actuaries ad Istitute of Actuaries Aim The aim of the Cotigecies subject is to provide a groudig i the mathematical

More information

Parametric Density Estimation: Maximum Likelihood Estimation

Parametric Density Estimation: Maximum Likelihood Estimation Parametric Desity stimatio: Maimum Likelihood stimatio C6 Today Itroductio to desity estimatio Maimum Likelihood stimatio Itroducto Bayesia Decisio Theory i previous lectures tells us how to desig a optimal

More information

1 Random Variables and Key Statistics

1 Random Variables and Key Statistics Review of Statistics 1 Radom Variables ad Key Statistics Radom Variable: A radom variable is a variable that takes o differet umerical values from a sample space determied by chace (probability distributio,

More information

Online appendices from Counterparty Risk and Credit Value Adjustment a continuing challenge for global financial markets by Jon Gregory

Online appendices from Counterparty Risk and Credit Value Adjustment a continuing challenge for global financial markets by Jon Gregory Olie appedices from Couterparty Risk ad Credit Value Adjustmet a APPENDIX 8A: Formulas for EE, PFE ad EPE for a ormal distributio Cosider a ormal distributio with mea (expected future value) ad stadard

More information

The Limit of a Sequence (Brief Summary) 1

The Limit of a Sequence (Brief Summary) 1 The Limit of a Sequece (Brief Summary). Defiitio. A real umber L is a it of a sequece of real umbers if every ope iterval cotaiig L cotais all but a fiite umber of terms of the sequece. 2. Claim. A sequece

More information

Hopscotch and Explicit difference method for solving Black-Scholes PDE

Hopscotch and Explicit difference method for solving Black-Scholes PDE Mälardale iversity Fiacial Egieerig Program Aalytical Fiace Semiar Report Hopscotch ad Explicit differece method for solvig Blac-Scholes PDE Istructor: Ja Röma Team members: A Gog HaiLog Zhao Hog Cui 0

More information

Stochastic Processes and their Applications in Financial Pricing

Stochastic Processes and their Applications in Financial Pricing Stochastic Processes ad their Applicatios i Fiacial Pricig Adrew Shi Jue 3, 1 Cotets 1 Itroductio Termiology.1 Fiacial.............................................. Stochastics............................................

More information

Anomaly Correction by Optimal Trading Frequency

Anomaly Correction by Optimal Trading Frequency Aomaly Correctio by Optimal Tradig Frequecy Yiqiao Yi Columbia Uiversity September 9, 206 Abstract Uder the assumptio that security prices follow radom walk, we look at price versus differet movig averages.

More information

4.5 Generalized likelihood ratio test

4.5 Generalized likelihood ratio test 4.5 Geeralized likelihood ratio test A assumptio that is used i the Athlete Biological Passport is that haemoglobi varies equally i all athletes. We wish to test this assumptio o a sample of k athletes.

More information

Overlapping Generations

Overlapping Generations Eco. 53a all 996 C. Sims. troductio Overlappig Geeratios We wat to study how asset markets allow idividuals, motivated by the eed to provide icome for their retiremet years, to fiace capital accumulatio

More information

Summary. Recap. Last Lecture. .1 If you know MLE of θ, can you also know MLE of τ(θ) for any function τ?

Summary. Recap. Last Lecture. .1 If you know MLE of θ, can you also know MLE of τ(θ) for any function τ? Last Lecture Biostatistics 60 - Statistical Iferece Lecture Cramer-Rao Theorem Hyu Mi Kag February 9th, 03 If you kow MLE of, ca you also kow MLE of τ() for ay fuctio τ? What are plausible ways to compare

More information

Exam 1 Spring 2015 Statistics for Applications 3/5/2015

Exam 1 Spring 2015 Statistics for Applications 3/5/2015 8.443 Exam Sprig 05 Statistics for Applicatios 3/5/05. Log Normal Distributio: A radom variable X follows a Logormal(θ, σ ) distributio if l(x) follows a Normal(θ, σ ) distributio. For the ormal radom

More information

1 The Black-Scholes model

1 The Black-Scholes model The Blac-Scholes model. The model setup I the simplest versio of the Blac-Scholes model the are two assets: a ris-less asset ba accout or bod)withpriceprocessbt) at timet, adarisyasset stoc) withpriceprocess

More information

Bayes Estimator for Coefficient of Variation and Inverse Coefficient of Variation for the Normal Distribution

Bayes Estimator for Coefficient of Variation and Inverse Coefficient of Variation for the Normal Distribution Iteratioal Joural of Statistics ad Systems ISSN 0973-675 Volume, Number 4 (07, pp. 7-73 Research Idia Publicatios http://www.ripublicatio.com Bayes Estimator for Coefficiet of Variatio ad Iverse Coefficiet

More information

1. Suppose X is a variable that follows the normal distribution with known standard deviation σ = 0.3 but unknown mean µ.

1. Suppose X is a variable that follows the normal distribution with known standard deviation σ = 0.3 but unknown mean µ. Chapter 9 Exercises Suppose X is a variable that follows the ormal distributio with kow stadard deviatio σ = 03 but ukow mea µ (a) Costruct a 95% cofidece iterval for µ if a radom sample of = 6 observatios

More information

BOUNDS FOR TAIL PROBABILITIES OF MARTINGALES USING SKEWNESS AND KURTOSIS. January 2008

BOUNDS FOR TAIL PROBABILITIES OF MARTINGALES USING SKEWNESS AND KURTOSIS. January 2008 BOUNDS FOR TAIL PROBABILITIES OF MARTINGALES USING SKEWNESS AND KURTOSIS V. Betkus 1,2 ad T. Juškevičius 1 Jauary 2008 Abstract. Let M = X 1 + + X be a sum of idepedet radom variables such that X k 1,

More information

A Bayesian perspective on estimating mean, variance, and standard-deviation from data

A Bayesian perspective on estimating mean, variance, and standard-deviation from data Brigham Youg Uiversity BYU ScholarsArchive All Faculty Publicatios 006--05 A Bayesia perspective o estimatig mea, variace, ad stadard-deviatio from data Travis E. Oliphat Follow this ad additioal works

More information

ASYMPTOTIC MEAN SQUARE ERRORS OF VARIANCE ESTIMATORS FOR U-STATISTICS AND THEIR EDGEWORTH EXPANSIONS

ASYMPTOTIC MEAN SQUARE ERRORS OF VARIANCE ESTIMATORS FOR U-STATISTICS AND THEIR EDGEWORTH EXPANSIONS J. Japa Statist. Soc. Vol. 8 No. 1 1998 1 19 ASYMPTOTIC MEAN SQUARE ERRORS OF VARIANCE ESTIMATORS FOR U-STATISTICS AND THEIR EDGEWORTH EXPANSIONS Yoshihiko Maesoo* This paper studies variace estimators

More information

We analyze the computational problem of estimating financial risk in a nested simulation. In this approach,

We analyze the computational problem of estimating financial risk in a nested simulation. In this approach, MANAGEMENT SCIENCE Vol. 57, No. 6, Jue 2011, pp. 1172 1194 iss 0025-1909 eiss 1526-5501 11 5706 1172 doi 10.1287/msc.1110.1330 2011 INFORMS Efficiet Risk Estimatio via Nested Sequetial Simulatio Mark Broadie

More information

CHAPTER 8 Estimating with Confidence

CHAPTER 8 Estimating with Confidence CHAPTER 8 Estimatig with Cofidece 8.2 Estimatig a Populatio Proportio The Practice of Statistics, 5th Editio Stares, Tabor, Yates, Moore Bedford Freema Worth Publishers Estimatig a Populatio Proportio

More information

Institute of Actuaries of India Subject CT5 General Insurance, Life and Health Contingencies

Institute of Actuaries of India Subject CT5 General Insurance, Life and Health Contingencies Istitute of Actuaries of Idia Subject CT5 Geeral Isurace, Life ad Health Cotigecies For 2017 Examiatios Aim The aim of the Cotigecies subject is to provide a groudig i the mathematical techiques which

More information

Bootstrapping high-frequency jump tests

Bootstrapping high-frequency jump tests Bootstrappig high-frequecy jump tests Prosper Dovoo Departmet of Ecoomics, Cocordia Uiversity Sílvia Goçalves Departmet of Ecoomics, Uiversity of Wester Otario Ulrich Houyo CREATES, Departmet of Ecoomics

More information

arxiv: v3 [math.st] 3 May 2016

arxiv: v3 [math.st] 3 May 2016 Bump detectio i heterogeeous Gaussia regressio arxiv:54.739v3 [math.st] 3 May 6 Farida Eikeeva farida.eikeeva@math.uiv-poitiers.fr Laboratoire de Mathématiques et Applicatios, Uiversité de Poitiers, Frace

More information

AY Term 2 Mock Examination

AY Term 2 Mock Examination AY 206-7 Term 2 Mock Examiatio Date / Start Time Course Group Istructor 24 March 207 / 2 PM to 3:00 PM QF302 Ivestmet ad Fiacial Data Aalysis G Christopher Tig INSTRUCTIONS TO STUDENTS. This mock examiatio

More information

INTERVAL GAMES. and player 2 selects 1, then player 2 would give player 1 a payoff of, 1) = 0.

INTERVAL GAMES. and player 2 selects 1, then player 2 would give player 1 a payoff of, 1) = 0. INTERVAL GAMES ANTHONY MENDES Let I ad I 2 be itervals of real umbers. A iterval game is played i this way: player secretly selects x I ad player 2 secretly ad idepedetly selects y I 2. After x ad y are

More information

SUPPLEMENTAL MATERIAL

SUPPLEMENTAL MATERIAL A SULEMENTAL MATERIAL Theorem (Expert pseudo-regret upper boud. Let us cosider a istace of the I-SG problem ad apply the FL algorithm, where each possible profile A is a expert ad receives, at roud, a

More information

CHANGE POINT TREND ANALYSIS OF GNI PER CAPITA IN SELECTED EUROPEAN COUNTRIES AND ISRAEL

CHANGE POINT TREND ANALYSIS OF GNI PER CAPITA IN SELECTED EUROPEAN COUNTRIES AND ISRAEL The 9 th Iteratioal Days of Statistics ad Ecoomics, Prague, September 0-, 05 CHANGE POINT TREND ANALYSIS OF GNI PER CAPITA IN SELECTED EUROPEAN COUNTRIES AND ISRAEL Lia Alatawa Yossi Yacu Gregory Gurevich

More information

43. A 000 par value 5-year bod with 8.0% semiaual coupos was bought to yield 7.5% covertible semiaually. Determie the amout of premium amortized i the 6 th coupo paymet. (A).00 (B).08 (C).5 (D).5 (E).34

More information

Limits of sequences. Contents 1. Introduction 2 2. Some notation for sequences The behaviour of infinite sequences 3

Limits of sequences. Contents 1. Introduction 2 2. Some notation for sequences The behaviour of infinite sequences 3 Limits of sequeces I this uit, we recall what is meat by a simple sequece, ad itroduce ifiite sequeces. We explai what it meas for two sequeces to be the same, ad what is meat by the -th term of a sequece.

More information

EQUIVALENCE OF FLOATING AND FIXED STRIKE ASIAN AND LOOKBACK OPTIONS

EQUIVALENCE OF FLOATING AND FIXED STRIKE ASIAN AND LOOKBACK OPTIONS EQUIVALENCE OF FLOATING AND FIXED STIKE ASIAN AND LOOKBACK OPTIONS ENST EBELEIN AND ANTONIS PAPAPANTOLEON Abstract. We prove a symmetry relatioship betwee floatig-strike ad fixed-strike Asia optios for

More information

Department of Mathematics, S.R.K.R. Engineering College, Bhimavaram, A.P., India 2

Department of Mathematics, S.R.K.R. Engineering College, Bhimavaram, A.P., India 2 Skewess Corrected Cotrol charts for two Iverted Models R. Subba Rao* 1, Pushpa Latha Mamidi 2, M.S. Ravi Kumar 3 1 Departmet of Mathematics, S.R.K.R. Egieerig College, Bhimavaram, A.P., Idia 2 Departmet

More information

Kernel Density Estimation. Let X be a random variable with continuous distribution F (x) and density f(x) = d

Kernel Density Estimation. Let X be a random variable with continuous distribution F (x) and density f(x) = d Kerel Desity Estimatio Let X be a radom variable wit cotiuous distributio F (x) ad desity f(x) = d dx F (x). Te goal is to estimate f(x). Wile F (x) ca be estimated by te EDF ˆF (x), we caot set ˆf(x)

More information

Math 312, Intro. to Real Analysis: Homework #4 Solutions

Math 312, Intro. to Real Analysis: Homework #4 Solutions Math 3, Itro. to Real Aalysis: Homework #4 Solutios Stephe G. Simpso Moday, March, 009 The assigmet cosists of Exercises 0.6, 0.8, 0.0,.,.3,.6,.0,.,. i the Ross textbook. Each problem couts 0 poits. 0.6.

More information

The material in this chapter is motivated by Experiment 9.

The material in this chapter is motivated by Experiment 9. Chapter 5 Optimal Auctios The material i this chapter is motivated by Experimet 9. We wish to aalyze the decisio of a seller who sets a reserve price whe auctioig off a item to a group of bidders. We begi

More information

Probability and statistics

Probability and statistics 4 Probability ad statistics Basic deitios Statistics is a mathematical disciplie that allows us to uderstad pheomea shaped by may evets that we caot keep track of. Sice we miss iformatio to predict the

More information

Solutions to Problem Sheet 1

Solutions to Problem Sheet 1 Solutios to Problem Sheet ) Use Theorem.4 to prove that p log for all real x 3. This is a versio of Theorem.4 with the iteger N replaced by the real x. Hit Give x 3 let N = [x], the largest iteger x. The,

More information

Sampling Distributions & Estimators

Sampling Distributions & Estimators API-209 TF Sessio 2 Teddy Svoroos September 18, 2015 Samplig Distributios & Estimators I. Estimators The Importace of Samplig Radomly Three Properties of Estimators 1. Ubiased 2. Cosistet 3. Efficiet I

More information

Topic-7. Large Sample Estimation

Topic-7. Large Sample Estimation Topic-7 Large Sample Estimatio TYPES OF INFERENCE Ò Estimatio: É Estimatig or predictig the value of the parameter É What is (are) the most likely values of m or p? Ò Hypothesis Testig: É Decidig about

More information

r i = a i + b i f b i = Cov[r i, f] The only parameters to be estimated for this model are a i 's, b i 's, σe 2 i

r i = a i + b i f b i = Cov[r i, f] The only parameters to be estimated for this model are a i 's, b i 's, σe 2 i The iformatio required by the mea-variace approach is substatial whe the umber of assets is large; there are mea values, variaces, ad )/2 covariaces - a total of 2 + )/2 parameters. Sigle-factor model:

More information

Bootstrapping high-frequency jump tests

Bootstrapping high-frequency jump tests Bootstrappig high-frequecy jump tests Prosper Dovoo Departmet of Ecoomics, Cocordia Uiversity Sílvia Goçalves Departmet of Ecoomics, McGill Uiversity Ulrich Houyo Departmet of Ecoomics, Uiversity at Albay,

More information

Estimation of a noisy subordinated Brownian Motion via two-scale power variations

Estimation of a noisy subordinated Brownian Motion via two-scale power variations Estimatio of a oisy subordiated Browia Motio via two-scale power variatios Jose E. Figueroa-López ad Kiseop Lee Departmet of Statistics Purdue Uiversity W. Lafayette, IN 7907, USA figueroa@purdue.edu Departmet

More information

Fourier Transform in L p (R) Spaces, p 1

Fourier Transform in L p (R) Spaces, p 1 Ge. Math. Notes, Vol. 3, No., March 20, pp.4-25 ISSN 229-784; Copyright c ICSS Publicatio, 200 www.i-csrs.org Available free olie at http://www.gema.i Fourier Trasform i L p () Spaces, p Devedra Kumar

More information

Inferential Statistics and Probability a Holistic Approach. Inference Process. Inference Process. Chapter 8 Slides. Maurice Geraghty,

Inferential Statistics and Probability a Holistic Approach. Inference Process. Inference Process. Chapter 8 Slides. Maurice Geraghty, Iferetial Statistics ad Probability a Holistic Approach Chapter 8 Poit Estimatio ad Cofidece Itervals This Course Material by Maurice Geraghty is licesed uder a Creative Commos Attributio-ShareAlike 4.0

More information

A New Constructive Proof of Graham's Theorem and More New Classes of Functionally Complete Functions

A New Constructive Proof of Graham's Theorem and More New Classes of Functionally Complete Functions A New Costructive Proof of Graham's Theorem ad More New Classes of Fuctioally Complete Fuctios Azhou Yag, Ph.D. Zhu-qi Lu, Ph.D. Abstract A -valued two-variable truth fuctio is called fuctioally complete,

More information

Chapter 8: Estimation of Mean & Proportion. Introduction

Chapter 8: Estimation of Mean & Proportion. Introduction Chapter 8: Estimatio of Mea & Proportio 8.1 Estimatio, Poit Estimate, ad Iterval Estimate 8.2 Estimatio of a Populatio Mea: σ Kow 8.3 Estimatio of a Populatio Mea: σ Not Kow 8.4 Estimatio of a Populatio

More information

(Hypothetical) Negative Probabilities Can Speed Up Uncertainty Propagation Algorithms

(Hypothetical) Negative Probabilities Can Speed Up Uncertainty Propagation Algorithms Uiversity of Texas at El Paso DigitalCommos@UTEP Departmetal Techical Reports (CS) Departmet of Computer Sciece 2-2017 (Hypothetical) Negative Probabilities Ca Speed Up Ucertaity Propagatio Algorithms

More information

Binomial Model. Stock Price Dynamics. The Key Idea Riskless Hedge

Binomial Model. Stock Price Dynamics. The Key Idea Riskless Hedge Biomial Model Stock Price Dyamics The value of a optio at maturity depeds o the price of the uderlyig stock at maturity. The value of the optio today depeds o the expected value of the optio at maturity

More information

FOUNDATION ACTED COURSE (FAC)

FOUNDATION ACTED COURSE (FAC) FOUNDATION ACTED COURSE (FAC) What is the Foudatio ActEd Course (FAC)? FAC is desiged to help studets improve their mathematical skills i preparatio for the Core Techical subjects. It is a referece documet

More information

Model checks for the volatility under microstructure noise

Model checks for the volatility under microstructure noise Model checks for the volatility uder microstructure oise Mathias Vetter ad Holger Dette Ruhr-Uiversität Bochum Fakultät für Mathematik 4478 Bochum Germay email: mathias.vetter@rub.de; holger.dette@rub.de

More information

These characteristics are expressed in terms of statistical properties which are estimated from the sample data.

These characteristics are expressed in terms of statistical properties which are estimated from the sample data. 0. Key Statistical Measures of Data Four pricipal features which characterize a set of observatios o a radom variable are: (i) the cetral tedecy or the value aroud which all other values are buched, (ii)

More information

Control Charts for Mean under Shrinkage Technique

Control Charts for Mean under Shrinkage Technique Helderma Verlag Ecoomic Quality Cotrol ISSN 0940-5151 Vol 24 (2009), No. 2, 255 261 Cotrol Charts for Mea uder Shrikage Techique J. R. Sigh ad Mujahida Sayyed Abstract: I this paper a attempt is made to

More information

NOTES ON ESTIMATION AND CONFIDENCE INTERVALS. 1. Estimation

NOTES ON ESTIMATION AND CONFIDENCE INTERVALS. 1. Estimation NOTES ON ESTIMATION AND CONFIDENCE INTERVALS MICHAEL N. KATEHAKIS 1. Estimatio Estimatio is a brach of statistics that deals with estimatig the values of parameters of a uderlyig distributio based o observed/empirical

More information

MODIFICATION OF HOLT S MODEL EXEMPLIFIED BY THE TRANSPORT OF GOODS BY INLAND WATERWAYS TRANSPORT

MODIFICATION OF HOLT S MODEL EXEMPLIFIED BY THE TRANSPORT OF GOODS BY INLAND WATERWAYS TRANSPORT The publicatio appeared i Szoste R.: Modificatio of Holt s model exemplified by the trasport of goods by ilad waterways trasport, Publishig House of Rzeszow Uiversity of Techology No. 85, Maagemet ad Maretig

More information

Non-Inferiority Logrank Tests

Non-Inferiority Logrank Tests Chapter 706 No-Iferiority Lograk Tests Itroductio This module computes the sample size ad power for o-iferiority tests uder the assumptio of proportioal hazards. Accrual time ad follow-up time are icluded

More information

Basic formula for confidence intervals. Formulas for estimating population variance Normal Uniform Proportion

Basic formula for confidence intervals. Formulas for estimating population variance Normal Uniform Proportion Basic formula for the Chi-square test (Observed - Expected ) Expected Basic formula for cofidece itervals sˆ x ± Z ' Sample size adjustmet for fiite populatio (N * ) (N + - 1) Formulas for estimatig populatio

More information

Monopoly vs. Competition in Light of Extraction Norms. Abstract

Monopoly vs. Competition in Light of Extraction Norms. Abstract Moopoly vs. Competitio i Light of Extractio Norms By Arkadi Koziashvili, Shmuel Nitza ad Yossef Tobol Abstract This ote demostrates that whether the market is competitive or moopolistic eed ot be the result

More information

A DOUBLE INCREMENTAL AGGREGATED GRADIENT METHOD WITH LINEAR CONVERGENCE RATE FOR LARGE-SCALE OPTIMIZATION

A DOUBLE INCREMENTAL AGGREGATED GRADIENT METHOD WITH LINEAR CONVERGENCE RATE FOR LARGE-SCALE OPTIMIZATION A DOUBLE INCREMENTAL AGGREGATED GRADIENT METHOD WITH LINEAR CONVERGENCE RATE FOR LARGE-SCALE OPTIMIZATION Arya Mokhtari, Mert Gürbüzbalaba, ad Alejadro Ribeiro Departmet of Electrical ad Systems Egieerig,

More information

Simulation Efficiency and an Introduction to Variance Reduction Methods

Simulation Efficiency and an Introduction to Variance Reduction Methods Mote Carlo Simulatio: IEOR E4703 Columbia Uiversity c 2017 by Marti Haugh Simulatio Efficiecy ad a Itroductio to Variace Reductio Methods I these otes we discuss the efficiecy of a Mote-Carlo estimator.

More information

Estimation of generalized Pareto distribution

Estimation of generalized Pareto distribution Estimatio of geeralized Pareto distributio Joa Del Castillo, Jalila Daoudi To cite this versio: Joa Del Castillo, Jalila Daoudi. Estimatio of geeralized Pareto distributio. Statistics ad Probability Letters,

More information

Estimation of integrated volatility of volatility with applications to goodness-of-fit testing

Estimation of integrated volatility of volatility with applications to goodness-of-fit testing Beroulli 21(4, 215, 2393 2418 DOI: 1.315/14-BEJ648 arxiv:126.5761v2 [math.st] 29 Sep 215 Estimatio of itegrated volatility of volatility with applicatios to goodess-of-fit testig MATHIAS VETTER 1 Fakultät

More information

Twitter: @Owe134866 www.mathsfreeresourcelibrary.com Prior Kowledge Check 1) State whether each variable is qualitative or quatitative: a) Car colour Qualitative b) Miles travelled by a cyclist c) Favourite

More information

We learned: $100 cash today is preferred over $100 a year from now

We learned: $100 cash today is preferred over $100 a year from now Recap from Last Week Time Value of Moey We leared: $ cash today is preferred over $ a year from ow there is time value of moey i the form of willigess of baks, busiesses, ad people to pay iterest for its

More information

Today: Finish Chapter 9 (Sections 9.6 to 9.8 and 9.9 Lesson 3)

Today: Finish Chapter 9 (Sections 9.6 to 9.8 and 9.9 Lesson 3) Today: Fiish Chapter 9 (Sectios 9.6 to 9.8 ad 9.9 Lesso 3) ANNOUNCEMENTS: Quiz #7 begis after class today, eds Moday at 3pm. Quiz #8 will begi ext Friday ad ed at 10am Moday (day of fial). There will be

More information

REVISIT OF STOCHASTIC MESH METHOD FOR PRICING AMERICAN OPTIONS. Guangwu Liu L. Jeff Hong

REVISIT OF STOCHASTIC MESH METHOD FOR PRICING AMERICAN OPTIONS. Guangwu Liu L. Jeff Hong Proceedigs of the 2008 Witer Simulatio Coferece S. J. Maso, R. R. Hill, L. Möch, O. Rose, T. Jefferso, J. W. Fowler eds. REVISIT OF STOCHASTIC MESH METHOD FOR PRICING AMERICAN OPTIONS Guagwu Liu L. Jeff

More information