Pure Jump Lévy Processes and Self-decomposability in Financial Modelling

Size: px
Start display at page:

Download "Pure Jump Lévy Processes and Self-decomposability in Financial Modelling"

Transcription

1 Pure Jump Lévy Processes and Self-decomposabiliy in Financial Modelling Ömer Önalan Faculy of Business Adminisraion and Economics, Marmara Universiy, Beykoz/Isanbul/ Turkey omeronalan@marmara.edu.r

2 Ouline Inroducion Levy Processes Self-decomposable Laws Ornsein-Uhlenbeck Processes Normal Inverse Gaussian Disribuion Modelling of Reurns Applicaion of Real Daa Conclusion

3 1. Inroducion The basic models of modern finance are based on he assumpion of normaliy of asse reurns. Bu a lo of empirical sudies are shown ha he assumpion which observaions are normally disribued is a poor approximaion in he real world. Because he reurns have feaures such as jumps, semi-heavy ails and asymmery. In radiional diffusion models, price movemens are very small in shor period of ime. Bu in real markes, prices may be show big jumps in shor ime periods. 3

4 Inroducion When price process model include he jumps, he perfec hedging is imposable. In his case, marke paricipans can no be hedge risks by using only underlying asses. For hese reasons, diffusion models used in finance is no a sufficien model. A good model mus be allow for disconinuiies and jumps in price process. Lévy processes is a valuable ool in financial modeling, hey provide a good fi wih reel daa. 4

5 Inroducion In probabiliy heory, a Lévy process, named afer he French mahemaician Paul Levy, is any coninuous-ime sochasic process ha sars a 0, admis cadlag modificaion and has "saionary independen incremens. The mos well-known examples are he Wiener process and he Poisson process. Lévy process is a simple Markov process wih jumps ha allow us o capure asse reurns wihou he necessiy of inroducing exreme parameer values 5

6 Inroducion Lévy models are no adequaely fiing implied volailiy surfaces of equiy opions across boh srike and mauriy. The incremens of addiive process presen a flexible model. These processes may be obained wih self-decomposable disribuions. 6

7 Inroducion A law in class of self-decomposable laws can be decomposable ino he sum of a scaled down version hemselves and an independen erm. The class of self-decomposable disribuions is obained as a limi laws of a sequence ha independen and suiably normalized 7

8 Inroducion The properies of he reurn disribuions depend on lengh of reurn inerval. Log reurns is aken monhly are reasonably represened by a normal disribuion. If one be dealing wih ick daa, hen reurn disribuions may have heavy ails. The aim of hese paper, o review o pure jump Lévy process wih self-decomposable disribuions in financial modeling and esing for he presence of a Brownian moion componen and discriminaing beween finie or infinie variaion 8

9 . LÉVY PROCESSES ( ) 0 Definiion.1 A cadlag sochasic process on a filered probabiliy spaces ( Ω, I, ( I ) 0, P) Lévy processes if he following condiions holds = is defined is called a i) 0 = 0 ii) has independen incremens of he pas i.e. independen of { } u ; u s for 0 s s iii) has saionary incremens i.e s has he same disribuion wih s 9

10 Lévy Processes iv) is sochasically coninuous, ( ε ) 0 lim P = ε > 0 k + k Process has cadlag pahs, i.e, is a.s. righ-coninuous wih lef limis. 10

11 Basic Examples Brownian moion B ~ N( μ, σ ) Poisson Process ~ Pois( λ ) λ ( 0, ) N 11

12 Independen and Saionary Incremens A coninuous-ime sochasic process assigns a random variable o each poin 0 in ime. I is a random funcion of. The incremens of such a process are he differences s beween is values a differen imes < s. The incremens of a process is called independen if incremens s and u v are independen random variables,whenever he wo ime inervals do no overlap. The incremens is called saionary, if he probabiliy disribuion of any incremen s depends only on he lengh s of he ime inerval; in he oher words, incremens wih equally long ime inervals are idenically disribued. 1

13 Independen and Saionary Incremens Lévy processes can be viewed as coninuous ime random walks. Lévy process form building blocks for boh Markov processes and semi-maringales There is one o one relaionship beween Lévy process and infiniely divisible disribuions. 13

14 Independen and Saionary Incremens ( ) 0 Theorem.1 Le be a Lévy process. Then has = an infiniely divisible disribuion F for every. Conversely if is an infiniely divisible disribuion, F hen here exis a = ( ) Lévy process, such ha disribuion of is 0 1 given F ( ) 0 The disribuion of a Lévy process = is compleely deermined by any of is marginal disribuions 14

15 Independen and Saionary Incremens The probabiliy disribuions of he incremens of any Lévy process are infiniely divisible There is a Lévy process for each infiniely divisible probabiliy disribuion. A coninuous process wih independen incremens is Gaussian. A linear combinaion of independen Levy processes is also a Levy process. 15

16 .1 Infiniely Divisibiliy ϕ ( u) Le be a characerisic funcion of a random variable. We say ha he law of random variable is infiniely divisible, if we can always find anoher characerisics funcion. u such ha, ( u ) [ ϕ ( )] n n u ϕ = ϕ n ( ) ( ) If ~ ϕ u is infiniely divisible, hen for all n Ν here exis i.i.d random variables ( 1 n ) ( 1 n ) 1,..., n 16

17 Infiniely Divisibiliy ( ) All disribued as u and such ha, ϕ n d = ( 1 n ) ( 1 n ) n In he oher words, random variable always is decomposable ino he sum of an arbirary finie number of i.i.d. random variables 17

18 . Jump Diffusion Model These models have jumps and random evoluion beween he jump imes, where W sandard Brownian moion, compound Poisson process. = μ + σ W + N i = 1 = Y N i= 1 [ ] iu u σ + ( iux exp iu μ λ e 1) F ( dx) E e i Y i 18

19 .3 The Probabilisic Properies of Lévy Processes Le be a Lévy process, we consider following characerisic funcion, ( u) = exp( ) ϕ iu [ ] For ime inerval 0, Δ = i i 1 = by he assumpion of independen and saionary incremens; n = ( ) ( ) 1 0 n n 1 d = n n 19

20 The Probabilisic Properies of Lévy Processes ϕ ( u) = E exp ( iu ) = E exp ( iun ) = ( E exp ( iu ) n n n for n =, ( ) ( ( )) ψ ( u u E exp iu e ) ϕ = = 1 ψ ( u ) = log E exp( iu ) 1 0

21 The Probabilisic Properies of Lévy Processes The every Lévy process can be represened in he following form, = μ + σ W + Z where Z is a jump process wih infiniely many jumps. 1

22 Lévy-Khinchine Formula ( ) o The characerisic funcion of a Lévy process = is given by he, Lévy- Khinchine formula, E e σ u = exp iub + ( ) iux e 1 iux1 v( dx) iu 1 { x } ( ) where, b,σ,ν called generaing riple, 0 ν does no have mass on

23 Lévy-Khinchine Formula ν ({ 0 }) = 0 and saisfies he following inegrabiliy condiion; R min ( x,1) ν ( dx ) < 3

24 Lévy-Io represenaion Lévy-Io represenaion describes he pah srucure of a Lévy process. A Lévy process can be represened in he following way, = μ + σw + Δ s1 s ε 0 s s s { } { } { } ( ) Δ 1 + lim Δ s1 ε Δ 1 x1 ε x 1 ν dx b ( W ) 0 0, σ 0, is a sandard Brownian moion. 4

25 Δ s = s s is called Lévy measure of Lévy measure : denoes he jumps a imes. ν { } ( ) Thus, ν dx is he inensiy of jumps of size x. We assumed ha he pahs of Lévy process is defined over a finie inervals 0,. [ ] [ ] As a consequence he sum of he jumps in ime inerval 0, wih absolue jump size bigger han 1 is a finie sum for each pah. 5

26 Lévy measure ( A) = E( card{ s [ 0,1] : Δ 0, Δ A} ) ν A Β( R) In he oher words, ν A is he average number of jumps of process [ ] ( ) in ime inerval 0,1 whose sizes fall in A s In general in a Lévy process, he frequency of he big jumps deermines exisence of momens of process. The fine srucure of he pahs of he process can be read of he frequency he small jumps. s 6

27 .4 Finie Aciviy A sochasic process has finie aciviy if all pahs of process have a finie number of jumps on every finie inerval i.e., ν ( R) ν ( dx) < = R If all pahs of process has an infinie number of jumps on every finie inerval i.e., ν ( R) = Lévy process has infinie aciviy. i.e. An infinie mass accumulaes around he origin. In a such case, small jumps of occur infiniely. 7

28 .5 Finie Variaion ( ) 0 Le be a Lévy process. = = 0 dx σ { } ( ) i) If and x 1 x ν < and process have finie variaion i.e., Δ s 1 < min( x,1) ν ( dx) < s i f and only if R ii) If σ 0 or hen process 1 dx { x 1 } x ν ( ) = have infinie variaion. 8

29 Subordinaors A Lévy process Z is called a subordinaor, if i has a.s. nondecreasing pahs. such a process canno have a diffusion componen. Given a Lévy process and an independen subordinaor Z, he process Y = Z is called a subordinaed Levy process. 9

30 3. SELF-DECOMPOSABLE LAWS There is a very close connecion he laws of class L and selfdecomposable laws. For hese reason, firs ime, we describe self-decomposable laws and relaed. 3.1 Laws of Class L The name of class L firs ime used by Khinchine, any random variable in class L is infiniely divisible ([30],Theorem 9.3), [6]. The infiniely divisible laws are he limi laws of riangle arrays, where arrays of independen random variables which individually negligible. Infinie divisibiliy is preserved under affine ransformaions. 30

31 ( ; n = 1,,... ) Y n 3.1 Laws of Class L Le is a sequence of independen random variables and denoes heir sum. Suppose ha, here exis cenering consans and scaling consans > 0 such ha he disribuion of b S + converges o he disribuion of some random n n a n b n n S n = Y i i= 1 a n R variable 31

32 3.1 Laws of Class L Then we say ha, random variable is a member of class L. As explained above, we shorly can say ha, if a random variable has same disribuion as he limi of some sequence of normalized sums of independen random variables, random variable has a disribuion of class L. 3

33 Self-decomposabiliy Definiion3.1 (Self-decomposabiliy) We suppose ha ϕ( u) is he characerisic funcion of a law. We say ha his law is self-decomposable and wrie if for c 0, 1 and ( u) ( cu) ϕ ( u) ϕ =ϕ ϕ SD ( ) u R c ϕ ( u) ( u) where is a characerisic funcion. is uniquely c deermined. For example normal case, ϕ c 33

34 For example normal case, Self-decomposabiliy ϕ ( u) = exp ( c) u c 1 We can resae his definiion for random variables as follows, 1 The disribuion of a random variable is self-decomposable, for any consan c ( 0,1) we can a find independen random variable such ha, c 34

35 Self-decomposabiliy d where = denoes equaliy in disribuion d = c + According o ([30], heorem 15.3) a random variable has a disribuion of class L if and only if he law of he is self-decomposable. The class of self-decomposable disribuions is a proper subse of infiniely divisible disribuions. A random variable or is probabiliy μ or probabiliy densiy f is self-decomposable,if he corresponding characeriic funcion is in he class SD. c 35

36 Self-decomposabiliy Self-decomposable laws arise as marginal laws in auoregressive ime series models c = 1 + ε The Lévy measure of he self-decomposable laws is absoluely coninuous wih following densiy form, ( ) ν dx = k( x) dx x 36

37 Self-decomposabiliy SD is closed under affine mappings, i.e., for all reals a and b if ϕ SD hen, e ibu ϕ a u SD ( ) SD If here exiss a unique Lévy process Z such ha d = e 0 s dz s 37

38 Self-decomposabiliy ϕ and ψ and Z1 Le denoe he characerisic funcions of respecively, hen ψ o dv () = ψ ( v) v log ϕ log () = exp ( logϕ() ), 0, ψ ( 0) = 1 38

39 Self-decomposabiliy ( x) k (, x) where, increasing for and decreasing for x,. ( ) The densiy of self-decomposable disribuions is unimodal. ( ) 0 = Le be a Lévy process. ( ) ( ) The process 1 self-decomposable if and only if selfdecomposable for every > 0. 39

40 Self-decomposabiliy The characerisic funcion of self-decomposable laws have following form, E u σ = R ( ) k ( x ) [ ] iu iux e exp iuμ + e 1 iux1 { x 1} < dx x μ R ( x) 0 Where,, and ( ) ( ) min 1, dx < k R x k x x 40

41 Self-decomposabiliy A self-decomposable random variable is he value a uni ime of some pure jump Lévy processes which sample pahs have bounded variaion. When he levy densiy inegraes x in he region x < 1, for μ = x < 1 x k ( x ) dx x 41

42 Self-decomposabiliy The characerisic funcion of, E [ ] ( ) iu iux k( x) e 1 exp e = x dx We can consider ha, he reurns are he sum of a suiable number of approximaely independen random variables. Furhermore he reurn disribuion is a limi disribuion. Self-decomposable disribuions can be consider as candidae for he uni period disribuion of asse reurns 4

43 Self-decomposabiliy For example, he Laplace (double exponenial) random variable has he probabiliy densiy, x ( x) = e, x R Is characerisic funcion is equal o, f 1 ϕ SD ϕ ( ) ( iux u = exp e 1) = exp u e x ( ivx e 1) 0 x e dx x dx dv v 43

44 3. Self-Decomposable Laws and Addiive Processes Addiive processes are obained from Lévy processes by relaxing he condiion of saionary of incremens A levy process is a addiive process wih requiremen ha he incremens mus be saionary. Addiive processes have markovian propery hence we can obain complee law of process, if we know he laws of incremens The rajecories of he addiive process can make jumps. A law is self-decomposable if and only if i is he law of uni ime of an addiive processes. 44

45 3.3 Self-similariy and Self-decomposabiliy ( ) 0 = A sochasic process is called self-similar for any given, if we can find such ha, a 0 d b R ( ) ( ) a 0 = b 0 Where, b can be expressed as b H = a, H 0 In he oher words, we say ha one sochasic process is selfsimilar such ha, he change an in ime scale can be compensaed by a corresponding change in he spaces scale. 45

46 Self-similariy and Self-decomposabiliy The H is called Hurs exponen. The connecion beween selfdecomposable laws and self-similar addiive process is given by [31]. A law is self-decomposable if and only if i is he law a uni ime of a self-similar addiive process. ( ) ϕ u Le be a characerisic funcion of a law. hen i can ake a characerisic funcion of Lévy process as follow, 46

47 Self-similariy and Self-decomposabiliy ( u) [ ϕ( u) ] N ϕ = Where, N is ime scale, if ϕ u is a infinie divisible, ϕ ( u) is a characerisic funcion. ( ) Now we describe a new funcion as follows, ψ k, H ( θ ) H H k h = ϕ θ ϕ θ N N 1 47

48 Self-similariy and Self-decomposabiliy I is a characerisic funcion if and only if is selfdecomposable ([30],p.99), ([7],p.1884). The saionary process and he self-similar process are relaed by using firs Lamperi represenaion 48

49 Self-similariy and Self-decomposabiliy Proposiion 3.1: Le = he a new process defined by ( ) R be a saionary process, H Z = log, > 0, Z0 = 0andH > is a self-similar process wih index H>0. 0 Conversely if he process Z is self-similar wih H>0, hen process defined by is a saionary process ( ) ( H ) R = e Z e R 49

50 Self-similariy and Self-decomposabiliy ( ) 0 In general a Lévy process = may be consruced from he a H-self-similar process in accordance wih Z = e s 1 H dz s 50

51 4. ORNSTEIN-UHLENBECK PROCESSES ( ) 0 = Ornsein-Uhlenbeck process (OU) saisfies following differenial equaion, d = λ d + dz λ = ( Z ) 0 where > 0 and Z is a Lévy process. Homogenous Lévy process Z [ ( )] < E log 1+ Z 1 has ha propery 51

52 Z Ornsein-Uhlenbeck Processes If is non-gaussian Lévy process, above differenial equaion has a unique soluion. I has following form, = e λ 0 + e λ ( s ) dz s 0 where λ R and 0 is iniial sae. The Ornsein-Uhlenbeck process is a Markov process and possesses a saionary regime under a mild log-momen condiion. 5

53 Ornsein-Uhlenbeck Processes ν 0 Le be Lévy measure such ha saisfies loginegrabiliy condiiion, x 1 log x ν 0 ( dx) < If a Lévy process saisfies above condiion, hen i has a selfdecomposable disribuion. 53

54 Ornsein-Uhlenbeck Processes O U process has he following, moving average represenaion: λ () ( s = e ) dz O U process has following correlaion funcion; s Corr [ ] λ s = e s, 54

55 Ornsein-Uhlenbeck Processes O U process can be defined as a Lamperi represenaion of he Brownian moion Z ( ) λ ( λ = e Z e ) 55

56 Ornsein-Uhlenbeck Processes ( ) 0 If an = OU process is saionary, i s characerisic funcion has following form, ϕ ( ) ( λ u = ϕ u e ) ϕ ( u) This denoes ha, he marginal disribuion of is selfdecomposable. ϕ c ( ) ( ( ) ( λ u = exp κ u κ u e ) c ( ) 56

57 Ornsein-Uhlenbeck Processes ( u) κ ( ) Where, is cumulan of and denoed by as Le Z be he background driving Lévy process of OU processes and ( u) = logϕ( u) Where v Lévy measure of Z process. κ log { x > } x ν ( dx) < 57

58 Ornsein-Uhlenbeck Processes Then he law of converge owards a selfdecomposable law as Characerisic funcion of his law is given by, ζ u = 0 ( ) ( λ s u exp φ u e ) φ 1 is he characerisic exponen of process ζ ds 58

59 Ornsein-Uhlenbeck Processes We can say ha he limi disribuion of an OU process is self-decomposable. The disribuion of a random variable is selfdecomposable if d = e 1 0 s dz s 1 ζ 59

60 Sochasic Volailiy In Black-Scholes model,volailiy is sandar deviaion of he reurn.volailiy is unobserved and has o be esimaed from pas prices. Volailiy is no consan, i can be ake as sochasic. For example, in BNS model,reurn process, d = [ ] μ + β σ d + σ dw Volailiy process may be as an OU process, σ dσ = λ σ d + dz The process is defined as a saionary markov process λ 60

61 Sochasic Volailiy Where Z is an increasing Levy process and W an independen Brownian moion ( ) 0 The price process has coninuous sample pahs. ( ) 0 The volailiy process σ is a non-gaussian Ornsein- Uhlenbeck process ( σ ) 0 The volailiy process is no direcly observable Finally,sochasic volailiy (SV) models are used o capure he impac of ime-varying volailiy on financial markes. 61

62 5. NORMAL INVERSE GAUSSIAN DISTRIBUTION Normal Inverse Gaussian Disribuion was inroduced by [5]. This disribuion is used in [9], [8], [9] o model equiy reurns. NIG disribuion has semi-heavy ails and a special case wih λ = 1 parameer in he class of Generalized Hyperbolic Disribuions (GH). GH class are infiniely divisible and self-decomposable. 6

63 Normal Inverse Gaussian Disribuion ( ) The probabiliy densiy funcion of he NIG α, β, μ, δ is defined by as follows, f NIG α δ ( x, α, β, μ, δ ) = exp δ α β + β ( x μ) π. K 1 α δ δ + + ( x μ) ( x μ) wih parameers, (, ) < < δ > 0, μ and 0 β 63

64 Normal Inverse Gaussian Disribuion K ν The funcion is modified Bessel funcion of hird kind wih index. ν K ν 1 = u ν 1 ( x) u exp x u du ( x > 0) 0 K 1 π x x ( x) =. e K ν ( x) = K ( x) ν 64

65 Normal Inverse Gaussian Disribuion The characerisic funcion of NIG disribuion, ϕ NIG = + iuμ ( u) exp δ α β δ α ( β + iu) The normal Inverse Gaussian disribuion has semi-heavy ails,i.e. f NIG ( x) ~ cons x 3 e ( α x + β x) as x 65

66 66 Normal Inverse Gaussian Disribuion The cenral momens of a random variable are [ ] 1 β α β δ μ + = E [ ] γ α δ β α β δ β α δ = + = Var

67 67 Normal Inverse Gaussian Disribuion [ ] 1 3 β α δ α β = Skew [ ] β α δ α β + + = Kurosis

68 Normal Inverse Gaussian Disribuion The Lévy measure of NIG has following form, α δ ( dx) =. exp{ β x} K ( α x )dx ν NIG 1 π x The probabiliy densiy funcion of NIG disribuion, f NIG α ( ) ( ) α + β x ~ exp δ α β e π x δ 3 x as x ± 68

69 5.1 Parameer Esimaion of NIG Disribuion Momen Esimaors These esimaes can be used as parameer esimaes which obained in real daa. Le m i, i = 1,,3,4 be sample momen of real daa. ˆ ρ = 3 m 3 m m ˆ μ ˆ β ˆ δ ˆ ρ = m 1 ˆ ( β = m m ρ ) ˆ 69

70 Parameer Esimaion of NIG Disribuion ˆ δ = m ˆ β ˆ ρ 3 + ˆ ρ ˆ ˆ α = ˆ ρ + β 70

71 71 Parameer Esimaion of NIG Disribuion ( ) 1 + = α β α β δ μ E ( ) 3 1 = α β α δ Var ( ) = α β α β α Skew ( ) = α β α β α Kur

72 Parameer Esimaion of NIG Disribuion β 3 ( Skew ( )) α = Kur ( ) 1 + β 4 α Skew 3 ( ) Ku( ) 5 7

73 5. The Simulaion of Normal Inverse Gaussian Disribuions The algorihm can be sae as following [9], [9], [10]. Generae from Y from N(0,1) Generae from Z from IG δ, α β Se = μ + β Z + Z Y ( ) Z The simulaion of Inverse Gaussian disribuion Sample Y from N(0,1) and se V = Y Se K η V 1 = η + 4η δ V + η δ δ V 73

74 74 Parameer Esimaion of NIG Disribuion wih Se β α δ η = + + < + = K U K U K K Z η η η η η

75 Parameer Esimaion of NIG Disribuion The simulaion of V Generae U1 and U uniform random variables N N 1 [ lnu ] cos ( π U ) 1 = 1 1 [ lnu ] sin ( π U ) = 1 Boh of N1 and N have sandard normal disribuion are independen. 75

76 6. THE MODELLING OF RETURNS We can consider a general exponenial Levy process model for sock prices, S = S exp ( ) 0 wih a Levy process ( : 0) The disribuion of a Levy process is uniquely deermined by any of is one dimensional marginal disribuion So we may use he disribuion of 1 The disribuion of 1 is infiniely divisible and is characerisic funcion is given by Levy-Khinchine formula. 76

77 6. THE MODELLING OF RETURNS An infiniely divisible disribuion can be esimaed from he real daa which available for he paricular asse. The log reurns of he model have independen and saionary incremens along ime inervals of lengh 1. The pah properies of he levy process carry over o S. If is a pure jump levy process, hen S is alsa a pure jump process. 77

78 6. THE MODELLING OF RETURNS We consider log price changes, = ln S+ Δ ln S reflecs he muliplicaive characer of price changes. When we use daily price daa, ypically will have he value 1. In generally, he properies of reurn disribuions depends on he lengh of he reurn inerval. Δ Δ 78

79 Δ The Modelling of Reurns For long as a resul of cenral limi heorem, he reurns can be described wih a Gaussian disribuion bu for ick daa, he normal disribuion can no be suiable In his case, as a model Generalized hyperbolic disribuions and i s subse can be used These disribuions are infiniely divisible and selfdecomposable Levy markes give good fi o reurn disribuions. 79

80 The Modelling of Reurns A non-parameric hreshold echnique is proposed o es he inegraed variance which based on discree observed prices by [13]. We can deermine wheher jump ype process is suiable or no as a model for price process, using wih his echnique. 80

81 ( ) o The Modelling of Reurns is he Lévy process generaed by he normal inverse Gaussian disribuion ha was fied o he real daa. The incremens of price process a long ime inerval 1, are disribued according o he NIG disribuion. +1 S The securiy price denoes he presen price of securiy The risk one day a head is calculaed as follow, 81

82 The Modelling of Reurns Var ( ) [ ( )] S S = S Var exp = S Var [ exp( )] 1 1 The disribuion of is he fied NIG 8

83 The Modelling of Reurns Levy processes are semimaringales. As a consequence, he following sochasic inegrals, 0 f () s d and A() s d s o are well definie for L 83

84 The Modelling of Reurns In classic Black-Scholes model, price process is soluion of below sochasic differenial equaion ds = S μ d + S σ dw, S 0 > 0 The soluion, S 1 exp μ d = S 0 σ + σ W 84

85 6.1 Tes Saisics We wan o make inference on o unobserved volailiy process from he observed price process, Blumenhal-Geoor, index is described as follows, p α inf p 0, x ν ( dx) < α [ 0,] = { x 1} This index measures aciviies of small jumps of Lévy process. For jump ype Lévy process, α = 1 For example, Normal inverse Gaussian moion has infiniely variaion and 85

86 Tes Saisics An infinie aciviy process wih Blumenhal-Geoor, index α < 1 has pahs wih finie variaion while henhesamplepahshaveinfinie variaion α > 1 An esimaor for he inegraed variance given by 0 σ ds s This esimaor is consisen in he presence of jumps and characerize is asympoic behavior via a cenral limi herem. V denoes hreshold esimaor for inegraed variance Under he very general assumpions, realised m-power variaion converge o quadraic variaion. 86

87 The esimaion of quadraic variaion(qv) δ Le denoes a ime periods beween observaions QV process is defined as following, [ δ ] ( δ ) = δ i δ ( i 1) lim δ 0 i= 1 ( ) ( ) σ ds δ We used daily reurns o measure incremens of he quadraic variaion p 0 s 87

88 6.1.1 Tes for he Presence of a Coninuous Maringale Componen of a Levy Model Tes procedure as follows; Choice a coefficien We choose a hreshold β [ 0.5,1 ] β r ( k) = k Se Δ i Y = Δ i + σ k Z i ( 0,1) Where as T k, Z i ~ N n = k 0 σ T is measured by as an annually, is known sandard deviaion, denoes log-reurns. Δ i 88

89 Tes for he Presence of a Coninuous Maringale Componen of a Lévy Model Null Hypohesis: H 0 : σ 0 Tes saisics is calculaed as following Vˆ k n Δ i= 1 ( iy ) 1 {( Δ Y ) r( k )} = i ˆ Qk = ( Δ iy ) 1 {( Δ Y ) r( k )} 1 3k n i= 1 4 i T k = ˆ k σ V k Qˆ n. k k N ( 0,1) 89

90 Tes for he Presence of a Coninuous Maringale Componen of a Lévy Model { T k > 1.96 } α If is near 0.05, we accep null hypohesis. P Under he alernaive hypohesis : σ 0, above es saisic diverges o +. { } α 05 1 If P > 1.96 >> 0., we rejec he null hypohesis T k H 90

91 6.1. The es wheher he jump componen has finie variaion Null Hypohesis : H : α 1 0 < Δ i M ˆ = Δ 1{ ( ) ( )} k Z, Z ~ N ( 0,1 ) i +σ Δ > r k i i i Vˆ k n Δ i= 1 ( ) i Mˆ 1 ( Δ M ) r( k ) = i Q ( ) k = ΔiM 1 ( Δ M ) r( k) { } ˆ 1 3k n i= 1 ˆ 4 { } i 91

92 The es wheher he jump componen has finie variaion T α k = ˆ k V α n. k kqˆ k N ( 0,1) { } α If P > 1.96 is near 0.05, we accep null. T k 9

93 6.1.3 The es of Presence Jumps The variance of a Lévy processes is esimaed as following form, S = n 1 i= 0 ( ) i+ 1 [ ] Where each i is a division of inerval 0, for each. and = n. Δ n Ν i = = Δ i i 1 i n i 93

94 94 The Tes of Presence Jumps ( )( ) = + = n i i i i i S = Δ Δ = n i n n i i n S

95 The Tes of Presence Jumps H : Null hypohesis: is a coninuous Lévy process 0 Z n S n S K S K 4 4 = K = ( π ) 0. 5 Z n > z α If, we rejec he null hypohesis z α where is a chosen quanile. 95

96 7. Applicaion o Real daa We consider he ime series of NIKKEI 5 in Japan reurns from o and ISE Compound 100 index in Turkey reurns from o , 75 We use daily daa and k = 1 5 = 0, and r ( k) = k We divided reurn series such ha each group have 500 daily reurn for NIKKEI 5 index Similarly, ISE Compound 100 index reurns divided sub groups which include 300 daily reurn. 96

97 Applicaion o Real daa Index Level (JPY) Year FIGURE 1. Hisory of NIKKEI 5 Index in Japan

98 Applicaion o Real daa 15% 10% 5% 0% -5% -10% -15% Reurns (daily coninuosly compounded rae %) 005 Year FIGURE. Hisory of NIKKEI 5 Index reurns in Japan

99 Applicaion o Real daa 0,5 NIKKEI 5 0, Normal disribuion PDF Daily mean: 0.01% Daily sd. Deviaion: 1.3% 0,15 0,1 0, % -1% -10% -8% -6% -4% -% 0% % 4% 6% 8% 10% 1% 14% 16% FIGURE 3. PDF of Daily Reurns of Nikkei 5 from 198/01/04 o 005/09/30 99

100 Applicaion o Real daa 15% 10% 5% 0% -5% % -15% FIGURE 4. Hisory of ISE Composie 100 Index reurn in Turkey 100

101 Applicaion o Real daa Close FIGURE 5. Hisory of ISE Composie 100 Index in Turkey

102 Applicaion o Real daa 14,00% 1,00% 10,00% 8,00% Normal Frequency 6,00% 4,00%,00% 0,00% -7% -6% -5% -4% -3% -% -1% 0% 1% % 3% 4% 5% 6% 7% FIGURE 6. Frequency of ISE Composie 1000 Index in Turkey 10

103 Applicaion o Real daa ISE ,095 0, , , ,07818 Vˆk Qˆ k 0, ,1095 0, , , T k ( j) 0, , ,5644 9, , Table 1. The esing for he presence of a Brownian componen for ISE Composie 100 index price process 103

104 Applicaion o Real daa ISE Vˆk 0, , , ,03769 Qˆ k T k ( j) 0, , ,89E-06 0,018 4,334696, , ,95618 Table. The esing finie variaion of jump componen of ISE Composie 100 index price process 104

105 Applicaion o Real daa NIKKEI 5 Vˆk Qˆ k T k ( j) , ,3476 0, , , ,7848 0, ,0776 0, , , , ,6515 5, ,68 16, , ,37838 Table 3. The esing for he presence of a Brownian componen for NIKKEI 5 Index reurns 105

106 Applicaion o Real daa NIKKEI 5 Vˆk Qˆ k T k ( j) , ,0595 0, ,0498 0, , ,09E-09 0, , , , ,39E-08 0, , , , , ,01413 Table 4. The esing finie variaion of jump componen of NIKKEI 5 Index reurns 106

107 Applicaion o Real daa INDE Z n NIKKIE 5 179,0877 ISE Composie ,6736 Table 5. The esing he presence of jumps under he hypohesis ha price process is coninuous Lévy process. 107

108 Applicaion o Real daa Index NIKKEI 5 ISE Composie 100 αˆ βˆ μˆ δˆ 0, , , , , ,0078 0, ,00038 Table 6. Esimaed parameers of he normal inverse Gaussian disribuion 108

109 8. CONCLUSION The laws of self-decomposable disribuions class are be consiued as a limi laws of Lévy driven Ornsein- Uhlenbeck process These laws can be relaed wih general addiive process. The Lévy processes has a flexible srucure o model in real phenomena in finance such as heavy-ails, jumps and volailiy smile In his sudy, we focus on Levy process and i s incremens laws. Self-decomposable laws are sub class of infiniely divisible laws 109

110 CONCLUSION The family of self-decomposable laws always has boh of selfsimilariy and saionary of incremens. We esed wheher he pure jump Levy process is a suiable model or no for financial reurn series For his es, we used a hreshold esimaor of quadraic variaion by proposed by Con and Mancini We applied his es saisics o wo inernaional index reurn series. 110

111 CONCLUSION Our empirical analysis show ha, Nikkei 5 index model mus include Brownian moion componen The jumps of Nikkei 5 reurn index have infinie variaion ISE 100 reurn index model mus include Brownian moion componen The jumps of ISE 100 reurn index have infinie variaion 111

112 CONCLUSION A self-decomposable random variable has he same disribuion of as a scaled version of iself and an independen residual random variable. Self-decomposable processes are Levy processes wih jump arrival raes ha are decreasing in he jump size A self-decomposable random variable also has a disribuion of class L which means i can moivaed as a limi law wih more general scaling han he Gaussian limi law. 11

113 CONCLUSION This means ha self-decomposable processes can be moivaed as limi laws where he independen infleunces being summed are of differen orders of magniude Thus hey are appropriae building blocks for sochasic processes used o model financial markes. Finally, we can say ha, a good fi model for real reurn daa have o include boh of a Brownian componen and a infinie variaion componen. 113

6. Stochastic calculus with jump processes

6. Stochastic calculus with jump processes A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio

More information

Sample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen

Sample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen Sample Auocorrelaions for Financial Time Series Models Richard A. Davis Colorado Sae Universiy Thomas Mikosch Universiy of Copenhagen Ouline Characerisics of some financial ime series IBM reurns NZ-USA

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Chapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull

Chapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull Chaper 14 Wiener Processes and Iô s Lemma Copyrigh John C. Hull 014 1 Sochasic Processes! Describes he way in which a variable such as a sock price, exchange rae or ineres rae changes hrough ime! Incorporaes

More information

Richard A. Davis Colorado State University Bojan Basrak Eurandom Thomas Mikosch University of Groningen

Richard A. Davis Colorado State University Bojan Basrak Eurandom Thomas Mikosch University of Groningen Mulivariae Regular Variaion wih Applicaion o Financial Time Series Models Richard A. Davis Colorado Sae Universiy Bojan Basrak Eurandom Thomas Mikosch Universiy of Groningen Ouline + Characerisics of some

More information

Empirical Process Theory

Empirical Process Theory Empirical Process heory 4.384 ime Series Analysis, Fall 27 Reciaion by Paul Schrimpf Supplemenary o lecures given by Anna Mikusheva Ocober 7, 28 Reciaion 7 Empirical Process heory Le x be a real-valued

More information

CH Sean Han QF, NTHU, Taiwan BFS2010. (Joint work with T.-Y. Chen and W.-H. Liu)

CH Sean Han QF, NTHU, Taiwan BFS2010. (Joint work with T.-Y. Chen and W.-H. Liu) CH Sean Han QF, NTHU, Taiwan BFS2010 (Join work wih T.-Y. Chen and W.-H. Liu) Risk Managemen in Pracice: Value a Risk (VaR) / Condiional Value a Risk (CVaR) Volailiy Esimaion: Correced Fourier Transform

More information

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Stochastic Model for Cancer Cell Growth through Single Forward Mutation Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com

More information

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1 Vecorauoregressive Model and Coinegraion Analysis Par V Time Series Analysis Dr. Sevap Kesel 1 Vecorauoregression Vecor auoregression (VAR) is an economeric model used o capure he evoluion and he inerdependencies

More information

On a Fractional Stochastic Landau-Ginzburg Equation

On a Fractional Stochastic Landau-Ginzburg Equation Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS

A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS Theory of Sochasic Processes Vol. 14 3), no. 2, 28, pp. 139 144 UDC 519.21 JOSEP LLUÍS SOLÉ AND FREDERIC UTZET A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS An explici procedure

More information

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1 Chaper 5 Heerocedasic Models Inroducion o ime series (2008) 1 Chaper 5. Conens. 5.1. The ARCH model. 5.2. The GARCH model. 5.3. The exponenial GARCH model. 5.4. The CHARMA model. 5.5. Random coefficien

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

in Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology

in Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology Risk and Saey in Engineering Pro. Dr. Michael Havbro Faber ETH Zurich, Swizerland Conens o Today's Lecure Inroducion o ime varian reliabiliy analysis The Poisson process The ormal process Assessmen o he

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility Saisics 441 (Fall 214) November 19, 21, 214 Prof Michael Kozdron Lecure #31, 32: The Ornsein-Uhlenbeck Process as a Model of Volailiy The Ornsein-Uhlenbeck process is a di usion process ha was inroduced

More information

Unit Root Time Series. Univariate random walk

Unit Root Time Series. Univariate random walk Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he

More information

OBJECTIVES OF TIME SERIES ANALYSIS

OBJECTIVES OF TIME SERIES ANALYSIS OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging

More information

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0.

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0. Advanced Financial Models Example shee 3 - Michaelmas 217 Michael Tehranchi Problem 1. Le f : [, R be a coninuous (non-random funcion and W a Brownian moion, and le σ 2 = f(s 2 ds and assume σ 2

More information

Lecture 4: Processes with independent increments

Lecture 4: Processes with independent increments Lecure 4: Processes wih independen incremens 1. A Wienner process 1.1 Definiion of a Wienner process 1.2 Reflecion principle 1.3 Exponenial Brownian moion 1.4 Exchange of measure (Girsanov heorem) 1.5

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If

More information

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND Asymmery and Leverage in Condiional Volailiy Models Michael McAleer WORKING PAPER

More information

The Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite

The Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite American Journal of Operaions Research, 08, 8, 8-9 hp://wwwscirporg/journal/ajor ISSN Online: 60-8849 ISSN Prin: 60-8830 The Opimal Sopping Time for Selling an Asse When I Is Uncerain Wheher he Price Process

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Stochastic Modelling in Finance - Solutions to sheet 8

Stochastic Modelling in Finance - Solutions to sheet 8 Sochasic Modelling in Finance - Soluions o shee 8 8.1 The price of a defaulable asse can be modeled as ds S = µ d + σ dw dn where µ, σ are consans, (W ) is a sandard Brownian moion and (N ) is a one jump

More information

Regular Variation and Financial Time Series Models

Regular Variation and Financial Time Series Models Regular Variaion and Financial Time Series Models Richard A. Davis Colorado Sae Universiy www.sa.colosae.edu/~rdavis Thomas Mikosch Universiy of Copenhagen Bojan Basrak Eurandom Ouline Characerisics of

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

Multivariate Generalized Ornstein-Uhlenbeck Processes

Multivariate Generalized Ornstein-Uhlenbeck Processes Mulivariae Generalized Ornsein-Uhlenbeck Processes Ania Behme and Alexander Lindner Absrac De Haan and Karandikar [12] inroduced generalized Ornsein Uhlenbeck processes as one-dimensional processes (V

More information

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance. 1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard

More information

Properties of Autocorrelated Processes Economics 30331

Properties of Autocorrelated Processes Economics 30331 Properies of Auocorrelaed Processes Economics 3033 Bill Evans Fall 05 Suppose we have ime series daa series labeled as where =,,3, T (he final period) Some examples are he dail closing price of he S&500,

More information

Representation of Stochastic Process by Means of Stochastic Integrals

Representation of Stochastic Process by Means of Stochastic Integrals Inernaional Journal of Mahemaics Research. ISSN 0976-5840 Volume 5, Number 4 (2013), pp. 385-397 Inernaional Research Publicaion House hp://www.irphouse.com Represenaion of Sochasic Process by Means of

More information

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based

More information

Reliability of Technical Systems

Reliability of Technical Systems eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy

More information

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,

More information

Estimation of Poses with Particle Filters

Estimation of Poses with Particle Filters Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

13.3 Term structure models

13.3 Term structure models 13.3 Term srucure models 13.3.1 Expecaions hypohesis model - Simples "model" a) shor rae b) expecaions o ge oher prices Resul: y () = 1 h +1 δ = φ( δ)+ε +1 f () = E (y +1) (1) =δ + φ( δ) f (3) = E (y +)

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

Distribution of Estimates

Distribution of Estimates Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion

More information

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate. Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011 2229-12 School and Workshop on Marke Microsrucure: Design, Efficiency and Saisical Regulariies 21-25 March 2011 Some mahemaical properies of order book models Frederic ABERGEL Ecole Cenrale Paris Grande

More information

Appendix to Creating Work Breaks From Available Idleness

Appendix to Creating Work Breaks From Available Idleness Appendix o Creaing Work Breaks From Available Idleness Xu Sun and Ward Whi Deparmen of Indusrial Engineering and Operaions Research, Columbia Universiy, New York, NY, 127; {xs2235,ww24}@columbia.edu Sepember

More information

Testing the Random Walk Model. i.i.d. ( ) r

Testing the Random Walk Model. i.i.d. ( ) r he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp

More information

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy Answer 4 of he following 5 quesions. 1. Consider a pure-exchange economy wih sochasic endowmens. The sae of he economy in period, 0,1,..., is he hisory of evens s ( s0, s1,..., s ). The iniial sae is given.

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

Backward stochastic dynamics on a filtered probability space

Backward stochastic dynamics on a filtered probability space Backward sochasic dynamics on a filered probabiliy space Gechun Liang Oxford-Man Insiue, Universiy of Oxford based on join work wih Terry Lyons and Zhongmin Qian Page 1 of 15 gliang@oxford-man.ox.ac.uk

More information

Solutions to Odd Number Exercises in Chapter 6

Solutions to Odd Number Exercises in Chapter 6 1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b

More information

Basic notions of probability theory (Part 2)

Basic notions of probability theory (Part 2) Basic noions of probabiliy heory (Par 2) Conens o Basic Definiions o Boolean Logic o Definiions of probabiliy o Probabiliy laws o Random variables o Probabiliy Disribuions Random variables Random variables

More information

Financial Econometrics Introduction to Realized Variance

Financial Econometrics Introduction to Realized Variance Financial Economerics Inroducion o Realized Variance Eric Zivo May 16, 2011 Ouline Inroducion Realized Variance Defined Quadraic Variaion and Realized Variance Asympoic Disribuion Theory for Realized Variance

More information

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In

More information

Lecture 6: Wiener Process

Lecture 6: Wiener Process Lecure 6: Wiener Process Eric Vanden-Eijnden Chapers 6, 7 and 8 offer a (very) brief inroducion o sochasic analysis. These lecures are based in par on a book projec wih Weinan E. A sandard reference for

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross.

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross. Homework (Sas 6, Winer 7 Due Tuesday April 8, in class Quesions are derived from problems in Sochasic Processes by S. Ross.. A sochasic process {X(, } is said o be saionary if X(,..., X( n has he same

More information

Lecture 2 April 04, 2018

Lecture 2 April 04, 2018 Sas 300C: Theory of Saisics Spring 208 Lecure 2 April 04, 208 Prof. Emmanuel Candes Scribe: Paulo Orensein; edied by Sephen Baes, XY Han Ouline Agenda: Global esing. Needle in a Haysack Problem 2. Threshold

More information

Comparison between the Discrete and Continuous Time Models

Comparison between the Discrete and Continuous Time Models Comparison beween e Discree and Coninuous Time Models D. Sulsky June 21, 2012 1 Discree o Coninuous Recall e discree ime model Î = AIS Ŝ = S Î. Tese equaions ell us ow e populaion canges from one day o

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information

HPCFinance research project 8

HPCFinance research project 8 HPCFinance research projec 8 Financial models, volailiy risk, and Bayesian algorihms Hanxue Yang Tampere Universiy of Technology March 14, 2016 Research projec 8 12/2012 11/2015, Tampere Universiy of Technology,

More information

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems Essenial Microeconomics -- 6.5: OPIMAL CONROL Consider he following class of opimizaion problems Max{ U( k, x) + U+ ( k+ ) k+ k F( k, x)}. { x, k+ } = In he language of conrol heory, he vecor k is he vecor

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion

More information

Computer Simulates the Effect of Internal Restriction on Residuals in Linear Regression Model with First-order Autoregressive Procedures

Computer Simulates the Effect of Internal Restriction on Residuals in Linear Regression Model with First-order Autoregressive Procedures MPRA Munich Personal RePEc Archive Compuer Simulaes he Effec of Inernal Resricion on Residuals in Linear Regression Model wih Firs-order Auoregressive Procedures Mei-Yu Lee Deparmen of Applied Finance,

More information

Generalized Snell envelope and BSDE With Two general Reflecting Barriers

Generalized Snell envelope and BSDE With Two general Reflecting Barriers 1/22 Generalized Snell envelope and BSDE Wih Two general Reflecing Barriers EL HASSAN ESSAKY Cadi ayyad Universiy Poly-disciplinary Faculy Safi Work in progress wih : M. Hassani and Y. Ouknine Iasi, July

More information

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix Wha Ties Reurn Volailiies o Price Valuaions and Fundamenals? On-Line Appendix Alexander David Haskayne School of Business, Universiy of Calgary Piero Veronesi Universiy of Chicago Booh School of Business,

More information

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

Recent Developments in the Unit Root Problem for Moving Averages

Recent Developments in the Unit Root Problem for Moving Averages Recen Developmens in he Uni Roo Problem for Moving Averages Richard A. Davis Colorado Sae Universiy Mei-Ching Chen Chaoyang Insiue of echnology homas Miosch Universiy of Groningen Non-inverible MA() Model

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H. ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models

More information

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Methodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.

Methodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size. Mehodology. Uni Roo Tess A ime series is inegraed when i has a mean revering propery and a finie variance. I is only emporarily ou of equilibrium and is called saionary in I(0). However a ime series ha

More information

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j = 1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

A proof of Ito's formula using a di Title formula. Author(s) Fujita, Takahiko; Kawanishi, Yasuhi. Studia scientiarum mathematicarum H Citation

A proof of Ito's formula using a di Title formula. Author(s) Fujita, Takahiko; Kawanishi, Yasuhi. Studia scientiarum mathematicarum H Citation A proof of Io's formula using a di Tile formula Auhor(s) Fujia, Takahiko; Kawanishi, Yasuhi Sudia scieniarum mahemaicarum H Ciaion 15-134 Issue 8-3 Dae Type Journal Aricle Tex Version auhor URL hp://hdl.handle.ne/186/15878

More information

Chapter 4. Location-Scale-Based Parametric Distributions. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University

Chapter 4. Location-Scale-Based Parametric Distributions. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Chaper 4 Locaion-Scale-Based Parameric Disribuions William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based on he auhors

More information

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS

More information

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange

More information

Regression with Time Series Data

Regression with Time Series Data Regression wih Time Series Daa y = β 0 + β 1 x 1 +...+ β k x k + u Serial Correlaion and Heeroskedasiciy Time Series - Serial Correlaion and Heeroskedasiciy 1 Serially Correlaed Errors: Consequences Wih

More information

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013 Mahemaical Theory and Modeling ISSN -580 (Paper) ISSN 5-05 (Online) Vol, No., 0 www.iise.org The ffec of Inverse Transformaion on he Uni Mean and Consan Variance Assumpions of a Muliplicaive rror Model

More information

I. Return Calculations (20 pts, 4 points each)

I. Return Calculations (20 pts, 4 points each) Universiy of Washingon Spring 015 Deparmen of Economics Eric Zivo Econ 44 Miderm Exam Soluions This is a closed book and closed noe exam. However, you are allowed one page of noes (8.5 by 11 or A4 double-sided)

More information

Elements of Stochastic Processes Lecture II Hamid R. Rabiee

Elements of Stochastic Processes Lecture II Hamid R. Rabiee Sochasic Processes Elemens of Sochasic Processes Lecure II Hamid R. Rabiee Overview Reading Assignmen Chaper 9 of exbook Furher Resources MIT Open Course Ware S. Karlin and H. M. Taylor, A Firs Course

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338

More information

Stochastic Modelling of Electricity and Related Markets: Chapter 3

Stochastic Modelling of Electricity and Related Markets: Chapter 3 Sochasic Modelling of Elecriciy and Relaed Markes: Chaper 3 Fred Espen Benh, Jūraė Šalyė Benh and Seen Koekebakker Presener: Tony Ware Universiy of Calgary Ocober 14, 2009 The Schwarz model S() = S(0)

More information