Semi-parametric estimation of non-stationary Pickands functions
|
|
- MargaretMargaret Bennett
- 5 years ago
- Views:
Transcription
1 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 Geneva, Switzerland 2 Faculty of Business and Economics, University of Lausanne, Switzerland 3 Laboratoire des Sciences du Climat et l'environnement, France NCAR, Boulder 2016 Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
2 Motivation Modelling multivariate extremes: univariate marginal distributions dependence in extremes (extremal dependence) One should be able to include covariates effects in the model: climate change, seasonality,... univariate case: regression methods have been proposed to model the parameters of extreme values distributions how can we relate the extremal dependence to a set of covariates? Aim: provide a very flexible estimation of a covariate-dependent extremal dependence. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
3 Overview 1 Multivariate extremes framework 2 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures Corrected Pickands dependence function estimator 3 Simulation Study No covariate dependence Covariate-dependent Pickands function 4 Application 5 Conclusion 6 References Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
4 Multivariate extremes framework 1 Multivariate extremes framework 2 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures Corrected Pickands dependence function estimator 3 Simulation Study No covariate dependence Covariate-dependent Pickands function 4 Application 5 Conclusion 6 References Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
5 Multivariate extremes framework Multivariate Extreme Value Theory Assume Z = (Z 1,..., Z d ) T has a multivariate max-stable distribution with unit-fréchet margins. Then { ( ) } ωi pr(z z) = exp { V (z)} = exp dh(ω 1,..., ω d ), S d 1 d max i=1 where S d 1 = { (ω 1,..., ω d ) R d + : ω ω d = 1 } and H is a finite measure that satisfies S d 1 ω i dh(ω 1,..., ω d ) = 1 for each i. Exponent function V (z) = z i ( ) A(ω), z 1 z d where ω = (z 1 / z,..., z d / z ) T S d 1. A(ω) is called the Pickands dependence function. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
6 Multivariate extremes framework Pickands Dependence function Necessary conditions on A( ) C1) A(ω) is continuous and convex C2) For any ω S d 1 max (ω 1,..., ω d ) A(ω) 1 C3) A(0,..., 0) = 1 and A(e i ) = 1 with e i = (0,..., 0, 1, 0,..., 0), for i = 1,..., d 1 A(ω) ω Many regularised non-parametric estimators have been proposed. But what if A(ω) varies with a set of covariates X? Benefit from the well-established framework of Generalized Additive Models. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
7 Multivariate extremes framework Generalized Additive Models 1 A model for Y conditional on X R q E (Y X = x) = g { u β + } g a link function K h k (t k ) (u 1,..., u p ) and (t 1,..., t K ) subsets of (x 1,..., x q ) k=1 h k C 2 (A k ) with A k R admits β R p a vector of parameters h k are smooth functions a finite (m k -)dim basis parametrized by h k = (h 1,, h mk ) R m k a quadratic penalty representation A k h k (t)2 dt = h T k S kh k 1 Textbooks: [Hastie and Tibshirani, 1990, Green and Silverman, 1993] Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
8 Estimation of non-stationary Pickands functions 1 Multivariate extremes framework 2 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures Corrected Pickands dependence function estimator 3 Simulation Study No covariate dependence Covariate-dependent Pickands function 4 Application 5 Conclusion 6 References Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
9 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures GAM for max-stable dependence structures Let Z(x) = (Z 1 (x),..., Z d (x)) T denote a d-dimensional max-stable process with margins transformed using a regression model to be unit Fréchet, i.e., F (z) = exp( 1/z) a covariate-varying extremal dependence, i.e., A(ω x) Define the max-projection of Z(x) as [ ] Y (ω) x = max F λ 1 {Z 1 (x)},..., F λ d {Z d (x)} with λ i = ω i ( 1 ω ω d ) 1 and ω = (ω1,..., ω d ) T S d 1 with non-zero entries. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
10 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures PROPOSITION The random variable Y (ω) x is Beta distributed with parameter Aω(x), i.e., pr(y (ω) y x) = y A ω(x), for 0 y 1. In particular, E [Y (ω) x] = g 1 {Aω(x)} where the concave link function equals g(u) = u 1 u with 0 < u < 1. The Pickands function Aω(x) is the transformed mean of Y (ω) x and can be related to a set of predictors using GAMs. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
11 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures A non-stationary extremal dependence The generalized additive model for Aω(x) A ω (x; θ) = u β + K h k (t k ) k=1 (u 1,..., u p ) and (t 1,..., t K ) subsets of {x 1,..., x q } β R p a vector of parameters h k C 2 parametrized by h k = (h 1,..., h mk ) R m k admits a quadratic penalty representation θ = (β, h 1,..., h K ) Θ R d with d = p + K k=1 m k Estimation procedure: Such models are fitted by penalized likelihood maximization. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
12 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures Penalized log-likelihood The penalized log-likelihood estimator for Aω(x) We define the penalized log-likelihood as l p (θ, γ; ω) = lω(θ) 1 2 K γ k h k (t k ) 2 dt k = lω(θ) 1 K γ k h k S k h k A k 2 k=1 and the penalized maximum log-likelihood estimator as θ n := argmax θ Θ l p (θ, γ; ω). k=1 γ = (γ 1,..., γ K ) R K + Observations at hand {y i (ω), x i } n i=1, then n lω(θ) = [{Aω(x i ; θ) 1} log y i (ω) + log Aω(x i ; θ)] i=1 Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
13 Estimation of non-stationary Pickands functions Corrected Pickands dependence function estimator Corrected Pickands function estimator Up to this point, the estimation procedure is valid for any d 2 and results in a covariate-varying Pickands function estimator Ã(ω x) := Aω(x; θ n ) Ã(ω x) ω = 0.98 ω = 0.8 ω = 0.6 ω = 0.4 ω = 0.2 ω = 0.01 x Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
14 Estimation of non-stationary Pickands functions Corrected Pickands dependence function estimator Corrected Pickands function estimator Up to this point, the estimation procedure is valid for any d 2 and results in a covariate-varying Pickands function estimator Ã(ω x) := Aω(x; θ n ) Ã(ω x) ω = 0.98 ω = 0.8 ω = 0.6 ω = 0.4 ω = 0.2 ω = 0.01 x But, Ã( x) is not a valid Pickands function! Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
15 Estimation of non-stationary Pickands functions Corrected Pickands dependence function estimator Corrected Pickands function estimator Regularization procedure when d = 2: 1 compute à {(ω i, 1 ω i ) x} = Ã(ω i x), for ω 1,..., ωñ in [0, 1] 2 apply the constrained median smoothing approach (cobs) 2 ĝ γ,l (ω x) = min g ñ i=1 Ã(ω i x) g(ω i x) + γ max g (ω x), ω with conditions C1-C2-C3 on the smoothing spline g(ω x) 3 define Â(ω x) := ĝ γ,l (ω x) Confidence intervals: Use the resampling cases bootstrap 3 : only assumes independence. Remark: Confidence intervals are not guaranteed to be convex. 2 [Ng and Maechler, 2007] 3 [Davison and Hinkley, 1997, Section 6.2.4] Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
16 Simulation Study 1 Multivariate extremes framework 2 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures Corrected Pickands dependence function estimator 3 Simulation Study No covariate dependence Covariate-dependent Pickands function 4 Application 5 Conclusion 6 References Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
17 Simulation Study No covariate dependence No covariate dependence The aim is to evaluate the impact of the max-projection step along with the constrained smoothing procedure on the estimation of the extremal dependence. The first stage of the estimation procedure boils down to maximising the log-likelihood l (A ω ) = n (A ω 1) log {y i (ω)} + log A ω. i=1 The constrained median smoothing approach is applied to the resulting MLEs for different values of ω [0, 1]. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
18 Simulation Study No covariate dependence True dependence model: Dirichlet dependence model with β = 3 and different values of α n = 1000 observations ñ = 800 regularly spaced values of ω in [0, 1] Asymptotic 95% confidence intervals α = 0.02 α = 0.2 α = 10 A(ω) True MLE 95% CI Pickands CFG A(ω) True MLE 95% CI Pickands CFG A(ω) True MLE 95% CI Pickands CFG ω ω ω Our estimator is competitive with the [Pickands, 1981] and the [Capéraà et al., 1997] estimators. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
19 Simulation Study Covariate-dependent Pickands function Logistic dependence with two covariates True dependence model: bivariate logistic dependence model with a covariate-dependent parameter α(t, x) = 1 {x= cata } α 0 sin(α 1 t) + 1 {x= catb } { β0 sin(β 1 t 2 ) + β 3 }, where t is a continuous covariate discretized into 12 equally spaced values, and x a factor covariate with two levels R = 1000 bootstrap samples n = 400 observations for each distinct pair of covariates ñ = 800 regularly spaced values of ω in [0, 1] Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
20 Simulation Study Covariate-dependent Pickands function A(ω) t ω A(ω) t ω A(ω) t ω A(ω) t ω A(ω) t ω A(ω) t ω Figure: Estimates of the Pickands function for the logistic model with weak, mild and strong dependence for cata (top) and catb (bottom) and their 95% bootstrap confidence intervals. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
21 Application 1 Multivariate extremes framework 2 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures Corrected Pickands dependence function estimator 3 Simulation Study No covariate dependence Covariate-dependent Pickands function 4 Application 5 Conclusion 6 References Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
22 Application Application Data description The data represents daily minimum and maximum temperature for long term stations from the U.S. COOP station network 4, in degrees Celcius. The temperatures were extracted for stations from 46 states and from to Aim: We are interested in the extremal dependence between daily maxima and daily minima temperatures, i.e., between hot days and warm nights. More specifically, how did this extremal dependence evolve through time? 4 [Climate and Global Dynamics Division, 2010] Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
23 Application Data analysis: Colorado We use daily maximum and minimum temperatures over all the 151 stations in Colorado. We model marginally monthly maxima of each time series. Temperature ( C) Temperature ( C) Time (years) Time (years) Figure: Monthly maxima of daily maximum temperatures (top) and daily minimum temperatures (bottom). Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
24 Application Deal with non-stationarity marginally µ max(t, m) = σ max(t, m) = 12 j=1 12 j=1 ξ max(t, m) = ξ maxt (µ 0,j + µ 1,j t) 1 {m=j} σ 0,j 1 {m=j} µ min (t, m) = σ min (t, m) = 12 j=1 12 j=1 ξ min (t, m) = ξ min (µ 0,j + µ 1,j t) 1 {m=j} σ 0,j 1 {m=j} where t and m are the year and month of the observed extreme temperature, respectively. Transform the data marginally to a stationary common scale (standardized Gumbel, unit Fréchet, etc). Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
25 Application Modelling time-varying Pickands dependence function 1.00 Pickands function ω ) 1970 ars 0.8 ye 1960 e( tim Figure: Time-dependent Pickands dependence function. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
26 Application Modelling time-varying Pickands dependence function 1.00 Pickands function ω ) 1970 ars 0.8 ye 1960 e( tim What happened in 1970? [Meehl et al., 2009]: Since 1970, less extreme low temperatures than expected (warmer nights). Figure: Time-dependent Pickands dependence function. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
27 Application Covariate-dependent extremal coefficient estimates Extremal coefficient Colorado Time (years) Extremal coefficient California Time (years) Extremal coefficient Alabama Time (years) Figure: Extremal coefficient estimates (smoothed) with 95% bootstrap confidence intervals (R=300). Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
28 Application Question: How does the extremal dependence evolve through time and altitude? Extremal coefficient Time (years) Altitude (m) Figure: Time and altitude dependent extremal coefficient. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
29 Conclusion 1 Multivariate extremes framework 2 Estimation of non-stationary Pickands functions GAM for max-stable dependence structures Corrected Pickands dependence function estimator 3 Simulation Study No covariate dependence Covariate-dependent Pickands function 4 Application 5 Conclusion 6 References Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
30 Conclusion Conclusion Very flexible framework for the estimation of a valid non-stationary Pickands dependence function in the bivariate case. The estimation procedure is valid in both the stationary and non-stationary cases. Extensions/limitations: Extension of the regularisation procedure to d > 2. Depending on how fine the discretization of the unit simplex is, computational times can be substantial! Application: What other covariates can we include in our model? mean temperatures? ideas? Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
31 References References Capéraà, P., Fougères, A.-L., and Genest, C. (1997). Thanks for your attention! A nonparametric estimation procedure for bivariate extreme value copulas. BIOMETRIKA, 84: Climate and Global Dynamics Division, National Center for Atmospheric Research, U. C. f. A. R. (2010). Daily minimum and maximum temperature and precipitation for long term stations from the u.s. coop data. accessed 7 april Davison, A. and Hinkley, D. (1997). Bootstrap Methods and Their Application. Cambridge Series in Statistical and Probabilistic Mathematics. Cambridge University Press. Green, P. and Silverman, B. (1993). Nonparametric Regression and Generalized Linear Models: A roughness penalty approach. Chapman & Hall/CRC Monographs on Statistics & Applied Probability. Taylor & Francis. Hastie, T. J. and Tibshirani, R. J. (1990). Generalized additive models. London: Chapman & Hall. Meehl, G. A., Tebaldi, C., Walton, G., Easterling, D., and McDaniel, L. (2009). Relative increase of record high maximum temperatures compared to record low minimum temperatures in the u.s. Geophysical Research Letters, 36. Ng, P. and Maechler, M. (2007). A fast and efficient implementation of qualitatively constrained quantile smoothing splines. Statistical Modelling, 7: Pickands, J. (1981). Multivariate extreme value distributions. volume 49, pages , Proceedings of the 43rd Session of the International Statistical Institute, Vol.2 (Buenos Aires, 1981). Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
32 References Marginal QQ-plots Empirical quantiles Theoretical quantiles Empirical quantiles Theoretical quantiles Figure: Gumbel marginal QQ-plots for monthly maxima of daily maximum (left) and daily minimum (right) temperatures. Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
33 References [Meehl et al., 2009] performed a study on the same dataset where the evolution through time of the ratio of annual numbers of record high maximum temperature and record low minimum temperature was of interest. Up to 1970, almost no warming. Since 1970, less extreme low temperatures than expected (warmer nights). Linda Mhalla (UNIGE) Non-stationary Pickands functions Boulder CO, / 31
Generalized 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 informationNonparametric Estimation of the Dependence Function for a Multivariate Extreme Value Distribution
Nonparametric Estimation of the Dependence Function for a Multivariate Extreme Value Distribution p. /2 Nonparametric Estimation of the Dependence Function for a Multivariate Extreme Value Distribution
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 211 http://stat.epfl.ch Multivariate Extremes 19 Componentwise maxima.................................................. 194 Standardization........................................................
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 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 informationAccounting for extreme-value dependence in multivariate data
Accounting for extreme-value dependence in multivariate data 38th ASTIN Colloquium Manchester, July 15, 2008 Outline 1. Dependence modeling through copulas 2. Rank-based inference 3. Extreme-value dependence
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 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 informationApplied Time Series Topics
Applied Time Series Topics Ivan Medovikov Brock University April 16, 2013 Ivan Medovikov, Brock University Applied Time Series Topics 1/34 Overview 1. Non-stationary data and consequences 2. Trends and
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 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 informationMultivariate Non-Normally Distributed Random Variables
Multivariate Non-Normally Distributed Random Variables An Introduction to the Copula Approach Workgroup seminar on climate dynamics Meteorological Institute at the University of Bonn 18 January 2008, Bonn
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 informationConditional Least Squares and Copulae in Claims Reserving for a Single Line of Business
Conditional Least Squares and Copulae in Claims Reserving for a Single Line of Business Michal Pešta Charles University in Prague Faculty of Mathematics and Physics Ostap Okhrin Dresden University of Technology
More informationRegularization in Cox Frailty Models
Regularization in Cox Frailty Models Andreas Groll 1, Trevor Hastie 2, Gerhard Tutz 3 1 Ludwig-Maximilians-Universität Munich, Department of Mathematics, Theresienstraße 39, 80333 Munich, Germany 2 University
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 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 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 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 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 informationGenerating gridded fields of extreme precipitation for large domains with a Bayesian hierarchical model
Generating gridded fields of extreme precipitation for large domains with a Bayesian hierarchical model Cameron Bracken Department of Civil, Environmental and Architectural Engineering University of Colorado
More informationPENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA
PENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA Kasun Rathnayake ; A/Prof Jun Ma Department of Statistics Faculty of Science and Engineering Macquarie University
More informationDiscussion of the paper Inference for Semiparametric Models: Some Questions and an Answer by Bickel and Kwon
Discussion of the paper Inference for Semiparametric Models: Some Questions and an Answer by Bickel and Kwon Jianqing Fan Department of Statistics Chinese University of Hong Kong AND Department of Statistics
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 informationSome Theories about Backfitting Algorithm for Varying Coefficient Partially Linear Model
Some Theories about Backfitting Algorithm for Varying Coefficient Partially Linear Model 1. Introduction Varying-coefficient partially linear model (Zhang, Lee, and Song, 2002; Xia, Zhang, and Tong, 2004;
More informationPreliminaries The bootstrap Bias reduction Hypothesis tests Regression Confidence intervals Time series Final remark. Bootstrap inference
1 / 171 Bootstrap inference Francisco Cribari-Neto Departamento de Estatística Universidade Federal de Pernambuco Recife / PE, Brazil email: cribari@gmail.com October 2013 2 / 171 Unpaid advertisement
More informationAn Introduction to Nonstationary Time Series Analysis
An Introduction to Analysis Ting Zhang 1 tingz@bu.edu Department of Mathematics and Statistics Boston University August 15, 2016 Boston University/Keio University Workshop 2016 A Presentation Friendly
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 informationIntroduction to Maximum Likelihood Estimation
Introduction to Maximum Likelihood Estimation Eric Zivot July 26, 2012 The Likelihood Function Let 1 be an iid sample with pdf ( ; ) where is a ( 1) vector of parameters that characterize ( ; ) Example:
More informationFinancial Econometrics and Volatility Models Copulas
Financial Econometrics and Volatility Models Copulas Eric Zivot Updated: May 10, 2010 Reading MFTS, chapter 19 FMUND, chapters 6 and 7 Introduction Capturing co-movement between financial asset returns
More informationFORECAST VERIFICATION OF EXTREMES: USE OF EXTREME VALUE THEORY
1 FORECAST VERIFICATION OF EXTREMES: USE OF EXTREME VALUE THEORY Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA Email: rwk@ucar.edu Web
More informationStatistical Models for Monitoring and Regulating Ground-level Ozone. Abstract
Statistical Models for Monitoring and Regulating Ground-level Ozone Eric Gilleland 1 and Douglas Nychka 2 Abstract The application of statistical techniques to environmental problems often involves a tradeoff
More informationWeather generators for studying climate change
Weather generators for studying climate change Assessing climate impacts Generating Weather (WGEN) Conditional models for precip Douglas Nychka, Sarah Streett Geophysical Statistics Project, National Center
More informationOPTIMISATION CHALLENGES IN MODERN STATISTICS. Co-authors: Y. Chen, M. Cule, R. Gramacy, M. Yuan
OPTIMISATION CHALLENGES IN MODERN STATISTICS Co-authors: Y. Chen, M. Cule, R. Gramacy, M. Yuan How do optimisation problems arise in Statistics? Let X 1,...,X n be independent and identically distributed
More informationGENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS
GENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA Joint work with Eva
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 29 http://stat.epfl.ch Multivariate Extremes 184 Europe............................................................. 185 Componentwise maxima..................................................
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 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 informationWhen is a copula constant? A test for changing relationships
When is a copula constant? A test for changing relationships Fabio Busetti and Andrew Harvey Bank of Italy and University of Cambridge November 2007 usetti and Harvey (Bank of Italy and University of Cambridge)
More informationPreliminaries The bootstrap Bias reduction Hypothesis tests Regression Confidence intervals Time series Final remark. Bootstrap inference
1 / 172 Bootstrap inference Francisco Cribari-Neto Departamento de Estatística Universidade Federal de Pernambuco Recife / PE, Brazil email: cribari@gmail.com October 2014 2 / 172 Unpaid advertisement
More informationInvestigation of an Automated Approach to Threshold Selection for Generalized Pareto
Investigation of an Automated Approach to Threshold Selection for Generalized Pareto Kate R. Saunders Supervisors: Peter Taylor & David Karoly University of Melbourne April 8, 2015 Outline 1 Extreme Value
More informationRelative increase of record high maximum temperatures compared to record low minimum temperatures in the U.S.
GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L23701, doi:10.1029/2009gl040736, 2009 Relative increase of record high maximum temperatures compared to record low minimum temperatures in the U.S. Gerald A. Meehl,
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 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 informationEXTREMAL MODELS AND ENVIRONMENTAL APPLICATIONS. Rick Katz
1 EXTREMAL MODELS AND ENVIRONMENTAL APPLICATIONS Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home page: www.isse.ucar.edu/hp_rick/
More informationLecture 4 September 15
IFT 6269: Probabilistic Graphical Models Fall 2017 Lecture 4 September 15 Lecturer: Simon Lacoste-Julien Scribe: Philippe Brouillard & Tristan Deleu 4.1 Maximum Likelihood principle Given a parametric
More informationA Brief Introduction to Copulas
A Brief Introduction to Copulas Speaker: Hua, Lei February 24, 2009 Department of Statistics University of British Columbia Outline Introduction Definition Properties Archimedean Copulas Constructing Copulas
More informationNonparametric inference in hidden Markov and related models
Nonparametric inference in hidden Markov and related models Roland Langrock, Bielefeld University Roland Langrock Bielefeld University 1 / 47 Introduction and motivation Roland Langrock Bielefeld University
More informationSpatial Statistics 2013, S2.2 6 th June Institute for Geoinformatics University of Münster.
Spatial Statistics 2013, S2.2 6 th June 2013 Institute for Geoinformatics University of Münster http://ifgi.uni-muenster.de/graeler Vine Vine 1 Spatial/spatio-temporal data Typically, spatial/spatio-temporal
More informationHigh Dimensional Empirical Likelihood for Generalized Estimating Equations with Dependent Data
High Dimensional Empirical Likelihood for Generalized Estimating Equations with Dependent Data Song Xi CHEN Guanghua School of Management and Center for Statistical Science, Peking University Department
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 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 informationFlexible Spatio-temporal smoothing with array methods
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session IPS046) p.849 Flexible Spatio-temporal smoothing with array methods Dae-Jin Lee CSIRO, Mathematics, Informatics and
More informationPhysically-Based Statistical Models of Extremes arising from Extratropical Cyclones
Lancaster University STOR603: PhD Proposal Physically-Based Statistical Models of Extremes arising from Extratropical Cyclones Author: Paul Sharkey Supervisors: Jonathan Tawn Jenny Wadsworth Simon Brown
More informationGENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS
GENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA Joint work with Eva
More informationIssues in estimating past climate change at local scale.
14th EMS Annual Meeting, Prague Oct. 6th 2014 Issues in estimating past climate change at local scale. Case study: the recent warming (1959-2009) over France. Lola Corre (Division de la Climatologie, Météo-France)
More information19. RECORD-BREAKING HEAT IN NORTHWEST CHINA IN JULY 2015: ANALYSIS OF THE SEVERITY AND UNDERLYING CAUSES
19. RECORD-BREAKING HEAT IN NORTHWEST CHINA IN JULY 2015: ANALYSIS OF THE SEVERITY AND UNDERLYING CAUSES Chiyuan Miao, Qiaohong Sun, Dongxian Kong, and Qingyun Duan The record-breaking heat over northwest
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 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 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 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 informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 3 Linear
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 informationProbability Distributions and Estimation of Ali-Mikhail-Haq Copula
Applied Mathematical Sciences, Vol. 4, 2010, no. 14, 657-666 Probability Distributions and Estimation of Ali-Mikhail-Haq Copula Pranesh Kumar Mathematics Department University of Northern British Columbia
More informationEstimation of cumulative distribution function with spline functions
INTERNATIONAL JOURNAL OF ECONOMICS AND STATISTICS Volume 5, 017 Estimation of cumulative distribution function with functions Akhlitdin Nizamitdinov, Aladdin Shamilov Abstract The estimation of the cumulative
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 informationUNIVERSITÄT POTSDAM Institut für Mathematik
UNIVERSITÄT POTSDAM Institut für Mathematik Testing the Acceleration Function in Life Time Models Hannelore Liero Matthias Liero Mathematische Statistik und Wahrscheinlichkeitstheorie Universität Potsdam
More informationLeast Absolute Shrinkage is Equivalent to Quadratic Penalization
Least Absolute Shrinkage is Equivalent to Quadratic Penalization Yves Grandvalet Heudiasyc, UMR CNRS 6599, Université de Technologie de Compiègne, BP 20.529, 60205 Compiègne Cedex, France Yves.Grandvalet@hds.utc.fr
More informationBivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories
Entropy 2012, 14, 1784-1812; doi:10.3390/e14091784 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Bivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories Lan Zhang
More informationDS-GA 1002 Lecture notes 12 Fall Linear regression
DS-GA Lecture notes 1 Fall 16 1 Linear models Linear regression In statistics, regression consists of learning a function relating a certain quantity of interest y, the response or dependent variable,
More informationOverview of Extreme Value Analysis (EVA)
Overview of Extreme Value Analysis (EVA) Brian Reich North Carolina State University July 26, 2016 Rossbypalooza Chicago, IL Brian Reich Overview of Extreme Value Analysis (EVA) 1 / 24 Importance of extremes
More informationClimate Change: the Uncertainty of Certainty
Climate Change: the Uncertainty of Certainty Reinhard Furrer, UZH JSS, Geneva Oct. 30, 2009 Collaboration with: Stephan Sain - NCAR Reto Knutti - ETHZ Claudia Tebaldi - Climate Central Ryan Ford, Doug
More informationApplications of Tail Dependence II: Investigating the Pineapple Express. Dan Cooley Grant Weller Department of Statistics Colorado State University
Applications of Tail Dependence II: Investigating the Pineapple Express Dan Cooley Grant Weller Department of Statistics Colorado State University Joint work with: Steve Sain, Melissa Bukovsky, Linda Mearns,
More informationA comparison of U.S. precipitation extremes under two climate change scenarios
A comparison of U.S. precipitation extremes under two climate change scenarios Miranda Fix 1 Dan Cooley 1 Steve Sain 2 Claudia Tebaldi 3 1 Department of Statistics, Colorado State University 2 The Climate
More informationInversion Base Height. Daggot Pressure Gradient Visibility (miles)
Stanford University June 2, 1998 Bayesian Backtting: 1 Bayesian Backtting Trevor Hastie Stanford University Rob Tibshirani University of Toronto Email: trevor@stat.stanford.edu Ftp: stat.stanford.edu:
More informationQuantifying Stochastic Model Errors via Robust Optimization
Quantifying Stochastic Model Errors via Robust Optimization IPAM Workshop on Uncertainty Quantification for Multiscale Stochastic Systems and Applications Jan 19, 2016 Henry Lam Industrial & Operations
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 informationBayes methods for categorical data. April 25, 2017
Bayes methods for categorical data April 25, 2017 Motivation for joint probability models Increasing interest in high-dimensional data in broad applications Focus may be on prediction, variable selection,
More informationCross-covariance Functions for Tangent Vector Fields on the Sphere
Cross-covariance Functions for Tangent Vector Fields on the Sphere Minjie Fan 1 Tomoko Matsuo 2 1 Department of Statistics University of California, Davis 2 Cooperative Institute for Research in Environmental
More informationTail Approximation of Value-at-Risk under Multivariate Regular Variation
Tail Approximation of Value-at-Risk under Multivariate Regular Variation Yannan Sun Haijun Li July 00 Abstract This paper presents a general tail approximation method for evaluating the Valueat-Risk of
More informationA Goodness-of-fit Test for Copulas
A Goodness-of-fit Test for Copulas Artem Prokhorov August 2008 Abstract A new goodness-of-fit test for copulas is proposed. It is based on restrictions on certain elements of the information matrix and
More informationMultivariate verification: Motivation, Complexity, Examples
Multivariate verification: Motivation, Complexity, Examples A.Hense, A. Röpnack, J. Keune, R. Glowienka-Hense, S. Stolzenberger, H. Weinert Berlin, May, 5th 2017 Motivations for MV verification Data assimilation
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 informationLecture Quantitative Finance Spring Term 2015
on bivariate Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 07: April 2, 2015 1 / 54 Outline on bivariate 1 2 bivariate 3 Distribution 4 5 6 7 8 Comments and conclusions
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 24 Paper 153 A Note on Empirical Likelihood Inference of Residual Life Regression Ying Qing Chen Yichuan
More informationChanges in frequency of heat and cold in the United States temperature record
Changes in frequency of heat and cold in the United States temperature record Summary: Have heat records increased and record low temperatures decreased as suggested by the IPCC greenhouse driven models?
More informationCalibration of the Empirical Likelihood Method for a Vector Mean
Calibration of the Empirical Likelihood Method for a Vector Mean Sarah C. Emerson and Art B. Owen Department of Statistics Stanford University Abstract The empirical likelihood method is a versatile approach
More informationof the 7 stations. In case the number of daily ozone maxima in a month is less than 15, the corresponding monthly mean was not computed, being treated
Spatial Trends and Spatial Extremes in South Korean Ozone Seokhoon Yun University of Suwon, Department of Applied Statistics Suwon, Kyonggi-do 445-74 South Korea syun@mail.suwon.ac.kr Richard L. Smith
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 informationReduced Overdispersion in Stochastic Weather Generators for Statistical Downscaling of Seasonal Forecasts and Climate Change Scenarios
Reduced Overdispersion in Stochastic Weather Generators for Statistical Downscaling of Seasonal Forecasts and Climate Change Scenarios Yongku Kim Institute for Mathematics Applied to Geosciences National
More informationOn block bootstrapping areal data Introduction
On block bootstrapping areal data Nicholas Nagle Department of Geography University of Colorado UCB 260 Boulder, CO 80309-0260 Telephone: 303-492-4794 Email: nicholas.nagle@colorado.edu Introduction Inference
More informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
More informationRISK AND EXTREMES: ASSESSING THE PROBABILITIES OF VERY RARE EVENTS
RISK AND EXTREMES: ASSESSING THE PROBABILITIES OF VERY RARE EVENTS Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu
More informationLecture 14: Variable Selection - Beyond LASSO
Fall, 2017 Extension of LASSO To achieve oracle properties, L q penalty with 0 < q < 1, SCAD penalty (Fan and Li 2001; Zhang et al. 2007). Adaptive LASSO (Zou 2006; Zhang and Lu 2007; Wang et al. 2007)
More informationTesting Restrictions and Comparing Models
Econ. 513, Time Series Econometrics Fall 00 Chris Sims Testing Restrictions and Comparing Models 1. THE PROBLEM We consider here the problem of comparing two parametric models for the data X, defined by
More informationMax-stable processes: Theory and Inference
1/39 Max-stable processes: Theory and Inference Mathieu Ribatet Institute of Mathematics, EPFL Laboratory of Environmental Fluid Mechanics, EPFL joint work with Anthony Davison and Simone Padoan 2/39 Outline
More informationLINEAR CLASSIFICATION, PERCEPTRON, LOGISTIC REGRESSION, SVC, NAÏVE BAYES. Supervised Learning
LINEAR CLASSIFICATION, PERCEPTRON, LOGISTIC REGRESSION, SVC, NAÏVE BAYES Supervised Learning Linear vs non linear classifiers In K-NN we saw an example of a non-linear classifier: the decision boundary
More informationRegression, Ridge Regression, Lasso
Regression, Ridge Regression, Lasso Fabio G. Cozman - fgcozman@usp.br October 2, 2018 A general definition Regression studies the relationship between a response variable Y and covariates X 1,..., X n.
More informationModelling spatially-dependent non-stationary extremes with application to hurricane-induced wave heights
Modelling spatially-dependent non-stationary extremes with application to hurricane-induced wave heights Paul Northrop and Philip Jonathan March 5, 2010 Abstract In environmental applications it is found
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 information