Consistent non-parametric Bayesian estimation for a time-inhomogeneous Brownian motion Gugushvili, S.; Spreij, P.J.C.

Size: px
Start display at page:

Download "Consistent non-parametric Bayesian estimation for a time-inhomogeneous Brownian motion Gugushvili, S.; Spreij, P.J.C."

Transcription

1 UvA-DARE Digital Academic Repository) Cosistet o-parametric Bayesia estimatio for a time-ihomogeeous Browia motio Gugushvili, S.; Spreij, P.J.C. Published i: ESAIM-Probability ad Statistics DOI:.5/ps/2339 Lik to publicatio Citatio for published versio APA): Gugushvili, S., & Spreij, P. 24). Cosistet o-parametric Bayesia estimatio for a time-ihomogeeous Browia motio. ESAIM-Probability ad Statistics, 8, DOI:.5/ps/2339 Geeral rights It is ot permitted to dowload or to forward/distribute the text or part of it without the coset of the authors) ad/or copyright holders), other tha for strictly persoal, idividual use, uless the work is uder a ope cotet licese like Creative Commos). Disclaimer/Complaits regulatios If you believe that digital publicatio of certai material ifriges ay of your rights or privacy) iterests, please let the Library kow, statig your reasos. I case of a legitimate complait, the Library will make the material iaccessible ad/or remove it from the website. Please Ask the Library: or a letter to: Library of the Uiversity of Amsterdam, Secretariat, Sigel 425, 2 WP Amsterdam, The Netherlads. You will be cotacted as soo as possible. UvA-DARE is a service provided by the library of the Uiversity of Amsterdam Dowload date: 23 Ja 29

2 ESAIM: PS 8 24) DOI:.5/ps/2339 ESAIM: Probability ad Statistics CONSISTENT NON-PARAMETRIC BAYESIAN ESTIMATION FOR A TIME-INHOMOGENEOUS BROWNIAN MOTION Shota Gugushvili ad Peter Spreij 2 Abstract. We establish posterior cosistecy for o-parametric Bayesia estimatio of the dispersio coefficiet of a time-ihomogeeous Browia motio. Mathematics Subject Classificatio. 62G2, 62M5. Received October 8, 22. Revised April 23, 23.. Itroductio Cosider a simple liear stochastic differetial equatio dx t = σt)dw t, X = x, t [, ],.) where W is a Browia motio o some give probability space ad the iitial coditio x ad the square itegrable dispersio coefficiet σ are determiistic. We iterpret equatio.) as a short-had otatio for the itegral equatio t X t = x + σs)dw s, t [, ], where the itegral is the Wieer itegral of σ with respect to the Browia motio W. The process X is thus a time-ihomogeeous Browia motio. The fuctio σ ca be viewed as a sigal trasmitted through a oisy chael, where the oise modelled by the Browia motio) is multiplicative. Note that X is a Gaussia process with mea mt) =x ad covariace ρs, t) = s t. By P σ we will deote the law of the solutio X to.). Assume for simplicity that x = ad deote t i, = i/, i =,...,. Suppose that correspodig to the true dispersio coefficiet σ = σ, oe has a sample X,,...,, from the process X at his disposal. Assumig that σ belogs to some o-parametric class X of dispersio coefficiets, our goal is to estimate σ. This problem for a similar model was treated i [5, 9, 9] usig a frequetist approach. However, a o-parametric Keywords ad phrases. Dispersio coefficiet, o-parametric Bayesia estimatio, posterior cosistecy, time-ihomogeeous browia motio. The research of the first author was supported by The Netherlads Orgaisatio for Scietific Research NWO). Mathematical Istitute, Leide Uiversity, P.O. Box 952, 23 RA Leide, The Netherlads. shota.gugushvili@math.leideuiv.l 2 Korteweg-de Vries Istitute for Mathematics, Uiversity of Amsterdam, P.O. Box 94248, 9 GE Amsterdam, The Netherlads. spreij@uva.l Article published by EDP Scieces c EDP Scieces, SMAI 24

3 BAYESIAN ESTIMATION FOR A TIME-INHOMOGENEOUS BROWNIAN MOTION 333 Bayesia approach to estimatio of σ is also possible. The likelihood correspodig to the observatios X is give by L σ) = 2π ψ X t i, X ti, t i, σ 2 u)du,.2) where ψu) =exp u 2 /2). For a prior Π o X, Bayes formula yields the posterior measure ΠΣ X t,...,x, )= Σ L σ)πdσ) X L σ)πdσ) of ay measurable set Σ X. I the Bayesia paradigm, the posterior ecodes all the iformatio required for iferetial purposes. Oce the posterior is available, oe ca proceed to computatio of other quatities of iterest i Bayesia statistics, such as Bayes poit estimates, Bayes factors ad so o. It has bee recogised sice log that Bayesia procedures should be theoretically grouded through establishig posterior cosistecy, see e.g. [3]. I our cotext posterior cosistecy will mea that for every eighbourhood U σ of σ i a suitable topology) ΠUσ c X t,,...,x t, ) Pσ.3) P σ as the otatio ξ ξ i.3) ad below stads for covergece of a sequece of radom variables ξ to a radom variable ξ i P σ -probability). I other words, a cosistet Bayesia procedure asymptotically puts posterior mass equal to oe o every fixed eighbourhood of the true parameter. This is similar to the study of cosistecy of frequetist estimators. A method that does ot appear to work i the idealised settig whe a ifiite amout of data is available formalised by assumig that the sample size ) should also be uattractive i the fiite sample settig. Hece the importace of a study of posterior cosistecy. The situatio is typically quite subtle i the ifiite-dimesioal Bayesia settig: it is kow that a careless choice of the prior might reder a Bayes procedure icosistet. For a itroductio to cosistecy issues i Bayesia o-parametric statistics see e.g. [6, 25]. Our task i this work is to establish.3) uder suitable assumptios o the class of dispersio coefficiets σ ad the prior Π. Asymptotic properties of Bayesia procedures i estimatio problems for stochastic differetial equatios have bee already cosidered uder various setups i [8, 2 4, 7], primarily i the cotext of oparametric Bayesia estimatio of the drift coefficiet of a stochastic differetial equatio. Computatioal approaches to o-parametric Bayesia iferece for stochastic differetial equatios were studied i [, 5]. Coveiet overviews of the available results are give i [6, 26]. However, i the above works dealig with Bayesia asymptotics it is assumed that either a cotiuous record of observatios is available o the solutio to a stochastic differetial equatio, or that the solutio is observed at equispaced time poits Δ, 2Δ,...,Δ, with asymptotics treated i the latter case uder the assumptio that Δ is idepedet of ad. Our problem, o the other had, requires a differet approach due to a differet samplig scheme ad the fact that ergodicity of the solutio to a stochastic differetial equatio, that played a promiet role i most of the previous works o o-parametric Bayesia approach to statistical iferece for stochastic differetial equatios, is irrelevat i our case. Although the setup we cosider looks simple, to the best of our kowledge our work is the first oe to treat a iferece problem for a stochastic differetial equatio i the so called high-frequecy data case whe Δ = Δ as usig a o-parametric Bayesia approach. The highfrequecy data settig is particularly relevat i fiacial mathematics, where asset prices are ofte modelled through stochastic differetial equatios ad where huge amouts of observatios o them separated by very short time istaces are available. Perhaps the most iterestig feature of the preset work is the method of proof of posterior cosistecy, which differs i certai respects from the curretly used techiques. See Sectio 4 for a discussio. Also the simplicity of our model should ot ecessarily be cosidered a disadvatage: ideed,

4 334 S. GUGUSHVILI AND P. SPREIJ the model is somewhat similar to the Gaussia white oise model see e.g. Chap.7,Sect.4i[]),which,as is kow, has triggered some importat developmets i mathematical statistics. The paper is orgaised as follows: i the ext sectio we formulate our mai result dealig with posterior cosistecy for o-parametric estimatio of the dispersio coefficiet. Sice posterior cosistecy is closely liked with properties of a prior Π, i Sectio 3 we provide a example of a reasoable prior satisfyig the assumptios made i Sectio 2. Sectio 4 cotais a brief discussio o the obtaied result. The proof of our mai theorem is deferred util Sectio 5, while the Appedix cotais two techical lemmas used i Sectio Results The o-parametric class of dispersio coefficiets we will be lookig at is give i the followig defiitio. Defiitio 2.. Let X be the collectio of dispersio coefficiets σ : [, ] [κ, K], such that σ X is Lipschitz with Lipschitz costat M. Here <κ<k< ad <M< are three costats idepedet of the particular σ X. Remark 2.2. Note that for a costat σ we have P σ = P σ. A positivity assumptio o σ X i Defiitio 2. ca hece be viewed as a simple ad atural idetifiability requiremet. Strict positivity assumptio σ κ > allows oe to escape techical complicatios whe maipulatig the likelihood.2) this coditio has already appeared e.g. i [9]), while the upper boud σt) K, t [, ], restricts the size of the o-parametric class X ad is reasoable i light of Defiitio 2.3 give below. Fially, Lipschitz cotiuity of σ comes i hady at various stages of the proof of posterior cosistecy. The otio of posterior cosistecy depeds o a topology o X. Defiitio 2.3. The topology T o X is the topology iduced by the L 2 -orm 2. We ow formalise the cocept of posterior cosistecy. Defiitio 2.4. Let the prior Π be defied o X. We say that posterior cosistecy holds, if for ay fixed σ X ad every eighbourhood U σ of σ i the topology T from Defiitio 2.3 we have as. We summarise our assumptios. ΠUσ c X t,...,x, ) Pσ Assumptio 2.5. Assume that a) the model.) isgivewithx =adσ X, where X is defied i Defiitio 2., b) σ X deotes the true dispersio coefficiet, c) a discrete-time sample {X } from the solutio to.) correspodig to σ is available, where t i, = i/, i =,...,. Let V σ,ε = {σ X : σ σ <ε}, where deotes the L -orm. The followig is the mai result of the paper. Theorem 2.6. Uder Assumptio 2.5 posterior cosistecy as i Defiitio 2.4 holds, provided the prior Π o X satisfies ΠV σ,ε ) > 2.) for ay ε> ad at ay σ X.

5 BAYESIAN ESTIMATION FOR A TIME-INHOMOGENEOUS BROWNIAN MOTION Example of a prior I this sectio we provide a example of a prior satisfyig coditio 2.). Fix <κ<k< ad take a fixed Lipschitz cotiuous fuctio f : R [,K κ] with Lipschitz costat N>adsetσt) =κ + t fhs))ds, where h :[, ] R rages over the set of Hölder cotiuous fuctios of order β, /2) o [, ] for some fixed β. The each σ maps the iterval [, ] ito the iterval [κ, K] adσ is also Lipschitz with Lipschitz costat K, because t σt) σs) = fhu))du K t s. s We take the collectio of these fuctios σ as the collectio X from Assumptio 2.5a). We will ow costruct apriorπox. Let W = W t ) t be a stadard Browia motio over the time iterval [, ] ad let Z be a stadard ormal radom variable idepedet of W. Defie the Browia motio W =W t ) t iitialised at Z by W t = Z + W t ad itroduce the process Y =Y t ) t, where t [,] Y t = κ + t fw s )ds. Our prior Π o X will be the law of the process Y. We have to check that the prior Π satisfies 2.). To that ed take a fixed σ t) =κ + t fh s))ds ad let w be a geeric realisatio of the process W, so that y t = κ + t fw s)ds is the correspodig geeric realisatio of the process Y. We have ΠV σ,ε) =Πy : y σ <ε) Π w : w h < ε ), N because t sup [fw s ) fh s))]ds fw) fh ) N w h. By Lemma 5.3 i [22], Π w : w h < ε ) exp N if g: g h <ε/2n) ) 2 g 2 H Π W < ε ) 2N Here H deotes the Reproducig Kerel Hilbert Space RKHS) of the process W,g H, while H is the RKHS orm see [22] for a detailed treatmet of these cocepts with a view towards o-parametric Bayesia statistics). I our case H cosists of absolutely cotiuous fuctios g :[, ] R, such that g 2 <, while the RKHS orm is give by g H = g 2 ) + g 2 2, see page 446 i [2]. Note that the set {g H : g h <ε/2n)} is ot empty, as h ca be approximated arbitrarily closely i the L -orm by the covolutio h k b of h with asmoothkerelk b ) =/b)k /b) cf. p. 446 i [2]; we assume that k 2 < ad b ). Furthermore, Π W <ε/2n)) >, because W has a strictly positive desity. Coditio 2.) easily follows. I case oe is iterested i a smoother class of dispersio coefficiets σ tha what we have just costructed, oe ca simply take i the above costructio of the prior Π a smoother, say β times differetiable fuctio f, ad replace the Browia motio W with a Riema Liouville process R =R t ) t with Hurst parameter β, R t = β Z k t k + k= t t s) β /2 d W t, where Z k s are stadard ormal radom variables, W is a stadard Browia motio ad Z,Z,...,Z β, W are idepedet. See Sectio 4.2 i [2] for more iformatio o the Riema Liouville processes. Argumets similar to the oes give above yield that i this case as well 2.) is satisfied.

6 336 S. GUGUSHVILI AND P. SPREIJ 4. Discussio I the preset work we established posterior cosistecy for a statistical model obtaied from a simple liear stochastic differetial equatio. Geeral techiques for provig posterior cosistecy for a wide rage of statistical models are by ow well-developed. I the i.i.d. settig, broadly speakig, two mai approaches exist i the literature: a classical approach as epitomised e.g. by [, 8] we combie these two papers ito oe category, because they i some sese make use of assumptios of similar type, although their actual assertios are differet), ad a martigale approach developed more recetly i [23,24]. The first approach was exteded to the settig of idepedet o-idetically distributed observatios i [2], see i particular Theorem A. there. The secod approach was exteded to the case of discretely observed Markov processes i [7]. A geeral theorem for posterior cosistecy i [2] makes two requiremets: firstly, the prior must put sufficiet mass i arbitrarily small eighbourhoods of the true parameter i a appropriate topology), ad secodly, a sequece of sieves icreasig sequece of subsets of the parameter set) guarateeig existece of certai expoetially cosistet tests has to be exhibited; see page 56 i [2] for additioal details. Although i our settig the observatios X,i =,...,, are ot idepedet, the icremets X X ti, are, ad it appears coceivable that Theorem A. i [2] could be used to establish posterior cosistecy i our model as well. However, we opted for a differet approach, see the proof of our posterior cosistecy result, Theorem 2.6. A similarity shared by Theorem A. i [2] adtheorem2.6 is that both theorems require that the prior puts sufficiet mass i arbitrarily small eighbourhoods of the true parameter i appropriate topologies). A differece is that due to the special structure of our model we do ot eed to make ay referece to tests ad sieves, but ca establish posterior cosistecy by directly maipulatig the posterior; see the Proof of Theorem 2.6 for details. I this sese our approach to provig posterior cosistecy appears to be more direct ad more elemetary tha the oe that would employ Theorem A. i [2]. Neither do we make referece to etropy argumets as doe e.g. i []. As far as the martigale approach to posterior cosistecy for ergodic Markov processes is cocered, we ca be brief here: ergodicity is irrelevat i our settig ad i fact our special samplig scheme seems to make geeralisatio or modificatio of the argumets from [7], [23, 24] impossible. Next a brief remark o coditio 2.) o the prior Π is i order. Although it is formulated i terms of eighbourhoods i the L -orm, the assertio retured by Theorem 2.6 employs the topology iduced by the L 2 -orm. A discrepacy betwee orms used is however ot ucommo i posterior cosistecy results. See for istace []. 5. Proofs ProofofTheorem2.6. Let U σ be a arbitrary, but fixed eighbourhood of σ i the topology T ad let Ũσ,ε = {σ X : σ σ 2 <ε}. There exists ε>, such that Ũσ,ε U σ, ad hece U c σ Ũ c σ,ε. Fix such a ε. I order to prove the theorem, it thus suffices to show that as. Write ΠŨ σ c X,ε t,...,x t, ) Pσ 5.) L ΠŨ σ c Ũσ,ε X t,...,x t, )= c σ)πdσ) R,ε Ũσ X L = c σ)πdσ),ε σ)πdσ) X R, 5.2) σ)πdσ) where R σ) =L σ)/l σ ) deotes the likelihood ratio. We will separately boud the umerator ad deomiator o the right-had side of the last equality i 5.2) we will use the otatio D for the deomiator ad N for the umerator) ad the combie the bouds to establish 5.). As we will see, the left-had side of 5.) i fact decays expoetially fast to zero. Note that whe establishig posterior cosistecy, [] ad[24] also treat the umerator ad deomiator i the expressio for the posterior separately, but similarity of our approach to the oe i those papers largely eds here.

7 BAYESIAN ESTIMATION FOR A TIME-INHOMOGENEOUS BROWNIAN MOTION 337 Let S σ) = log R σ). The D = X exps σ))πdσ). Now S σ) = t log i, σu)du 2 2 X t i, X ti, ) 2 2 = T, σ)+t 2, σ), X X ti,)2 σ 2 u)du with obvious defiitios of T, σ) adt 2, σ). Let ε > be a costat with its value to be chose appropriately later o. Sice D V σ, ε R σ)πdσ), Lemmas A. ad A.2 from the Appedix ad formula A.) givethat with probability tedig to oe, D V σ, ε exp 4K ) κ 2 ε Πdσ) =exp 4K ) κ 2 ε ΠV σ, ε). By assumptio 2.), ΠV σ, ε) >. Fix a costat β =5K ε/κ 2. The for all large eough, exp 4K ) κ 2 ε ΠV σ, ε) e β. As a cosequece, with probability tedig to oe, D e β 5.3) as. This is our required lower boud for the deomiator D. Usig similar techiques, we will ext treat the umerator N. Firstly, ote that by elemetary argumets oe ca show that for a arbitrary fixed costat C>there exists aother costat c>, such that the iequality log + y) y cy 2, <y C holds oe ca take c =[2C +)] ). Hece t log i, σ 2u)du [σ 2 u) σ2 u)]du c [σ 2 u) σ2 u)]du for some costat c idepedet of σ X,iad. Therefore, after a simple, but legthy computatio employig Assumptio 2.5a), cf. the proof of Lemma A., T, σ) [σu) 2 σ 2 u)]du 2 c t i, [σu) 2 σ 2 2 u)]du 2 = σu) 2 σ 2 u) 2 σ 2 du c σu) 2 σ 2 u)) 2 ) u) 2 σ 4 du + O, u) where the remaider term is of order uiformly i σ X. Hece S σ) c σu) 2 σ 2 u)) 2 2 σ 4 du 5.4) u) + T 2, σ)+ 2 + O σ 2u) σ2 u) σ 2 du 5.5) u) ). 5.6) 2

8 338 S. GUGUSHVILI AND P. SPREIJ To boud from above the term o the right-had side of iequality 5.4), use the fact that c 2 σ 2 u) σ 2 u)) 2 σ 4 u) du 2κ2 c K 4 ε2 for σ Ũ σ c. Furthermore, by Lemma A.2, uiformly i σ Ũ c,ε σ,ε ad with probability tedig to oe as, the term 5.5) is smaller tha ay positive umber fixed beforehad. So is the term 5.6). Therefore, uiformly i σ Ũ σ c,ε ad with probability tedig to oe as, S σ) κ2 c K 4 ε2, say. Thus with probability tedig to oe as, ) N = exps σ))πdσ) exp κ2 c K 4 ε2. 5.7) Ũ c σ,ε This fiishes boudig from above the umerator N. We ow combie bouds 5.3) ad5.7) to coclude that with probability tedig to oe as, κ ΠŨ σ c 2 X c,ε t,...,x t, ) exp K 4 ε2 5K ) ) κ 2 ε. Pickig ε small eough, so that κ 2 c K 4 ε2 5K ε >, κ2 implies 5.) ad completes the proof of the theorem. Appedix A Lemma A.. Uder the same assumptios as i Theorem 2.6, fort, σ) as i the Proof of Theorem 2.6, all σ V σ, ε simulateously ad for large eough, T, σ) 2K ε/κ 2. Proof. The elemetary iequality gives that y log + y), y + y t log i, σu)du 2 > [σu) 2 σ 2 u)]du σ 2u)du Next, employig Assumptio 2.5a) ad c), by a simple computatio oe ca show that [σ 2u) σ2 u)]du σ 2u)du = σ2 ) σ 2 ) ) σ 2t + O, i,) where the remaider term is of order uiformly i σ X. Therefore, T, σ) 2 σ 2 ) σ 2 ) σ 2 ) + O )

9 BAYESIAN ESTIMATION FOR A TIME-INHOMOGENEOUS BROWNIAN MOTION 339 By aother simple computatio, σt 2 i, ) σ 2 ) σ 2) = σu) 2 σ 2 ) u) σ 2u) du + O, where the remaider term is of order uiformly i σ X. For σ V σ, ε we have uder Assumptio 2.5a) that σu) 2 σ 2 u) σ 2u) du 2K κ 2 ε. A.) This implies the statemet of the lemma. Lemma A.2. Deote Q σ) = T 2,σ)+ σu) 2 σ 2 u) 2 σ 2 du u), where T 2, σ) is defied i the Proof of Theorem 2.6. The uder the same assumptios as i Theorem 2.6, sup σ X Q σ) Pσ as. Furthermore, for ay fixed ε >, for all σ V σ, ε simultaeously, with probability tedig to oe as,t 2, σ) 2K ε/κ 2. Proof. The first statemet of the lemma will be derived from a applicatio of Theorem 8.4 from [2]. I particular, viewig Q as a process o X with bouded sample paths, we will show that it coverges i distributio to a zero process o X. The first statemet of the lemma will the be a cosequece of equivalece of covergece i distributio ad i probability for costat limits. Note that i order to circumvet possible o)-measurability issues, outer probability is employed i the formulatio of Theorem 8.4 from [2] see Sect. 8.2 i [2] for more iformatio o outer probability). Sice o such problems will arise i our settig, we ca istead directly work uder probability P σ. Ideed, the summads i T 2, σ) are of the form F i, σ)x X ti, ) 2,withF i, the obvious fuctioal of σ. Hece takig the supremum over σ does ot affect the measurability property of T 2, σ). I order to apply Theorem 8.4 from [2], we eed to verify its coditios. I our settig they reduce to the followig oes: firstly, margial vectors of Q must coverge i distributio to zero vectors, i.e. Q σ ),...,Q σ l )) Dσ,...,), }{{} l N. A.2) l Secodly, the tightess coditio must be satisfied: for arbitrary costats η>adξ>, oe must be able to fid a partitio of X ito fiitely may X,...,X l, such that ) lim sup P σ sup sup Q σ ) Q σ 2 ) ξ η. A.3) k l σ,σ 2 X k Deote F i, = σx tj,,j =,...,i) ad χ i, σ) = 2 X t i, X ti, ) 2 [σ 2u) σ2 u)]du σ 2u)du t i, Also let E σ be the expectatio operator with respect to measure P σ. Note that E σ χ i, σ) F i, )= 2 [σ 2u) σ2 u)]du = 2 σ 2u) σ2 u) σ 2 du + O u) ),

10 34 S. GUGUSHVILI AND P. SPREIJ where the remaider term is of order uiformly i σ X. Furthermore, Assumptio 2.5a) yields that E σ χ 2 i, σ) F i,) = where the order boud is uiform i σ X. It follows that [σ 2 u) σ2 u)]du E σ χ 2 i, σ) F i,). 2 ) = O 2, Lemma 9 i [4] the implies that Q σ) Pσ. This verifies A.2). We will ow check A.3). Fix ξ ad η i A.3). By a legthy, but simple computatio employig Assumptio 2.5a) ad the triagle iequality, Q σ ) Q σ 2 ) K κ 4 σ σ 2 X X ti,) 2 + K3 κ 4 σ σ 2. A.4) By the Arzelà Ascoli theorem, uder Assumptio 2.5a) the family X is totally bouded for the supremum metric. By defiitio this meas that for every ζ>there exists a fiite set X X, such that for ay σ X there is some σ X with σ σ <ζ/2. This ad the triagle iequality imply existece of a fiite partitio X,...,X l of X, such that sup sup σ σ 2 <ζ. A.5) k l σ,σ 2 X k Furthermore, by the defiitio of the quadratic variatio of the process X, X X ti,) 2 Pσ σ 2 u)du. A.6) Combiatio of A.4) A.6) yields A.3) forζ small eough, ad cosequetly the first statemet of the lemma too. The secod statemet of the lemma is a cosequece of the first oe, the fact that σ V σ, ε, a aalogue of iequality A.), σ 2u) σ2 u) 2 σ 2 du u) K κ 2 ε, ad of a simple rearragemet T 2, σ) = { T 2, σ)+ 2 This completes the proof of the lemma. σu) 2 σ 2 } u) σ 2 du u) 2 σu) 2 σ 2 u) σ 2 du. u) Ackowledgemets. The authors would like to thak the referees for a careful readig of the mauscript ad useful suggestios. Refereces [] A. Barro, M.J. Schervish ad L. Wasserma, The cosistecy of posterior distributios i oparametric problems. A. Statist ) [2] N. Choudhuri, S. Ghosal ad A. Roy, Bayesia estimatio of the spectral desity of a time series. J. Amer. Statist. Assoc ) 5 59.

11 BAYESIAN ESTIMATION FOR A TIME-INHOMOGENEOUS BROWNIAN MOTION 34 [3] P. Diacois ad D. Freedma, O the cosistecy of Bayes estimates. With a discussio ad a rejoider by the authors. A. Statist ) 67. [4] V. Geo-Catalot ad J. Jacod, O the estimatio of the diffusio coefficiet for multi-dimesioal diffusio processes. A. Ist. Heri Poicaré Probab. Statist ) 9 5. [5] V. Geo-Catalot, C. Laredo ad D. Picard, Noparametric estimatio of the diffusio coefficiet by wavelets methods. Scad. J. Statist ) [6] S. Ghosal, J.K. Ghosh ad R.V. Ramamoorthi, Cosistecy issues i Bayesia oparametrics. Asymptotics, Noparametrics, ad Time Series. Vol. 58 of Textbooks Moogr. Dekker, New York 999) [7] S. Ghosal ad Y. Tag, Bayesia cosistecy for Markov processes. Sakhyā 68 26) [8] S. Gugushvili ad P. Spreij, No-parametric Bayesia drift estimatio for stochastic differetial equatios 22). Preprit arxiv: [math.st]. [9] M. Hoffma, Miimax estimatio of the diffusio coefficiet through irregular sampligs. Statist. Probab. Lett ) 24. [] I.A. Ibragimov ad R.Z. Has miskiĭ, Asimptoticheskaya teoriya otseivaiya [Asymptotic Theory of Estimatio] Russia). Nauka, Moscow 979). [] F. va der Meule, M. Schauer ad H. va Zate, Reversible jump MCMC for oparametric drift estimatio for diffusio processes. Comput. Statist. Data Aal. 7 24) Available o [2] F.H. va der Meule, A.W. va der Vaart ad J.H. va Zate, Covergece rates of posterior distributios for Browia semimartigale models. Beroulli 2 26) [3] F. va der Meule ad H. va Zate, Cosistet oparametric Bayesia estimatio for discretely observed scalar diffusios. Beroulli 9 23) [4] L. Pazar ad H. va Zate, Noparametric Bayesia iferece for ergodic diffusios. J. Statist. Pla. Iferece 39 29) [5] O. Papaspiliopoulos, Y. Poker, G.O. Roberts ad A.M. Stuart, Noparametric estimatio of diffusios: a differetial equatios approach. Biometrika 99 22) [6] G.A. Pavliotis, Y. Poker ad A.M. Stuart, Parameter estimatio for multiscale diffusios: a overview. Statistical Methods for Stochastic Differetial Equatios. Vol. 24 of Moogr. Statist. Appl. Probab. CRC Press, Boca Rato, FL 22) [7] Y. Poker, A.M. Stuart ad J.H. va Zate. Posterior cosistecy via precisio operators for oparametric drift estimatio i SDEs. Stoch. Process. Appl ) [8] L. Schwartz, O Bayes procedures. Z. Wahrscheilichkeitstheorie ud Verw. Gebiete 4 965) 26. [9] P. Soulier, Noparametric estimatio of the diffusio coefficiet of a diffusio process. Stochastic Aal. Appl ) [2] A.W. va der Vaart, Asymptotic Statistics. Vol. 3 of Cambr. Ser. Stat. Probab. Math. Cambridge Uiversity Press, Cambridge 998). [2] A.W. va der Vaart ad J.H. va Zate, Rates of cotractio of posterior distributios based o Gaussia process priors. A. Statist a) [22] A.W. va der Vaart ad J.H. va Zate, Reproducig kerel Hilbert spaces of Gaussia priors. Pushig the Limits of Cotemporary Statistics: Cotributios i Hoor of Jayata K. Ghosh. Vol. 3 of Ist. Math. Stat. Collect. Ist. Math. Statist., Beachwood, OH 28) [23] S. Walker, O sufficiet coditios for Bayesia cosistecy. Biometrika 9 23) [24] S. Walker, New approaches to Bayesia cosistecy. A. Statist ) [25] L. Wasserma, Asymptotic properties of oparametric Bayesia procedures. Practical Noparametric ad Semiparametric Bayesia Statistics. Vol. 33 of Lect. Notes Statist. Spriger, New York 998) [26] H. va Zate, Noparametric Bayesia methods for oe-dimesioal diffusio models. Math. Biosci. 23). Available o

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

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

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

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

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

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

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

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

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

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

Sequences and Series

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

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

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

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

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

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

Research Article The Probability That a Measurement Falls within a Range of n Standard Deviations from an Estimate of the Mean

Research Article The Probability That a Measurement Falls within a Range of n Standard Deviations from an Estimate of the Mean Iteratioal Scholarly Research Network ISRN Applied Mathematics Volume 0, Article ID 70806, 8 pages doi:0.540/0/70806 Research Article The Probability That a Measuremet Falls withi a Rage of Stadard Deviatios

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Estimation of Population Variance Utilizing Auxiliary Information

Estimation of Population Variance Utilizing Auxiliary Information Iteratioal Joural of Statistics ad Systems ISSN 0973-675 Volume 1, Number (017), pp. 303-309 Research Idia Publicatios http://www.ripublicatio.com Estimatio of Populatio Variace Utilizig Auxiliary Iformatio

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

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

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

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

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

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

NORMALIZATION OF BEURLING GENERALIZED PRIMES WITH RIEMANN HYPOTHESIS

NORMALIZATION OF BEURLING GENERALIZED PRIMES WITH RIEMANN HYPOTHESIS Aales Uiv. Sci. Budapest., Sect. Comp. 39 2013) 459 469 NORMALIZATION OF BEURLING GENERALIZED PRIMES WITH RIEMANN HYPOTHESIS We-Bi Zhag Chug Ma Pig) Guagzhou, People s Republic of Chia) Dedicated to Professor

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

STAT 135 Solutions to Homework 3: 30 points

STAT 135 Solutions to Homework 3: 30 points STAT 35 Solutios to Homework 3: 30 poits Sprig 205 The objective of this Problem Set is to study the Stei Pheomeo 955. Suppose that θ θ, θ 2,..., θ cosists of ukow parameters, with 3. We wish to estimate

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

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

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

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

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

On Regret and Options - A Game Theoretic Approach for Option Pricing

On Regret and Options - A Game Theoretic Approach for Option Pricing O Regret ad Optios - A Game Theoretic Approach for Optio Pricig Peter M. DeMarzo, Ila Kremer ad Yishay Masour Staford Graduate School of Busiess ad Tel Aviv Uiversity October, 005 This Revisio: 9/7/05

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

On the Composition of 2-Prover Commitments and Multi-Round Relativistic Commitments. CWI Amsterdam

On the Composition of 2-Prover Commitments and Multi-Round Relativistic Commitments. CWI Amsterdam O the Compositio of 2-Prover Commitmets ad Multi-Roud Relativistic Commitmets Serge Fehr Max Filliger CWI Amsterdam Warig My talk: maily classical Still iterestig for you: topic (relativistic commitmets)

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

ESTIMATING THE VOLATILITY OCCUPATION TIME VIA REGULARIZED LAPLACE INVERSION

ESTIMATING THE VOLATILITY OCCUPATION TIME VIA REGULARIZED LAPLACE INVERSION Ecoometric Theory, 32, 216, 1253 1288. doi:1.117/s266466615171 ESTIMATING THE VOLATILITY OCCUPATION TIME VIA REGULARIZED LAPLACE INVERSION JIA LI Duke Uiversity VIKTOR TODOROV Northwester Uiversity GEORGE

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

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

Section Mathematical Induction and Section Strong Induction and Well-Ordering

Section Mathematical Induction and Section Strong Induction and Well-Ordering Sectio 4.1 - Mathematical Iductio ad Sectio 4. - Strog Iductio ad Well-Orderig A very special rule of iferece! Defiitio: A set S is well ordered if every subset has a least elemet. Note: [0, 1] is ot well

More information

Research Article The Average Lower Connectivity of Graphs

Research Article The Average Lower Connectivity of Graphs Applied Mathematics, Article ID 807834, 4 pages http://dx.doi.org/10.1155/2014/807834 Research Article The Average Lower Coectivity of Graphs Ersi Asla Turgutlu Vocatioal Traiig School, Celal Bayar Uiversity,

More information

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models Optimally Thresholded Realized Power Variatios for Lévy Jump Diffusio Models José E. Figueroa-López Jeffrey Nise February 5, 1 Abstract: Thresholded Realized Power Variatios (TPV are oe of the most popular

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

AN INTRODUCTION TO STOCHASTIC CALCULUS AND BLACK-SCHOLES OPTION PRICING

AN INTRODUCTION TO STOCHASTIC CALCULUS AND BLACK-SCHOLES OPTION PRICING AN INRODUCION O OCHAIC CALCULU AND BLACK-CHOLE OPION PRICING ZACHRY WANG Abstract. his paper is a expositio of the mathematics behid the Black- choles model of pricig a Europea optio. After we briefly

More information

Marking Estimation of Petri Nets based on Partial Observation

Marking Estimation of Petri Nets based on Partial Observation Markig Estimatio of Petri Nets based o Partial Observatio Alessadro Giua ( ), Jorge Júlvez ( ) 1, Carla Seatzu ( ) ( ) Dip. di Igegeria Elettrica ed Elettroica, Uiversità di Cagliari, Italy {giua,seatzu}@diee.uica.it

More information

A STOCHASTIC GROWTH PRICE MODEL USING A BIRTH AND DEATH DIFFUSION GROWTH RATE PROCESS WITH EXTERNAL JUMP PROCESS *

A STOCHASTIC GROWTH PRICE MODEL USING A BIRTH AND DEATH DIFFUSION GROWTH RATE PROCESS WITH EXTERNAL JUMP PROCESS * Page345 ISBN: 978 0 9943656 75; ISSN: 05-6033 Year: 017, Volume: 3, Issue: 1 A STOCHASTIC GROWTH PRICE MODEL USING A BIRTH AND DEATH DIFFUSION GROWTH RATE PROCESS WITH EXTERNAL JUMP PROCESS * Basel M.

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

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

An Improved Estimator of Population Variance using known Coefficient of Variation

An Improved Estimator of Population Variance using known Coefficient of Variation J. Stat. Appl. Pro. Lett. 4, No. 1, 11-16 (017) 11 Joural of Statistics Applicatios & Probability Letters A Iteratioal Joural http://dx.doi.org/10.18576/jsapl/04010 A Improved Estimator of Populatio Variace

More information

5 Decision Theory: Basic Concepts

5 Decision Theory: Basic Concepts 5 Decisio Theory: Basic Cocepts Poit estimatio of a ukow parameter is geerally cosidered the most basic iferece problem. Speakig geerically, if θ is some ukow parameter takig values i a suitable parameter

More information

Calculation of the Annual Equivalent Rate (AER)

Calculation of the Annual Equivalent Rate (AER) Appedix to Code of Coduct for the Advertisig of Iterest Bearig Accouts. (31/1/0) Calculatio of the Aual Equivalet Rate (AER) a) The most geeral case of the calculatio is the rate of iterest which, if applied

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

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