Max stable Processes & Random Fields: Representations, Models, and Prediction
|
|
- Candice Lyons
- 5 years ago
- Views:
Transcription
1 Max stable Processes & Random Fields: Representations, Models, and Prediction Stilian Stoev University of Michigan, Ann Arbor March 2, 2011 Based on joint works with Yizao Wang and Murad S. Taqqu.
2 1 Preliminaries 2 Representations 3 Further Examples 4 Prediction for Max stable Random Fields Future work 6 References 7 Appendix
3 Preliminaries
4 Max stable processes. Definition and Motivation. A stochastic process (field) X = {X t } t T is max stable if for all n N, { fdd {X t } t T = a n max X (i) t + b n 1 i n }t T, for some a n > 0, b n R, where X (i), i = 1, 2, are independent copies of X.
5 Max stable processes. Definition and Motivation. A stochastic process (field) X = {X t } t T is max stable if for all n N, { fdd {X t } t T = a n max X (i) t + b n 1 i n }t T, for some a n > 0, b n R, where X (i), i = 1, 2, are independent copies of X. Fact: If for i.i.d. processes ξ (i) = {ξ (i) t } t T, i = 1, 2,, we have n a n i=1 ξ t (i) + b n fdd X t, as n, for some a n > 0, b n R, then the limit X = {X t } t T is max stable.
6 Extreme Value Distributions The following is the counterpart to the CLT for maxima: Theorem (Fisher Tippett (1928) & Gnedenko (1943)) Suppose that X n, n = 1, 2,... are iid random variables and a n max X i + b n a n 1 i n 1 i n X i + b n d Z, as n, for some a n > 0, b n R. Then, G(x) = P{Z x} is one of three types of extreme value distributions (EVD): e x α, x > 0(α > 0) (Fréchet) G(x) = e e x, x R (Gumbel) e ( x)α, x < 0(α > 0) (negative Fréchet) The type of the limit is determined by the tails of the X i s. Max stable processes have extreme value marginals.
7 α Fréchet max stable processes Z is said to be α Fréchet, if P(Z x) = exp{ σ α /x α }, x > 0, with some scale parameter σ > 0.
8 α Fréchet max stable processes Z is said to be α Fréchet, if P(Z x) = exp{ σ α /x α }, x > 0, with some scale parameter σ > 0. Heavy tails: P(Z > x) = 1 e σα /x α σ α /x α, as x.
9 α Fréchet max stable processes Z is said to be α Fréchet, if P(Z x) = exp{ σ α /x α }, x > 0, with some scale parameter σ > 0. Heavy tails: P(Z > x) = 1 e σα /x α σ α /x α, as x. The process X = {X t } t T is called α Fréchet, if the max linear combination n a i X ti is α Fréchet i=1 for all a i > 0, t i T, 1 i n. Note: Recall the definition of Gaussian processes.
10 α Fréchet max stable processes Z is said to be α Fréchet, if P(Z x) = exp{ σ α /x α }, x > 0, with some scale parameter σ > 0. Heavy tails: P(Z > x) = 1 e σα /x α σ α /x α, as x. The process X = {X t } t T is called α Fréchet, if the max linear combination n a i X ti is α Fréchet i=1 for all a i > 0, t i T, 1 i n. Note: Recall the definition of Gaussian processes. Fact: A process X with α Fréchet marginals is max stable if and only if it is an α Fréchet process. No generality is lost if one focuses on the class of α Fréchet process max stable processes.
11 Asymptotics in a Picture: The Chocolate Hills Where?
12 Asymptotics in a Picture: The Chocolate Hills Where? Bohol Island, Philippines.
13 Representations
14 The de Haan s spectral represenation Theorem (de Haan (1984)) Let X = {X t } t R be a stochastically continuous α Fréchet max stable process. Then, there exist functions f t (u) L α +([0, 1], du), such that f.d.d. {X t } t R = { f t (U i ) } i 1 Γ 1/α i where 0 Γ 1 Γ 2 is a standard Poisson point process on (0, ) and the U i s are independent Uniform(0, 1) random variables, independent of the Γ i s. t R,
15 Extremal integral representations Theorem (S. & Taqqu (2006); S. (2008)) Let X = {X t } t T be an α Fréchet max stable process (separable in probability). Then, {X t } t X fdd = { e S } f t (u)m α (du), t T for some deterministic f t (u) L α +(S, µ). Here M α is an α Fréchet random sup measure.
16 Extremal integral representations Theorem (S. & Taqqu (2006); S. (2008)) Let X = {X t } t T be an α Fréchet max stable process (separable in probability). Then, {X t } t X fdd = { e S } f t (u)m α (du), t T for some deterministic f t (u) L α +(S, µ). Here M α is an α Fréchet random sup measure. The sup measure: Given a measure space (S, S, µ), M α is said to be an α Fréchet sup measure with control measure µ, if: independently scattered: M α (A i ) s are independent, for disjoint A i s. α Fréchet: P(M α (A) x) = exp{ µ(a)x α }, x > 0. σ sup-additive: For disjoint measurable A i s: M α ( i=1a i ) = sup M α (A i ), i N almost surely.
17 Extremal Integrals For a simple and non negative function f (u) := n i=1 a i1 Ai (u), (A i S), with disjoint A i s, we set e S f (u)m α (du) def = n a i M α (A i ), i=1
18 Extremal Integrals For a simple and non negative function f (u) := n i=1 a i1 Ai (u), (A i S), with disjoint A i s, we set Properties: S e S f (u)m α (du) def = i=1 n a i M α (A i ), i=1 α Fréchet: For all x > 0, we have ( e ) { n P fdm α x = exp ai α µ(a i )x α} { = exp f α L α (µ) /x α}.
19 Extremal Integrals For a simple and non negative function f (u) := n i=1 a i1 Ai (u), (A i S), with disjoint A i s, we set Properties: S e S f (u)m α (du) def = i=1 n a i M α (A i ), i=1 α Fréchet: For all x > 0, we have ( e ) { n P fdm α x = exp ai α µ(a i )x α} { = exp f α L α (µ) /x α}. max linearity: For all a, b > 0 and f, g L α +(µ), e S (af bg)dm α = a e S fdm α b e S gdm α.
20 Extremal Integrals For a simple and non negative function f (u) := n i=1 a i1 Ai (u), (A i S), with disjoint A i s, we set Properties: S e S f (u)m α (du) def = i=1 n a i M α (A i ), i=1 α Fréchet: For all x > 0, we have ( e ) { n P fdm α x = exp ai α µ(a i )x α} { = exp f α L α (µ) /x α}. max linearity: For all a, b > 0 and f, g L α +(µ), e S (af bg)dm α = a e S fdm α b e S gdm α. e (convergence) S f P ndm α ξ if and only if S f n α f α 0, for some f L α e +(µ).thus, extends by continuity in P to all f L α e +(µ) and ξ := S fdm α.
21 Back to max stable processes For any collection of non negative f t L α +(µ), the process X t := is α Fréchet (max stable). e S f t dm α, (t T ) { ( f ti (u) ) αµ(du) } P(X ti x i, 1 i n) = exp max, (x i > 0). S 1 i n x i Follows: from max linearity of the extremal integral.
22 Back to max stable processes For any collection of non negative f t L α +(µ), the process X t := is α Fréchet (max stable). e S f t dm α, (t T ) { ( f ti (u) ) αµ(du) } P(X ti x i, 1 i n) = exp max, (x i > 0). S 1 i n x i Examples: stationary α Fréchet processes/fields. (moving maxima) With (S, S, µ) (R, B, Leb), X t := e f (t x)m α (dx), R t R (mixed moving maxima) (S, S, µ) (R d V, B R d V, Leb ν), X t := e R V f ( t x, v)m α (d x, dv), t R d. Known as M3 process in Zhang & Smith (2004)
23 A Basic Example: The Extremal Process Take (S, S, µ) (R, B, Leb) and let X t = e dm α (dx) M α ([0, t)), t > 0. [0,t)
24 A Basic Example: The Extremal Process Take (S, S, µ) (R, B, Leb) and let X t = e dm α (dx) M α ([0, t)), t > 0. [0,t) For all t > 0, we have X t = M α ([0, t)) d = t 1/α Z, where Z is standard α Fréchet: P(Z x) = e x α, x > 0.
25 A Basic Example: The Extremal Process Take (S, S, µ) (R, B, Leb) and let X t = e dm α (dx) M α ([0, t)), t > 0. [0,t) For all t > 0, we have X t = M α ([0, t)) d = t 1/α Z, where Z is standard α Fréchet: P(Z x) = e x α, x > 0. The process {X t } t 0 has independent max increments: (X t1,x t2,,x tn ) d = (t 1/α 1 Z 1,t 1/α 1 Z 1 (t 2 t 1 ) 1/α Z 2,,t 1/α 1 Z 1 (t n t n 1 ) 1/α Z n ), where the Z i s are iid standard α Fréchet. Can be viewed as an analog of the Brownian motion.
26 α Fréchet Extremal Processes: Sample Paths 16 Five paths of the α Frechet Extremal Process: α = X(t) t
27 Further Examples
28 Smith Random Fields Let Σ be a positive definite d d symmetric matrix and set φ Σ ( x) := 1 { (2π) d/2 Σ 1/2 exp for the Normal density with covariance Σ. Define the random field X t := i 1 x T Σ 1 x }, x R d 2 φ Σ (t U i ) Γ 1/α, t R d, i where {(Γ i, U i )} is a Poisson point process on (0, ) R d with the Lebesgue intensity. Then {X t } t R d is a moving average α Fréchet random field: {X t } t R d f.d.d. = { e R d φ Σ (t x)m α (dx) } t R d.
29 Smith Random Fields: Sample Surface 2D Smith random field α = 3, ρ = Y X
30 Brown Resnick Processes Let w t (ω ) be a standard Brownian motion supported on (Ω, F, P ). Consider a random α Fréchet sup measure M α (dω ) supported on (Ω, F, P). Define the α Fréchet process: X t := e Ω e wt(ω ) α t /2 M α (dω ), t R.
31 Brown Resnick Processes Let w t (ω ) be a standard Brownian motion supported on (Ω, F, P ). Consider a random α Fréchet sup measure M α (dω ) supported on (Ω, F, P). Define the α Fréchet process: X t := e Ω e wt(ω ) α t /2 M α (dω ), t R. Surprisingly, the process {X t } t R is stationary. It is called a Brown Resnick stationary α Fréchet process.
32 Brown Resnick Processes Let w t (ω ) be a standard Brownian motion supported on (Ω, F, P ). Consider a random α Fréchet sup measure M α (dω ) supported on (Ω, F, P). Define the α Fréchet process: X t := e Ω e wt(ω ) α t /2 M α (dω ), t R. Surprisingly, the process {X t } t R is stationary. It is called a Brown Resnick stationary α Fréchet process. More generally, if {w t } t R d is a zero-mean Gaussian random field with stationary increments, then X t := e wt(ω ) ασ2 e t /2 M α (dω ), t R d, Ω is a stationary α Fréchet random field. For more details, see Kabluchko, Schlather & de Haan (2009).
33 Prediction
34 Prediction for Max stable Random Fields Given is a max stable random field: X t := Problem: We observe X t1,, X tn. e S f t (u)m α (du), t T. (1)
35 Prediction for Max stable Random Fields Given is a max stable random field: X t := Problem: We observe X t1,, X tn. Goal: Predict X s1,, X sm. e S f t (u)m α (du), t T. (1)
36 Prediction for Max stable Random Fields Given is a max stable random field: X t := e S f t (u)m α (du), t T. (1) Problem: We observe X t1,, X tn. Goal: Predict X s1,, X sm. Need: To compute functionals of the conditional distribution (X s1,, X sm ) (X t1,, X tn ).
37 Prediction for Max stable Random Fields Given is a max stable random field: X t := Problem: We observe X t1,, X tn. Goal: Predict X s1,, X sm. e S f t (u)m α (du), t T. (1) Need: To compute functionals of the conditional distribution Challenges: (X s1,, X sm ) (X t1,, X tn ). The c.d.f. s are theoretically available from (1), but the densities are not.
38 Prediction for Max stable Random Fields Given is a max stable random field: X t := Problem: We observe X t1,, X tn. Goal: Predict X s1,, X sm. e S f t (u)m α (du), t T. (1) Need: To compute functionals of the conditional distribution Challenges: (X s1,, X sm ) (X t1,, X tn ). The c.d.f. s are theoretically available from (1), but the densities are not. The current state-of-the-art: only bi variate densities are known in closed form, for just a couple of models.
39 The Max Linear Model Consider the model X = A Z,
40 The Max Linear Model Consider the model X = A Z, that is p X i = a ij Z j, (1 i n) j=1 where A = (a ij ) n p, X = (Xi ) n i=1, Z = (Zj ) p j=1.
41 The Max Linear Model Consider the model X = A Z, that is p X i = a ij Z j, (1 i n) j=1 where A = (a ij ) n p, X = (Xi ) n i=1, Z = (Zj ) p j=1. If the Z j s are i.i.d. α Fréchet, then X is max stable.
42 The Max Linear Model Consider the model X = A Z, that is p X i = a ij Z j, (1 i n) j=1 where A = (a ij ) n p, X = (Xi ) n i=1, Z = (Zj ) p j=1. If the Z j s are i.i.d. α Fréchet, then X is max stable. In fact, X i = e [0,1] p f i dm α, with f i (u) = p 1/α a ij 1 [(j 1)/p,j/p) (u) j=1
43 The Max Linear Model Consider the model X = A Z, that is p X i = a ij Z j, (1 i n) j=1 where A = (a ij ) n p, X = (Xi ) n i=1, Z = (Zj ) p j=1. If the Z j s are i.i.d. α Fréchet, then X is max stable. In fact, X i = e [0,1] p f i dm α, with f i (u) = p 1/α a ij 1 [(j 1)/p,j/p) (u) j=1 Max linear models with large p can approximate, arbitrary e max stable models: X t = [0,1] f tdm α.
44 Another Representation: Spectral Measure For all max stable (α Fréchet) random vectors X : P( X { ( x) = exp S n 1 + max 1 i n w ) i αγ(d } w), ( x 0) x i where S n 1 + = { w 0 : w = max 1 i n w i = 1}, and Γ is unique finite measure on S n 1 + the spectral measure of X.
45 Another Representation: Spectral Measure For all max stable (α Fréchet) random vectors X : P( X { ( x) = exp S n 1 + max 1 i n w ) i αγ(d } w), ( x 0) x i where S n 1 + = { w 0 : w = max 1 i n w i = 1}, and Γ is unique finite measure on S n 1 + the spectral measure of X. For the max-linear model, Γ is discrete: Γ(d w) = where a j = (a ij ) n i=1 Rn +. p a j α δ aj / a j (d w), j=1
46 Another Representation: Spectral Measure For all max stable (α Fréchet) random vectors X : P( X { ( x) = exp S n 1 + max 1 i n w ) i αγ(d } w), ( x 0) x i where S n 1 + = { w 0 : w = max 1 i n w i = 1}, and Γ is unique finite measure on S n 1 + the spectral measure of X. For the max-linear model, Γ is discrete: Γ(d w) = where a j = (a ij ) n i=1 Rn +. p a j α δ aj / a j (d w), j=1 Hence the max linear models are called spectrally discrete.
47 Examples Spectrally discrete random field: Let φ t (j) 0 be deterministic functions and p X t := φ j (t)z j t T. j=1 For any t 1,, t n, and s 1,, s m, we have X = A Z, with a ij := φ ti (j) and X = (X ti ) n i=1. and Y = B Z, with b ij := φ si (j) and Y = (X si ) m i=1. Moving maxima: Let φ(u) L α +(R d, du), then X t := e R d φ(t u)m α (du) P φ(t u j )Z j. j [ M,M] d
48 Prediction: a Computational Solution Problem: Consider the max linear model X = A Z, where A = (a ij ) n p, a ij > 0, and the Z j s are independent α Fréchet. Goal: sample from the conditional distribution of Z X. Plug in: Y = B Z and obtain samples from the conditional distribution Y X : (X s1,, X sm ) (X t1,, X tn ). }{{}}{{} = Y = X
49 Prediction: a Computational Solution Problem: Consider the max linear model X = A Z, where A = (a ij ) n p, a ij > 0, and the Z j s are independent α Fréchet. Goal: sample from the conditional distribution of Z X. Plug in: Y = B Z and obtain samples from the conditional distribution Y X : (X s1,, X sm ) (X t1,, X tn ). }{{}}{{} = Y Idea: If Y 1,, Y N are independent samples from Y X = x, then by the LLN, for P X almost all x: 1 N h( Y N ( i ) a.s. E h( Y ) σ( X ) ), as N. j=1 = X
50 Benefits & Caveats Benefits:
51 Benefits & Caveats Benefits: Approximations to the optimal predictor for h( Z) given X in the mean square sense.
52 Benefits & Caveats Benefits: Approximations to the optimal predictor for h( Z) given X in the mean square sense. Monte Carlo approximations to conditional medians, quantiles other optimal predictors and prediction intervals.
53 Benefits & Caveats Benefits: Approximations to the optimal predictor for h( Z) given X in the mean square sense. Monte Carlo approximations to conditional medians, quantiles other optimal predictors and prediction intervals.
54 Benefits & Caveats Benefits: Approximations to the optimal predictor for h( Z) given X in the mean square sense. Monte Carlo approximations to conditional medians, quantiles other optimal predictors and prediction intervals.
55 Benefits & Caveats Benefits: Approximations to the optimal predictor for h( Z) given X in the mean square sense. Monte Carlo approximations to conditional medians, quantiles other optimal predictors and prediction intervals.
56 Benefits & Caveats Benefits: Approximations to the optimal predictor for h( Z) given X in the mean square sense. Monte Carlo approximations to conditional medians, quantiles other optimal predictors and prediction intervals.
57 Benefits & Caveats Benefits: Caveats: Approximations to the optimal predictor for h( Z) given X in the mean square sense. Monte Carlo approximations to conditional medians, quantiles other optimal predictors and prediction intervals. Need to sample efficiently from the conditional distribution. Must be able to handle large p = dim(z)!
58 Basic Observations Suppose x = A Z. (2) Observations: For all 1 j n: x i 0 Z j ẑ j ( x) min ( where 1 1 i n a 0 = ) ij For some j s the upper bounds are attained, so that (2) holds! Hitting scenario: such a set J = J( x) {1,, p}, that yields (2). Z j = ẑ j ( x), j J and Z j < ẑ j ( x), j J,
59 The Simple Solution Theorem (Wang & S. (2010)) The regular conditional probability p Z X ( x) equals: p Z X (E x) = p J ( x)ν J (E x), J hitting scenario where p J ( x) = w J / K w K for some weights w j ν J (E x) = δẑj (π j (E)) P(Z j π j (E) Z j ẑ j ). j J }{{} active constraints j J } {{ } inequality constraints
60 The Simple Solution Theorem (Wang & S. (2010)) The regular conditional probability p Z X ( x) equals: p Z X (E x) = p J ( x)ν J (E x), J hitting scenario where p J ( x) = w J / K w K for some weights w j ν J (E x) = δẑj (π j (E)) P(Z j π j (E) Z j ẑ j ). j J j J The weights: w J = 0, if J > r( x) and if J = r( x): w J = j J ẑjf Zj (ẑ j ) j J F Z j (ẑ j ). Notes: Here r( x) is the minimal number of equality constraints: r( x) = min J hitting scenario for x J. The result applies to independent Z j s with p.d.f. s f Z s and c.d.f. s F Z s.
61 Example Consider the toy model X 1 X 2 = X Z 1 Z 2 Z 3.
62 Example Consider the toy model X 1 X 2 = X Z 1 Z 2 Z 3. We have X 1 = Z 1, X 2 = Z 1 Z 2, X 3 = Z 1 Z 2 Z 3.
63 Example Consider the toy model X 1 X 2 = X Z 1 Z 2 Z 3. We have X 1 = Z 1, X 2 = Z 1 Z 2, X 3 = Z 1 Z 2 Z 3. If x = (2 2 2) T, then the hitting scenarios are: J 1 = {1}, J 2 = {1, 2}, J 3 = {1, 3}, J 4 = {1, 2, 3}.
64 Example Consider the toy model X 1 X 2 = X Z 1 Z 2 Z 3. We have X 1 = Z 1, X 2 = Z 1 Z 2, X 3 = Z 1 Z 2 Z 3. If x = (2 2 2) T, then the hitting scenarios are: J 1 = {1}, J 2 = {1, 2}, J 3 = {1, 3}, J 4 = {1, 2, 3}. Thus, r = r( x) = 1, and the conditional probability p Z X (dz x) = δ 2 (dz 1 )P(Z 2 dz 2 Z 2 < 2)P(Z 3 dz 3 Z 3 < 2) has only one component.
65 The Efficient Solution. Implementation Big Problem: Finding all hitting scenarios with brute force is NP hard (as a function of n and p)! Harder than the set covering problem.
66 The Efficient Solution. Implementation Big Problem: Finding all hitting scenarios with brute force is NP hard (as a function of n and p)! Harder than the set covering problem. Big Solution: Yizao found an equivalent representation of the regular conditional probability. Theorem (Wang & S. (2010)) For P X almost all x, we have p Z X (E x) = r( x) k=1 ν(k) (E x), with ν (k) (E x) = j J (k) ( x) p (k) j ( x)ν (k) j (E x), p (k) j ( x) = where w (k) j ( x) = ẑ j f Zj (ẑ j ) l J (k) \{j} F Z l (ẑ l ), ν (k) j (E x) = δẑj (π j (E)) l J (k) \{j} ( x) i J (k) ( x) w (k) i ( x), w (k) j P(Z l π l (E) Z l ẑ l ).
67 Some ingredients in the proof Hitting matrix: Given x = A Z, construct { 0, if aij ẑ H = (h ij ) n p, where h ij = j < x i 1, if a ij ẑ j = x i Intuition: want to pick r = r( x) columns of H (i.e., j s) that cover all rows.
68 Some ingredients in the proof Hitting matrix: Given x = A Z, construct { 0, if aij ẑ H = (h ij ) n p, where h ij = j < x i 1, if a ij ẑ j = x i Intuition: want to pick r = r( x) columns of H (i.e., j s) that cover all rows. Key facts: With P X probability one, the structure of H = H( x) is nice! One can decompose {1,, p} into r = r( x) disjoint classes: J(k), k = 1,, r. One and only one j in each J (k) J (k), k = 1,, r is active. For more details, please see Wang & S. (2010).
69 An application: Max AR(q) Davis & Resnick (1989): Max AR(q) are the stationary solutions to: X t = φ 1 X t 1 φ m X t m Z t, t Z, (3) with φ i 0, 1 i q and i.i.d. 1 Fréchet Z t s.
70 An application: Max AR(q) Davis & Resnick (1989): Max AR(q) are the stationary solutions to: X t = φ 1 X t 1 φ m X t m Z t, t Z, (3) with φ i 0, 1 i q and i.i.d. 1 Fréchet Z t s. If φ = q i=1 φ i < 1, then (3) has a unique solution. X t = ψ j Z t j. j=0
71 An application: Max AR(q) Davis & Resnick (1989): Max AR(q) are the stationary solutions to: X t = φ 1 X t 1 φ m X t m Z t, t Z, (3) with φ i 0, 1 i q and i.i.d. 1 Fréchet Z t s. If φ = q i=1 φ i < 1, then (3) has a unique solution. X t = ψ j Z t j. j=0 ψ j s decay exponentially: by truncation, we get a max linear model X t P Xt = M ψ j Z t j. j=0
72 An application: Max AR(q) Davis & Resnick (1989): Max AR(q) are the stationary solutions to: X t = φ 1 X t 1 φ m X t m Z t, t Z, (3) with φ i 0, 1 i q and i.i.d. 1 Fréchet Z t s. If φ = q i=1 φ i < 1, then (3) has a unique solution. X t = ψ j Z t j. j=0 ψ j s decay exponentially: by truncation, we get a max linear model X t P Xt = M ψ j Z t j. Max AR(q): we get the projection predictors by recursively solving j=0 X t+h = φ 1 Xt+h 1 φ q Xt+h q.
73 Prediction for Max AR Using i.i.d. conditional samples, we get: the conditional median, 9% predition intervals. estimates of P{X t+h X t+h X s, s t}. X Prediction of MARMA(3,0) processes MARMA process Conditional 9% quantile Conditional median Projection predictor Lag (h) Coverage Width Proj Pred t
74 12 Moving Maxima Random Fields Smith moving maxima random field model where φ(t) = e X t = φ(t u)mα(du), R2 β 1 β 2 2π 1 ρ 2 exp { [β2 1 t2 1 2ρβ 1β 2 t 1 t 2 + β 2 2 t2 2 ] 2(1 ρ 2 ) Figure on the left: 4 conditional samples with β 1 = β 2 = 1, ρ = 0, given 7 observed values (all equal to ). }
75 The Conditional Median Conditional expected median of the Smith model Parameters:ρ=0, β 1=1, β 2=1
76 An application: Rainfall in Bourgogne, France The model: Discretization of the 2 component moving maxima (Smith models) where X t = Intuition: e φ β1 (t u)m α (1) (du) R 2 e φ β2 (t u)m α (2) (du), R 2 φ(u) φ β (u 1, u 2 ) = 1 2πβ 2 e (u2 1 +u2 2 )/2β2, β > 0. Need 2 components to represent large and small scale effects. For simplicity, use isotropic components (ρ = 0 and equal scales).
77 A Realization of the Two component Smith Model A Realization of a two component Smith model
78 Observations: X ti, for i = 1,, 146 stations maxima over 1 years of daily rainfall Lambert Longitude Lambert Latitude Conditional Median
79 Rainfall in Bourgogne, France Observations: X ti, for i = 1,, 146 stations maxima over 1 years of daily rainfall Lambert Longitude Conditional Median Model fitting via Cross Validation Do a grid search over (β 1, β 2 ): 1. Condition on largest stations. 2. Compute p values for the rest Check if uniform. If not: 4. Change (β 1, β 2 ), goto 1.
80 Future work General estimation methodology?? Hard but recent exciting work via partial likelihood (sandwich) methods available. Construction and estimation of max stable models area of interest. Spatio temporal models many environmental applications. Computer networks extremal dependence in delay/loads/congestions. Prediction for the spectrally continuous case?? Important but hard!
81 Thank you!
82 Some References Brown, B. M. & Resnick, S. I. (1977), Extreme values of independent stochastic processes, J. Appl. Probability 14(4), de Haan, L. (1984), A spectral representation for max stable processes, Annals of Probability, 12(4): Kabluchko & Schlather (2009) Ergodic properties of max infinitely divisible processes, Stoch. Process. Appl. Smith, R.L. (1990) Max-Stable Processes and Spatial Extremes. Unpubished manuscript. S. (2008) On the ergodicity and mixing of max-stable processes,stoch. Process. Appl. Wang (2010) maxlinear an R package. yizwang/software/maxlinear. Wang & S. (2010) On the structure and representations of max stable processes, Adv. Appl. Probab. 42(3): Wang & S. (2010) Conditional sampling for spectrally discrete max stable random fields. Adv. Appl. Probab. To appear. Zhang & Smith (2004) The behavior of multivariate maxima of moving maxima processes, Journal of Applied Probability, 41(4):
83 Appendix: Some Derivations
84 The idea behind de Haan s spectral representation { P{X tk x k, k = 1,, n} = P { =: P i 1 g(u i ) Γ 1/α i i 1 1 Γ 1/α i } 1 = P{N A = } = e µ(a), } max f t k (U i )/x k 1 1 k n where µ is the intensity measure of the Poisson point process N = {(U i, Γ i )} and A = {(u, x) : g(u) x 1/α } c {(u, x) : g(u) α > x}. We have (by Fubini) that µ(a) = 1 Thus, { P{X tk x k, k = 1,, n} = exp {g α (u)>x}dudx = 1 0 g α (u)du. ( ) 1 αdu } max 0 1 k n f tk (u)/x k.
85 Stationarity of the standard Brown Resnick process For simplicity α 1 and let 0 t 1 t n, τ > 0: P{X tk +τ x k, 1 k n} { } = exp max x 1 Ω 1 k n k ewt k +τ (ω ) (t k +τ)/2 P (dω ) { ( = exp E P max x 1 1 k n k e(wt k +τ wτ tk/2) e (wτ τ/2)) } { ( ) = exp E P max x 1 1 k n k e(wt k t k/2) } = P{X tk x k, 1 k n}. Since {w t } t 0 has stationary and independent increments and since E P e wt = e t /2.
Max stable processes: representations, ergodic properties and some statistical applications
Max stable processes: representations, ergodic properties and some statistical applications Stilian Stoev University of Michigan, Ann Arbor Oberwolfach March 21, 2008 The focus Let X = {X t } t R be a
More informationConditional Sampling for Max Stable Random Fields
Conditional Sampling for Max Stable Random Fields Yizao Wang Department of Statistics, the University of Michigan April 30th, 0 th GSPC at Duke University, Durham, North Carolina Joint work with Stilian
More informationSimulation of Max Stable Processes
Simulation of Max Stable Processes Whitney Huang Department of Statistics Purdue University November 6, 2013 1 / 30 Outline 1 Max-Stable Processes 2 Unconditional Simulations 3 Conditional Simulations
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationOn the estimation of the heavy tail exponent in time series using the max spectrum. Stilian A. Stoev
On the estimation of the heavy tail exponent in time series using the max spectrum Stilian A. Stoev (sstoev@umich.edu) University of Michigan, Ann Arbor, U.S.A. JSM, Salt Lake City, 007 joint work with:
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationExtremogram and ex-periodogram for heavy-tailed time series
Extremogram and ex-periodogram for heavy-tailed time series 1 Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia) and Yuwei Zhao (Ulm) 1 Zagreb, June 6, 2014 1 2 Extremal
More informationCRPS M-ESTIMATION FOR MAX-STABLE MODELS WORKING MANUSCRIPT DO NOT DISTRIBUTE. By Robert A. Yuen and Stilian Stoev University of Michigan
arxiv: math.pr/ CRPS M-ESTIMATION FOR MAX-STABLE MODELS WORKING MANUSCRIPT DO NOT DISTRIBUTE By Robert A. Yuen and Stilian Stoev University of Michigan Max-stable random fields provide canonical models
More informationExtremogram and Ex-Periodogram for heavy-tailed time series
Extremogram and Ex-Periodogram for heavy-tailed time series 1 Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia) and Yuwei Zhao (Ulm) 1 Jussieu, April 9, 2014 1 2 Extremal
More informationAsymptotics and Simulation of Heavy-Tailed Processes
Asymptotics and Simulation of Heavy-Tailed Processes Department of Mathematics Stockholm, Sweden Workshop on Heavy-tailed Distributions and Extreme Value Theory ISI Kolkata January 14-17, 2013 Outline
More informationThe largest eigenvalues of the sample covariance matrix. in the heavy-tail case
The largest eigenvalues of the sample covariance matrix 1 in the heavy-tail case Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia NY), Johannes Heiny (Aarhus University)
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationEXACT SIMULATION OF BROWN-RESNICK RANDOM FIELDS
EXACT SIMULATION OF BROWN-RESNICK RANDOM FIELDS A.B. DIEKER AND T. MIKOSCH Abstract. We propose an exact simulation method for Brown-Resnick random fields, building on new representations for these stationary
More informationarxiv:math/ v1 [math.st] 16 May 2006
The Annals of Statistics 006 Vol 34 No 46 68 DOI: 04/009053605000000886 c Institute of Mathematical Statistics 006 arxiv:math/0605436v [mathst] 6 May 006 SPATIAL EXTREMES: MODELS FOR THE STATIONARY CASE
More informationRichard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC
EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu AMS Committee on Probability and Statistics
More informationStationary particle Systems
Stationary particle Systems Kaspar Stucki (joint work with Ilya Molchanov) University of Bern 5.9 2011 Introduction Poisson process Definition Let Λ be a Radon measure on R p. A Poisson process Π is a
More informationThe extremal elliptical model: Theoretical properties and statistical inference
1/25 The extremal elliptical model: Theoretical properties and statistical inference Thomas OPITZ Supervisors: Jean-Noel Bacro, Pierre Ribereau Institute of Mathematics and Modeling in Montpellier (I3M)
More informationThe Behavior of Multivariate Maxima of Moving Maxima Processes
The Behavior of Multivariate Maxima of Moving Maxima Processes Zhengjun Zhang Department of Mathematics Washington University Saint Louis, MO 6313-4899 USA Richard L. Smith Department of Statistics University
More informationBayesian Modelling of Extreme Rainfall Data
Bayesian Modelling of Extreme Rainfall Data Elizabeth Smith A thesis submitted for the degree of Doctor of Philosophy at the University of Newcastle upon Tyne September 2005 UNIVERSITY OF NEWCASTLE Bayesian
More informationRegular Variation and Extreme Events for Stochastic Processes
1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for
More informationIntroduction to Algorithmic Trading Strategies Lecture 10
Introduction to Algorithmic Trading Strategies Lecture 10 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationBeyond the color of the noise: what is memory in random phenomena?
Beyond the color of the noise: what is memory in random phenomena? Gennady Samorodnitsky Cornell University September 19, 2014 Randomness means lack of pattern or predictability in events according to
More informationSequential Monte Carlo Samplers for Applications in High Dimensions
Sequential Monte Carlo Samplers for Applications in High Dimensions Alexandros Beskos National University of Singapore KAUST, 26th February 2014 Joint work with: Dan Crisan, Ajay Jasra, Nik Kantas, Alex
More informationNegative Association, Ordering and Convergence of Resampling Methods
Negative Association, Ordering and Convergence of Resampling Methods Nicolas Chopin ENSAE, Paristech (Joint work with Mathieu Gerber and Nick Whiteley, University of Bristol) Resampling schemes: Informal
More informationSpatial extreme value theory and properties of max-stable processes Poitiers, November 8-10, 2012
Spatial extreme value theory and properties of max-stable processes Poitiers, November 8-10, 2012 November 8, 2012 15:00 Clement Dombry Habilitation thesis defense (in french) 17:00 Snack buet November
More informationAsymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½
University of Pennsylvania ScholarlyCommons Statistics Papers Wharton Faculty Research 1998 Asymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½ Lawrence D. Brown University
More informationMultivariate Generalized Ornstein-Uhlenbeck Processes
Multivariate Generalized Ornstein-Uhlenbeck Processes Anita Behme TU München Alexander Lindner TU Braunschweig 7th International Conference on Lévy Processes: Theory and Applications Wroclaw, July 15 19,
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2 Biostatistics,
More informationEstimation of spatial max-stable models using threshold exceedances
Estimation of spatial max-stable models using threshold exceedances arxiv:1205.1107v1 [stat.ap] 5 May 2012 Jean-Noel Bacro I3M, Université Montpellier II and Carlo Gaetan DAIS, Università Ca Foscari -
More informationMonte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan
Monte-Carlo MMD-MA, Université Paris-Dauphine Xiaolu Tan tan@ceremade.dauphine.fr Septembre 2015 Contents 1 Introduction 1 1.1 The principle.................................. 1 1.2 The error analysis
More informationPractical conditions on Markov chains for weak convergence of tail empirical processes
Practical conditions on Markov chains for weak convergence of tail empirical processes Olivier Wintenberger University of Copenhagen and Paris VI Joint work with Rafa l Kulik and Philippe Soulier Toronto,
More informationWeak quenched limiting distributions of a one-dimensional random walk in a random environment
Weak quenched limiting distributions of a one-dimensional random walk in a random environment Jonathon Peterson Cornell University Department of Mathematics Joint work with Gennady Samorodnitsky September
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry
More informationLARGE DEVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILED DEPENDENT RANDOM VECTORS*
LARGE EVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILE EPENENT RANOM VECTORS* Adam Jakubowski Alexander V. Nagaev Alexander Zaigraev Nicholas Copernicus University Faculty of Mathematics and Computer Science
More informationSTOCHASTIC GEOMETRY BIOIMAGING
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING 2018 www.csgb.dk RESEARCH REPORT Anders Rønn-Nielsen and Eva B. Vedel Jensen Central limit theorem for mean and variogram estimators in Lévy based
More informationMAXIMA OF LONG MEMORY STATIONARY SYMMETRIC α-stable PROCESSES, AND SELF-SIMILAR PROCESSES WITH STATIONARY MAX-INCREMENTS. 1.
MAXIMA OF LONG MEMORY STATIONARY SYMMETRIC -STABLE PROCESSES, AND SELF-SIMILAR PROCESSES WITH STATIONARY MAX-INCREMENTS TAKASHI OWADA AND GENNADY SAMORODNITSKY Abstract. We derive a functional limit theorem
More informationENSC327 Communications Systems 19: Random Processes. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 19: Random Processes Jie Liang School of Engineering Science Simon Fraser University 1 Outline Random processes Stationary random processes Autocorrelation of random processes
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationEvgeny Spodarev WIAS, Berlin. Limit theorems for excursion sets of stationary random fields
Evgeny Spodarev 23.01.2013 WIAS, Berlin Limit theorems for excursion sets of stationary random fields page 2 LT for excursion sets of stationary random fields Overview 23.01.2013 Overview Motivation Excursion
More informationChapter 5. Chapter 5 sections
1 / 43 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationHierarchical Modeling for Univariate Spatial Data
Hierarchical Modeling for Univariate Spatial Data Geography 890, Hierarchical Bayesian Models for Environmental Spatial Data Analysis February 15, 2011 1 Spatial Domain 2 Geography 890 Spatial Domain This
More informationStatistical inference on Lévy processes
Alberto Coca Cabrero University of Cambridge - CCA Supervisors: Dr. Richard Nickl and Professor L.C.G.Rogers Funded by Fundación Mutua Madrileña and EPSRC MASDOC/CCA student workshop 2013 26th March Outline
More informationSTA205 Probability: Week 8 R. Wolpert
INFINITE COIN-TOSS AND THE LAWS OF LARGE NUMBERS The traditional interpretation of the probability of an event E is its asymptotic frequency: the limit as n of the fraction of n repeated, similar, and
More informationNonparametric estimation of extreme risks from heavy-tailed distributions
Nonparametric estimation of extreme risks from heavy-tailed distributions Laurent GARDES joint work with Jonathan EL METHNI & Stéphane GIRARD December 2013 1 Introduction to risk measures 2 Introduction
More informationScaling properties of Poisson germ-grain models
Scaling properties of Poisson germ-grain models Ingemar Kaj Department of Mathematics Uppsala University ikaj@math.uu.se Reisensburg Workshop, Feb 2007 Purpose In spatial stochastic structure, investigate
More informationA Conditional Approach to Modeling Multivariate Extremes
A Approach to ing Multivariate Extremes By Heffernan & Tawn Department of Statistics Purdue University s April 30, 2014 Outline s s Multivariate Extremes s A central aim of multivariate extremes is trying
More informationGaussian Random Fields: Excursion Probabilities
Gaussian Random Fields: Excursion Probabilities Yimin Xiao Michigan State University Lecture 5 Excursion Probabilities 1 Some classical results on excursion probabilities A large deviation result Upper
More informationTail process and its role in limit theorems Bojan Basrak, University of Zagreb
Tail process and its role in limit theorems Bojan Basrak, University of Zagreb The Fields Institute Toronto, May 2016 based on the joint work (in progress) with Philippe Soulier, Azra Tafro, Hrvoje Planinić
More informationGaussian, Markov and stationary processes
Gaussian, Markov and stationary processes Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ November
More informationHierarchical Modelling for Univariate Spatial Data
Hierarchical Modelling for Univariate Spatial Data Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department
More information1 Simulating normal (Gaussian) rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions
Copyright c 2007 by Karl Sigman 1 Simulating normal Gaussian rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions Fundamental to many applications
More informationON EXTREME VALUE ANALYSIS OF A SPATIAL PROCESS
REVSTAT Statistical Journal Volume 6, Number 1, March 008, 71 81 ON EXTREME VALUE ANALYSIS OF A SPATIAL PROCESS Authors: Laurens de Haan Erasmus University Rotterdam and University Lisbon, The Netherlands
More informationNonlinear Time Series Modeling
Nonlinear Time Series Modeling Part II: Time Series Models in Finance Richard A. Davis Colorado State University (http://www.stat.colostate.edu/~rdavis/lectures) MaPhySto Workshop Copenhagen September
More informationStructural Reliability
Structural Reliability Thuong Van DANG May 28, 2018 1 / 41 2 / 41 Introduction to Structural Reliability Concept of Limit State and Reliability Review of Probability Theory First Order Second Moment Method
More informationEstimation of the long Memory parameter using an Infinite Source Poisson model applied to transmission rate measurements
of the long Memory parameter using an Infinite Source Poisson model applied to transmission rate measurements François Roueff Ecole Nat. Sup. des Télécommunications 46 rue Barrault, 75634 Paris cedex 13,
More informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More informationHandbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp
Handbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp Marcela Alfaro Córdoba August 25, 2016 NCSU Department of Statistics Continuous Parameter
More informationHardy-Stein identity and Square functions
Hardy-Stein identity and Square functions Daesung Kim (joint work with Rodrigo Bañuelos) Department of Mathematics Purdue University March 28, 217 Daesung Kim (Purdue) Hardy-Stein identity UIUC 217 1 /
More informationECE Lecture #10 Overview
ECE 450 - Lecture #0 Overview Introduction to Random Vectors CDF, PDF Mean Vector, Covariance Matrix Jointly Gaussian RV s: vector form of pdf Introduction to Random (or Stochastic) Processes Definitions
More informationExtreme Value Theory and Applications
Extreme Value Theory and Deauville - 04/10/2013 Extreme Value Theory and Introduction Asymptotic behavior of the Sum Extreme (from Latin exter, exterus, being on the outside) : Exceeding the ordinary,
More informationTail dependence in bivariate skew-normal and skew-t distributions
Tail dependence in bivariate skew-normal and skew-t distributions Paola Bortot Department of Statistical Sciences - University of Bologna paola.bortot@unibo.it Abstract: Quantifying dependence between
More informationCONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS
EVA IV, CONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS Jose Olmo Department of Economics City University, London (joint work with Jesús Gonzalo, Universidad Carlos III de Madrid) 4th Conference
More informationSupplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION. September 2017
Supplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION By Degui Li, Peter C. B. Phillips, and Jiti Gao September 017 COWLES FOUNDATION DISCUSSION PAPER NO.
More informationarxiv: v1 [math.pr] 3 Nov 2009
Ergodic Properties of Sum and Max Stable Stationary Random Fields via Null and Positive Group Actions Yizao Wang, Parthanil Roy 2 and Stilian A. Stoev University of Michigan, Michigan State University
More informationSingle Equation Linear GMM with Serially Correlated Moment Conditions
Single Equation Linear GMM with Serially Correlated Moment Conditions Eric Zivot October 28, 2009 Univariate Time Series Let {y t } be an ergodic-stationary time series with E[y t ]=μ and var(y t )
More informationLaplace Distribution Based Stochastic Models
Laplace Distribution Based Stochastic Models Anastassia Baxevani Department of Mathematics and Statistics University of Cyprus results based on work of/with Podgorski, Wegener, Aberg, Bolin, Walin July
More informationSpring 2012 Math 541B Exam 1
Spring 2012 Math 541B Exam 1 1. A sample of size n is drawn without replacement from an urn containing N balls, m of which are red and N m are black; the balls are otherwise indistinguishable. Let X denote
More informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
More informationSample Spaces, Random Variables
Sample Spaces, Random Variables Moulinath Banerjee University of Michigan August 3, 22 Probabilities In talking about probabilities, the fundamental object is Ω, the sample space. (elements) in Ω are denoted
More informationThe Convergence Rate for the Normal Approximation of Extreme Sums
The Convergence Rate for the Normal Approximation of Extreme Sums Yongcheng Qi University of Minnesota Duluth WCNA 2008, Orlando, July 2-9, 2008 This talk is based on a joint work with Professor Shihong
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationComputational statistics
Computational statistics Markov Chain Monte Carlo methods Thierry Denœux March 2017 Thierry Denœux Computational statistics March 2017 1 / 71 Contents of this chapter When a target density f can be evaluated
More informationLogarithmic scaling of planar random walk s local times
Logarithmic scaling of planar random walk s local times Péter Nándori * and Zeyu Shen ** * Department of Mathematics, University of Maryland ** Courant Institute, New York University October 9, 2015 Abstract
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationThreshold estimation in marginal modelling of spatially-dependent non-stationary extremes
Threshold estimation in marginal modelling of spatially-dependent non-stationary extremes Philip Jonathan Shell Technology Centre Thornton, Chester philip.jonathan@shell.com Paul Northrop University College
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationTopics in fractional Brownian motion
Topics in fractional Brownian motion Esko Valkeila Spring School, Jena 25.3. 2011 We plan to discuss the following items during these lectures: Fractional Brownian motion and its properties. Topics in
More informationTwo applications of max stable distributions: Random sketches and Heavy tail exponent estimation
Two applications of max stable distributions: Random sketches and Heavy tail exponent estimation Stilian Stoev (sstoev@umich.edu) Department of Statistics University of Michigan Boston University Probability
More informationLectures on Stochastic Stability. Sergey FOSS. Heriot-Watt University. Lecture 4. Coupling and Harris Processes
Lectures on Stochastic Stability Sergey FOSS Heriot-Watt University Lecture 4 Coupling and Harris Processes 1 A simple example Consider a Markov chain X n in a countable state space S with transition probabilities
More informationReview. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda
Review DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Probability and statistics Probability: Framework for dealing with
More informationTail negative dependence and its applications for aggregate loss modeling
Tail negative dependence and its applications for aggregate loss modeling Lei Hua Division of Statistics Oct 20, 2014, ISU L. Hua (NIU) 1/35 1 Motivation 2 Tail order Elliptical copula Extreme value copula
More informationDoes k-th Moment Exist?
Does k-th Moment Exist? Hitomi, K. 1 and Y. Nishiyama 2 1 Kyoto Institute of Technology, Japan 2 Institute of Economic Research, Kyoto University, Japan Email: hitomi@kit.ac.jp Keywords: Existence of moments,
More informationOn the Estimation and Application of Max-Stable Processes
On the Estimation and Application of Max-Stable Processes Zhengjun Zhang Department of Statistics University of Wisconsin Madison, WI 53706, USA Co-author: Richard Smith EVA 2009, Fort Collins, CO Z. Zhang
More informationNonlinear time series
Based on the book by Fan/Yao: Nonlinear Time Series Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 27, 2009 Outline Characteristics of
More informationEstimation de mesures de risques à partir des L p -quantiles
1/ 42 Estimation de mesures de risques à partir des L p -quantiles extrêmes Stéphane GIRARD (Inria Grenoble Rhône-Alpes) collaboration avec Abdelaati DAOUIA (Toulouse School of Economics), & Gilles STUPFLER
More informationOn the Fisher Bingham Distribution
On the Fisher Bingham Distribution BY A. Kume and S.G Walker Institute of Mathematics, Statistics and Actuarial Science, University of Kent Canterbury, CT2 7NF,UK A.Kume@kent.ac.uk and S.G.Walker@kent.ac.uk
More informationHeavy Tailed Time Series with Extremal Independence
Heavy Tailed Time Series with Extremal Independence Rafa l Kulik and Philippe Soulier Conference in honour of Prof. Herold Dehling Bochum January 16, 2015 Rafa l Kulik and Philippe Soulier Regular variation
More informationEXTREMAL QUANTILES OF MAXIMUMS FOR STATIONARY SEQUENCES WITH PSEUDO-STATIONARY TREND WITH APPLICATIONS IN ELECTRICITY CONSUMPTION ALEXANDR V.
MONTENEGRIN STATIONARY JOURNAL TREND WITH OF ECONOMICS, APPLICATIONS Vol. IN 9, ELECTRICITY No. 4 (December CONSUMPTION 2013), 53-63 53 EXTREMAL QUANTILES OF MAXIMUMS FOR STATIONARY SEQUENCES WITH PSEUDO-STATIONARY
More informationLarge sample distribution for fully functional periodicity tests
Large sample distribution for fully functional periodicity tests Siegfried Hörmann Institute for Statistics Graz University of Technology Based on joint work with Piotr Kokoszka (Colorado State) and Gilles
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationRandom Process Lecture 1. Fundamentals of Probability
Random Process Lecture 1. Fundamentals of Probability Husheng Li Min Kao Department of Electrical Engineering and Computer Science University of Tennessee, Knoxville Spring, 2016 1/43 Outline 2/43 1 Syllabus
More informationModule 9: Stationary Processes
Module 9: Stationary Processes Lecture 1 Stationary Processes 1 Introduction A stationary process is a stochastic process whose joint probability distribution does not change when shifted in time or space.
More informationSystem Identification, Lecture 4
System Identification, Lecture 4 Kristiaan Pelckmans (IT/UU, 2338) Course code: 1RT880, Report code: 61800 - Spring 2012 F, FRI Uppsala University, Information Technology 30 Januari 2012 SI-2012 K. Pelckmans
More information1: PROBABILITY REVIEW
1: PROBABILITY REVIEW Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 1: Probability Review 1 / 56 Outline We will review the following
More informationChapter 3 - Estimation by direct maximization of the likelihood
Chapter 3 - Estimation by direct maximization of the likelihood 02433 - Hidden Markov Models Martin Wæver Pedersen, Henrik Madsen Course week 3 MWP, compiled June 7, 2011 Recall: Recursive scheme for the
More informationGARCH processes probabilistic properties (Part 1)
GARCH processes probabilistic properties (Part 1) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/
More informationVine Copulas. Spatial Copula Workshop 2014 September 22, Institute for Geoinformatics University of Münster.
Spatial Workshop 2014 September 22, 2014 Institute for Geoinformatics University of Münster http://ifgi.uni-muenster.de/graeler 1 spatio-temporal data Typically, spatio-temporal data is given at a set
More informationSome Results on the Ergodicity of Adaptive MCMC Algorithms
Some Results on the Ergodicity of Adaptive MCMC Algorithms Omar Khalil Supervisor: Jeffrey Rosenthal September 2, 2011 1 Contents 1 Andrieu-Moulines 4 2 Roberts-Rosenthal 7 3 Atchadé and Fort 8 4 Relationship
More informationSystem Identification, Lecture 4
System Identification, Lecture 4 Kristiaan Pelckmans (IT/UU, 2338) Course code: 1RT880, Report code: 61800 - Spring 2016 F, FRI Uppsala University, Information Technology 13 April 2016 SI-2016 K. Pelckmans
More information1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1].
1 Introduction The content of these notes is also covered by chapter 3 section B of [1]. Diffusion equation and central limit theorem Consider a sequence {ξ i } i=1 i.i.d. ξ i = d ξ with ξ : Ω { Dx, 0,
More information