Bayesian model mergings for multivariate extremes Application to regional predetermination of oods with incomplete data
|
|
- Cecil Sims
- 6 years ago
- Views:
Transcription
1 Bayesian model mergings for multivariate extremes Application to regional predetermination of oods with incomplete data Anne Sabourin PhD defense University of Lyon 1 Committee: Anthony Davison, Johan Segers (Rapporteurs) Anne-Laure Fougères (Director), Philippe Naveau (Co-director), Clémentine Prieur, Stéphane Robin, Éric Sauquet. September 24 th,
2 Multivariate extreme values Risk management: Largest events, largest losses Hydrology: `ood predetermination'. Return levels (extreme quantiles) Return periods (1 / probability of occurrence ) digs, dams, land use plans. Simultaneous occurrence of rare events can be catastrophic Multivariate extremes: Probability of jointly extreme events? 2
3 Censored Multivariate extremes: oods in the `Gardons' joint work with Benjamin Renard Daily streamow data at 4 neighbouring sites : St Jean du Gard, Mialet, Anduze, Alès. Joint distributions of extremes? probability of simultaneous oods. Recent, `clean' series very short Historical data from archives, depending on `perception thresholds' for oods (Earliest: 1604). censored data Gard river Neppel et al. (2010) How to use all dierent kinds of data? 3
4 Multivariate extremes for regional analysis in hydrology Many sites, many parameters for marginal distributions, short observation period. `Regional analysis': replace time with space. Assume some parameters constant over the region and use extreme data from all sites. Independence between extremes at neighbouring sites? Dependence structure? Idea: use multivariate extreme value models 4
5 Outline Multivariate extremes and model uncertainty Bayesian model averaging (`Mélange de modèles') Dirichlet mixture model (`Modèle de mélange'): a re-parametrization Historical, censored data in the Dirichlet model 5
6 Multivariate extremes Random vectors Y = (Y 1,..., Y d, ) ; Y j 0 Margins: Y j F j, 1 j d (Generalized Pareto above large thresholds) Standardization ( unit Fréchet margins) X j = 1/ log [F j (Y j )] ; P(X j x) = e 1/x, 1 j d Joint behaviour of extremes: distribution of X above large thresholds? P(X A X A 0 )? (A A 0, 0 / A 0 ), A 0 `far from the origin'. X2 X A u 2 A 0 : Extremal region u 1 X1 6
7 Polar decomposition and angular measure Polar coordinates: R = d j=1 X j (L 1 norm) ; W = X R. W simplex S d = {w : w j 0, j w j = 1}. Angular probability measure: H(B) = P(W B) (B S d ). X 2nd component w B x S d W R B st component 7
8 Fundamental Result de Haan, Resnick, 70's, 80's Radial homogeneity (regular variation) P(R > r, W B R r 0 ) r 0 r 0 r H(B) (r = c r 0, c > 1) Above large radial thresholds, R is independent from W H (+ margins) entirely determines the joint distribution X1 X2 x w x X1 X2 x w x One condition only for genuine H: moments constraint w dh(w) = ( 1 d,..., 1 d ). Center of mass at the center of the simplex. Few constraints: non parametric family! 8
9 Estimating the angular measure: non parametric problem Non parametric estimation (empirical likelihood, Einmahl et al., 2001, Einmahl, Segers, 2009, Guillotte et al, 2011.) No explicit expression for asymptotic variance, Bayesian inference with d = 2 only. Restriction to parametric family: Gumbel, logistic, pairwise Beta... Coles & Tawn, 91, Cooley et al., 2010, Ballani & Schlather, 2011 : How to take into account model uncertainty? 9
10 Outline Multivariate extremes and model uncertainty Bayesian model averaging (`Mélange de modèles') Dirichlet mixture model (`Modèle de mélange'): a re-parametrization Historical, censored data in the Dirichlet model 10
11 BMA: Averaging estimates from dierent models Disjoint union of several parametric models Sabourin, Naveau, Fougères, 2013 (Extremes) Package R: `BMAmevt', available on CRAN repositories
12 Outline Multivariate extremes and model uncertainty Bayesian model averaging (`Mélange de modèles') Dirichlet mixture model (`Modèle de mélange'): a re-parametrization Historical, censored data in the Dirichlet model 12
13 Dirichlet distribution w S d, diri(w µ, ν) = Γ(ν) d i=1 Γ(νµ i) d i=1 w νµ i 1 i. µ S d : location parameter (point on the simplex): `center'; ν > 0 : concentration parameter. w w3 w1 13
14 Dirichlet distribution w S d, diri(w µ, ν) = Γ(ν) d i=1 Γ(νµ i) d i=1 w νµ i 1 i. µ S d : location parameter (point on the simplex): `center'; ν > 0 : concentration parameter. w w3 w1 13
15 Dirichlet mixture model Boldi, Davison, 2007 µ = µ,1:k, ν = ν 1:k, p = p 1:k, ψ = (µ, p, ν), k h ψ (w) = p m diri(w µ, m, ν m ) m=1 Moments constraint on (µ, p): k p m µ.,m = ( 1 d,..., 1 d ). m=1 Weakly dense family (k N) in the space of admissible angular measures 14
16 Bayesian inference and censored data Two issues : (i) parameters constraints (ii) censorship (i) Bayesian framework: MCMC methods to sample the posterior distribution. Constraints Sampling issues for d > 2. Re-parametrization: No more constraint, tting is manageable for d = 5: Sabourin, Naveau, 2013 (ii) Censoring: data points but segments or boxes in R d. Intervals overlapping threshold: extreme or not? Likelihood: density dr dh(w) integrated over boxes. r 2 Sabourin, under review ; Sabourin, Renard, in preparation 15
17 Bayesian inference and censored data Two issues : (i) parameters constraints (ii) censorship (i) Bayesian framework: MCMC methods to sample the posterior distribution. Constraints Sampling issues for d > 2. Re-parametrization: No more constraint, tting is manageable for d = 5: Sabourin, Naveau, 2013 (ii) Censoring: data points but segments or boxes in R d. Intervals overlapping threshold: extreme or not? Likelihood: density dr dh(w) integrated over boxes. r 2 Sabourin, under review ; Sabourin, Renard, in preparation 15
18 Re-parametrization: intermediate variables (γ 1,..., γ k 1), partial barycenters ex: k = 4 γ 0 γ m : Barycenter of kernels 'following µ.,m : µ.,m+1,..., µ.,k. γ m = ( ) 1 p j p j µ.,j j>m j>m 16
19 γ 1 on a line segment: eccentricity parameter e 1 (0, 1). ex: k = 4 I1 γ0 γ1 µ1 Draw (µ,1 S d, e 1 (0, 1)) γ 1 dened by γ 0 γ 1 γ 0 I 1 = e 1 ; p 1 = γ 0 γ 1 µ,1 γ 1. 17
20 γ 2 on a line segment: eccentricity parameter e 2 (0, 1). ex: k = 4 I1 µ 2 γ1 γ2 I2 γ0 µ 1 Draw (µ,2, e 2 ) γ 2 : γ 1 γ 2 γ 1 I 2 = e 2 p 2 18
21 Last density kernel = last center µ,k. ex: k = 4 I3 µ4 I1 µ 2 γ1 γ2 I2 γ0 µ 3 µ 1 Draw (µ,3, e 3 ) γ 3 p 3, µ.,4 = γ 3. p 4 19
22 Summary I3 µ4 I1 µ 2 γ1 γ2 I2 γ0 µ 3 µ 1 Given (µ.,1:k 1, e 1:k 1 ), One obtains (µ.,1:k, p 1:k ). The density h may thus be parametrized by θ = (µ.,1:k 1, e 1:k 1, ν 1:k ). 20
23 Bayesian model Unconstrained parameter space : union of product spaces (`rectangles') } Θ = Θ k ; Θ k = {(S d ) k 1 [0, 1) k 1 (0, ] k 1 k=1 Inference: Gibbs + Reversible-jumps. Restriction (numerical convenience) : k 15, ν < ν max, etc... `Reasonable' prior `at' and rotation invariant. Balanced weight and uniformly scattered centers. 21
24 MCMC sampling: Metropolis-within-Gibbs, reversible jumps. Three transition types for the Markov chain: Classical (Gibbs): one µ.,m, e m or a ν m is modied. Trans-dimensional (Green, 1995): One component (µ.,k, e k, ν k+1 ) is added or deleted. `Shue': Indices permutation of the original mixture: Re-allocating mass from old components to new ones. 22
25 Results: model's and algorithm's consistency Ergodicity: The generated MC is φ-irréducible, aperiodic and admits π n (= posterior W 1:n ) as invariant distribution. Consequence: g C b (Θ), 1 T T g(θ t ) E πn (g). Key point: π n is invariant under the `shue' moves. Posterior consistency for π n under `weak conditions' 2, π-a.s., U weakly open containing θ 0, t=1 π n (U) n 1. Consequence: E πn (g) g(θ 0). n Key: The Euclidian topology is ner than the Kullback topology in this model. 2 If the prior grants some mass to every Euclidian neighbourhood of Θ and if θ 0 is in the Kullback-Leibler closure of Θ 23
26 Convergence checking (simulated data, d = 5, k = 4) θ summarized by a scalar quantity (integrating the DM density against a test function) Original algorithm standard tests: Re-parametrized version Stationarity (Heidelberger & Welch, 83) variance ratio (inter/intra chains, Gelman, 92) 24
27 Estimating H in dimension 3, simulated data MCMC: T 2 = ; T 1 = w w3 w1 25
28 Dimension 5, simulated data Angular measure density for one pair (T 2 = , T 1 = ) X2/(X2+X5) 0 1 X2/(X2+X5) Gelman ratio: Original version: 2.18 ; Re-parametrized: Credibility sets (posterior quantiles): wider. 26
29 Bayesian inference and censored data Two issues : (i) parameters constraints (ii) censorship (i) Bayesian framework: MCMC methods to sample the posterior distribution. Constraints Sampling issues for d > 2. Re-parametrization: No more constraint, tting is manageable for d = 5: Sabourin, Naveau, 2013 (ii) Censoring: data points but segments or boxes in R d. Intervals overlapping threshold: extreme or not? Likelihood: density dr dh(w) integrated over boxes. r 2 Sabourin, under review ; Sabourin, Renard, in preparation 27
30 Bayesian inference and censored data Two issues : (i) parameters constraints (ii) censorship (i) Bayesian framework: MCMC methods to sample the posterior distribution. Constraints Sampling issues for d > 2. Re-parametrization: No more constraint, tting is manageable for d = 5: Sabourin, Naveau, 2013 (ii) Censoring: data points but segments or boxes in R d. Intervals overlapping threshold: extreme or not? Likelihood: density dr dh(w) integrated over boxes. r 2 Sabourin, under review ; Sabourin, Renard, in preparation 27
31 Outline Multivariate extremes and model uncertainty Bayesian model averaging (`Mélange de modèles') Dirichlet mixture model (`Modèle de mélange'): a re-parametrization Historical, censored data in the Dirichlet model 28
32 Censored data: univariate and pairwise plots Univariate time series: St Jean Mialet Débit (m3/s) u Année Anduze Débit (m3/s) u Année Ales Débit (m3/s) u Année Débit (m3/s) u Année 29
33 Censored data: univariate and pairwise plots Bivariate plots: Streamflow at Mialet (m3/s) Streamflow at Anduze (m3/s) Streamflow at Ales (m3/s) Streamflow at StJean (m3/s) Streamflow at StJean (m3/s) Streamflow at StJean (m3/s) 29
34 Using censored data: wishes and reality Take into account as many data as possible Censored likelihood, integration problems Information transfer from well gauged to poorly gauged sites using the dependence structure Estimate together marginal parameters + dependence 30
35 Data overlapping threshold and Poisson model How to include the rectangles overlapping threshold in the likelihood? {( t n, X t n ), 1 t n } Poisson Process (Leb λ) on [0, 1] A u,n λ: ` exponent measure', with Dirichlet Mixture angular component dλ dr dw (r, w) = d r 2 h(w). Overlapping events appear in Poisson likelihood as [ { P N ( t 2 n t 1 n ) 1 } ] n A i = 0 = exp [ (t 2 t 1 )λ(a i )] 31
36 `Censored' likelihood: model density integrated over boxes Ledford & Tawn, 1996: partially extreme data censored at threshold, GEV models Explicit expression for censored likelihood. Here: idem + natural censoring Poisson model No closed form expression for integrated likelihood. Two terms without closed form: Censored regions A i overlapping threshold: exp { (t 2 t 1 )λ(a i )} Classical censoring above threshold censored region dλ dx. 32
37 Data augmentation One more Gibbs step, no more numerical integration. Objective: sample [θ Obs] likelihood (censored obs) Additional variables (replace missing data component): Z Full conditionals [Z i Z j j, θ, Obs], [θ Z, Obs],... explicit (Thanks Dirichlet): Gibbs sampling. Sample [z, θ Obs] + (augmented distribution) on Θ Z. 33
38 Censored regions above threshold Censored region dλ dx dx j1:j r : Generate missing components under univariate conditional distributions Z j 1:r [X missing X obs, θ] x2 Augmentation data Z j = [X censored X observed, θ] Censored interval u 2 /n Extremal region u 1 /n x1 Dirichlet Explicit univariate conditionals Exact sampling of censored data on censored interval 34
39 Censored regions overlapping threshold augmentation Poisson process N i on E i A i. e (t 2,i t1,i )λ(a i ) + Functional ϕ(n i ) X2 Augmentation Z i = PP(τ.λ) on E i φ(#{points in Ai}) U' 1 E i A i u 2 /n Censored region u 1 /n U' 2 X1 [z, θ Obs] }{{}... [N i ]ϕ(n i ) density terms, prior, augmented missing components 35
40 Simulated data (Dirichlet, d = 4, k = 3 components), same censoring as real data Pairwise plot and angular measure density (true/ posterior predictive) S h S X3/( X3 + X4 ) 36
41 Simulated data (Dirichlet, d = 4, k = 3 components), same censoring as real data Marginal quantile curves: better in joint model. S3 discharge : dependent independent true return period (years, log scale) 36
42 Angular predictive density for Gardons data h St Jean/(St Jean+Mialet ) h St Jean/(St Jean+Anduze ) h St Jean/(St Jean+Ales ) h Mialet/(Mialet+Anduze ) h Mialet/(Mialet+Ales ) h Anduze/(Anduze+Ales ) 37
43 Conditional exceedance probability P( Mialet > y St Jean > 300 ) P( Anduze > y St Jean > 300 ) P( Ales > y St Jean > 300 ) Condit. exceed. proba empirical probability of an excess empirical error posterior predictive posterior quantiles Condit. exceed. proba Condit. exceed. proba threshold (y) threshold (y) threshold (y) P( Anduze > y Mialet > 320 ) P( Ales > y Mialet > 320 ) P( Ales > y Anduze > 520 ) Condit. exceed. proba Condit. exceed. proba Condit. exceed. proba threshold (y) threshold (y) threshold (y) 38
44 Conclusion Building Bayesian multivariate models for excesses: Dirichlet mixture family: `non' parametric, Bayesian inference possible up to re-parametrization Censoring data augmenting (Dirichlet conditioning properies) Two packages R: DiriXtremes, MCMC algorithm for Dirichlet mixtures, DiriCens, implementation with censored data. High dimensional sample space (GCM grid, spatial elds)? Impose reasonable structure (sparse) on Dirichlet parameters Dirichlet Process? Challenges : Discrete random measure continuous framework 39
45 Bibliographie I M.-O. Boldi and A. C. Davison. A mixture model for multivariate extremes. JRSS: Series B (Statistical Methodology), 69(2):217229, Coles, SG and Tawn, JA Modeling extreme multivariate events JR Statist. Soc. B, 53:377392, 1991 Gómez, G., Calle, M. L., and Oller, R. Frequentist and bayesian approaches for interval-censored data. Statistical Papers, 45(2):139173, Hosking, J.R.M. and Wallis, J.R.. Regional frequency analysis: an approach based on L-moments Cambridge University Press, Ledford, A. and Tawn, J. (1996). Statistics for near independence in multivariate extreme values. Biometrika, 83(1): Neppel, L., Renard, B., Lang, M., Ayral, P., Coeur, D., Gaume, E., Jacob, N., Payrastre, O., Pobanz, K., and Vinet, F. (2010). Flood frequency analysis using historical data: accounting for random and systematic errors. Hydrological Sciences JournalJournal des Sciences Hydrologiques, 55(2): Resnick, S. (1987). Extreme values, regular variation, and point processes, volume 4 of Applied Probability. A Series of the Applied Probability Trust. Springer-Verlag, New York. 40
46 Bibliographie II Sabourin, A., Naveau, P. and Fougères, A-L. (2013) Bayesian model averaging for multivariate extremes. Extremes, 16(3) Sabourin, A., Naveau, P. (2013) Bayesian Dirichlet mixture model for multivariate extremes: a re-parametrization. Computation. Stat and Data Analysis Schnedler, W. (2005). Likelihood estimation for censored random vectors. Econometric Reviews, 24(2): Tanner, M. and Wong, W. (1987). The calculation of posterior distributions by data augmentation. Journal of the American Statistical Association, 82(398): Van Dyk, D. and Meng, X. (2001). The art of data augmentation. Journal of Computational and Graphical Statistics, 10(1):
47 Outline Multivariate extremes and model uncertainty Bayesian model averaging (`Mélange de modèles') Dirichlet mixture model (`Modèle de mélange'): a re-parametrization Historical, censored data in the Dirichlet model 42
48 Bayesian Model Averaging: reducing model uncertainty. Parametric framework: arbitrary restriction, dierent models can yield dierent estimates. First option: Fight! Choose one model (Information criterions: BIC/ AIC /AICC) 43
49 Bayesian Model Averaging: reducing model uncertainty. Parametric framework: arbitrary restriction, dierent models can yield dierent estimates. First option: Fight! Choose one model (Information criterions: BIC/ AIC /AICC) H's family M 2 M 1 43
50 Bayesian Model Averaging: reducing model uncertainty. Parametric framework: arbitrary restriction, dierent models can yield dierent estimates. BMA = averaging predictions based on posterior model weights H's family M 2 Models' average M 1 Already widely studied and used in several contexts ( weather forecast... ). Hoeting et al. (99), Madigan & Raftery (94), Raftery et al. (05) 43
51 BMA: principle J statistical models M (1),... M (J), with parametrization Θ j, 1 j J and priors π j dened on Θ j 44
52 BMA: principle J statistical models M (1),... M (J), with parametrization Θ j, 1 j J and priors π j dened on Θ j BMA model = disjoint union: Θ = J 1 Θ j, with prior p on index set {1,..., J}: p(m j ) = `prior marginal model weight' for M j prior on Θ: π( J 1 B j) = J 1 p(m j) π j (B j ) (B j Θ j ) 44
53 BMA: principle J statistical models M (1),... M (J), with parametrization Θ j, 1 j J and priors π j dened on Θ j BMA model = disjoint union: Θ = J 1 Θ j, with prior p on index set {1,..., J}: p(m j ) = `prior marginal model weight' for M j prior on Θ: π( J 1 B j) = J 1 p(m j) π j (B j ) (B j Θ j ) posterior (conditioning on data X ) = weighted average J π( B j X ) = 1 posterior inm j J {}}{ p(m j X ) π j (B j X ) }{{} 1 posterior marginal model weight Key: posterior weights. (Laplace approx or standard MC?) p(m j X ) p(m j ) Likelihood(X θ) dπ j (θ) Θ j 44
54 BMA for multivariate extremes Sabourin, Naveau, Fougères (2013) Background: univariate EVD's of dierent types Stephenson & Tawn (04) or multivariate, asymptotically dependent/ independent EVD's Apputhurai & Stephenson (10) Our approach: averaging angular measure models, with angular data W. F 1,..., F J max-stable distributions j p jf j not max-stable. H 1,... H J angular measures (moments constraint) j p j H j is a valid angular measure! (linearity) 45
55 BMA for multivariate extremes Sabourin, Naveau, Fougères (2013) Background: univariate EVD's of dierent types Stephenson & Tawn (04) or multivariate, asymptotically dependent/ independent EVD's Apputhurai & Stephenson (10) Our approach: averaging angular measure models, with angular data W. F 1,..., F J max-stable distributions j p jf j not max-stable. H 1,... H J angular measures (moments constraint) j p j H j is a valid angular measure! (linearity) 45
56 Implementing and scoring BMA for Multivariate extremes Does the BMA framework perform signicantly better than selecting models based on AIC? Yes, in terms of logarithmic score for the predictive density (Kullback-Leibler divergence to the truth) Madigan & Raftery (94) In average over the union model, w.r.t prior! Simulation study : Evaluation via proper scoring rules (Logarithmic + probability of failure regions) 2 models of same dimension: Pairwise-Beta / Nested asymmetric logistic 100 data sets simulated from another model Results: The BMA framework performs consistently (for all scores), slightly (1/20 to 1/100), better than model selection. 46
57 Implementing and scoring BMA for Multivariate extremes Does the BMA framework perform signicantly better than selecting models based on AIC? Yes, in terms of logarithmic score for the predictive density (Kullback-Leibler divergence to the truth) Madigan & Raftery (94) In average over the union model, w.r.t prior! Simulation study : Evaluation via proper scoring rules (Logarithmic + probability of failure regions) 2 models of same dimension: Pairwise-Beta / Nested asymmetric logistic 100 data sets simulated from another model Results: The BMA framework performs consistently (for all scores), slightly (1/20 to 1/100), better than model selection. 46
58 Implementing and scoring BMA for Multivariate extremes Does the BMA framework perform signicantly better than selecting models based on AIC? Yes, in terms of logarithmic score for the predictive density (Kullback-Leibler divergence to the truth) Madigan & Raftery (94) In average over the union model, w.r.t prior! Simulation study : Evaluation via proper scoring rules (Logarithmic + probability of failure regions) 2 models of same dimension: Pairwise-Beta / Nested asymmetric logistic 100 data sets simulated from another model Results: The BMA framework performs consistently (for all scores), slightly (1/20 to 1/100), better than model selection. 46
59 Implementing and scoring BMA for Multivariate extremes Does the BMA framework perform signicantly better than selecting models based on AIC? Yes, in terms of logarithmic score for the predictive density (Kullback-Leibler divergence to the truth) Madigan & Raftery (94) In average over the union model, w.r.t prior! Simulation study : Evaluation via proper scoring rules (Logarithmic + probability of failure regions) 2 models of same dimension: Pairwise-Beta / Nested asymmetric logistic 100 data sets simulated from another model Results: The BMA framework performs consistently (for all scores), slightly (1/20 to 1/100), better than model selection. 46
60 Discussion BMA vs selection: Moderate gain for large sample size : Posterior concentration on `asymptotic carrier regions' = points (parameters) of minimal KL divergence from truth BMA : simple if several models have already been tted (`only' compute posterior weights) Way out: Mixture models for increased dimension of the parameter space. (product) 47
Non parametric modeling of multivariate extremes with Dirichlet mixtures
Non parametric modeling of multivariate extremes with Dirichlet mixtures Anne Sabourin Institut Mines-Télécom, Télécom ParisTech, CNRS-LTCI Joint work with Philippe Naveau (LSCE, Saclay), Anne-Laure Fougères
More informationBayesian nonparametrics for multivariate extremes including censored data. EVT 2013, Vimeiro. Anne Sabourin. September 10, 2013
Bayesian nonparametrics for multivariate extremes including censored data Anne Sabourin PhD advisors: Anne-Laure Fougères (Lyon 1), Philippe Naveau (LSCE, Saclay). Joint work with Benjamin Renard, IRSTEA,
More informationCOMBINING REGIONAL ESTIMATION AND HISTORICAL FLOODS: A MULTIVARIATE SEMI-PARAMETRIC PEAKS-OVER-THRESHOLD MODEL WITH CENSORED DATA
COMBINING REGIONAL ESTIMATION AND HISTORICAL FLOODS: A MULTIVARIATE SEMI-PARAMETRIC PEAKS-OVER-THRESHOLD MODEL WITH CENSORED DATA A. SABOURIN AND B. RENARD Multivariate extremes; censored data; semi-parametric
More informationarxiv: v1 [stat.ap] 28 Nov 2014
COMBINING REGIONAL ESTIMATION AND HISTORICAL FLOODS: A MULTIVARIATE SEMI-PARAMETRIC PEAKS-OVER-THRESHOLD MODEL WITH CENSORED DATA Anne Sabourin 1, Benjamin Renard 2 arxiv:1411.7782v1 [stat.ap] 28 Nov 2014
More informationBayesian Model Averaging for Multivariate Extreme Values
Bayesian Model Averaging for Multivariate Extreme Values Philippe Naveau naveau@lsce.ipsl.fr Laboratoire des Sciences du Climat et l Environnement (LSCE) Gif-sur-Yvette, France joint work with A. Sabourin
More informationCombining regional estimation and historical floods: a multivariate semi-parametric peaks-over-threshold model with censored data
Combining regional estimation and historical floods: a multivariate semi-parametric peaks-over-threshold model with censored data Anne Sabourin, Benjamin Renard To cite this version: Anne Sabourin, Benjamin
More informationOverview of Extreme Value Theory. Dr. Sawsan Hilal space
Overview of Extreme Value Theory Dr. Sawsan Hilal space Maths Department - University of Bahrain space November 2010 Outline Part-1: Univariate Extremes Motivation Threshold Exceedances Part-2: Bivariate
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 informationWorkshop Copulas and Extremes
Workshop Copulas and Extremes Organization : Marianne Clausel & Stéphane Girard November, 19-20th, 2013 Tuesday, 9:00 9:30, Welcome coffee. Tuesday, 9:30 12:00, Session I, chair Anne-Catherine Favre. Some
More informationBayesian Inference for Clustered Extremes
Newcastle University, Newcastle-upon-Tyne, U.K. lee.fawcett@ncl.ac.uk 20th TIES Conference: Bologna, Italy, July 2009 Structure of this talk 1. Motivation and background 2. Review of existing methods Limitations/difficulties
More informationNew Classes of Multivariate Survival Functions
Xiao Qin 2 Richard L. Smith 2 Ruoen Ren School of Economics and Management Beihang University Beijing, China 2 Department of Statistics and Operations Research University of North Carolina Chapel Hill,
More informationBayesian inference for multivariate extreme value distributions
Bayesian inference for multivariate extreme value distributions Sebastian Engelke Clément Dombry, Marco Oesting Toronto, Fields Institute, May 4th, 2016 Main motivation For a parametric model Z F θ of
More informationSparse Representation of Multivariate Extremes with Applications to Anomaly Ranking
Sparse Representation of Multivariate Extremes with Applications to Anomaly Ranking Nicolas Goix Anne Sabourin Stéphan Clémençon Abstract Capturing the dependence structure of multivariate extreme events
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 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 informationZwiers FW and Kharin VV Changes in the extremes of the climate simulated by CCC GCM2 under CO 2 doubling. J. Climate 11:
Statistical Analysis of EXTREMES in GEOPHYSICS Zwiers FW and Kharin VV. 1998. Changes in the extremes of the climate simulated by CCC GCM2 under CO 2 doubling. J. Climate 11:2200 2222. http://www.ral.ucar.edu/staff/ericg/readinggroup.html
More informationESTIMATING BIVARIATE TAIL
Elena DI BERNARDINO b joint work with Clémentine PRIEUR a and Véronique MAUME-DESCHAMPS b a LJK, Université Joseph Fourier, Grenoble 1 b Laboratoire SAF, ISFA, Université Lyon 1 Framework Goal: estimating
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning MCMC and Non-Parametric Bayes Mark Schmidt University of British Columbia Winter 2016 Admin I went through project proposals: Some of you got a message on Piazza. No news is
More informationBivariate generalized Pareto distribution
Bivariate generalized Pareto distribution in practice Eötvös Loránd University, Budapest, Hungary Minisymposium on Uncertainty Modelling 27 September 2011, CSASC 2011, Krems, Austria Outline Short summary
More informationSome conditional extremes of a Markov chain
Some conditional extremes of a Markov chain Seminar at Edinburgh University, November 2005 Adam Butler, Biomathematics & Statistics Scotland Jonathan Tawn, Lancaster University Acknowledgements: Janet
More informationApproximate Bayesian computation for spatial extremes via open-faced sandwich adjustment
Approximate Bayesian computation for spatial extremes via open-faced sandwich adjustment Ben Shaby SAMSI August 3, 2010 Ben Shaby (SAMSI) OFS adjustment August 3, 2010 1 / 29 Outline 1 Introduction 2 Spatial
More informationExtreme value statistics: from one dimension to many. Lecture 1: one dimension Lecture 2: many dimensions
Extreme value statistics: from one dimension to many Lecture 1: one dimension Lecture 2: many dimensions The challenge for extreme value statistics right now: to go from 1 or 2 dimensions to 50 or more
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 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 informationAPPLICATION OF EXTREMAL THEORY TO THE PRECIPITATION SERIES IN NORTHERN MORAVIA
APPLICATION OF EXTREMAL THEORY TO THE PRECIPITATION SERIES IN NORTHERN MORAVIA DANIELA JARUŠKOVÁ Department of Mathematics, Czech Technical University, Prague; jarus@mat.fsv.cvut.cz 1. Introduction The
More informationBayesian Modeling of Conditional Distributions
Bayesian Modeling of Conditional Distributions John Geweke University of Iowa Indiana University Department of Economics February 27, 2007 Outline Motivation Model description Methods of inference Earnings
More informationExtreme Precipitation: An Application Modeling N-Year Return Levels at the Station Level
Extreme Precipitation: An Application Modeling N-Year Return Levels at the Station Level Presented by: Elizabeth Shamseldin Joint work with: Richard Smith, Doug Nychka, Steve Sain, Dan Cooley Statistics
More informationInvariant HPD credible sets and MAP estimators
Bayesian Analysis (007), Number 4, pp. 681 69 Invariant HPD credible sets and MAP estimators Pierre Druilhet and Jean-Michel Marin Abstract. MAP estimators and HPD credible sets are often criticized in
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
More informationData. Climate model data from CMIP3
Data Observational data from CRU (Climate Research Unit, University of East Anglia, UK) monthly averages on 5 o x5 o grid boxes, aggregated to JJA average anomalies over Europe: spatial averages over 10
More informationStatistics for extreme & sparse data
Statistics for extreme & sparse data University of Bath December 6, 2018 Plan 1 2 3 4 5 6 The Problem Climate Change = Bad! 4 key problems Volcanic eruptions/catastrophic event prediction. Windstorms
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods By Oleg Makhnin 1 Introduction a b c M = d e f g h i 0 f(x)dx 1.1 Motivation 1.1.1 Just here Supresses numbering 1.1.2 After this 1.2 Literature 2 Method 2.1 New math As
More informationApproximating the Conditional Density Given Large Observed Values via a Multivariate Extremes Framework, with Application to Environmental Data
Approximating the Conditional Density Given Large Observed Values via a Multivariate Extremes Framework, with Application to Environmental Data Daniel Cooley 1 Richard A. Davis 2 Philippe Naveau 3 1 Department
More informationThe Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations
The Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations John R. Michael, Significance, Inc. and William R. Schucany, Southern Methodist University The mixture
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 informationSTAT 425: Introduction to Bayesian Analysis
STAT 425: Introduction to Bayesian Analysis Marina Vannucci Rice University, USA Fall 2017 Marina Vannucci (Rice University, USA) Bayesian Analysis (Part 2) Fall 2017 1 / 19 Part 2: Markov chain Monte
More informationThe Bayesian Choice. Christian P. Robert. From Decision-Theoretic Foundations to Computational Implementation. Second Edition.
Christian P. Robert The Bayesian Choice From Decision-Theoretic Foundations to Computational Implementation Second Edition With 23 Illustrations ^Springer" Contents Preface to the Second Edition Preface
More informationMarkov Chain Monte Carlo (MCMC)
Markov Chain Monte Carlo (MCMC Dependent Sampling Suppose we wish to sample from a density π, and we can evaluate π as a function but have no means to directly generate a sample. Rejection sampling can
More informationCTDL-Positive Stable Frailty Model
CTDL-Positive Stable Frailty Model M. Blagojevic 1, G. MacKenzie 2 1 Department of Mathematics, Keele University, Staffordshire ST5 5BG,UK and 2 Centre of Biostatistics, University of Limerick, Ireland
More informationUniversity of Toronto Department of Statistics
Norm Comparisons for Data Augmentation by James P. Hobert Department of Statistics University of Florida and Jeffrey S. Rosenthal Department of Statistics University of Toronto Technical Report No. 0704
More informationMultivariate generalized Pareto distributions
Multivariate generalized Pareto distributions Holger Rootzén and Nader Tajvidi Abstract Statistical inference for extremes has been a subject of intensive research during the past couple of decades. One
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 informationMULTIVARIATE EXTREMES AND RISK
MULTIVARIATE EXTREMES AND RISK Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu Interface 2008 RISK: Reality Durham,
More informationBernstein polynomial angular densities of multivariate extreme value distributions
Bernstein polynomial angular densities of multivariate extreme value distributions Timothy E. Hanson Department of Statistics, University of South Carolina, Columbia, SC 29208, USA email: hansont@stat.sc.edu
More informationOn the occurrence times of componentwise maxima and bias in likelihood inference for multivariate max-stable distributions
On the occurrence times of componentwise maxima and bias in likelihood inference for multivariate max-stable distributions J. L. Wadsworth Department of Mathematics and Statistics, Fylde College, Lancaster
More information16 : Approximate Inference: Markov Chain Monte Carlo
10-708: Probabilistic Graphical Models 10-708, Spring 2017 16 : Approximate Inference: Markov Chain Monte Carlo Lecturer: Eric P. Xing Scribes: Yuan Yang, Chao-Ming Yen 1 Introduction As the target distribution
More informationBayesian Point Process Modeling for Extreme Value Analysis, with an Application to Systemic Risk Assessment in Correlated Financial Markets
Bayesian Point Process Modeling for Extreme Value Analysis, with an Application to Systemic Risk Assessment in Correlated Financial Markets Athanasios Kottas Department of Applied Mathematics and Statistics,
More informationSTATISTICAL INFERENCE WITH DATA AUGMENTATION AND PARAMETER EXPANSION
STATISTICAL INFERENCE WITH arxiv:1512.00847v1 [math.st] 2 Dec 2015 DATA AUGMENTATION AND PARAMETER EXPANSION Yannis G. Yatracos Faculty of Communication and Media Studies Cyprus University of Technology
More informationMULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES
MULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES Philip Jonathan, Kevin Ewans, David Randell, Yanyun Wu philip.jonathan@shell.com www.lancs.ac.uk/ jonathan Wave Hindcasting & Forecasting
More informationLecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE. Rick Katz
1 Lecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home
More informationMachine Learning using Bayesian Approaches
Machine Learning using Bayesian Approaches Sargur N. Srihari University at Buffalo, State University of New York 1 Outline 1. Progress in ML and PR 2. Fully Bayesian Approach 1. Probability theory Bayes
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 211 http://stat.epfl.ch Multivariate Extremes 19 Componentwise maxima.................................................. 194 Standardization........................................................
More informationA Fully Nonparametric Modeling Approach to. BNP Binary Regression
A Fully Nonparametric Modeling Approach to Binary Regression Maria Department of Applied Mathematics and Statistics University of California, Santa Cruz SBIES, April 27-28, 2012 Outline 1 2 3 Simulation
More informationRonald Christensen. University of New Mexico. Albuquerque, New Mexico. Wesley Johnson. University of California, Irvine. Irvine, California
Texts in Statistical Science Bayesian Ideas and Data Analysis An Introduction for Scientists and Statisticians Ronald Christensen University of New Mexico Albuquerque, New Mexico Wesley Johnson University
More informationEstimating Bivariate Tail: a copula based approach
Estimating Bivariate Tail: a copula based approach Elena Di Bernardino, Université Lyon 1 - ISFA, Institut de Science Financiere et d'assurances - AST&Risk (ANR Project) Joint work with Véronique Maume-Deschamps
More informationTwo practical tools for rainfall weather generators
Two practical tools for rainfall weather generators Philippe Naveau naveau@lsce.ipsl.fr Laboratoire des Sciences du Climat et l Environnement (LSCE) Gif-sur-Yvette, France FP7-ACQWA, GIS-PEPER, MIRACLE
More informationA Bayesian perspective on GMM and IV
A Bayesian perspective on GMM and IV Christopher A. Sims Princeton University sims@princeton.edu November 26, 2013 What is a Bayesian perspective? A Bayesian perspective on scientific reporting views all
More informationPenalized Loss functions for Bayesian Model Choice
Penalized Loss functions for Bayesian Model Choice Martyn International Agency for Research on Cancer Lyon, France 13 November 2009 The pure approach For a Bayesian purist, all uncertainty is represented
More informationModels for Spatial Extremes. Dan Cooley Department of Statistics Colorado State University. Work supported in part by NSF-DMS
Models for Spatial Extremes Dan Cooley Department of Statistics Colorado State University Work supported in part by NSF-DMS 0905315 Outline 1. Introduction Statistics for Extremes Climate and Weather 2.
More informationSpatial Extremes in Atmospheric Problems
Spatial Extremes in Atmospheric Problems Eric Gilleland Research Applications Laboratory (RAL) National Center for Atmospheric Research (NCAR), Boulder, Colorado, U.S.A. http://www.ral.ucar.edu/staff/ericg
More informationModelling Multivariate Peaks-over-Thresholds using Generalized Pareto Distributions
Modelling Multivariate Peaks-over-Thresholds using Generalized Pareto Distributions Anna Kiriliouk 1 Holger Rootzén 2 Johan Segers 1 Jennifer L. Wadsworth 3 1 Université catholique de Louvain (BE) 2 Chalmers
More informationGeneralized additive modelling of hydrological sample extremes
Generalized additive modelling of hydrological sample extremes Valérie Chavez-Demoulin 1 Joint work with A.C. Davison (EPFL) and Marius Hofert (ETHZ) 1 Faculty of Business and Economics, University of
More informationModels and estimation.
Bivariate generalized Pareto distribution practice: Models and estimation. Eötvös Loránd University, Budapest, Hungary 7 June 2011, ASMDA Conference, Rome, Italy Problem How can we properly estimate the
More informationA spatio-temporal model for extreme precipitation simulated by a climate model
A spatio-temporal model for extreme precipitation simulated by a climate model Jonathan Jalbert Postdoctoral fellow at McGill University, Montréal Anne-Catherine Favre, Claude Bélisle and Jean-François
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 informationStat260: Bayesian Modeling and Inference Lecture Date: February 10th, Jeffreys priors. exp 1 ) p 2
Stat260: Bayesian Modeling and Inference Lecture Date: February 10th, 2010 Jeffreys priors Lecturer: Michael I. Jordan Scribe: Timothy Hunter 1 Priors for the multivariate Gaussian Consider a multivariate
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 informationOn Reparametrization and the Gibbs Sampler
On Reparametrization and the Gibbs Sampler Jorge Carlos Román Department of Mathematics Vanderbilt University James P. Hobert Department of Statistics University of Florida March 2014 Brett Presnell Department
More informationIndex. Pagenumbersfollowedbyf indicate figures; pagenumbersfollowedbyt indicate tables.
Index Pagenumbersfollowedbyf indicate figures; pagenumbersfollowedbyt indicate tables. Adaptive rejection metropolis sampling (ARMS), 98 Adaptive shrinkage, 132 Advanced Photo System (APS), 255 Aggregation
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 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 informationObnoxious lateness humor
Obnoxious lateness humor 1 Using Bayesian Model Averaging For Addressing Model Uncertainty in Environmental Risk Assessment Louise Ryan and Melissa Whitney Department of Biostatistics Harvard School of
More informationOverall Objective Priors
Overall Objective Priors Jim Berger, Jose Bernardo and Dongchu Sun Duke University, University of Valencia and University of Missouri Recent advances in statistical inference: theory and case studies University
More informationFinancial Econometrics and Volatility Models Extreme Value Theory
Financial Econometrics and Volatility Models Extreme Value Theory Eric Zivot May 3, 2010 1 Lecture Outline Modeling Maxima and Worst Cases The Generalized Extreme Value Distribution Modeling Extremes Over
More informationBayesian Regularization
Bayesian Regularization Aad van der Vaart Vrije Universiteit Amsterdam International Congress of Mathematicians Hyderabad, August 2010 Contents Introduction Abstract result Gaussian process priors Co-authors
More informationarxiv: v1 [stat.me] 16 Apr 2012
A Euclidean Likelihood Estimator for Bivariate Tail Dependence Miguel de Carvalho, Boris Oumow, Johan Segers, and Michał Warchoł arxiv:124.3524v1 [stat.me] 16 Apr 212 September 19, 218 Abstract The spectral
More informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationSharp statistical tools Statistics for extremes
Sharp statistical tools Statistics for extremes Georg Lindgren Lund University October 18, 2012 SARMA Background Motivation We want to predict outside the range of observations Sums, averages and proportions
More informationBayesian Multivariate Extreme Value Thresholding for Environmental Hazards
Bayesian Multivariate Extreme Value Thresholding for Environmental Hazards D. Lupton K. Abayomi M. Lacer School of Industrial and Systems Engineering Georgia Institute of Technology Institute for Operations
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationBayesian Inference for the Multivariate Normal
Bayesian Inference for the Multivariate Normal Will Penny Wellcome Trust Centre for Neuroimaging, University College, London WC1N 3BG, UK. November 28, 2014 Abstract Bayesian inference for the multivariate
More informationContents. Part I: Fundamentals of Bayesian Inference 1
Contents Preface xiii Part I: Fundamentals of Bayesian Inference 1 1 Probability and inference 3 1.1 The three steps of Bayesian data analysis 3 1.2 General notation for statistical inference 4 1.3 Bayesian
More informationPREPRINT 2005:38. Multivariate Generalized Pareto Distributions HOLGER ROOTZÉN NADER TAJVIDI
PREPRINT 2005:38 Multivariate Generalized Pareto Distributions HOLGER ROOTZÉN NADER TAJVIDI Department of Mathematical Sciences Division of Mathematical Statistics CHALMERS UNIVERSITY OF TECHNOLOGY GÖTEBORG
More informationPractical Bayesian Quantile Regression. Keming Yu University of Plymouth, UK
Practical Bayesian Quantile Regression Keming Yu University of Plymouth, UK (kyu@plymouth.ac.uk) A brief summary of some recent work of us (Keming Yu, Rana Moyeed and Julian Stander). Summary We develops
More informationA Search and Jump Algorithm for Markov Chain Monte Carlo Sampling. Christopher Jennison. Adriana Ibrahim. Seminar at University of Kuwait
A Search and Jump Algorithm for Markov Chain Monte Carlo Sampling Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj Adriana Ibrahim Institute
More informationSemi-parametric estimation of non-stationary Pickands functions
Semi-parametric estimation of non-stationary Pickands functions Linda Mhalla 1 Joint work with: Valérie Chavez-Demoulin 2 and Philippe Naveau 3 1 Geneva School of Economics and Management, University of
More informationPseudo-marginal MCMC methods for inference in latent variable models
Pseudo-marginal MCMC methods for inference in latent variable models Arnaud Doucet Department of Statistics, Oxford University Joint work with George Deligiannidis (Oxford) & Mike Pitt (Kings) MCQMC, 19/08/2016
More informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
More informationExpectation Propagation for Approximate Bayesian Inference
Expectation Propagation for Approximate Bayesian Inference José Miguel Hernández Lobato Universidad Autónoma de Madrid, Computer Science Department February 5, 2007 1/ 24 Bayesian Inference Inference Given
More informationStat 451 Lecture Notes Markov Chain Monte Carlo. Ryan Martin UIC
Stat 451 Lecture Notes 07 12 Markov Chain Monte Carlo Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapters 8 9 in Givens & Hoeting, Chapters 25 27 in Lange 2 Updated: April 4, 2016 1 / 42 Outline
More informationBayesian nonparametric estimation of finite population quantities in absence of design information on nonsampled units
Bayesian nonparametric estimation of finite population quantities in absence of design information on nonsampled units Sahar Z Zangeneh Robert W. Keener Roderick J.A. Little Abstract In Probability proportional
More informationEstimation and Maintenance of Measurement Rates for Multiple Extended Target Tracking
FUSION 2012, Singapore 118) Estimation and Maintenance of Measurement Rates for Multiple Extended Target Tracking Karl Granström*, Umut Orguner** *Division of Automatic Control Department of Electrical
More informationStatistical Methods. Missing Data snijders/sm.htm. Tom A.B. Snijders. November, University of Oxford 1 / 23
1 / 23 Statistical Methods Missing Data http://www.stats.ox.ac.uk/ snijders/sm.htm Tom A.B. Snijders University of Oxford November, 2011 2 / 23 Literature: Joseph L. Schafer and John W. Graham, Missing
More informationBayesian Learning. HT2015: SC4 Statistical Data Mining and Machine Learning. Maximum Likelihood Principle. The Bayesian Learning Framework
HT5: SC4 Statistical Data Mining and Machine Learning Dino Sejdinovic Department of Statistics Oxford http://www.stats.ox.ac.uk/~sejdinov/sdmml.html Maximum Likelihood Principle A generative model for
More informationConstruction of an Informative Hierarchical Prior Distribution: Application to Electricity Load Forecasting
Construction of an Informative Hierarchical Prior Distribution: Application to Electricity Load Forecasting Anne Philippe Laboratoire de Mathématiques Jean Leray Université de Nantes Workshop EDF-INRIA,
More informationStat 535 C - Statistical Computing & Monte Carlo Methods. Arnaud Doucet.
Stat 535 C - Statistical Computing & Monte Carlo Methods Arnaud Doucet Email: arnaud@cs.ubc.ca 1 1.1 Outline Introduction to Markov chain Monte Carlo The Gibbs Sampler Examples Overview of the Lecture
More informationBayesian Methods for Estimating the Reliability of Complex Systems Using Heterogeneous Multilevel Information
Statistics Preprints Statistics 8-2010 Bayesian Methods for Estimating the Reliability of Complex Systems Using Heterogeneous Multilevel Information Jiqiang Guo Iowa State University, jqguo@iastate.edu
More informationOn the spatial dependence of extreme ocean storm seas
Accepted by Ocean Engineering, August 2017 On the spatial dependence of extreme ocean storm seas Emma Ross a, Monika Kereszturi b, Mirrelijn van Nee c, David Randell d, Philip Jonathan a, a Shell Projects
More informationContents. Preface to Second Edition Preface to First Edition Abbreviations PART I PRINCIPLES OF STATISTICAL THINKING AND ANALYSIS 1
Contents Preface to Second Edition Preface to First Edition Abbreviations xv xvii xix PART I PRINCIPLES OF STATISTICAL THINKING AND ANALYSIS 1 1 The Role of Statistical Methods in Modern Industry and Services
More informationFitting Narrow Emission Lines in X-ray Spectra
Outline Fitting Narrow Emission Lines in X-ray Spectra Taeyoung Park Department of Statistics, University of Pittsburgh October 11, 2007 Outline of Presentation Outline This talk has three components:
More information