An Introduction to Wavelets with Applications in Environmental Science
|
|
- Warren Bates
- 5 years ago
- Views:
Transcription
1 An Introduction to Wavelets with Applications in Environmental Science Don Percival Applied Physics Lab, University of Washington Data Analysis Products Division, MathSoft overheads for talk available at joint work with: B. Whitcher, S. Byers & P. Guttorp (UW) S. Sardy (UW, MathSoft, EPFL) P. Craigmile (UW)
2 Overview of Talk overview of discrete wavelet transform (DWT) wavelet coefficients and their interpretation DWT as a time series decorrelator three uses for wavelets 1. testing for variance changes 2. bootstrapping auto/cross-correlation estimates 3. estimating d for stationary/nonstationary fractional difference processes with trend warning: results for 2& 3 are preliminary
3 Overview of DWT let X =[X 0,X 1,...,X N 1 ] T be observed time series (for convenience, assume N integer multiple of 2 J 0) let W be N N orthonormal DWT matrix (more precisely: partial DWT of level J 0 ) W = WX is vector of DWT coefficients can partition W as follows: W = W 1. W J0 V J0 W j contains N j = N/2 j wavelet coefficients related to changes of averages at scale τ j =2 j 1 (τ j is jth dyadic scale) related to times spaced 2 j units apart V J0 contains N J0 = N/2 J 0 scaling coefficients related to averages at scale λ J0 =2 J 0 related to times spaced 2 J 0 units apart
4 Example: DWT of FDP X t called fractional difference process (FDP) if it has a spectral density function (SDF) given by S X (f) = where σ 2 > 0 and 1 2 d<1 2 σ 2 2 sin(πf) 2d, note: for small f, have S X (f) C/ f 2d ; i.e., 1/f type process if d = 0, FDP is white noise if 0 <d< 1 2, FDP stationary with long memory can extend definition to d 1 2 nonstationary 1/f type process also called FARIMA(0,d,0) process example: DWT of FDP, d =0.4
5 Two Consequences of Orthonormality multiresolution analysis (MRA) X = W T W = J 0 j=1 W T j W j + V T J 0 V J0 J 0 j=1 (W j partitions W commensurate with W j ) scale-based additive decomposition D j s & S J0 called details & smooth example: Nile River minimum flood levels analysis of variance: ˆσ 2 X 1 N N 1 t=0 (X t X) 2 = 1 N J 0 j=1 D j + S J0 W j 2 + V J0 2 X 2 scale-based decomposition (cf. frequency-based) can define wavelet variance νx(τ 2 j ) for FDP, can deduce d from log/log plots since ν 2 X(τ j ) Cτ 2d 1 j example: Nile River minimum flood levels
6 DWT in Terms of Filters filter X 0,X 1,...,X N 1 to obtain 2 j/2 W j,t L j 1 l=0 h j,l X t l mod N, t =0, 1,...,N 1 where h j,l is jth level wavelet filter note: circular filtering subsample to obtain wavelet coefficients: W j,t =2 j/2 W j,2 j (t+1) 1, t =0, 1,...,N j 1, where W j,t is tth element of W j examples: Haar, D(4), C(6) & LA(8) wavelet filters 1 jth wavelet filter is band-pass with pass-band [, 1 ] 2 j+1 2 j note: jth scale related to interval of frequencies similarly, scaling filters yield V J0 examples: Haar, D(4), C(6) & LA(8) scaling filters J 0 th scaling filter is low-pass with pass-band [0, 1 2 J 0 +1 ]
7 Wavelets as Whitening Filters recall DWT of FDP for d =0.4 since FDP is stationary process, W j is also (ignoring terms influenced by circularity) can compute SDFs for each W j see figure DWT acts as whitening filter requires SDF of X to be flat over pass-band 1 [, 1 ] 2 j+1 2 j if not true, can use wavelet packet transform (DWPT) used by Flandrin, Tewfik & Kim, Wornell, McCoy & Walden three examples built on whitening property 1. testing for variance changes 2. bootstrapping auto/cross-correlation estimates 3. estimating d for stationary/nonstationary fractional difference processes with trend whitening property should help with other problems
8 Homogeneity of Variance: I claim: DWT approximately decorrelates FDPs implication: W j should resemble white noise (ignoring coefficients influenced by circularity) cov {W j,t,w j,t } 0 when t t var {W j,t } should not vary with t (homogeneity of variance) can test for homogeneity of variance using W j suppose Y 0,...,Y N 1 independent normal RVs with E{Y t } = 0 and var {Y t } = σ 2 t want to test null hypothesis H 0 : σ 2 0 = σ 2 1 = = σ 2 N 1 can test H 0 versus a variety of alternatives, e.g., H 1 : σ 2 0 = = σ 2 k σ 2 k+1 = = σ 2 N 1 using normalized cumulative sum of squares
9 Homogeneity of Variance: II D + to define test statistic D, start with P k and then compute max k +1 0 k N 2 k j=0 Yj 2 N 1 j=0 Yj 2, k =0,...,N 2 N 1 P k & D max 0 k N 2 from which we form D max (D +,D ) can reject H 0 if observed D is too large can quantify too large by considering distribution of D under H 0 need to find critical value x α such that P[D x α ]=α for, e.g., α =0.01, 0.05 or 0.1 P k once determined, can perform α level test of H 0 : compute D statistic from data Y 0,...,Y N 1 reject H 0 at level α if D x α k N 1
10 Homogeneity of Variance: III can determine critical values x α in two ways Monte Carlo simulations large sample approximation to distribution of D: P[(N/2) 1/2 D x] 1+2 ( 1) l e 2l2 x 2 l=1 (reasonable approximation for N 128) idea: given time series X, compute D using W j,t =2 j/2 W j,2 j (t+1) 1, (L 2) j t N 2 1, j where L is length of j = 1 level wavelet filter and 2 j/2 W j,t L j 1 l=0 h j,l X t l mod N results in level by level tests above formulation allows for general N (i.e., N need not be multiple of 2 J 0) Q: is DWT decorrelation of FDPs good enough? Fig. 9 says yes!
11 Example: Nile River Minima recall MRA & wavelet variance plots application of homogeneity of variance test: scale D x 0.1 x 0.05 x year years years years if H 0 rejected, use nondecimated DWT to detect change point: compute rotated cumulative variance curve look for time of largest excursion from 0 Fig. 10: change point detection 720 AD for level j =1 722 AD for level j =2 agrees well with mosque construction in 715 AD interpretation differs from Beran & Terrin (1996)
12 Wavelet-Based Bootstrapping: I Davison & Hinkley, 1998, Chapter 8, discusses bootstrapping in context of time series analysis whitening allows wavelet-based bootstrapping for certain statistics (but not all) first example: lag 1 autocorrelation estimate ˆr 1 N 2 t=0 X t X t+1 N 1 t=0 X 2 t idea: to get standard error of ˆr 1, compute DWT of X sample with replacement from W j to form W (b) j (do same with V J0 ) synthesize X (b) using W (b) j s & V (b) J 0 compute ˆr (b) 1 for X (b) repeat until computer gets tired use standard error of ˆr (b) 1 s for X (b) s to assess standard error of ˆr 1 for X
13 Wavelet-Based Bootstrapping: II to test scheme, did Monte Carlo study involving AR(1) process: X t =0.9X t 1 + ɛ t MA(1) process: X t = ɛ t + ɛ t 1 FDP with d =0.45 average ˆr (b) 1 s have negligible bias comparison of standard errors, N = 128: LA(8) DWT LA(8) DWPT true AR(1) ± ± MA(1) ± ± FDP ± ± comparison of standard errors, N = 1024: LA(8) DWT LA(8) DWPT true AR(1) ± ± MA(1) ± ± FDP ± ± handles both short & long memory models
14 Wavelet-Based Bootstrapping: III second example: cross-correlation estimate ˆr (XY ) 0 N 1 ( N 1 t=0 X 2 t t=0 X t Y t N 1 t=0 Yt 2 ) 1/2 to assess null hypothesis r (XY ) 0 =0, separately generate X (b) & Y (b) bootstrapped ˆr (XY ) 0 should reflect variability in ˆr (XY ) 0 under null example: cross-correlation ˆr (XY ) 0 = 0.26 between maximum annual snow-pack level at Mt. Rainier Pacific decadel oscillation index can reject null at critical level p =0.01 lots of questions to be addressed! what is class of applicable statistics?
15 Estimation for FDPs: I extension of work by Wornell, McCoy & Walden problem: estimate d from time series U t such that where U t = T t + X t T t r j=0 a j t j is polynomial trend X t is FDP process, but can have d 1 2 DWT wavelet filter of width L has embedded differencing operation of order L/2 if L 2 r + 1, reduces polynomial trend to 0 can partition DWT coefficients as W = W s + W b + W w where W s has scaling coefficients and 0s elsewhere W s has boundary-dependent wavelet coefficients W w has boundary-independent wavelet coefficients
16 Estimation for FDPs: II since U = W T W, can write U = W T (W s + W b )+W T W w T + X example: microbolometric camera data can use values in W w to form likelihood: L(d, σ 2 ɛ ) J 0 N j j=1 t=1 where σ 2 ɛ is innovations variance; 1 ( 2πσ 2 j ) 1/2 e W 2 j,t+l j 1/(2σ2 j ) σj 2 1/2 H 1/2 j(f) df ; 2 sin(πf) 2d and H j ( ) is squared gain for {h j,l } σ 2 ɛ works well in Monte Carlo simulations for Hansen Lebedeff series, get ˆd. =0.40 ± 0.08 lots of questions to be addressed! what are properties of T? how to assess significance of T? (extension of Brillinger, 1994)
Wavelet Methods for Time Series Analysis. Part IV: Wavelet-Based Decorrelation of Time Series
Wavelet Methods for Time Series Analysis Part IV: Wavelet-Based Decorrelation of Time Series DWT well-suited for decorrelating certain time series, including ones generated from a fractionally differenced
More informationWavelet Methods for Time Series Analysis. Motivating Question
Wavelet Methods for Time Series Analysis Part VII: Wavelet-Based Bootstrapping start with some background on bootstrapping and its rationale describe adjustments to the bootstrap that allow it to work
More informationWavelet Methods for Time Series Analysis. Part IX: Wavelet-Based Bootstrapping
Wavelet Methods for Time Series Analysis Part IX: Wavelet-Based Bootstrapping start with some background on bootstrapping and its rationale describe adjustments to the bootstrap that allow it to work with
More informationTesting for Homogeneity of Variance in Time Series: Long Memory, Wavelets and the Nile River
Testing for Homogeneity of Variance in Time Series: Long Memory, Wavelets and the Nile River Brandon Whitcher Geophysical Statistics Project National Center for Atmospheric Research Boulder, Colorado Peter
More informationTesting for homogeneity of variance in time series: Long memory, wavelets, and the Nile River
WATER RESOURCES RESEARCH, VOL. 38, NO. 5, 1054, 10.109/001WR000509, 00 Testing for homogeneity of variance in time series: Long memory, wavelets, and the Nile River B. Whitcher, 1 S. D. Byers, P. Guttorp,
More informationA Tutorial on Stochastic Models and Statistical Analysis for Frequency Stability Measurements
A Tutorial on Stochastic Models and Statistical Analysis for Frequency Stability Measurements Don Percival Applied Physics Lab, University of Washington, Seattle overheads for talk available at http://staff.washington.edu/dbp/talks.html
More informationWavestrapping Time Series: Adaptive Wavelet-Based Bootstrapping
Wavestrapping Time Series: Adaptive Wavelet-Based Bootstrapping D. B. Percival, S. Sardy and A. C. Davison 1 Introduction Suppose we observe a time series that can be regarded as a realization of a portion
More informationWavelet-Based Bootstrapping for Non-Gaussian Time Series. Don Percival
Wavelet-Based Bootstrapping for Non-Gaussian Time Series Don Percival Applied Physics Laboratory Department of Statistics University of Washington Seattle, Washington, USA overheads for talk available
More informationWavelet Methods for Time Series Analysis
Wavelet Methods for Time Series Analysis Donald B. Percival UNIVERSITY OF WASHINGTON, SEATTLE Andrew T. Walden IMPERIAL COLLEGE OF SCIENCE, TECHNOLOGY AND MEDICINE, LONDON CAMBRIDGE UNIVERSITY PRESS Contents
More informationUsing wavelet tools to estimate and assess trends in atmospheric data
NRCSE Using wavelet tools to estimate and assess trends in atmospheric data Wavelets Fourier analysis uses big waves Wavelets are small waves Requirements for wavelets Integrate to zero Square integrate
More informationWavelet Analysis of Discrete Time Series
Wavelet Analysis of Discrete Time Series Andrew T. Walden Abstract. We give a brief review of some of the wavelet-based techniques currently available for the analysis of arbitrary-length discrete time
More informationA Tutorial on Stochastic Models and Statistical Analysis for Frequency Stability Measurements
A Tutorial on Stochastic Models and Statistical Analysis for Frequency Stability Measurements Don Percival Applied Physics Lab, University of Washington, Seattle overheads for talk available at http://staffwashingtonedu/dbp/talkshtml
More informationTesting for Homogeneity of Variance in Time Series: Long Memory, Wavelets and the Nile River NRCSE. T e c h n i c a l R e p o r t S e r i e s
Testing for Homogeneity of Variance in Time Series: Long Memory, Wavelets and the Nile River Brandon Whitcher Peter Guttorp Simon D Byers Donald B Percival NRCSE T e c h n i c a l R e p o r t S e r i e
More informationWavelet Methods for Time Series Analysis. Part II: Wavelet-Based Statistical Analysis of Time Series
Wavelet Methods for Time Series Analysis Part II: Wavelet-Based Statistical Analysis of Time Series topics to covered: wavelet variance (analysis phase of MODWT) wavelet-based signal extraction (synthesis
More informationDefining the Discrete Wavelet Transform (DWT)
Defining the Discrete Wavelet Transform (DWT) can formulate DWT via elegant pyramid algorithm defines W for non-haar wavelets (consistent with Haar) computes W = WX using O(N) multiplications brute force
More informationWavelet-based estimation for Gaussian time series and spatio-temporal processes
Wavelet-based estimation for Gaussian time series and spatio-temporal processes Dissertation Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School
More informationWavelet Methods for Time Series Analysis. Quick Comparison of the MODWT to the DWT
Wavelet Methods for Time Series Analysis Part IV: MODWT and Examples of DWT/MODWT Analysis MODWT stands for maximal overlap discrete wavelet transform (pronounced mod WT ) transforms very similar to the
More informationWavelet Methods for Time Series Analysis. Quick Comparison of the MODWT to the DWT
Wavelet Methods for Time Series Analysis Part III: MODWT and Examples of DWT/MODWT Analysis MODWT stands for maximal overlap discrete wavelet transform (pronounced mod WT ) transforms very similar to the
More informationAsymptotic Decorrelation of Between-Scale Wavelet Coefficients
Asymptotic Decorrelation of Between-Scale Wavelet Coefficients Peter F. Craigmile 1 and Donald B. Percival 2 1 Department of Statistics, The Ohio State University, 44 Cockins Hall, 1958 Neil Avenue, Columbus.
More informationNontechnical introduction to wavelets Continuous wavelet transforms Fourier versus wavelets - examples
Nontechnical introduction to wavelets Continuous wavelet transforms Fourier versus wavelets - examples A major part of economic time series analysis is done in the time or frequency domain separately.
More informationExamples of DWT & MODWT Analysis: Overview
Examples of DWT & MODWT Analysis: Overview look at DWT analysis of electrocardiogram (ECG) data discuss potential alignment problems with the DWT and how they are alleviated with the MODWT look at MODWT
More informationWavelet domain test for long range dependence in the presence of a trend
Wavelet domain test for long range dependence in the presence of a trend Agnieszka Jach and Piotr Kokoszka Utah State University July 24, 27 Abstract We propose a test to distinguish a weakly dependent
More informationMultiscale Detection and Location of Multiple Variance Changes in the Presence of Long Memory Brandon Whitcher EURANDOM Peter
Multiscale Detection and Location of Multiple Variance Changes in the Presence of Long Memory Brandon Whitcher EURANDOM whitcher@eurandom.tue.nl Peter Guttorp University of Washington peter@stat.washington.edu
More informationMultiresolution Models of Time Series
Multiresolution Models of Time Series Andrea Tamoni (Bocconi University ) 2011 Tamoni Multiresolution Models of Time Series 1/ 16 General Framework Time-scale decomposition General Framework Begin with
More informationDifferencing Revisited: I ARIMA(p,d,q) processes predicated on notion of dth order differencing of a time series {X t }: for d = 1 and 2, have X t
Differencing Revisited: I ARIMA(p,d,q) processes predicated on notion of dth order differencing of a time series {X t }: for d = 1 and 2, have X t 2 X t def in general = (1 B)X t = X t X t 1 def = ( X
More informationSerial Correlation Tests with Wavelets
Serial Correlation Tests with Wavelets Ramazan Gençay April 010, This version: November 011 Abstract This paper offers two new statistical tests for serial correlation with better power properties. The
More informationInference. Undercoverage of Wavelet-Based Resampling Confidence Intervals
Communications in Statistics Simulation and Computation, 37: 1307 1315, 2008 Copyright Taylor & Francis Group, LLC ISSN: 0361-0918 print/1532-4141 online DOI: 10.1080/03610910802050902 Inference Undercoverage
More informationFGN 1.0 PPL 1.0 FD
FGN PPL FD f f Figure SDFs for FGN, PPL and FD processes (top to bottom rows, respectively) on both linear/log and log/log aes (left- and right-hand columns, respectively) Each SDF S X ( ) is normalized
More informationDecline of Arctic Sea-Ice Thickness as Evidenced by Submarine Measurements. Don Percival
Decline of Arctic Sea-Ice Thickness as Evidenced by Submarine Measurements Don Percival Applied Physics Laboratory Department of Statistics University of Washington Seattle, Washington, USA http://facultywashingtonedu/dbp
More informationClassic Time Series Analysis
Classic Time Series Analysis Concepts and Definitions Let Y be a random number with PDF f Y t ~f,t Define t =E[Y t ] m(t) is known as the trend Define the autocovariance t, s =COV [Y t,y s ] =E[ Y t t
More informationQuantitative Finance II Lecture 10
Quantitative Finance II Lecture 10 Wavelets III - Applications Lukas Vacha IES FSV UK May 2016 Outline Discrete wavelet transformations: MODWT Wavelet coherence: daily financial data FTSE, DAX, PX Wavelet
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationWavelet-Generalized Least Squares: A New BLU Estimator of Linear Regression Models with 1/f Errors
NeuroImage 15, 217 232 (2002) doi:10.1006/nimg.2001.0955, available online at http://www.idealibrary.com on Wavelet-Generalized Least Squares: A New BLU Estimator of Linear Regression Models with 1/f Errors
More information[y i α βx i ] 2 (2) Q = i=1
Least squares fits This section has no probability in it. There are no random variables. We are given n points (x i, y i ) and want to find the equation of the line that best fits them. We take the equation
More informationDetection of Anomalous Observations Using Wavelet Analysis in Frequency Domain
International Journal of Scientific and Research Publications, Volume 7, Issue 8, August 2017 76 Detection of Anomalous Observations Using Wavelet Analysis in Frequency Domain Aideyan D.O. Dept of Mathematical
More informationEcon 424 Time Series Concepts
Econ 424 Time Series Concepts Eric Zivot January 20 2015 Time Series Processes Stochastic (Random) Process { 1 2 +1 } = { } = sequence of random variables indexed by time Observed time series of length
More informationEXAMINATION OF RELATIONSHIPS BETWEEN TIME SERIES BY WAVELETS: ILLUSTRATION WITH CZECH STOCK AND INDUSTRIAL PRODUCTION DATA
EXAMINATION OF RELATIONSHIPS BETWEEN TIME SERIES BY WAVELETS: ILLUSTRATION WITH CZECH STOCK AND INDUSTRIAL PRODUCTION DATA MILAN BAŠTA University of Economics, Prague, Faculty of Informatics and Statistics,
More informationStatistical signal processing
Statistical signal processing Short overview of the fundamentals Outline Random variables Random processes Stationarity Ergodicity Spectral analysis Random variable and processes Intuition: A random variable
More informationCointegrated VARIMA models: specification and. simulation
Cointegrated VARIMA models: specification and simulation José L. Gallego and Carlos Díaz Universidad de Cantabria. Abstract In this note we show how specify cointegrated vector autoregressive-moving average
More informationCalculating wavelet variance associated with discrete wavelet transform (DWT) 1
Calculating wavelet variance associated with discrete wavelet transform (DWT) 1 Xuejun Dong a and Haitao Li b a North Dakota State University Central Grasslands Research Extension Center Streeter, ND 58483
More informationTrend Assessment of Correlated Data
Trend Assessment of Correlated Data Malaak Kallache, Henning Rust and Jürgen Kropp Potsdam Institute for Climate Impact Research P.O. Box 60 12 03 14412 Potsdam, Germany Abstract The assessment of trends
More informationSimple and Multiple Linear Regression
Sta. 113 Chapter 12 and 13 of Devore March 12, 2010 Table of contents 1 Simple Linear Regression 2 Model Simple Linear Regression A simple linear regression model is given by Y = β 0 + β 1 x + ɛ where
More informationEconometrics 2, Class 1
Econometrics 2, Class Problem Set #2 September 9, 25 Remember! Send an email to let me know that you are following these classes: paul.sharp@econ.ku.dk That way I can contact you e.g. if I need to cancel
More informationWavelet Analysis of Variance for Time Series with Missing Values. Technical Report no. 535 Department of Statistics, University of Washington
Wavelet Analysis of Variance for Time Series with Missing Values Debashis Mondal 1 and Donald B. Percival Technical Report no. 535 Department of Statistics, University of Washington Abstract The wavelet
More informationWavelet Method for Locally Stationary Seasonal Long Memory Processes
Wavelet Method for Locally Stationary Seasonal Long Memory Processes Dominique Guegan, Zhiping Lu To cite this version: Dominique Guegan, Zhiping Lu. Wavelet Method for Locally Stationary Seasonal Long
More informationWavelet methods and null models for spatial pattern analysis
Wavelet methods and null models for spatial pattern analysis Pavel Dodonov Part of my PhD thesis, by the Federal University of São Carlos (São Carlos, SP, Brazil), supervised by Dr Dalva M. Silva-Matos
More informationSTAT 520: Forecasting and Time Series. David B. Hitchcock University of South Carolina Department of Statistics
David B. University of South Carolina Department of Statistics What are Time Series Data? Time series data are collected sequentially over time. Some common examples include: 1. Meteorological data (temperatures,
More informationWavelet-Based Test for Time Series Non-Stationarity 1
STATISTIKA 2015 95 (4) Wavelet-Based Test for Time Series Non-Stationarity 1 Milan Bašta 2 University of Economics, Prague, Czech Republic Abstract In the present paper, we propose a wavelet-based hypothesis
More informationProf. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis
Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation
More informationWavelet methods for the detection of anomalies and their application to network traffic analysis
Wavelet methods for the detection of anomalies and their application to network traffic analysis D.W. Kwon, K. Ko, M. Vannucci, A.L.N. Reddy, and S. Kim March 24, 2005 Abstract Here we develop an integrated
More informationFE570 Financial Markets and Trading. Stevens Institute of Technology
FE570 Financial Markets and Trading Lecture 5. Linear Time Series Analysis and Its Applications (Ref. Joel Hasbrouck - Empirical Market Microstructure ) Steve Yang Stevens Institute of Technology 9/25/2012
More informationNonstationary Time Series:
Nonstationary Time Series: Unit Roots Egon Zakrajšek Division of Monetary Affairs Federal Reserve Board Summer School in Financial Mathematics Faculty of Mathematics & Physics University of Ljubljana September
More informationSTAT 536: Genetic Statistics
STAT 536: Genetic Statistics Tests for Hardy Weinberg Equilibrium Karin S. Dorman Department of Statistics Iowa State University September 7, 2006 Statistical Hypothesis Testing Identify a hypothesis,
More informationWeek 5 Quantitative Analysis of Financial Markets Characterizing Cycles
Week 5 Quantitative Analysis of Financial Markets Characterizing Cycles Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036
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 informationAN INTRODUCTION TO SPECTRAL ANALYSIS AND WAVELETS
International Workshop on Advanced Mathematical Tools in Metrology Torino, October 1993 AN INTRODUCTION TO SPECTRAL ANALYSIS AND WAVELETS Donald B. Percival Applied Physics Laboratory, HN-10 University
More informationComputationally efficient banding of large covariance matrices for ordered data and connections to banding the inverse Cholesky factor
Computationally efficient banding of large covariance matrices for ordered data and connections to banding the inverse Cholesky factor Y. Wang M. J. Daniels wang.yanpin@scrippshealth.org mjdaniels@austin.utexas.edu
More informationReview of Concepts from Fourier & Filtering Theory. Fourier theory for finite sequences. convolution/filtering of infinite sequences filter cascades
Review of Concepts from Fourier & Filtering Theory precise definition of DWT requires a few basic concepts from Fourier analysis and theory of linear filters will start with discussion/review of: basic
More informationTMA4285 December 2015 Time series models, solution.
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 5 TMA4285 December 205 Time series models, solution. Problem a) (i) The slow decay of the ACF of z t suggest that
More informationMA Advanced Econometrics: Applying Least Squares to Time Series
MA Advanced Econometrics: Applying Least Squares to Time Series Karl Whelan School of Economics, UCD February 15, 2011 Karl Whelan (UCD) Time Series February 15, 2011 1 / 24 Part I Time Series: Standard
More informationTopic 4 Unit Roots. Gerald P. Dwyer. February Clemson University
Topic 4 Unit Roots Gerald P. Dwyer Clemson University February 2016 Outline 1 Unit Roots Introduction Trend and Difference Stationary Autocorrelations of Series That Have Deterministic or Stochastic Trends
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 informationMA 575 Linear Models: Cedric E. Ginestet, Boston University Non-parametric Inference, Polynomial Regression Week 9, Lecture 2
MA 575 Linear Models: Cedric E. Ginestet, Boston University Non-parametric Inference, Polynomial Regression Week 9, Lecture 2 1 Bootstrapped Bias and CIs Given a multiple regression model with mean and
More informationECON/FIN 250: Forecasting in Finance and Economics: Section 7: Unit Roots & Dickey-Fuller Tests
ECON/FIN 250: Forecasting in Finance and Economics: Section 7: Unit Roots & Dickey-Fuller Tests Patrick Herb Brandeis University Spring 2016 Patrick Herb (Brandeis University) Unit Root Tests ECON/FIN
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 informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 9 Jakub Mućk Econometrics of Panel Data Meeting # 9 1 / 22 Outline 1 Time series analysis Stationarity Unit Root Tests for Nonstationarity 2 Panel Unit Root
More informationCHANGE DETECTION IN TIME SERIES
CHANGE DETECTION IN TIME SERIES Edit Gombay TIES - 2008 University of British Columbia, Kelowna June 8-13, 2008 Outline Introduction Results Examples References Introduction sunspot.year 0 50 100 150 1700
More informationWavelet-based confidence intervals for the self-similarity parameter
Journal of Statistical Computation and Simulation Vol. 78, No. 12, December 2008, 1179 1198 Wavelet-based confidence intervals for the self-similarity parameter AGNIESZKA JACH and PIOTR KOKOSZKA* Department
More informationLinear models and their mathematical foundations: Simple linear regression
Linear models and their mathematical foundations: Simple linear regression Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/21 Introduction
More informationMULTI-SCALE TESTS FOR SERIAL CORRELATION. October 2013, First Version: April 2010
MULTI-SCALE TESTS FOR SERIAL CORRELATION RAMAZAN GENÇAY AND DANIELE SIGNORI October 2013, First Version: April 2010 Abstract This paper introduces a new family of portmanteau tests for serial correlation.
More informationBooth School of Business, University of Chicago Business 41914, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Midterm
Booth School of Business, University of Chicago Business 41914, Spring Quarter 017, Mr Ruey S Tsay Solutions to Midterm Problem A: (51 points; 3 points per question) Answer briefly the following questions
More informationPractical Statistics
Practical Statistics Lecture 1 (Nov. 9): - Correlation - Hypothesis Testing Lecture 2 (Nov. 16): - Error Estimation - Bayesian Analysis - Rejecting Outliers Lecture 3 (Nov. 18) - Monte Carlo Modeling -
More informationLecture 5: Unit Roots, Cointegration and Error Correction Models The Spurious Regression Problem
Lecture 5: Unit Roots, Cointegration and Error Correction Models The Spurious Regression Problem Prof. Massimo Guidolin 20192 Financial Econometrics Winter/Spring 2018 Overview Stochastic vs. deterministic
More informationMEI Exam Review. June 7, 2002
MEI Exam Review June 7, 2002 1 Final Exam Revision Notes 1.1 Random Rules and Formulas Linear transformations of random variables. f y (Y ) = f x (X) dx. dg Inverse Proof. (AB)(AB) 1 = I. (B 1 A 1 )(AB)(AB)
More information7. MULTIVARATE STATIONARY PROCESSES
7. MULTIVARATE STATIONARY PROCESSES 1 1 Some Preliminary Definitions and Concepts Random Vector: A vector X = (X 1,..., X n ) whose components are scalar-valued random variables on the same probability
More informationNonstationary time series models
13 November, 2009 Goals Trends in economic data. Alternative models of time series trends: deterministic trend, and stochastic trend. Comparison of deterministic and stochastic trend models The statistical
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 informationIt is easily seen that in general a linear combination of y t and x t is I(1). However, in particular cases, it can be I(0), i.e. stationary.
6. COINTEGRATION 1 1 Cointegration 1.1 Definitions I(1) variables. z t = (y t x t ) is I(1) (integrated of order 1) if it is not stationary but its first difference z t is stationary. It is easily seen
More informationTime Series 4. Robert Almgren. Oct. 5, 2009
Time Series 4 Robert Almgren Oct. 5, 2009 1 Nonstationarity How should you model a process that has drift? ARMA models are intrinsically stationary, that is, they are mean-reverting: when the value of
More informationSingle Equation Linear GMM with Serially Correlated Moment Conditions
Single Equation Linear GMM with Serially Correlated Moment Conditions Eric Zivot November 2, 2011 Univariate Time Series Let {y t } be an ergodic-stationary time series with E[y t ]=μ and var(y t )
More informationFinal Exam November 24, Problem-1: Consider random walk with drift plus a linear time trend: ( t
Problem-1: Consider random walk with drift plus a linear time trend: y t = c + y t 1 + δ t + ϵ t, (1) where {ϵ t } is white noise with E[ϵ 2 t ] = σ 2 >, and y is a non-stochastic initial value. (a) Show
More informationUtility of Correlation Functions
Utility of Correlation Functions 1. As a means for estimating power spectra (e.g. a correlator + WK theorem). 2. For establishing characteristic time scales in time series (width of the ACF or ACV). 3.
More informationAsymptotically sufficient statistics in nonparametric regression experiments with correlated noise
Vol. 0 0000 1 0 Asymptotically sufficient statistics in nonparametric regression experiments with correlated noise Andrew V Carter University of California, Santa Barbara Santa Barbara, CA 93106-3110 e-mail:
More informationStratified Market Equilibria: An Application of Maximal Overlap Discrete Wavelet Analysis
Empirical Economics Review 7(1): (March 017) ISSN -9736 Stratified Market Equilibria: An Application of Maximal Overlap Discrete Wavelet Analysis Wen-Den Chen Department of Economics, Tunghai University
More informationEconometric Methods for Panel Data
Based on the books by Baltagi: Econometric Analysis of Panel Data and by Hsiao: Analysis of Panel Data Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies
More informationIntroduction to Wavelets and Wavelet Transforms
Introduction to Wavelets and Wavelet Transforms A Primer C. Sidney Burrus, Ramesh A. Gopinath, and Haitao Guo with additional material and programs by Jan E. Odegard and Ivan W. Selesnick Electrical and
More informationWavelet-Based Multiresolution Analysis of Wivenhoe Dam Water Temperatures
WATER RESOURCES RESEARCH, VOL.???, XXXX, DOI:0.029/, 2 Wavelet-Based Multiresolution Analysis of Wivenhoe Dam Water Temperatures D. B. Percival, S. M. Lennox, 2 Y.-G. Wang, 2,3 and R. E. Darnell 2 S. M.
More informationChristopher Dougherty London School of Economics and Political Science
Introduction to Econometrics FIFTH EDITION Christopher Dougherty London School of Economics and Political Science OXFORD UNIVERSITY PRESS Contents INTRODU CTION 1 Why study econometrics? 1 Aim of this
More information1 Linear Difference Equations
ARMA Handout Jialin Yu 1 Linear Difference Equations First order systems Let {ε t } t=1 denote an input sequence and {y t} t=1 sequence generated by denote an output y t = φy t 1 + ε t t = 1, 2,... with
More informationApplied Time. Series Analysis. Wayne A. Woodward. Henry L. Gray. Alan C. Elliott. Dallas, Texas, USA
Applied Time Series Analysis Wayne A. Woodward Southern Methodist University Dallas, Texas, USA Henry L. Gray Southern Methodist University Dallas, Texas, USA Alan C. Elliott University of Texas Southwestern
More information3 Time Series Regression
3 Time Series Regression 3.1 Modelling Trend Using Regression Random Walk 2 0 2 4 6 8 Random Walk 0 2 4 6 8 0 10 20 30 40 50 60 (a) Time 0 10 20 30 40 50 60 (b) Time Random Walk 8 6 4 2 0 Random Walk 0
More informationWavelets and Multiresolution Processing. Thinh Nguyen
Wavelets and Multiresolution Processing Thinh Nguyen Multiresolution Analysis (MRA) A scaling function is used to create a series of approximations of a function or image, each differing by a factor of
More informationECE 541 Stochastic Signals and Systems Problem Set 11 Solution
ECE 54 Stochastic Signals and Systems Problem Set Solution Problem Solutions : Yates and Goodman,..4..7.3.3.4.3.8.3 and.8.0 Problem..4 Solution Since E[Y (t] R Y (0, we use Theorem.(a to evaluate R Y (τ
More informationHydrological statistics for engineering design in a varying climate
EGS - AGU - EUG Joint Assembly Nice, France, 6- April 23 Session HS9/ Climate change impacts on the hydrological cycle, extremes, forecasting and implications on engineering design Hydrological statistics
More informationHeteroskedasticity; Step Changes; VARMA models; Likelihood ratio test statistic; Cusum statistic.
47 3!,57 Statistics and Econometrics Series 5 Febrary 24 Departamento de Estadística y Econometría Universidad Carlos III de Madrid Calle Madrid, 126 2893 Getafe (Spain) Fax (34) 91 624-98-49 VARIANCE
More informationTrending Models in the Data
April 13, 2009 Spurious regression I Before we proceed to test for unit root and trend-stationary models, we will examine the phenomena of spurious regression. The material in this lecture can be found
More informationLocally Self-Similar Processes and Their Wavelet Analysis
Locally Self-Similar Processes and Their Wavelet Analysis Joseph E. Cavanaugh Department of Statistics, University of Missouri, Columbia, MO 65211, USA Yazhen Wang Department of Statistics, University
More informationCh 2: Simple Linear Regression
Ch 2: Simple Linear Regression 1. Simple Linear Regression Model A simple regression model with a single regressor x is y = β 0 + β 1 x + ɛ, where we assume that the error ɛ is independent random component
More informationTesting for Panel Unit Roots in the Presence of an Unknown Structural Break and Cross-Sectional Dependency
Testing for Panel Unit Roots in the Presence of an Unknown Structural Break and Cross-Sectional Dependency A. Almasri, K. Månsson 2, P. Sjölander 3 and G. Shukur 4 Department of Economics and Statistics,
More informationEC408 Topics in Applied Econometrics. B Fingleton, Dept of Economics, Strathclyde University
EC408 Topics in Applied Econometrics B Fingleton, Dept of Economics, Strathclyde University Applied Econometrics What is spurious regression? How do we check for stochastic trends? Cointegration and Error
More information