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

Size: px
Start display at page:

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

Transcription

1 Beroulli 21(4, 215, DOI: 1.315/14-BEJ648 arxiv: v2 [math.st] 29 Sep 215 Estimatio of itegrated volatility of volatility with applicatios to goodess-of-fit testig MATHIAS VETTER 1 Fakultät für Mathematik, Ruhr-Uiversität Bochum, 4478 Bochum, Germay. mathias.vetter@rub.de I this paper, we are cocered with oparametric iferece o the volatility of volatility process i stochastic volatility models. We costruct several estimators for its itegrated versio i a high-frequecy settig, all based o icremets of spot volatility estimators. Some of those are positive by costructio, others are bias corrected i order to attai the optimal rate 1/4. Associated cetral limit theorems are prove which ca be widely used i practice, as they are the key to essetially all tools i model validatio for stochastic volatility models. As a illustratio we give a brief idea o a goodess-of-fit test i order to check for a certai parametric form of volatility of volatility. Keywords: cetral limit theorem; goodess-of-fit testig; high-frequecy observatios; model validatio; stable covergece; stochastic volatility model 1. Itroductio Nowadays, stochastic volatility models are stadard tools i the cotiuous-time modellig of fiacial time series. Typically, the uderlyig (log price process is assumed to follow a diffusio process of the form X t =X + µ s ds+ σ s dw s, (1.1 where µ ad σ ca be quite geeral stochastic processes themselves. A classical case is where the volatility σ 2 s =σ 2 (s,x s is a fuctio of time ad state a situatio referred to as the oe of a local volatility model. It has tured out i empirical fiace that such models do ot fit the data very well, as some stylised facts such as the leverage effect or volatility clusterig caot be explaied usig local volatility oly. Stochastic volatility models, however, are able to reproduce such features, as they bear a additioal source of radomess. I these models, the volatility process is a diffusio process itself, ad we This is a electroic reprit of the origial article published by the ISI/BS i Beroulli, 215, Vol. 21, No. 4, This reprit differs from the origial i pagiatio ad typographic detail c 215 ISI/BS

2 2 M. Vetter focus o a rather geeral situatio, amely σt 2 =σ2 + ν s ds+ β s dw s + η s dw s, (1.2 where ν, β ad η agai are suitable stochastic processes ad W is aother Browia motio, idepedet of W. This model obviously icludes the widely used special case of a volatility with oly oe drivig Browia motio, which is dσ 2 t =ν tdt+τ t dv t, where V ad W are joitly Browia with some correlatio ρ. Stochastic volatility models are typically parametric oes, ad probably the prime example amog those is the Hesto model of [14], give by X t =X + ( β σ2 s 2 ds+ σ s dw s, σt 2 =σ2 +κ (α σs 2 ds+ξ σ s dv s, for some parameters β,κ,α ad ξ, ad with Corr(W,V=ρ. Here, the volatility process follows a Cox Igersoll Ross model, that meas it is mea-revertig with mea α ad speed κ, ad both diffusio coefficiets are proportioal with parameter ξ. Particularly the latter property appears to be rather typical for stochastic volatility models, ad i this sese the Hesto model ca be regarded as prototypic. Popular alteratives are for example comig from the more geeral (but agai parametric class of (oe factor CEV models, where the diffusio coefficiet τ of σ 2 becomes a geeral power fuctio of σ, whereasthe drift part ofthe volatility remaisi priciple the same. See [21] for a survey. For this reaso, statistical iferece for stochastic volatility models has focused o parametric methods for most times, ad usually the authors provide tools for a specific class of models. However, oe is faced with two severe problems: First, it is i most cases impossible to assess the distributio of X (or its icremets, which makes stadard maximum likelihood theory uavailable. Secod, the volatility process σ 2 is ot observable, ad may statistical cocepts have i commo that they propose to reproduce the ukow volatility process from observed optio prices, typically by usig proxies based o implied volatility. A survey o early estimatio methods i this cotext ca be foud i [8]. Oe remarkable exceptio where stock price data oly is used is the paper of [7] who costruct a GMM estimator for the parameters of the Hesto model from icremets of realised variace. But also i a geeral settig with o specific model i mid, the focus has bee o parametric approaches. A early approach o parameter estimatio whe σ 2 is ergodic is the work of [12], optimal rates are discussed i [15] ad [13], ad a maximum likelihood approach based o proxies for the volatility ca be foud i[1]. Eve oparametric cocepts have bee used to idetify parameters of a stochastic volatility model; see, for example, [3] or [25]. Geuie oparametric iferece for stochastic volatility models has typically focused o fuctio estimatio. Both [24] ad [9] discuss techiques for the estimatio of f ad g,

3 Itegrated volatility of volatility 3 whethevolatilityprocesssatisfiesdσt 2 =f(σ2 t dt+g(σ2 t dv t.ithemoregeeralmodelfree cotext of(1.2, oly [4] ad [28] have discussed estimatio of fuctioals of volatility of volatility. While the latter focus o estimatio of a kid of leverage effect which ivolves the volatility of volatility process(es, the work of [4] provides a cosistet estimator for itegratedvolatilityofvolatility τ2 s ds i theoe-factorcase.theirapproachisispired by the asymptotic behaviour of realised variace, which states that the sum of squared icremets of σ 2 coverges i probability to the quatity of iterest. Sice σ 2 is ot observable, the authors use spot volatility estimators istead. We will pursue their approach ad discuss i detail the asymptotic behaviour of several estimators for itegrated volatility of volatility, all based o icremets of spot volatility estimators, thus usig observatios of X oly. It turs out that i order to attai the optimal rate of covergece i this cotext, it is ecessary to coduct a certai bias correctio which destroys positivity of the estimator a feature which is well kow from the related problem of volatility estimatio uder microstructure oise. Several stable cetral limit theorems are provided, ad by defiig appropriate estimators for the asymptotic (coditioal variace we obtai feasible versios as well. The latter results are of theoretical iterest o oe had, but are extremely importat from a applied poit of view as well, as they make model validatio for stochastic volatility models possible. Give the tremedous umber of such models with etirely differet qualitative behaviours, there is a lack of techiques that help decidig whether a certai model fits the data appropriately or ot. As a first approach to model validatio i this framework, we give a brief idea o how to do goodess-of-fit testig, but our method is by o meas limited to it. Related procedures ca be used to test for example, whether a Browia compoet or jumps are preset i the volatility process ad what i geeral the structure of the jump part is. Such problems have bee solved for the price process X i recet years (see [18] for a overview, ad i priciple the methods are all based o the estimatio of plai itegrated volatility σ2 sds ad further quatities, such as trucated versiosor bipower variatio. Usig our mai results, these cocepts ca be traslated to the stochastic volatility case by usig estimators for itegrated volatility of volatility istead, but usually with the slower rate of covergece 1/4. The paper is orgaised as follows: I Sectio 2, we itroduce our estimators ad state the cetral limit theorems, whereas Sectio 3 is o goodess-of-fit testig i stochastic volatility models. Some Mote Carlo results ca be foud i Sectio 4, followed by some cocludig remarks i Sectio 5. A overview o some proofs plus a couple of details ca be foud i the Appedix, whereas large parts of them have bee relegated to a supplemetary article [26]. 2. Mai results Let us start with some coditios o the processes ivolved. All of these are rather mild ad covered by a variety of(stochastic volatility models used. The oly major restrictio is that we will assume most processesto be cotiuous for a while ad oly discuss briefly

4 4 M. Vetter later how possible adjustmets i order to hadle jumps i price ad volatility could look like. Assumptio 2.1. Suppose that the process X is give by (1.1, where W is a stadard Browia motio ad the drift process µ is left cotiuous. We assume further that the volatility process σ 2 is a cotiuous Itô semimartigale itself, havig the represetatio (1.2. ν is assumed to be left cotiuous as well, whereas β satisfies the regularity coditio βs 2 =β2 + ω s ds+ ϑ (1 s dw s + ϑ (2 s dw s, (2.1 where ω is locally bouded ad each ϑ (l is left cotiuous,,2. A similar coditio is assumed to hold for η as well. Fially, all processes are defied o the same probability space (Ω,F,(F t t,p, ad all coefficiets are specified i such a way that σ 2 is almost surely positive ad that β 2 ad η 2 are either almost surely positive or vaishig idetically, respectively. As oted i the Itroductio, (1.2 covers a large class of volatility models used. For η we are essetially i the case of a local volatility model, whereas β s = ρτ s ad η s = 1 ρ 2 τ s for some process τ ad ρ ( 1,1 refers to the settig of the typical stochastic volatility models metioed before, i which both drivig Browia motios arecorrelatedwith ρ.themodeli(1.2isevemoreflexible,aditisstraight-forwardto extedallresultstothecaseofamulti-factormodeldrivebymorethatwoidepedet Browia motios as well. Our aim i the followig is to draw iferece o the itegrated volatility of volatility up to time t, which becomes (β2 s +η2 s ds i our cotext. Ay statistical iferece will be based o high-frequecy observatios of X, ad we assume that the data is recorded at equidistat times. Thus, without loss of geerality let the process be defied o the iterval [,1] ad observed at the time poits i/, i=,...,. Before we discuss several cocepts to assess itegrated volatility of volatility i detail, let us recall the priciples of estimatio of stadard itegrated volatility 1 σ2 sds. The usual estimator i the geeral model-free settig of (1.1 is realised volatility, give by RV t t = i=1 ( i X2, where we set i Z = Z i/ Z (i 1/ for ay process Z. This estimator is optimal i several respects, eve though Itô formula proves i/ i/ ( i X 2 = σsds+2 2 (X s X (i 1/ dx s (2.2 (i 1/ (i 1/ oly,fromwhichitissimpletoseethateachsquaredicremet( i X2 isolyoaverage equal to itegrated volatility over the correspodig time iterval, but ot cosistet for

5 Itegrated volatility of volatility 5 it. (Realised volatility, the sum of the squared icremets, however, is cosistet for the etire itegrated volatility, which is basically due to a martigale argumet. Our estimators for itegrated volatility of volatility will be based o a similar ituitio: Defie statistics via sums of icremets such that each summad is o average equal to itegrated volatility of volatility over the correspodig time iterval, but ot ecessarily cosistet. As before, oe would like to build those estimators upo icremets of σ 2. These are i geeral ot observable, so a proxy for them is eeded. Sice we are i a model-free world, a atural estimator for spot volatility σi/ 2 is give by ˆσ i/ 2 = k ( i+j k X2, i=,..., k, j=1 for some auxiliary (iteger-valued sequece k. See [2] or [25] for details o the asymptotic behaviour of this estimator. Itô formula agai gives ˆσ i/ 2 = k (i+j/ 2 (X s X (i+j 1/ dx s + (i+k/ σs 2 k (i+j 1/ k ds=:a i +B i, (2.3 j=1 so that ˆσ i/ 2 σ2 i/ cosists of two sources of error. From the proofs later o, we see that A i =O p( 1/k, whereas Bi σ2 i/ =O p( k /. Therefore, it appears atural to choose k to be of the order 1/2 i order to miimize the error of the spot volatility estimator (ad we will see later that this is ideed the best thig to do, but we will keep this sequece arbitrary i order to allow for other estimators as well. While the choice of the spot volatility estimators depeds o the auxiliary sequece k, we will itroduce a secod sequece of itegers l which govers the legth of the itervals over which icremets of ˆσ 2 are computed. Thus, the basic elemet of our fial estimators will be (ˆσ 2 (i+l / ˆσ2 i/ 2, which ca be decomposed as (ˆσ 2 (i+l / ˆσ2 i/ 2 =(A i+l A i 2 +(B i+l B i 2 +2(A i+l A i (B i+l B i. The average behaviour of the terms above is discussed i the followig lemma, ad it depeds crucially o the size of both k ad l. Lemma 2.2. Suppose that Assumptio 2.1 holds ad let E i [Z] deote coditioal expectatio of some variable Z with respect to F i/. Set also M = max(k,l ad m =mi(k,l. The we have E i [(A i+l A i 2 ]=4l (k M 1 σ 4 i/ (1+O p(m 1/2 1/2, E i [(B i+l B i 2 ]=l m (M m /3(k M 1 (β 2 i/ +η2 i/ (1+O p(m 1/2 1/2. The previous lemma gives us several hits o how to obtai a estimator for itegrated volatility of volatility via sums over (ˆσ 2 (i+l / ˆσ2 i/ 2. First, iformatio about β 2 i/ + i/

6 6 M. Vetter ηi/ 2 is cotaied i icremets over the B i oly. Therefore, it appears to be reasoable to choose k ad l later o i such a way that these terms are at least ot smaller tha the bias terms due to icremets of A i. Or i other words, the coditio becomes that Ck l for some geeric C >. Also, there are basically two ways to costruct a estimator. Either, pick k ad l such that the bias due to icremets of A i is egligible eve after dividig by the rate of covergece. This cocept will lead to the estimator t (k +l ˆT t = k M (l m (M m /3 1 (ˆσ (i+l 2 / ˆσ2 i/ 2 i= which is positive by costructio. As oted i the Itroductio, this is the kid of estimator [4] were lookig at. Alteratively, oe ca use a bias correctio ad subtract a estimator for the local quarticity σi/ 4. I this case oe loses positivity, but we will see later that the rate of covergece is much faster i this situatio. Let us pursue the first path for a momet, however. I order to uderstad what the rate of covergece for estimatio of itegrated volatility of volatility will be, the ext result is extremely helpful, as it gives the cetral limit theorem for the oracle estimator t (k +l Ŝt = k M (l m (M m /3 1 (Bi+l Bi 2 i= which depeds o the uobservable icremets of B i oly. All results i this sectio will be poitwise i t, eve though it is likely that fuctioal versios hold as well. Propositio 2.3. Suppose that Assumptio 2.1 holds ad that both k c α ad l d β hold for some α,β (,1 ad c,d>. Let also M ad m be defied as before. (a If α β, we have M ( Ŝ t (b For k =l we have M ( Ŝ t (βs 2 +ηsds 2 L (s 4/3 (βs 2 +ηsds 2 L (s 151/7 (β 2 s +η 2 sdw s. (β 2 s +η 2 sdw s. I both cases, W is a Browia motio defied o a extesio of the origial probability space ad idepedet of F ad the covergece i (2.6 is F-stable i law. Remark 2.4. It is obvious from Propositio 2.3 that the rate of covergece becomes faster the smaller M is chose. O the other had, the coditio Ck l forces M to be at least of the order 1/2. I this case, the rate of covergece i Propositio 2.3

7 Itegrated volatility of volatility 7 becomes 1/4, ad this rate is kow to be optimal for this statistical problem. Ideed, a related parametric settig has bee discussed i [15] a decade ago, ad it was show therei that this rate is optimal i the special case, where β vaishes idetically ad η is a fuctio of time ad state, kow up to a parameter θ. Our first mai theorem specifies coditios for a cetral limit theorem for ˆT t ad is a simple cosequece of Lemma 2.2 ad Propositio 2.3. Theorem 2.5. Suppose that all the assumptios of Propositio 2.3 hold true. If further 3/2 M 3/2 m 1 ad α β, the the stable cetral limit theorem t (ˆT t (βs 2 +η2 s ds L (s 4/3 (βs 2 +η2 s dw s. (2.4 holds true. M The optimal rate of covergece i this case is obtaied for the choice of M = O( 3/5+ε ad m = O( 3/5 ad approaches 1/5 for ε. This proves also that it is o restrictio to assume α β above. I order to obtai a estimator with the optimal rate of covergece, we choose l ad k to be both the same ad of the order 1/2, but as oted above we eed a bias correctio the. Therefore, we defie with a slight abuse of otatio ˆR t = t 2k i= ( 3 (ˆσ (i+k 2 2k / ˆσ2 i/ k 2 ˆσ i/ 4, (2.5 where ˆσ 4 i/ = 2 3k k j=1 i+j X 4 is i geeral differet from (ˆσ 2 i/ 2. Its asymptotic behaviour is discussed i the followig theorem. Theorem 2.6. Suppose that Assumptio 2.1 holds ad let k =c 1/2 +o( 1/4 for some c>. The ( t ˆR t (βs 2 +ηsds 2 L (s U t (2.6 k for all t>, where the limitig variable has the represetatio U t = α s dw s, α 2 s = 48 c 4σ8 s + 12 c 2σ4 s(β 2 s +η 2 s (β2 s +η 2 s 2. (2.7 Remark 2.7. The situatio ecoutered above has a iterestig coectio to the problem of elimiatig microstructure oise, as we face similar problems regardig optimal rates of covergece ad positivity of the estimators. Whereas the optimal rate of covergece for estimatig itegrated volatility i the oisy settig is 1/4, stadard estimators attaiig this rate are ot always positive. To esure positivity, oe typically

8 8 M. Vetter accepts a drop i the rate of covergece to 1/5 as well. See, for example, [6] for a thorough discussio i a geeral multivariate settig. Remark 2.8. Recetly, [19] discussed efficiet estimatio of g(σ2 s ds for geeralfuctios g. It tured out that Riema sums based o g(ˆσ i/ 2 ideed attai the optimal rate of covergece 1/2 i this cotext, but agai the choice of k affects the limitig distributio. The optimal k 1/2 leads to additioal bias terms i their settig, ad at least some of these ca be avoided by choosig k i a differet way. The limitig distributio i Theorem 2.5 ad Theorem 2.6 is mixed ormal, ad i order to obtai a feasible cetral limit theorem we have to itroduce cosistet estimators for the respective coditioal variaces. These are costructed usig the same ituitio as before, ad precisely we obtai the followig theorem. Theorem 2.9. ( a Uder the coditios of Theorem 2.5, we have t (k +l ˆQ 4k 2 t = M 2 9(l m (M m /3 2(ˆσ2 (i+l / ˆσ2 i/ 4 P 4 3 (β2 s +ηs 2 2 ds. i= (b I the situatio of Theorem 2.6, we have G (1 t, = 1 t k (ˆσ i/ 4 2 P i=1 t 2k G (2 t, = i=1 t 2k G (3 t, = i=1 σ 8 sds, ( 3 (ˆσ (i+k 2 2k / ˆσ2 i/ k 2 ˆσ i/ 4 ˆσ i/ 4 k 2 Therefore, as a cosequece (ˆσ 2 (i+k / ˆσ2 i/ 4 P P ˆP t = G(3 t, 486 k 2 35 G(2 t, k 4 35 G(1 P t, σ 4 s (β2 s +η2 s ds, ( 48 c 4σ8 s + 16 c 2σ4 s (β2 s +η2 s (β2 s +η2 s 2 ds. α 2 sds. Remark 2.1. Theorem 2.9 shows that a cosistet estimator for (β2 s +η 2 s 2 ds is for example, give by 3 4 G(3 t, 12 k 2 G (2 t, 362 k 4 G (1 t,, ad its proof suggests that a cetral limit theorem holds with the same rate of covergece as before. I geeral, it is quite likely that this methods provides estimates for arbitrary eve powers of itegrated volatility of volatility. A precise theory is left for future research.

9 Itegrated volatility of volatility 9 The properties of stable covergece guaratee that dividig by the square root of a cosistet estimator for the coditioal variace gives a feasible cetral limit theorem for the estimatio of itegrated volatility of volatility. See, for example, [23] for details. Therefore, the followig corollary ca be cocluded easily. Corollary (a Uder the assumptios of Theorem 2.5, we have for all t> M (ˆT t (β 2 s +η2 s ds (ˆQ t 1/2 L N(,1. (2.8 (b Uder the assumptios of Theorem 2.6, we have for all t> k ( ˆR t (βs 2 +ηsds 2 (ˆP t 1/2 N(,1. L (2.9 Remark So far we have oly discussed the case where both processes have cotiuous paths. Extesios to the situatio of additioal jumps i the price process seem to be possible, but are already quite ivolved. The followig observatio is useful: Wheever there is a jump withi the iterval [i/,(i+2k /], it appears squared ad blow up by 2 /k 2 withi (ˆσ2 (i+k / ˆσ2 i/ 2. This is a much larger order tha the usual k / i the cotiuous case. For this reaso,it appearsas if the trucatio method due to [22] ca be applied, ad a similar ituitio holds for the bias correctio as well. Note, however, that the raw statistics i this cotext are sums of squared icremets of X rather tha plai icremets of X as for the power variatios ecoutered i [22]. Therefore, the required techiques are differet tha the stadard oes i this area. The case of jumps i the volatility appears to be eve more complicated, as these come ito play via Bi+k Bi = (i+k/ (σs+k 2 k / σ2 s ds, i/ ad therefore the amout to which each icremet is affected by a jump depeds crucially o the time at which the jump occurs. Thus, plai trucatio might ot be sufficiet i this case ad a etirely differet estimator was ecessary. Both topics are left for future research. 3. Model checks for stochastic volatility models I this sectio, we propose a first approach to goodess-of-fit testig for stochastic volatility models. Assume we have represetatio (1.1 for the log price process X, whereas the volatility process satisfies dσ 2 t =ν t dt+τ t dv t as i typical SV models. There is still a lot of freedom i the modellig of σ 2, ad the various proposals i the literature typically differ i the represetatio of its diffusio part τ. As oted i the Itroductio, a quite geeral class of stochastic volatility models is give by the so-called CEV models, i

10 1 M. Vetter which τ 2 s =θ(σ2 s γ for some oegative γ ad a ukow parameter θ, ad the most popular amog these is the Hesto model from [14], correspodig to γ =1. I order to costruct a test whether a certai fuctioal relatioship betwee σ ad τ ispreset,weemployatechiquewhich wasalreadyusedi [1] or[27] whe dealigwith localvolatilitymodels.supposeweareiteresteditestigforτ 2 s =τ 2 (s,x s,σ 2 s,θ,where τ 2 is a give fuctio ad θ is some ukow (i geeral multidimesioal parameter. For simplicity, we will focus o the oe-dimesioal liear case oly, that is H :τ 2 s =θτ2 (s,x s,σ 2 s for all s [,1] (a.s. Extesios to the geeral case follow alog the lies of Sectio 5 i [27]. A test for the ull hypothesis will be based o the observatio that H is equivalet to N t = for all t [,1] (a.s., where the process N t is give by N t = θ mi =argmi θ (τ 2 s θ mi τ 2 (s,x s,σ 2 sds, 1 (τ 2 s θτ2 (s,x s,σ 2 s 2 ds. Assume that the fuctio τ 2 is bouded awayfrom zero. The a stadard argumet from Hilbert space theory shows that θ mi =D 1 C (ad therefore N t =R t B t D 1 C, where we have set R t = τ2 s ds ad B t = D = C = 1 1 τ 2 (s,x s,σ 2 s ds, τ 4 (s,x s,σ 2 sds, τ 2 sτ 2 (s,x s,σ 2 sds. To defie estimators let k as before ad recall (2.5. We set ˆτ 2 i/ =3(2k 1 (ˆσ 2 (i+k / ˆσ2 i/ 2 6k 2 ˆσ4 i/ (3.1 ad also ˆN t = ˆR t ˆB t (ˆD 1 Ĉ with ˆR t from the previous sectio, whereas we deote ˆB t = 1 ˆD = 1 t k i= k i= τ 2 ( i,x i/,ˆσ 2 i/ τ 4 ( i,x i/,ˆσ 2 i/,,

11 Itegrated volatility of volatility 11 Ĉ = 1 2k i= ˆτ 2 i/ τ2 ( i,x i/,ˆσ 2 i/ I the sequel, we will prove weak covergeceof ˆN t N t, up to a suitable ormalisatio. Theorem 2.6 suggests that /k is a reasoable choice, ad the followig claim proves that two of the estimators coverge at a faster speed, at least if we impose a additioal smoothess coditio o the fuctio τ 2. Lemma 3.1. Suppose that the fuctio τ 2 has cotiuous partial derivatives of secod order. The we have. ˆB t B t =o p ( 1/4, ˆD D=o p ( 1/4, the first result holdig uiformly i t [,1]. The above claim idicates that we have to focus o the terms ivolvig ˆτ 2 i/ oly, which is familiar groud due to the results of Sectio 2. We start with a propositio o the joit asymptotic behaviour of ˆR t ad Ĉ. Lemma 3.2. Let d be a iteger ad t 1,...,t d be arbitrary i [,1]. Set Σ t1,...,t d (s,x s,σ 2 s=α 2 sh t1,...,t d (s,x s,σ 2 sh t1,...,t d (s,x s,σ 2 s T with h t1,...,t d (s,x s,σs 2=(1 [,t 1],...,1 [,td ],τ 2 (s,x s,σs 2T ad α 2 s as i Theorem 2.6. Uder the previous assumptios we have the stable covergece (ˆR t k 1 R t1,..., ˆR t d R td,ĉ C T L (s 1 Σ 1/2 t 1,...,t d (s,x s,σ 2 sdw s, where W is a (d + 1-dimesioal stadard Browia motio defied o a extesio of the origial space ad idepedet of F. We are iterested i the asymptotics of the process A (t= /k ( ˆN t N t, ad the precedig lemma basically leads to its fiite dimesioal covergece. The etire result o weak covergece of A reads as follows. Theorem 3.3. Assume that the previous assumptios hold. The the process (A (t t [,1] coverges weakly to a mea zero process (A(t t [,1], which is Gaussia coditioally o F ad whose coditioal covariace equals the oe of the process where U U[,1], idepedet of F. {α U (1 {U t} B t D 1 τ 2 (U,X U,σ 2 U} t [,1]

12 12 M. Vetter As idicated before, covergece of the fiite dimesioal distributios is a direct cosequece of Lemma 3.2, usig the Delta method for stable covergece (see, e.g., [11]. Tightess follows from Theorem VI. 4.5 i [2] with a miimal amout of work. Recall that N t = for all t uder the ull hypothesis. Therefore Theorem 3.3 shows that a cosistet test is obtaied by rejectig the ull hypothesis for large values of a suitable fuctioal of the process { /k ˆN t } t [,1]. If we choose the Kolmogorov Smirov fuctioal K =sup t [,1] /k ˆN t for example, we have weak covergece uder the ull to sup t [,1] A t as a cosequece of Theorem 3.3. The distributio of the latter statistic is extremely difficult to assess, as it typically depeds o the etire process (X,σ 2. We therefore propose to obtai critical values via a simple bootstrap procedure, which will be itroduced i the ext sectio. To ed this sectio, we defie a appropriate estimator for the coditioal variace of A(t, which is give by s 2 t = α 2 s ds 2B td 1 α 2 s τ2 (s,x s,σs 2 ds+b2 t D 2 α 2 s τ4 (s,x s,σs 2 ds, due to Theorem 3.3. Empirical couterparts for B t ad D are obviously defied by the statistics ˆB t ad ˆD, whereastheorem 2.9 suggeststhat a local estimator for α 2 i/ is give by ( ˆα 2 i/ = k 2 28 (ˆσ2 (i+k / ˆσ2 i/ ˆτ2 i/ˆσ4 i/ k k j=1 i+jx 8. We obtai the followig result, which ca be prove i the same way as Theorem 2.9. Theorem 3.4. Let t be arbitrary ad set (ŝ t 2 = 1 t 2k i=1 + ˆB 2 t ˆD 21 The (ŝ t 2 is cosistet for s 2 t. ˆα 2 i/ 2ˆB t ˆD 1 1 t 2k i=1 t 2k i=1 ˆα 2 i/ τ2 ( i,x i/,ˆσ 2 i/ ˆα 2 i/ τ2 ( i,x i/,ˆσ 2 i/ As a cosequece, each statistic /k ˆN t /ŝ t coverges weakly to a ormal distributio. This result will be used to costruct a feasible bootstrap statistic i the followig. 4. Simulatio study. Let us start with a simulatio study cocerig the performace of the rate-optimal ˆR t as a estimator for itegrated volatility of volatility. Throughout this sectio, we will

13 Itegrated volatility of volatility 13 Table 1. Mea/variace ad simulated quatiles of the feasible test statistic (2.9 for ρ =. The last colum gives the relative amout of egative estimates Mea Variace Neg work with the Hesto model oly, ad the parameters are chose as follows: β =.3, κ=5, α=.2 ad ξ=.5. Furthermore, we set X = ad σ 2 =α. Note that the Feller coditio 2κα ξ 2 is satisfied, which esuresthat the process σ 2 is almost surelypositive asrequested.sodoesτ 2,aditisobviousthat (2.1holdsaswell.Thereforeallcoditios from Sectio 2 are satisfied. We discuss the fiite sample properties of ˆR t for differet choices of the correlatio parameter ρ ad the umber of observatios, ad for comparability oly we take to be a square umber ad k equal to 1/2 i all cases, so we have c=1. Theorem 2.6 suggests that such a medium size of c is reasoable for fiite samples, ad additioal results ot reported here also poit towards the fact that k should be chose close to 1/2. Fially, we set t=1. Tables 1 6 below are based o 1 simulatios. Table 1 shows the performace for ρ=, for which we see that it takes quite some time for the asymptotics to kick i. Apparet is a slight overestimatioof the lower tails of the distributio, which seems to origiate from the relatio of the estimators ˆR 1 ad G (3 1,. By costructio, i cases where ˆR 1 is uderestimatig the true quatity, it is typically the case that icremets of ˆσ 2 are relatively small. As these icremets occur i G (3 1, as well, most likely the asymptotic variace is uderestimated as well, which explais a too large egative stadardised statistic. The same effect is visible for the upper quatiles as well (but resultig i a overestimatio, ad this simple explaatio is supported by a Table 2. Mea/variace ad simulated quatiles of the feasible test statistic (2.9 for ρ =.2. The last colum gives the relative amout of egative estimates Mea Variace Neg

14 14 M. Vetter Table 3. Mea/variace ad simulated quatiles of the feasible test statistic (2.8 for k = 3/4, l = 1/2 ad ρ= Mea Variace detailed look at simulatio results ot reported here which reveal that the estimatio of the asymptotic variace is extremely accurate for moderate sizes of ˆR1 1 τ2 s ds, but becomes worse whe the deviatio is rather large. Similar coclusios ca be draw for the case of a moderately egative ρ=.2. We proceed with the fiite sample behaviour of the statistics ˆT t, for which we have a lot of freedom i choosig k ad l. However, i order for both M to be rather small ad the coditio 3/2 M 3/2 m 1 to be satisfied, we choose M = 3/4 ad m = 1/2, resultig i a rate of covergece of about 1/8. Also, we restrict ourselves to ρ=. As expected, the approximatio of the omial level is rather poor i this situatio, both whe reproducig mea/variace ad the quatiles i the tails. Empirically the results do ot improve for other choices of k ad l. Note from Table 3 ad Table 4 that results do ot differ very much whe choosig either k or l large, apart from the remarkable expectio of a larger l ad =1. But eve i this case, the results are ot better tha for the rate-optimal ˆR t, which is why we recommed to choose this oe rather tha ˆT t, eve though oly the latter estimator is esured to be positive. As a example for a applicatio i goodess-of-fit testig, we have costructed a test for a Hesto-like volatility structure via a bootstrap procedure as follows: Based o the Table 4. Mea/variace ad simulated quatiles of the feasible test statistic (2.8 for l = 3/4, k = 1/2 ad ρ= Mea Variace

15 Itegrated volatility of volatility 15 Table 5. Simulated level of the bootstrap test based o the stadardised Kolmogorov Smirov statistic Y observatio that for each t, /k ˆN t /ŝ t coverges weakly to a stadard ormal distributio if the ull is satisfied, it seems reasoable to reject the hypothesis for large values of the stadardised Kolmogorov Smirov statistic Y =sup i 2k /k ˆN i/ /ŝ i/. Sice its (asymptotic distributio is i geeral hard to assess, we used bootstrap quatilesistead,adpreciselywehavegeeratedbootstrapdatax (b i/,b=1,...,b,followig the equatio X t = σ s dw s, (σ t 2 = ˆα+ ˆκ(ˆα (σ s 2 ds+ ˆξ σ s dv s. Here, W ad V are idepedet Browia motios, ad we have idetified ˆα with the realised volatility of the origial data (which is a measure for the average volatility over [,1] ad defied ˆξ= ˆθ 1/2, sice both quatities coicide uder the ull. Fially, we have simply set ˆκ=5ˆθ/ˆα such that Feller s coditio is satisfied. Settig B =2, we have ru 5 simulatios each. Table 6. Simulated rejectio probabilities of the bootstrap test based o the stadardised Kolmogorov Smirov fuctioal statistic Y for various alteratives Alt γ = γ =

16 16 M. Vetter Table 5 shows that the simulated levels are rather close to the expected oes, irrespectively of. We have tested two alteratives from the class of CEV models, amely ad σt 2 =σ2 +κ (α σs 2 ds+v t σt 2 =σ2 +κ (α σs 2 ds+ κ σ 2 s dv s, correspodig to γ = ad γ =2, respectively, ad usig the parameters from above. We see from the simulatio results that the rejectio probabilities are much larger for the secod alterative tha for the first, which ca partially explaied from two observatios: First, the Vasicek model does ot satisfy the assumptios from the previous sectios sice the volatilitymay become egative(i which case it is set to zero; secod,our choiceof ˆκ is resposible for a large speed of mea reversio i the bootstrap algorithm which makes it difficult to distiguish betwee a Hesto-like volatility of volatility ad a costat oe. It is expected that the power improves for a etirely data-drive choice of ˆκ. 5. Coclusio I this paper, we have discussed a oparametric method to estimate the itegrated volatility of volatility process i stochastic volatility models. Our cocept is based o spot volatility estimators, ad just as for stadard realised volatility we use sums of squares of these spot volatility estimators to obtai a global estimator for itegrated volatility of volatility. Two classes of estimators have bee ivestigated oe cosistig of positive estimators with a slow rate of covergece, the other oe beig bias corrected but covergig at the optimal rate 1/4. I both cases, cetral limit theorems are provided, ad we also discuss briefly why a trucated versio could be useful whe there are additioal jumps i the price process. Give the variety of stochastic volatility models (i cotiuous time which are used to describe fiacial data, there is a severe lack i tools o model validatio. Our results fill this gap to a first extet, as we provide a bootstrap method for goodess-of-fit testig i such models which ivestigates whether a specific parametric model for volatility of volatility is appropriate give the data or ot. A rigorous proof that the proposed procedure keeps the asymptotic level ad is cosistet agaist a large class of alteratives has ot bee provided, however, ad is left for future research. A differet issue to take microstructure issues ito accout which are likely to be preset whe data is observed at high-frequecy. Agai it is promisig to combie filterig methods for oisy diffusios with the method proposed i this paper to obtai a estimator for itegrated volatility of volatility i such models as well, but the rate of covergece is expected to drop further. Precise statemets are beyod the scope of the paper as well.

17 Itegrated volatility of volatility 17 Appedix Note first that every left-cotiuous process is locally bouded, thus all processes appearig are. Secod, stadard localisatio procedures as i [5] or [17] allow us to assume that ay locally bouded process is actually bouded, ad that almost surely positive processes ca be regarded as bouded away from zero. Uiversal costats are deoted by C or C r, the latter if we wat to emphasise depedece o some additioal parameter r. Withi the mai corpus, we give the proof of Theorem 2.6 oly, which is the by far most complicated result of this work. Aalogues of Lemma 2.2 ad Propositio 2.3 for the special case of l =k are of course parts of it, ad it is ot difficult to geeralise the proofs i order for both claims to be covered as well. Therefore, these results are ot show explicitly. Let us start with a brief sketch of what we will be doig. I geeral, F- stable covergece of a sequece Z to some limitig variable Z defied o a extesio ( Ω, F, P of the origial space is equivalet to E[h(Z Y] Ẽ[h(ZY] (A.1 for ay bouded Lipschitz fuctio h ad ay bouded F-measurable Y. For details, see, for example, [2] ad related work. Suppose ow that there are additioal variables Z,p ad Z p (the latter defied o the same extesio as Z such that lim limsup E[ Z Z,p ] =, p (A.2 Z,p L (s Z p for all p, (A.3 lim Ẽ[ Z p Z ] =, p (A.4 hold. The the desired stable covergece Z L (s Z follows. Ideed, let ε >. The there exists a δ > such that x y <δ implies h(x h(y <ε. Thus we have E[h(Z Y] E[h(Z,p Y] C(E[h(Z h(z,p 1 { Z Z,p δ}]+e[ h(z h(z,p 1 { Z Z,p <δ}] C(P( Z Z,p δ+ε. We have lim p limsup E[h(Z Y] E[h(Z,p Y] = from Markov iequality, (A.2 ad as ε was arbitrary. lim p Ẽ[h(ZpY] Ẽ[h(ZY] = ca be show similarly usig (A.4, ad (A.3 is by defiitio equivalet to lim E[h(Z,p Y] Ẽ[h(Z py] =. Puttig the latter three claims together (plus the triagle iequality ad the fact that all three limitig coditios o p ad are actually the same gives (A.1. Our aim i this proof is to employ a certai blockig techique, which allows us to make use of a type of coditioal idepedece betwee the summads withi ˆR t. To this ed, we apply the above methodology, so we have to defie a appropriate double

18 18 M. Vetter sequece U,p t, which will correspod to a approximated versio of ˆR t where we sum over the big blocks oly. Some additioal otatio is ecessary. Let p N be arbitrary. We set a l (p=(l 1(p+2k, b l (p=a l (p+pk, c(p=j (p(p+2k +1, the first two for ay,...,j (p with J (p= t 2k /((p+2k. These umbers deped o as well, eve though it does ot show up i the otatio. We defie further Hi = i/ (i 1/ (W s W (i 1/ dw s. I order to exploit the afore-metioed coditioal idepedece, we eed approximatios for A i ad Bi from (2.3. For the sake of brevity, we will oly state the approximated icremets explicitly, which are give by à (i+k/ Ãi/ := k 2σa 2 k l (p/ (H i+j+k Hi+j j=1 = k σ 2 a l (p/ k j=1 (( i+k +jw 2 ( i+jw 2, (A.5 where the latter idetity is a cosequece of Itô formula, ad B (i+k/ B i/ := (i+k/ (β al (p/(w s+k/ W s +η al (p/(w s+k k / W s ds. i/ (A.6 These quatities are defied for i=a l (p,...,b l (p 1, thus over the big blocks. For later reasos, we itroduce similar approximatios over the small blocks. Set C (i+k/ C i/ = k σ 2 b l (p/ k j=1 (( i+k +j W2 ( i+j W2, D (i+k/ D i/ = (i+k/ (β bl (p/(w s+k/ W s +η bl (p/(w s+k k / W s dsds, i/ both for i=b l (p,...,a l+1 (p 1. The the followig claim holds, whose proof is postpoed to the supplemetal file [26]. Lemma A.1. We have E[ A (i+k/ A i/ (Ã(i+k / Ãi/ r ] C r (p 1 r/2, E[ B (i+k/ B i/ ( B (i+k/ B i/ r ] C r (p 1 r/2,

19 Itegrated volatility of volatility 19 as well as E[ A (i+k/ A i/ r ] C r r/4 ad E[ B (i+k/ B i/ r ] C r r/4 for every r >. The latter bouds hold also for the approximated versios, ad the same results are true for the approximatio via icremets of C ad D over the small blocks. Up to a differet stadardisatio, the role of Z,p i this proof will be played by U,p t = J (p U,p l, where b l (p 1 U,p 3 l = ((Ã(i+k 2k / Ãi/+( B (i+k/ B i/ 2 i=a l (p pk [ 6 ] k 2 σa 4 l (p/ +(β2 a l (p/ +η2 a l (p/ (A.7 ivolves quatities from the big blocks oly. The U,p l ca be show to be martigale differeces, ad the most ivolved part i the proof is to use Lemma A.1 to obtai [ t lim limsup E (ˆR t p k (βs 2 +η2 s ds U,p t ] =, (A.8 which is the aalogue of (A.2. Let us focus o the remaiig two steps as well. We set U p t = α(p s dw s, α(p 2 s = p ( 48p+d1 p+2 pc 4 σs p+d 2 pc 2 σs(β 4 s 2 +ηs p+d 3 (βs 2 +η 2 7p s 2 for certai uspecified costats d l,,2,3. I order to prove the stable covergece U,p t k L (s U p t (A.9 we use a well-kow result for triagular arrays of martigale differeces, which is due to Jacod [16]. I particular, the followig three coditios have to be checked. k J (p E a l (p [U,p k J (p E a l (p [(U,p 2 J (p k 2 E a l (p [(U,p P l 2 ] P l 4 ], l (N al+1 (p/ N al (p/], P α(p 2 s ds, (A.1 (A.11 (A.12

20 2 M. Vetter where N is ay compoet of (W,W or a bouded martigale orthogoal to both W ad W. The fial step lim p Ẽ U p t U t = is obvious. A.1. Proof of (A.8 For simplicity, we set η ad ϑ (2 from ow o, as otherwise the proof is exactly the same. I a brief first step, we replace ˆR t by a versio i which the ukow bias ad ot the estimator for it is subtracted, that is we itroduce t 2k Ut 3 = (ˆσ (i+k 2 2k / ˆσ2 i/ 2 6 c 2 i= σ 4 s ds β 2 s ds. Theorem 2.1 i [5] shows that itegrals over σ ca be estimated with rate 1/2, so the assumptio o k ad a stadard argumet regardig boudary terms prove that [ t ] E (ˆR t k βs U 2 ds t =o(1, uiformly i t. A simple cosequece of Lemma A.1 is that the remaider terms i U t are egligible, that is [ lim limsup E p k t 2k i=c(p 3 (ˆσ (i+k 2 2k / ˆσ2 i/ 2 6 ] t c 2 σs 4 ds βs 2 ds =, c(p/ c(p/ usig also boudedess of the processes o the right had side ad the defiitio of c(p. Therefore, we are left to show with Ũ,p t = lim limsup /k p E[ Ũ,p t U,p t ] = J (p ( bl (p 1 i=a l (p + a l+1 (p 1 i=b l (p (ˆσ 2 (i+k / ˆσ2 i/ 2 6 c(p/ c(p/ c 2 σs 4 ds βs 2 ds. (A.13 (A.14 For the itegrals withi (A.14, recall that these are replaced by approximated versios i U,p t. Therefore we have to show for example, [ J lim limsup (p bl (p/ E (β p s 2 βa 2 k l (p/ ]=. ds (A.15 a l (p/

21 Itegrated volatility of volatility 21 For its proof, recall (2.1. The result above follows from ad E E[ J (p bl (p/ s ( J(p bl (p/ s a l (p/ a l (p/ a l (p/ Of course, the similar claim a l (p/ [ lim limsup E p k ] ω r drds C pk 2 J (p ϑ (1 r dw r ds = E Cp 2 1. J (p bl (p/ a l (p/ ( 2 pk Cp 1/2 ( bl (p/ s a l (p/ a l (p/ 2 ϑ (1 r dw r ds (σ 4 s σ4 a l (p/ ds ]= (A.16 holds for the same reasos. We have further [ J lim limsup (p ( ] pk E p k k 2 1 σ 4 c 2 a l (p/ =, (A.17 which by boudedess of σ amouts to prove 3/4 (k 2 c2 =o(1, ad the latter is satisfiedbydefiitio ofk.note that aaloguesof(a.15, (A.16 ad(a.17 aresatisfied over the small blocks as well. The latter claims prove that we are left to show the approximatio over the big blocks, which is lim limsup E p k [ J (p b l (p 1 i=a l (p ad the egligibility of the small blocks, that is 3 2k (((A (i+k/ A i/ +(B (i+k/ B i/ 2 (A.18 ] ((Ã(i+k/ Ãi/+( B (i+k/ B i/ 2 =, [ lim limsup E p k J (p ( al+1 (p 1 i=b l (p 2k 3 2k (ˆσ 2 (i+k / ˆσ2 i/ 2 [ 6 k 2 σ(b 4 l (p/ +β2 b l (p/] ] = (A.19 to obtai (A.8.

22 22 M. Vetter To prove (A.18, the biomial theorem tells us that we ca discuss the approximatio for B, the oe for A ad the mixed part separately. Usig further x 2 y 2 =2y(x y+ (x y 2 adxx yy =(x yy +y(x y +(x y(x y,weseefromlemmaa.1ad the growth coditios that (A.18 follows from lim p limsup 4 r=1 E[ L(j,p ] = with L (1,p = k L (2,p = k L (3,p = k L (4,p = k J (p J (p J (p J (p b l (p 1 i=a l (p b l (p 1 i=a l (p b l (p 1 i=a l (p b l (p 1 i=a l (p 1 k ((B (i+k/ B i/ ( B (i+k/ B i/ ( B (i+k/ B i/, 1 k ((B (i+k/ B i/ ( B (i+k/ B i/ (Ã(i+k / Ãi/, 1 k ((A (i+k/ A i/ (Ã(i+k / Ãi/(Ã(i+k / Ãi/, 1 k ((A (i+k/ A i/ (Ã(i+k / Ãi/( B (i+k/ B i/. (A.2 (A.21 (A.22 (A.23 Proofs of these claims ca be foud i the supplemetary material [26]. Fially, to obtai (A.19, we compute the coditioal expectatio of the approximated icremets, ad we will do this for the à ad B terms oly. We have E a l (p [(Ã(i+k / Ãi/ 2 ] = 2 k 2 k σa 4 l (p/ j=1 = 4 k σ 4 a l (p/ E[(( i+k +jw 2 ( i+k W 2 2 ] (A.24 as well as E a l (p [( B (i+k/ B i/ 2 ] =2 2 k 2 βa 2 l (p/ (i+k/ s i/ i/ E[(W s+k/ W s (W r+k/ W r ]drds

23 Itegrated volatility of volatility 23 =2 2 k 2 βa 2 l (p/ = 2k 3 β2 a l (p/. (i+k/ s i/ i/ ( r+k / s drds The expectatio of the mixed part is zero. Obviously, we have E b l (p [( al+1 (p 1 i=b l (p 2k 3 2k (( C (i+k/ C i/ +( D (i+k/ D i/ 2 [ 6 k 2 σb 4 l (p/ +β2 b l (p/] ] = as well. (A.19 the follows from the fact that k J (p E [( al+1 (p 1 i=b l (p 3 2k (( C (i+k/ C i/ +( D (i+k/ D i/ 2 2k [ ] 2] 6 k 2 σb 4 l (p/ +β2 b l (p/ is bouded by a costat times p 1, usig Lemma A.1. A.2. Proof of (A.9 Let us check the coditios for stable covergece i this step, where particularly the proof of (A.1 is tedious. Write U,p l = 3 s=1 U,p,s l with b l (p 1 U,p,1 l = i=a l (p b l (p 1 U,p,2 l = i=a l (p 3 2k 3 2k ( (Ã(i+k/ Ãi/ 2 4 σa 4 k l (p/, ( ( B (i+k/ B i/ 2 2k 3 β2 a l (p/, b l (p 1 U,p,3 3 l = (Ã(i+k k / Ãi/( B (i+k/ B i/. i=a l (p We have see i the fial step above that these terms are ideed martigale differeces, ad it turs out that olythe (U,p,s l 2 terms areresposiblefor the coditioalvariace, whereas the remaiig mixed oes are of small order each. To summarize, the followig lemma holds which is prove i the supplemetary material [26].

24 24 M. Vetter Lemma A.2. We have E a l (p [(U,p,1 l 2 ]=σa 8 48p+d 1 l (p/ k 2 +O P (p 3/2, ( E a l (p [(U,p,2 l k ]=βa 4 2 l (p/ 7 p+d 2 2 +O P(p 3/2, E a l (p [(U,p,3 l 2 ]=σa 4 12p+d 3 l (p/ β2 a l (p/ +O P (p 3/2, for certai uspecified costats d m, m=1,2,3, as well as for each r s We use Lemma A.2 to obtai k J (p = pk E a l (p [(U,p l 2 ] ( 2 J (p k 4 ( 48+ d 1 p + k 2 E a l (p [U,p,r l U,p,s l ]=O P (p 3/2. σ 8 a l (p/ ( 12+ d ( σa 4 p l (p/ β2 a l (p/ d ( 3 1 βa 4 p l (p/ +O p, 1/2 thus (A.1 holds usig k c 1/2. Simpler to obtai is (A.11, as Lemma A.1 gives 2 J (p k 2 E a l (p [(U,p l 4 ] C 3 pk 3 p 4 2, which coverges to zero i the usual sese. Fially, oe ca prove E a l (p [ bl (p 1 i=a l (p 3 2k ((Ã(i+k / Ãi/ +( B (i+k/ B i/ 2 (N al+1 (p/ N al (p/ ] =, (A.25 where N is either W or W or whe N is a bouded martigale, orthogoal to (W,W.Focusothefirstcaseaddecompose((Ã(i+k / Ãi/+( B (i+k/ B i/ 2 via the biomial theorem. For the pure à ad the pure B term, the claim follows immediately from properties of the ormal distributio upo usig that σ al (p/ or β al (p/ are F al (p/ measurable. For the mixed term, oe has to use the special form

25 Itegrated volatility of volatility 25 of à (i+k/ Ãi/ as a differece of two sums, ad a symmetry argumet proves (A.25 i this case. For a orthogoal N, we use stadard calculus. By Itô formula, both (Ã(i+k / Ãi/ 2 ad ( B (i+k/ B i/ 2 are a measurable variable times the sum of a costat ad a stochastic itegral with respect to W ad W, respectively. Thus (A.25 holds. I the mixed case, we use itegratio by parts formula to reduce (Ã(i+k / Ãi/( B (i+k/ B i/ tothe sum ofacostat,adw- ad a dw -itegral. The the same argumet applies. Altogether, this gives (A.12. Ackowledgemets The author is grateful for fiacial support through the collaborative research ceter Statistical modelig of oliear dyamic processes (SFB 823 of the Germa Research Foudatio (DFG. Special thaks go to two aoymous referees for their valuable commets o earlier versios of this paper. Supplemetary Material Additioal proofs for claims made i the article (DOI: 1.315/14-BEJ648SUPP;.pdf. We provide several proofs for either theorems from the mai corpus or additioal steps discussed i the Appedix. Refereces [1] Aït-Sahalia, Y. ad Kimmel, R. (27. Maximum likelihood estimatio of stochastic volatility models. J. Fiacial Ecoomics [2] Alvarez, A., Paloup, F., Potier, M. ad Savy, N. (212. Estimatio of the istataeous volatility. Stat. Iferece Stoch. Process MR [3] Badi, F. ad Reò, R. (28. Noparametric stochastic volatility. Techical report. [4] Bardorff-Nielse, O. ad Veraart, A. (29. Stochastic volatility of volatility i cotiuous time. Techical report. [5] Bardorff-Nielse, O.E., Graverse, S.E., Jacod, J., Podolskij, M. ad Shephard, N. (26. A cetral limit theorem for realised power ad bipower variatios of cotiuous semimartigales. I From Stochastic Calculus to Mathematical Fiace Berli: Spriger. MR [6] Bardorff-Nielse, O.E., Hase, P.R., Lude, A. ad Shephard, N. (211. Multivariate realised kerels: Cosistet positive semi-defiite estimators of the covariatio of equity prices with oise ad o-sychroous tradig. J. Ecoometrics MR [7] Bollerslev, T. ad Zhou, H. (22. Estimatig stochastic volatility diffusio usig coditioal momets of itegrated volatility. J. Ecoometrics MR [8] Cherov, M. ad Ghysels, E. (2. Estimatio of stochastic volatility models for the purpose of optio pricig. I Computatioal Fiace 1999 (Y. Abu-Mostafa, B. LeBaro, A. Lo ad A. Weiged, eds Cambridge: MIT Press.

26 26 M. Vetter [9] Comte, F., Geo-Catalot, V. ad Rozeholc, Y. (21. Noparametric estimatio for a stochastic volatility model. Fiace Stoch MR [1] Dette, H. ad Podolskij, M. (28. Testig the parametric form of the volatility i cotiuous time diffusio models A stochastic process approach. J. Ecoometrics MR [11] Dette, H., Podolskij, M. ad Vetter, M. (26. Estimatio of itegrated volatility i cotiuous-time fiacial models with applicatios to goodess-of-fit testig. Scad. J. Stat MR [12] Geo-Catalot, V., Jeatheau, T. ad Laredo, C. (1999. Parameter estimatio for discretely observed stochastic volatility models. Beroulli MR [13] Gloter, A. (27. Efficiet estimatio of drift parameters i stochastic volatility models. Fiace Stoch MR [14] Hesto, S. (1993. A closed-form solutio for optios with stochastic volatility with applicatios to bods ad currecy optios. Rev. Fiacial Studies [15] Hoffma, M. (22. Rate of covergece for parametric estimatio i a stochastic volatility model. Stochastic Process. Appl MR [16] Jacod, J. (1997. O cotiuous coditioal Gaussia martigales ad stable covergece i law. I Sémiaire de Probabilités XXXI. Lecture Notes i Math Berli: Spriger. MR [17] Jacod, J. (28. Asymptotic properties of realized power variatios ad related fuctioals of semimartigales. Stochastic Process. Appl MR [18] Jacod, J. ad Protter, P. (212. Discretizatio of Processes. Stochastic Modellig ad Applied Probability 67. Heidelberg: Spriger. MR [19] Jacod, J. ad Rosebaum, M. (213. Quarticity ad other fuctioals of volatility: Efficiet estimatio. A. Statist MR [2] Jacod, J. ad Shiryaev, A.N. (23. Limit Theorems for Stochastic Processes, 2d ed. Grudlehre der Mathematische Wisseschafte [Fudametal Priciples of Mathematical Scieces] 288. Berli: Spriger. MR [21] Joes, C.S. (23. The dyamics of stochastic volatility: Evidece from uderlyig ad optios markets. J. Ecoometrics Frotiers of fiacial ecoometrics ad fiacial egieerig. MR22525 [22] Macii, C. (29. No-parametric threshold estimatio for models with stochastic diffusio coefficiet ad jumps. Scad. J. Stat MR [23] Podolskij, M. ad Vetter, M. (21. Uderstadig limit theorems for semimartigales: A short survey. Stat. Neerl MR [24] Reò, R. (26. Noparametric estimatio of stochastic volatility models. Ecoom. Lett MR [25] Vetter, M. (212. Estimatio of correlatio for cotiuous semimartigales. Scad. J. Stat MR3847 [26] Vetter, M. (214. Supplemet to Estimatio of itegrated volatility of volatility with applicatios to goodess-of-fit testig. DOI:1.315/14-BEJ648SUPP. [27] Vetter, M. ad Dette, H. (212. Model checks for the volatility uder microstructure oise. Beroulli MR [28] Wag, C.D. ad Myklad, P.A. (214. The estimatio of leverage effect with highfrequecy data. J. Amer. Statist. Assoc MR Received July 212 ad revised March 214

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

. (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

CENTRE FOR ECONOMETRIC ANALYSIS

CENTRE FOR ECONOMETRIC ANALYSIS CENTRE FOR ECONOMETRIC ANALYSIS CEA@Cass http://www.cass.city.ac.uk/cea/idex.html Cass Busiess School Faculty of Fiace 106 Buhill Row Lodo EC1Y 8TZ Testig for Oe-Factor Models versus Stochastic Volatility

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

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

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

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

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

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

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

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

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

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

ST 305: Exam 2 Fall 2014

ST 305: Exam 2 Fall 2014 ST 305: Exam Fall 014 By hadig i this completed exam, I state that I have either give or received assistace from aother perso durig the exam period. I have used o resources other tha the exam itself ad

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

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

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

EFFICIENT ESTIMATION OF INTEGRATED VOLATILITY FUNCTIONALS VIA MULTISCALE JACKKNIFE

EFFICIENT ESTIMATION OF INTEGRATED VOLATILITY FUNCTIONALS VIA MULTISCALE JACKKNIFE EFFICIENT ESTIMATION OF INTEGRATED VOLATILITY FUNCTIONALS VIA MULTISCALE JACKKNIFE By Jia Li ad Yuxiao Liu ad Dacheg Xiu Duke Uiversity, Uiversity of North Carolia at Chapel Hill ad Uiversity of Chicago

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

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

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

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

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

Sequences and Series

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

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

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

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

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

SEMIPARAMETRIC INFERENCE FOR INTEGRATED VOLATILITY FUNCTIONALS USING HIGH-FREQUENCY FINANCIAL DATA. Yunxiao Liu

SEMIPARAMETRIC INFERENCE FOR INTEGRATED VOLATILITY FUNCTIONALS USING HIGH-FREQUENCY FINANCIAL DATA. Yunxiao Liu SEMIPARAMETRIC INFERENCE FOR INTEGRATED VOLATILITY FUNCTIONALS USING HIGH-FREQUENCY FINANCIAL DATA Yuxiao Liu A dissertatio submitted to the faculty of the Uiversity of North Carolia at Chapel Hill i partial

More information

APPLICATION OF GEOMETRIC SEQUENCES AND SERIES: COMPOUND INTEREST AND ANNUITIES

APPLICATION OF GEOMETRIC SEQUENCES AND SERIES: COMPOUND INTEREST AND ANNUITIES APPLICATION OF GEOMETRIC SEQUENCES AND SERIES: COMPOUND INTEREST AND ANNUITIES Example: Brado s Problem Brado, who is ow sixtee, would like to be a poker champio some day. At the age of twety-oe, he would

More information

Research Paper Number From Discrete to Continuous Time Finance: Weak Convergence of the Financial Gain Process

Research Paper Number From Discrete to Continuous Time Finance: Weak Convergence of the Financial Gain Process Research Paper Number 197 From Discrete to Cotiuous Time Fiace: Weak Covergece of the Fiacial Gai Process Darrell Duffie ad Philip Protter November, 1988 Revised: September, 1991 Forthcomig: Mathematical

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

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

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

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

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

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

BASIC STATISTICS ECOE 1323

BASIC STATISTICS ECOE 1323 BASIC STATISTICS ECOE 33 SPRING 007 FINAL EXAM NAME: ID NUMBER: INSTRUCTIONS:. Write your ame ad studet ID.. You have hours 3. This eam must be your ow work etirely. You caot talk to or share iformatio

More information

Lecture 5 Point Es/mator and Sampling Distribu/on

Lecture 5 Point Es/mator and Sampling Distribu/on Lecture 5 Poit Es/mator ad Samplig Distribu/o Fall 03 Prof. Yao Xie, yao.xie@isye.gatech.edu H. Milto Stewart School of Idustrial Systems & Egieerig Georgia Tech Road map Poit Es/ma/o Cofidece Iterval

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

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

Models of Asset Pricing

Models of Asset Pricing 4 Appedix 1 to Chapter Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

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

Testing for Jumps: A Delta-Hedging Perspective

Testing for Jumps: A Delta-Hedging Perspective Testig for Jumps: A Delta-Hedgig Perspective Jia Li Priceto Uiversity Departmet of Ecoomics ad Bedheim Ceter for Fiace Priceto, NJ, 8544 This Versio: Jue 6, 211 Abstract We measure asset price jumps by

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

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

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

Appendix 1 to Chapter 5

Appendix 1 to Chapter 5 Appedix 1 to Chapter 5 Models of Asset Pricig I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy a asset, we are

More information

Models of Asset Pricing

Models of Asset Pricing APPENDIX 1 TO CHAPTER 4 Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

More information

Models of Asset Pricing

Models of Asset Pricing APPENDIX 1 TO CHAPTER4 Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

More information

A point estimate is the value of a statistic that estimates the value of a parameter.

A point estimate is the value of a statistic that estimates the value of a parameter. Chapter 9 Estimatig the Value of a Parameter Chapter 9.1 Estimatig a Populatio Proportio Objective A : Poit Estimate A poit estimate is the value of a statistic that estimates the value of a parameter.

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

Notes on Expected Revenue from Auctions

Notes on Expected Revenue from Auctions Notes o Epected Reveue from Auctios Professor Bergstrom These otes spell out some of the mathematical details about first ad secod price sealed bid auctios that were discussed i Thursday s lecture You

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

of Asset Pricing R e = expected return

of Asset Pricing R e = expected return Appedix 1 to Chapter 5 Models of Asset Pricig EXPECTED RETURN I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy

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

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

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

Valuation of options on discretely sampled variance: A general analytic approximation

Valuation of options on discretely sampled variance: A general analytic approximation Valuatio of optios o discretely sampled variace: A geeral aalytic approximatio Gabriel Drimus Walter Farkas, Elise Gourier 3 Previous versio: Jauary 3 This versio: July 4 Abstract The values of optios

More information

ON DIFFERENTIATION AND HARMONIC NUMBERS

ON DIFFERENTIATION AND HARMONIC NUMBERS ON DIFFERENTIATION AND HARMONIC NUMBERS ERIC MORTENSON Abstract. I a paper of Adrews ad Uchimura [AU, it is show how differetiatio applied to hypergeometric idetities produces formulas for harmoic ad q-harmoic

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

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

Estimation of a noisy subordinated Brownian Motion via two-scales power variations Estimatio of a oisy subordiated Browia Motio via two-scales power variatios José E. Figueroa-López Departmet of Mathematics, Washigto Uiversity i St. Louis ad Kiseop Lee Departmet of Statistics, Purdue

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

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

Confidence Intervals Introduction

Confidence Intervals Introduction Cofidece Itervals Itroductio A poit estimate provides o iformatio about the precisio ad reliability of estimatio. For example, the sample mea X is a poit estimate of the populatio mea μ but because of

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

ii. Interval estimation:

ii. Interval estimation: 1 Types of estimatio: i. Poit estimatio: Example (1) Cosider the sample observatios 17,3,5,1,18,6,16,10 X 8 X i i1 8 17 3 5 118 6 16 10 8 116 8 14.5 14.5 is a poit estimate for usig the estimator X ad

More information

CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Means and Proportions

CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Means and Proportions CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Meas ad Proportios Itroductio: I this chapter we wat to fid out the value of a parameter for a populatio. We do t kow the value of this parameter for the etire

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

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

1 Estimating the uncertainty attached to a sample mean: s 2 vs.

1 Estimating the uncertainty attached to a sample mean: s 2 vs. Political Sciece 100a/200a Fall 2001 Cofidece itervals ad hypothesis testig, Part I 1 1 Estimatig the ucertaity attached to a sample mea: s 2 vs. σ 2 Recall the problem of descriptive iferece: We wat to

More information

ISBN Copyright 2015 The Continental Press, Inc.

ISBN Copyright 2015 The Continental Press, Inc. TABLE OF CONTENTS Itroductio 3 Format of Books 4 Suggestios for Use 7 Aotated Aswer Key ad Extesio Activities 9 Reproducible Tool Set 183 ISBN 978-0-8454-7897-4 Copyright 2015 The Cotietal Press, Ic. Exceptig

More information

Time-Varying Periodicity in Intraday Volatility. Torben G. Andersen, Martin Thyrsgaard and Viktor Todorov. CREATES Research Paper

Time-Varying Periodicity in Intraday Volatility. Torben G. Andersen, Martin Thyrsgaard and Viktor Todorov. CREATES Research Paper ime-varyig Periodicity i Itraday Volatility orbe G. Aderse, Marti hyrsgaard ad Viktor odorov CREAES Research Paper 2018-5 Departmet of Ecoomics ad Busiess Ecoomics Aarhus Uiversity Fuglesags Allé 4 DK-8210

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

Confidence Intervals. CI for a population mean (σ is known and n > 30 or the variable is normally distributed in the.

Confidence Intervals. CI for a population mean (σ is known and n > 30 or the variable is normally distributed in the. Cofidece Itervals A cofidece iterval is a iterval whose purpose is to estimate a parameter (a umber that could, i theory, be calculated from the populatio, if measuremets were available for the whole populatio).

More information

of Asset Pricing APPENDIX 1 TO CHAPTER EXPECTED RETURN APPLICATION Expected Return

of Asset Pricing APPENDIX 1 TO CHAPTER EXPECTED RETURN APPLICATION Expected Return APPENDIX 1 TO CHAPTER 5 Models of Asset Pricig I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy a asset, we are

More information

Inverse Realized Laplace Transforms for Nonparametric Volatility Density Estimation in Jump-Diffusions

Inverse Realized Laplace Transforms for Nonparametric Volatility Density Estimation in Jump-Diffusions Iverse Realized Laplace Trasforms for Noparametric Volatility Desity Estimatio i Jump-Diffusios Viktor Todorov ad George Tauche April 212 Abstract We develop a oparametric estimator of the stochastic volatility

More information