Investigation of an Automated Approach to Threshold Selection for Generalized Pareto
|
|
- Thomasina Tiffany Hill
- 6 years ago
- Views:
Transcription
1 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
2 Outline 1 Extreme Value Theory 2
3 Problem What are the climate processes that drive extreme rainfall? (El Niño Southern Oscillation, Interdecadal Pacific Oscillation) How do these drivers differ at different timescales; sub-daily, daily, consecutive day totals?
4 Data Extreme Value Theory
5 Extreme Value Theory
6 Block Maxima Extreme Value Theory
7 Block Maxima Let X 1, X 2,..., X n be a sequence of i.i.d. random variables with distribution function F. Define M n = max{x 1, X 2,..., X n }. (X i might be daily rainfall observations and M 365 the annual maximum rainfall.) Pr(M n x) = Pr(X 1 x,..., X n x) = Pr(X 1 x) Pr(X n x) = F (x) n. As n, the distribution of the M n converges to a generalised extreme value distribution.
8 Generalized Extreme Value Theorem (Fisher-Tippett-Gnendenko) If there exists sequences of constants {a n > 0} and {b n } such that ( ) Mn b n Pr z G(z) as n a n for a non-degenerate distribution function G, then G is a member of the Generalized Extreme Value family [ ( )] 1 z µ G(z) = exp 1 + ξ ξ σ defined on {z : 1 + ξ(z µ)/σ > 0}, where < µ <, σ > 0 and < ξ <.
9 Leveraging more data Extreme Value Theory
10 Generalized Pareto Distribution Let X 1, X 2,..., X n be a sequence of iid random variables with marginal distribution function F. Pr{X > u + y X > u} = 1 F (u + y) 1 F (u) y > 0. If F satisfies Generalized Extreme Value Theorem then for a large enough threshold u, the distribution function of (X u) conditional on X > u is the GPD. Generalized Pareto Distribution - Picklands (1975) ( H(y) = ξỹ ) 1/ξ σ defined on {y : y > 0} and (1 + ξy/ σ > 0) where, σ = σ + ξ(u µ).
11 Dependence Rainfall observations are dependent Heavy rainfall yesterday effects the probability of heavy rain today Heavy rainfall a year ago doesn t Extreme Value Theory extends to stationary series with weak long range dependence However, for processes with short range dependence extremes occur in clusters
12 Clusters Extreme Value Theory
13 Dependent Series Let {X i } i 1 be a stationary series and {X i } i 1 be an independent series of variables with the same marginal distribution. Define M n = max{x 1,..., X n } and Mn = max{x1,..., X n }. Under suitable regularity conditions, { } (M Pr n b n ) z G(z), a n as n for normalizing sequences {a n > 0} and {b n }, where G is a non-degenerate distribution functions, if and only if { } (Mn b n ) Pr z G θ (z), a n for a constant θ such that 0 < θ 1.
14 Extremal Index θ = {Limiting mean cluster size} 1 (0, 1] θ = observations per cluster on average.
15
16 Select a threshold Decluster the data for independent observations
17 Declustering Blocks Partition the observation sequence into blocks of length, b Assume extreme observations within the same block belong to the same same cluster. Runs Specify a run length, K Assume extreme observations with an inter-exceedance time of less than K belong to the same cluster.
18 Intervals The limiting process of exceedance times is compound Poisson for stationary series (Hsing et al. 1988). Ferro and Segers (2003) showed the limiting distribution of inter-exceedance times is a mixture distribution with weight θ, T θ (t) = (1 θ)ɛ 0 + θ θ exp( θt), where ɛ 0 is a degenerate distribution, T θ is the distribution of arrival times of exceedances at threshold u. By equating moments a non-parametric estimator can be found for θ. The largest θ(n 1) inter-exceedance times can be interpreted as between cluster arrivals.
19 Select a threshold Decluster the data for independent observations
20 Mean Residual Life Plots For sufficiently high thresholds, as the threshold increases the expected exceedance above the threshold should grow linearly.
21 Parameter Stability Plots Parameter estimates of (modified) scale and shape parameters should be constant for the range of valid thresholds.
22 Alternative Set the threshold according to a high quantile of non-zero observations Eg. 90 th percentile. Is this an appropriate threshold? Is our model is misspecified? Suggested approach by Süveges and Davison et al. (2010) is to test the threshold, u, and run parameter, K pair for model misspecification.
23 Log-Likelihood Limiting distribution of inter-exceedance times: T θ (t) = (1 θ)ɛ 0 + θ 2 exp( θt), Log-Likelihood (strictly positive inter-exceedance times): N 1 i=1 log ( (1 θ) I(t i =0) (θ 2 exp(θt i ) I(t i >0) ) = N 1 i=1 [2I(t i > 0) log(θ) θt i ], where t i = NT i n, n is the total number of observations and N is the number of exceedances. However as n gets large our estimate, ˆθ, tends to 1 suggesting independence.
24 Log-Likelihood Adjustment of the inter-exceedance times using the run parameter K: Log-likelihood: l(θ; c i ) = N 1 i=1 c i = max{t i K, 0} [I(c i = 0) log(1 θ) + 2I(c i > 0) log(θ) θc i ] Approach used in Fukutome et al. (2014) and Süveges and Davison (2010). Test combinations of threshold, u, and run parameter, K, for misspecification of the likelihood function. Select the (u, K) pair that maximizes the number of independent clusters.
25 Model Misspecification If a parametric model is misspecified then there is no θ such that g = f (θ), where g is the true model and f is the misspecified parametric model. For a well specified model, the Fisher s information matrix, I (θ) = E{l (θ; c j } is equal to the variance of the score vector, J(θ) = Var{l (θ; c j )}. Test the hypothesis: D(θ) = J(θ) I (θ), where H 0 : D(θ) = 0 and H 1 : D(θ) 0.
26 Empirically: I N 1 (ˆθ) = 1 N 1 l (ˆθ; c j ) (N 1) j=1 J N 1 (ˆθ) = N 1 1 l (ˆθ; c j ) 2 (N 1) j=1 D N 1 (ˆθ) = J N 1 (ˆθ) I N 1 (ˆθ) V N 1 (ˆθ) = 1 (N 1) N 1 j=1 [ (dj (ˆθ; c j ) D N 1 (ˆθ)I N 1 (ˆθ) 1 l (ˆθ; c j ) ) ] 2 where V N 1 (ˆθ) is the sample variance of D N 1 (ˆθ).
27 Model Misspecification Theorem: (Information Matrix Test - Whyte 1982) If the assumed model l(θ; c i ) contains the true model for some θ = θ 0, then as n, (i) (N 1)D N 1 (ˆθ) w N(0, V (θ 0 )), (ii) V N 1 ( θ N 1 ˆ ) a.s. V (θ 0 ), and V N 1 (ˆθ) is non-singular for sufficiently large N, (iii) Then the Information Matrix Test statistic, (N 1)D N 1 (ˆθ) V N 1 (ˆθ) 1 D N 1 (ˆθ) is asymptotically χ 2 1 distributed.
28 Example: AR(2) Y i = 0.95Y i Y i 2 + Z i where Z i GP(1, 1/2) and n = simulations
29 Adjusting inter-exceedance times Common to assume stationarity by enforcing seasonal blocking. Collapse inter-exceedance times across seasonal blocks using the memoryless property of the exponential for fitting.
30 Results: Gatton, South East Queensland
31 Results: Oenpelli, Northern Territory
32 Summary Shown how to check if the threshold and run parameter selected violate the assumptions of the model Given confidence to threshold selection in the absence of a hard and fast rule and in the presence of subjectivity
33 References Ferro, C. and Segers, J. (2003). Inference for clusters of extreme values. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 65(2), pp Fukutome, S., Liniger, M. and Sveges, M. (2014). Automatic threshold and run parameter selection: a climatology for extreme hourly precipitation in Switzerland. Theoretical and Applied Climatology. Hsing, T., Hüsler, J. and Leadbetter, M. (1988). On the exceedance point process for a stationary sequence. Probability Theory and Related Fields, 78(1), pp Süveges, M. and Davison, A. (2010). Model misspecification in peaks over threshold analysis. The Annals of Applied Statistics, 4(1), pp White, H. (1982). Maximum Likelihood Estimation of Misspecified Models. Econometrica, 50(1), p.1.
34 ANZAPW 2015: Barossa Valley, South Australia This work has been supported by the ARC through the Laureate Fellowship FL Questions?
35 Results: Kalamia, Far North Queensland
36 Results: Yamba, New South Wales
Bayesian 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 informationEVA Tutorial #2 PEAKS OVER THRESHOLD APPROACH. Rick Katz
1 EVA Tutorial #2 PEAKS OVER THRESHOLD APPROACH Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home page: www.isse.ucar.edu/staff/katz/
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 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 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 informationSTATISTICAL METHODS FOR RELATING TEMPERATURE EXTREMES TO LARGE-SCALE METEOROLOGICAL PATTERNS. Rick Katz
1 STATISTICAL METHODS FOR RELATING TEMPERATURE EXTREMES TO LARGE-SCALE METEOROLOGICAL PATTERNS Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder,
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 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 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 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 informationHIERARCHICAL MODELS IN EXTREME VALUE THEORY
HIERARCHICAL MODELS IN EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research, University of North Carolina, Chapel Hill and Statistical and Applied Mathematical Sciences
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 informationWhat Can We Infer From Beyond The Data? The Statistics Behind The Analysis Of Risk Events In The Context Of Environmental Studies
What Can We Infer From Beyond The Data? The Statistics Behind The Analysis Of Risk Events In The Context Of Environmental Studies Sibusisiwe Khuluse, Sonali Das, Pravesh Debba, Chris Elphinstone Logistics
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 informationExtreme value modelling of rainfalls and
Universite de Paris Sud Master s Thesis Extreme value modelling of rainfalls and flows Author: Fan JIA Supervisor: Pr Elisabeth Gassiat Dr Elena Di Bernardino Pr Michel Bera A thesis submitted in fulfilment
More informationEmma Simpson. 6 September 2013
6 September 2013 Test What is? Beijing during periods of low and high air pollution Air pollution is composed of sulphur oxides, nitrogen oxides, carbon monoxide and particulates. Particulates are small
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 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 informationPeaks-Over-Threshold Modelling of Environmental Data
U.U.D.M. Project Report 2014:33 Peaks-Over-Threshold Modelling of Environmental Data Esther Bommier Examensarbete i matematik, 30 hp Handledare och examinator: Jesper Rydén September 2014 Department of
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 informationInference for clusters of extreme values
J. R. Statist. Soc. B (2003) 65, Part 2, pp. 545 556 Inference for clusters of extreme values Christopher A. T. Ferro University of Reading, UK and Johan Segers EURANDOM, Eindhoven, the Netherlands [Received
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 informationExtreme Event Modelling
Extreme Event Modelling Liwei Wu, SID: 52208712 Department of Mathematics City University of Hong Kong Supervisor: Dr. Xiang Zhou March 31, 2014 Contents 1 Introduction 4 2 Theory and Methods 5 2.1 Asymptotic
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 informationMFM Practitioner Module: Quantitiative Risk Management. John Dodson. October 14, 2015
MFM Practitioner Module: Quantitiative Risk Management October 14, 2015 The n-block maxima 1 is a random variable defined as M n max (X 1,..., X n ) for i.i.d. random variables X i with distribution function
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 informationStatistical Methods for Clusters of Extreme Values
Statistical Methods for Clusters of Extreme Values Christopher A. T. Ferro, B.Sc. M.Sc. Submitted for the degree of Doctor of Philosophy at Lancaster University, September 2003. I declare that the work
More informationStatistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach
Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Jae-Kwang Kim Department of Statistics, Iowa State University Outline 1 Introduction 2 Observed likelihood 3 Mean Score
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 informationOutline of GLMs. Definitions
Outline of GLMs Definitions This is a short outline of GLM details, adapted from the book Nonparametric Regression and Generalized Linear Models, by Green and Silverman. The responses Y i have density
More informationHigh-frequency data modelling using Hawkes processes
High-frequency data modelling using Hawkes processes Valérie Chavez-Demoulin 1 joint work J.A McGill 1 Faculty of Business and Economics, University of Lausanne, Switzerland Boulder, April 2016 Boulder,
More informationAssessing Dependence in Extreme Values
02/09/2016 1 Motivation Motivation 2 Comparison 3 Asymptotic Independence Component-wise Maxima Measures Estimation Limitations 4 Idea Results Motivation Given historical flood levels, how high should
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 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 informationHigh-frequency data modelling using Hawkes processes
Valérie Chavez-Demoulin joint work with High-frequency A.C. Davison data modelling and using A.J. Hawkes McNeil processes(2005), J.A EVT2013 McGill 1 /(201 High-frequency data modelling using Hawkes processes
More informationRegional Estimation from Spatially Dependent Data
Regional Estimation from Spatially Dependent Data R.L. Smith Department of Statistics University of North Carolina Chapel Hill, NC 27599-3260, USA December 4 1990 Summary Regional estimation methods are
More informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationRecall that in order to prove Theorem 8.8, we argued that under certain regularity conditions, the following facts are true under H 0 : 1 n
Chapter 9 Hypothesis Testing 9.1 Wald, Rao, and Likelihood Ratio Tests Suppose we wish to test H 0 : θ = θ 0 against H 1 : θ θ 0. The likelihood-based results of Chapter 8 give rise to several possible
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 information2.6.3 Generalized likelihood ratio tests
26 HYPOTHESIS TESTING 113 263 Generalized likelihood ratio tests When a UMP test does not exist, we usually use a generalized likelihood ratio test to verify H 0 : θ Θ against H 1 : θ Θ\Θ It can be used
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 informationExtreme Value Theory and Applications
Extreme Value Theory and Deauville - 04/10/2013 Extreme Value Theory and Introduction Asymptotic behavior of the Sum Extreme (from Latin exter, exterus, being on the outside) : Exceeding the ordinary,
More informationGraduate Econometrics I: Maximum Likelihood II
Graduate Econometrics I: Maximum Likelihood II Yves Dominicy Université libre de Bruxelles Solvay Brussels School of Economics and Management ECARES Yves Dominicy Graduate Econometrics I: Maximum Likelihood
More informationTheory of Maximum Likelihood Estimation. Konstantin Kashin
Gov 2001 Section 5: Theory of Maximum Likelihood Estimation Konstantin Kashin February 28, 2013 Outline Introduction Likelihood Examples of MLE Variance of MLE Asymptotic Properties What is Statistical
More informationUNIVERSITY OF CALGARY. Inference for Dependent Generalized Extreme Values. Jialin He A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES
UNIVERSITY OF CALGARY Inference for Dependent Generalized Extreme Values by Jialin He A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF
More informationPENULTIMATE APPROXIMATIONS FOR WEATHER AND CLIMATE EXTREMES. Rick Katz
PENULTIMATE APPROXIMATIONS FOR WEATHER AND CLIMATE EXTREMES Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA Email: rwk@ucar.edu Web site:
More informationFuture extreme precipitation events in the Southwestern US: climate change and natural modes of variability
Future extreme precipitation events in the Southwestern US: climate change and natural modes of variability Francina Dominguez Erick Rivera Fernandez Hsin-I Chang Christopher Castro AGU 2010 Fall Meeting
More informationEmpirical Likelihood
Empirical Likelihood Patrick Breheny September 20 Patrick Breheny STA 621: Nonparametric Statistics 1/15 Introduction Empirical likelihood We will discuss one final approach to constructing confidence
More informationBayesian Model Diagnostics and Checking
Earvin Balderama Quantitative Ecology Lab Department of Forestry and Environmental Resources North Carolina State University April 12, 2013 1 / 34 Introduction MCMCMC 2 / 34 Introduction MCMCMC Steps in
More informationRISK ANALYSIS AND EXTREMES
RISK ANALYSIS AND EXTREMES Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu Opening Workshop SAMSI program on
More informationJournal of Environmental Statistics
jes Journal of Environmental Statistics February 2010, Volume 1, Issue 3. http://www.jenvstat.org Exponentiated Gumbel Distribution for Estimation of Return Levels of Significant Wave Height Klara Persson
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 informationMISCELLANEOUS TOPICS RELATED TO LIKELIHOOD. Copyright c 2012 (Iowa State University) Statistics / 30
MISCELLANEOUS TOPICS RELATED TO LIKELIHOOD Copyright c 2012 (Iowa State University) Statistics 511 1 / 30 INFORMATION CRITERIA Akaike s Information criterion is given by AIC = 2l(ˆθ) + 2k, where l(ˆθ)
More informationClimate predictability beyond traditional climate models
Climate predictability beyond traditional climate models Rasmus E. Benestad & Abdelkader Mezghani Rasmus.benestad@met.no More heavy rain events? More heavy rain events? Heavy precipitation events with
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 informationThe Spatial Variation of the Maximum Possible Pollutant Concentration from Steady Sources
International Environmental Modelling and Software Society (iemss) 2010 International Congress on Environmental Modelling and Software Modelling for Environment s Sake, Fifth Biennial Meeting, Ottawa,
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 informationZambia. General Climate. Recent Climate Trends. UNDP Climate Change Country Profiles. Temperature. C. McSweeney 1, M. New 1,2 and G.
UNDP Climate Change Country Profiles Zambia C. McSweeney 1, M. New 1,2 and G. Lizcano 1 1. School of Geography and Environment, University of Oxford. 2. Tyndall Centre for Climate Change Research http://country-profiles.geog.ox.ac.uk
More informationRobust Backtesting Tests for Value-at-Risk Models
Robust Backtesting Tests for Value-at-Risk Models Jose Olmo City University London (joint work with Juan Carlos Escanciano, Indiana University) Far East and South Asia Meeting of the Econometric Society
More informationMath 576: Quantitative Risk Management
Math 576: Quantitative Risk Management Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 11 Haijun Li Math 576: Quantitative Risk Management Week 11 1 / 21 Outline 1
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 informationR.Garçon, F.Garavaglia, J.Gailhard, E.Paquet, F.Gottardi EDF-DTG
Homogeneous samples and reliability of probabilistic models : using an atmospheric circulation patterns sampling for a better estimation of extreme rainfall probability R.Garçon, F.Garavaglia, J.Gailhard,
More informationGraduate Econometrics I: Maximum Likelihood I
Graduate Econometrics I: Maximum Likelihood I Yves Dominicy Université libre de Bruxelles Solvay Brussels School of Economics and Management ECARES Yves Dominicy Graduate Econometrics I: Maximum Likelihood
More informationR&D Research Project: Scaling analysis of hydrometeorological time series data
R&D Research Project: Scaling analysis of hydrometeorological time series data Extreme Value Analysis considering Trends: Methodology and Application to Runoff Data of the River Danube Catchment M. Kallache,
More informationModelling Risk on Losses due to Water Spillage for Hydro Power Generation. A Verster and DJ de Waal
Modelling Risk on Losses due to Water Spillage for Hydro Power Generation. A Verster and DJ de Waal Department Mathematical Statistics and Actuarial Science University of the Free State Bloemfontein ABSTRACT
More informationReview of existing statistical methods for flood frequency estimation in Greece
EU COST Action ES0901: European Procedures for Flood Frequency Estimation (FloodFreq) 3 rd Management Committee Meeting, Prague, 28 29 October 2010 WG2: Assessment of statistical methods for flood frequency
More informationNormalising constants and maximum likelihood inference
Normalising constants and maximum likelihood inference Jakob G. Rasmussen Department of Mathematics Aalborg University Denmark March 9, 2011 1/14 Today Normalising constants Approximation of normalising
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationAccurate directional inference for vector parameters
Accurate directional inference for vector parameters Nancy Reid October 28, 2016 with Don Fraser, Nicola Sartori, Anthony Davison Parametric models and likelihood model f (y; θ), θ R p data y = (y 1,...,
More informationON THE TWO STEP THRESHOLD SELECTION FOR OVER-THRESHOLD MODELLING
ON THE TWO STEP THRESHOLD SELECTION FOR OVER-THRESHOLD MODELLING Pietro Bernardara (1,2), Franck Mazas (3), Jérôme Weiss (1,2), Marc Andreewsky (1), Xavier Kergadallan (4), Michel Benoît (1,2), Luc Hamm
More informationParameter Estimation
Parameter Estimation Chapters 13-15 Stat 477 - Loss Models Chapters 13-15 (Stat 477) Parameter Estimation Brian Hartman - BYU 1 / 23 Methods for parameter estimation Methods for parameter estimation Methods
More informationChange Point Analysis of Extreme Values
Change Point Analysis of Extreme Values TIES 2008 p. 1/? Change Point Analysis of Extreme Values Goedele Dierckx Economische Hogeschool Sint Aloysius, Brussels, Belgium e-mail: goedele.dierckx@hubrussel.be
More informationQuasi-likelihood Scan Statistics for Detection of
for Quasi-likelihood for Division of Biostatistics and Bioinformatics, National Health Research Institutes & Department of Mathematics, National Chung Cheng University 17 December 2011 1 / 25 Outline for
More informationAccounting for Choice of Measurement Scale in Extreme Value Modeling
Accounting for Choice of Measurement Scale in Extreme Value Modeling J. L. Wadsworth Department of Mathematics and Statistics Lancaster University, LA1 4YF, UK J. A. Tawn Department of Mathematics and
More information1 Degree distributions and data
1 Degree distributions and data A great deal of effort is often spent trying to identify what functional form best describes the degree distribution of a network, particularly the upper tail of that distribution.
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 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 informationAccurate directional inference for vector parameters
Accurate directional inference for vector parameters Nancy Reid February 26, 2016 with Don Fraser, Nicola Sartori, Anthony Davison Nancy Reid Accurate directional inference for vector parameters York 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 informationIT S TIME FOR AN UPDATE EXTREME WAVES AND DIRECTIONAL DISTRIBUTIONS ALONG THE NEW SOUTH WALES COASTLINE
IT S TIME FOR AN UPDATE EXTREME WAVES AND DIRECTIONAL DISTRIBUTIONS ALONG THE NEW SOUTH WALES COASTLINE M Glatz 1, M Fitzhenry 2, M Kulmar 1 1 Manly Hydraulics Laboratory, Department of Finance, Services
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 informationModelação de valores extremos e sua importância na
Modelação de valores extremos e sua importância na segurança e saúde Margarida Brito Departamento de Matemática FCUP (FCUP) Valores Extremos - DemSSO 1 / 12 Motivation Consider the following events Occurance
More information4.5.1 The use of 2 log Λ when θ is scalar
4.5. ASYMPTOTIC FORM OF THE G.L.R.T. 97 4.5.1 The use of 2 log Λ when θ is scalar Suppose we wish to test the hypothesis NH : θ = θ where θ is a given value against the alternative AH : θ θ on the basis
More informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 00 MODULE : Statistical Inference Time Allowed: Three Hours Candidates should answer FIVE questions. All questions carry equal marks. The
More informationDA Freedman Notes on the MLE Fall 2003
DA Freedman Notes on the MLE Fall 2003 The object here is to provide a sketch of the theory of the MLE. Rigorous presentations can be found in the references cited below. Calculus. Let f be a smooth, scalar
More informationStatistics. Lecture 2 August 7, 2000 Frank Porter Caltech. The Fundamentals; Point Estimation. Maximum Likelihood, Least Squares and All That
Statistics Lecture 2 August 7, 2000 Frank Porter Caltech The plan for these lectures: The Fundamentals; Point Estimation Maximum Likelihood, Least Squares and All That What is a Confidence Interval? Interval
More informationAutomated, Efficient, and Practical Extreme Value Analysis with Environmental Applications
Automated, Efficient, and Practical Extreme Value Analysis with Environmental Applications arxiv:1611.08261v1 [stat.me] 24 Nov 2016 Brian M. Bader, Ph.D. University of Connecticut, 2016 ABSTRACT Although
More informationNon stationary extremes
5 Non stationary extremes 5.1 Introduction In the context of environmental processes, it is common to observe non stationarity for example, due to different seasons having different climate patterns, or
More informationSTAT 135 Lab 5 Bootstrapping and Hypothesis Testing
STAT 135 Lab 5 Bootstrapping and Hypothesis Testing Rebecca Barter March 2, 2015 The Bootstrap Bootstrap Suppose that we are interested in estimating a parameter θ from some population with members x 1,...,
More information5.2 Annual maximum sea levels in Venice
Chapter 5 Non stationary extremes 5.1 Introduction In the context of environmental processes, it is common to observe non stationarity for example, due to different seasons having different climate patterns,
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 informationEstimating AR/MA models
September 17, 2009 Goals The likelihood estimation of AR/MA models AR(1) MA(1) Inference Model specification for a given dataset Why MLE? Traditional linear statistics is one methodology of estimating
More informationExtending clustered point process-based rainfall models to a non-stationary climate
Extending clustered point process-based rainfall models to a non-stationary climate Jo Kaczmarska 1, 2 Valerie Isham 2 Paul Northrop 2 1 Risk Management Solutions 2 Department of Statistical Science, University
More informationNew Bayesian methods for model comparison
Back to the future New Bayesian methods for model comparison Murray Aitkin murray.aitkin@unimelb.edu.au Department of Mathematics and Statistics The University of Melbourne Australia Bayesian Model Comparison
More informationAnalysing River Discharge Data for Flood Analysis
Analysing River Discharge Data for Flood Analysis This suggested response provides answers to the eleven questions posed on the web-site. Your report should have covered all of this material. (A) Produce
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 information1. Fisher Information
1. Fisher Information Let f(x θ) be a density function with the property that log f(x θ) is differentiable in θ throughout the open p-dimensional parameter set Θ R p ; then the score statistic (or score
More informationNew design rainfalls. Janice Green, Project Director IFD Revision Project, Bureau of Meteorology
New design rainfalls Janice Green, Project Director IFD Revision Project, Bureau of Meteorology Design Rainfalls Design Rainfalls Severe weather thresholds Flood forecasting assessing probability of rainfalls
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 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 information