Change point in trending regression
|
|
- Crystal Hill
- 5 years ago
- Views:
Transcription
1 Charles University, Prague ROBUST 2010 joint work with A. Aue, L.Horváth, J. Picek (Charles University, Prague)
2 1 (Charles University, Prague)
3 1 2 model-formulation (Charles University, Prague)
4 1 2 model-formulation 3 (Charles University, Prague)
5 1 2 model-formulation 3 4 (Charles University, Prague)
6 1 2 model-formulation (Charles University, Prague)
7 1 2 model-formulation (Charles University, Prague)
8 1 2 model-formulation (Charles University, Prague)
9 Procedures on stability of statistical models - structural breaks, disorder, stability, segmented, switching, change point problem, etc. Observations Y 1,..., Y n obtained at the ordered time points t 1 < < t n such that Y 1,..., Y k - model I Y k,..., Y n - model II k change point unknown The problem: to detect ( to test H 0 : no change & H 1 : there is a change), to identify k ( to estimate k ) to estimate the model before after the change. (Charles University, Prague)
10 Many variants multiple changes, abrupt changes, gradual changes, changes in various parameters, changes in distributions, independent observations, dependent observations. Construction of tests estimators various approaches as in most of the statistical problems. Theoretical problems theoretical problems. s (meteorology, climatology, hydrology or environmental studies, econometric time series, statistical quality control, etc.) (Charles University, Prague)
11 1 model-formulation 2 model-formulation (Charles University, Prague)
12 model-formulation model- formulation Y 1,..., Y n are observed at time points t 1 < < t n : Y i = x T i β + e i, i = 1..., k = x T i β + x T i δ + e i, i = k , n, e 1,..., e n rom errors i.i.d., zero mean, nonzero variance σ 2 finite E e i 2+ with some > 0 β, δ 0 parameters k... change point (Charles University, Prague)
13 model-formulation x 1,..., x n p-dim. design points (rom or nonrom): nontrending : 1 n k i=1 x ix T i k n C, k n (Charles University, Prague)
14 model-formulation x 1,..., x n p-dim. design points (rom or nonrom): nontrending : 1 n k i=1 x ix T i k n C, k n trending : x i = h(i/n), i = 1,..., n, h smooth nonconstant vector function Main problems: H 0 : no jump in & H 1 : at most one jump estimators of change points (Charles University, Prague)
15 model-formulation Test statistics T n = max p k<n p β k LSE of β based on Y 1,..., Y k {( βk β 0 k) T Σ 1 k ( βk β 0 k)} β 0 k LSE of β based on Y k+1,..., Y n 1 Σ k is an estimator of the variance matrix of β k β 0 k Equivalently T n = { max S T k p k<n p C 1 k C n(c 0 1 } k ) 1 S k σ n 2, S k = k i=1 x iê i, k = 1,..., n, ê i = Y i x T i β n, i = 1,..., n - residuals C k = k i=1 x ix T i, C 0 k = C n C k (Charles University, Prague)
16 model-formulation Test statistics T n = max p k<n p {( βk β 0 k) T Σ 1 k ( βk β 0 k)} β k LSE of β based on Y 1,..., Y k β 0 k LSE of β based on Y k+1,..., Y n 1 Σ k is an estimator of the variance matrix of β k β 0 k Equivalently T n = { max S T k p k<n p C 1 k C n(c 0 1 } k ) 1 S k σ n 2, S k = k i=1 x iê i, k = 1,..., n, ê i = Y i x T i β n, i = 1,..., n - residuals C k = k i=1 x ix T i, C 0 k = C n C k σ n 2 is an estimator of var e i = σ 2 (Charles University, Prague)
17 model-formulation Critical regions c n (α) critical value α level T n > c n (α) Approximation of the critical values: (i) limit distribution of T n under H 0 ; (ii) resampling methods (bootstrap) Estimator k of the change point k defined as such k it maximizes w.r.t. k { ( β k β 0 k) T Σ 1 k ( β k β 0 } k) (Charles University, Prague)
18 1 2 model-formulation (Charles University, Prague)
19 Test statistics: T n = max p k<n p T n (ε) = {( βk β 0 k) T Σ 1 k ( βk β 0 k)} { } max εn k (1 ε)n..., 0 < ε < 1/2 nontrending polynomial : a(y) = (2 log y) 1/2, lim P(a(log n)(t n) 1/2 t + b p (log n)) n = exp{ 2 exp{ t}}, t R 1, b p (y) = 2 log y + p 2 log log y log(2γ(p/2)), y > 1, x i = (1, i/n, (i/n) 2,..., (i/n) p 1 ) T, i = 1,..., n (Charles University, Prague)
20 nontrending { p T n (ε) d i=1 sup B2 i (t) } ε<t<1 ε t(1 t) {B j (t); t (0, 1)}, j = 1,..., p, independent Brownian bridges trending x i = h(i/n) S(t) = T n (ε) d t 0 sup S T (t)c(t)c 1 (1)C 0 (t)s(t) ε<t<1 ε h(x)db(x) C(t)C 1 (1) 1 with {B(x), x [0, 1]} being a Brownian bridge, C(t) = lim n 1 n C nt. 0 h(x)db(x), t [0, 1] (Charles University, Prague)
21 1 2 model-formulation (Charles University, Prague)
22 Assumptions (A.1) The sequence (e i : i 1) satisfy E[e i ] = 0, E[e 2 i ] = σ 2 > 0 (1) (A.2) There are independent stard Brownian motions (W 1,n (t): t 0) (W 2,n (t): t 0) such that 1 k max 1 k n/2 k 1/ν e i τw 1,n (k) = O P(1) (n ) (2) max n/2<k<n i=1 1 (n k) 1/ν n i=k+1 with some ν > 2 τ > 0. (Charles University, Prague) e i τw 2,n (n k) = O P(1) (n ) (3)
23 (A.3) The components of h(.) are continuous on [0, 1]. The matrices t 0 h(x)ht (x)dx 1 t h(x)ht (x)dx are regular for all t (t 0, 1 t 0 ) for all t 0 (0, 1/2). (A.4) There are p linearly independent p-dimensional vectors a 01,..., a 0p nonnegative 0 γ 01 <... < γ 0p such that lim sup t t γ 0p+1 h(t) p a 0l t γ 0l <. (4) l=1 (A.5) There are p linearly independent p-dimensional vectors a 11,..., a 1p nonnegative 0 γ 11 <... < γ 1p such that lim sup t 1 1 t γ 1p+1 h(t) p a 1l t γ 1l <. (5) l=1 (Charles University, Prague)
24 Particular cases: polynomial h(t) = (t γ 1,..., t γp ) T, t [0, 1], harmonic 0 γ 1 <... < γ p h(t) = (cos(2πtω 1 ), sin(2πtω 1 ),..., cos(2πtω p ), cos(2πtω p )) T, t [0, 1] ω 1,..., ω p known. independent e 1,..., e n : τ = σ. dependent e 1,..., e n : the estimator with flat top kernel used: τ 2 n = 1 n ê j = Y j x T j n (ê j e n ) q n n j=1 β n, e n = 1 n n i=1 êi j=1 n j w j i=1 (ê i e n )(ê i+j e n ) w j = 1I {1 j q n /2} + 2(1 j/q n )I {q n /2 < j q n } (Charles University, Prague)
25 1 2 model-formulation (Charles University, Prague)
26 Recall notation: T n = max p k<n p V k (τ 2 / τ 2 n ) V k = S T k C 1 k C n(c 0 k ) 1 S k /τ 2 S k = k i=1 x iê i, k = 1,..., n, ê i = Y i x T i β n, i = 1,..., n - residuals C k = k i=1 x ix T i, C 0 k = C n C k x i = h(i/n) (Charles University, Prague)
27 Main step of the proof The limit behavior of T n is the same as max(t n0, T n1 ) T n0 = T n1 = max p k<ns n V k max p n k<ns n V k for some s n 0 +, T n0 T n1 are asymptotically independent. T n0 has the same limit distribution as T 0 n0 = max Z T k p k<ns M 1 k Z k/τ 2 n Z k = k i=1 x ie i, M k = k i=1 x ix T i. T 0 n0 will not change if x i is replaced by x 0 i = Ax i with arbitrary nonsingular p p matrix A. max 1 k n/2 Z k τ k j=1 x j(w 1,n (j) W 1,n (j 1)) = o P (n κ ) for some κ > 0. (Charles University, Prague)
28 After some steps we get that Tn0 0 has the same distribution as 1 max k k Q k n k<ns n kτ QT k R 1 with Q k = k i=1 ((i/k)γ 01,..., (i/k) γ 0p ) T e i R k = ( 1 k k j=1 (j/k)γ 0v +γ 0r ) p v,r=1. Finally, we get that Tn0 0 has the same distribution as sup Q(t) T V(t) 1 Q(t) k n t ns n ( 1 t ) Q(t) = t γ x γ 0l dw n,1 (x); l = 1,..., p 0l+1/2 0 ( 1 t ) V(t) = t γ x γ 0l+γ 0l dx; l, l = 1,..., p 0l+γ 0l +1 0 After exponential transformation Q(t) T V(t) 1 Q(t) is a norm of stationary Gaussian process results of Albin, Piterbarg, Jaruskova etc. are applied. (Charles University, Prague)
29 1 2 model-formulation (Charles University, Prague)
30 Figure 1: Monthly air carrier traffic in the United Stated from January 1996 to March (Charles University, Prague)
31 Y j = 4+ 6j n + 3 2πj cos( 2 12n )+ 3 2π sin( 2 12n ) cos(2π 4n sin(2π 4n )+e j, j = 1,..., 200 e j either AR(1) or MA(1), normal distr. k = 100 either in the intercept or in one harmonic regressor 1000 repetitions critical value obtained through circular block bootstrap H (1) A H (1) A change in intercept, AR (1) or MA(1) change in one regressor, AR (1) or MA(1) (Charles University, Prague)
32 Figure 4: Time series plots of the processes under H (1) A with δ = 2 (upper panel) ROBUST H(2) A (lower 2010 joint work with A. Aue, L.Hor (Charles University, Prague)
33 Figure 2: Circular bootstrap distribution versus finite sample distribution for AR(1) innovations (left) MA(1) innovations (right) with scalings ϕ 2 = 2 (- - -), 1.5 ( - - -), 1.0 ( ), 0.5 (- - -). (Charles University, Prague)
34 Figure 3: Size-power curves for H (1) A (upper half) H(2) A (lower half) for the asymptotic modification (first third line) the the circular bootstrap (second fourth line) with AR(1) innovations (left) MA(1) innovations (right). (Charles University, Prague) 17
35 ROBUST joint work with A. Aue, L.Hor (Charles University, Prague)
36 ROBUST joint work with A. Aue, L.Hor (Charles University, Prague)
37 Data Monthly air traffic data model through the root of data, n = 159 Y j = β 0 + β 1 j/n + q ( l=1 β c l cos(2πω ellj/n) + β (s) l ) sin(2πω ell j/n) + e j q = 4, ω 1 = 2/160, ω 3 = 13/160, ω 3 = 40/160, ω 1 = 80/160 ω 2 annual cycle ω 3 quarterly cycle ω 4 two months cycle k = 69 (Charles University, Prague)
38 bwidth = Figure 6: periodogram (Charles University, (lower panel). Prague) Square root transformation of the monthly air carrier traffic data (upper panel) its
39 Figure 7: The fitted model based on the proposed data segmentation procedure (upper panel) based (Charles on auniversity, global fit (lower Prague) panel). The Change original point datainis trending shown both plots as dashed line.
Understanding Regressions with Observations Collected at High Frequency over Long Span
Understanding Regressions with Observations Collected at High Frequency over Long Span Yoosoon Chang Department of Economics, Indiana University Joon Y. Park Department of Economics, Indiana University
More informationTime Series and Forecasting Lecture 4 NonLinear Time Series
Time Series and Forecasting Lecture 4 NonLinear Time Series Bruce E. Hansen Summer School in Economics and Econometrics University of Crete July 23-27, 2012 Bruce Hansen (University of Wisconsin) Foundations
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationEcon 2148, fall 2017 Gaussian process priors, reproducing kernel Hilbert spaces, and Splines
Econ 2148, fall 2017 Gaussian process priors, reproducing kernel Hilbert spaces, and Splines Maximilian Kasy Department of Economics, Harvard University 1 / 37 Agenda 6 equivalent representations of the
More informationAbrupt change in mean avoiding variance estimation and block bootstrap. Barbora Peštová The Czech Academy of Sciences Institute of Computer Science
Abrupt change in mean avoiding variance estimation and block bootstrap Barbora Peštová The Czech Academy of Sciences Institute of Computer Science Structural changes 2 24 Motivation To know whether a change
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationRobustní monitorování stability v modelu CAPM
Robustní monitorování stability v modelu CAPM Ondřej Chochola, Marie Hušková, Zuzana Prášková (MFF UK) Josef Steinebach (University of Cologne) ROBUST 2012, Němčičky, 10.-14.9. 2012 Contents Introduction
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationStatistics & Data Sciences: First Year Prelim Exam May 2018
Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book
More informationCUSUM TEST FOR PARAMETER CHANGE IN TIME SERIES MODELS. Sangyeol Lee
CUSUM TEST FOR PARAMETER CHANGE IN TIME SERIES MODELS Sangyeol Lee 1 Contents 1. Introduction of the CUSUM test 2. Test for variance change in AR(p) model 3. Test for Parameter Change in Regression Models
More informationRegression of Time Series
Mahlerʼs Guide to Regression of Time Series CAS Exam S prepared by Howard C. Mahler, FCAS Copyright 2016 by Howard C. Mahler. Study Aid 2016F-S-9Supplement Howard Mahler hmahler@mac.com www.howardmahler.com/teaching
More informationUpper Bound for Intermediate Singular Values of Random Sub-Gaussian Matrices 1
Upper Bound for Intermediate Singular Values of Random Sub-Gaussian Matrices 1 Feng Wei 2 University of Michigan July 29, 2016 1 This presentation is based a project under the supervision of M. Rudelson.
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 informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationCointegrating Regressions with Messy Regressors: J. Isaac Miller
NASMES 2008 June 21, 2008 Carnegie Mellon U. Cointegrating Regressions with Messy Regressors: Missingness, Mixed Frequency, and Measurement Error J. Isaac Miller University of Missouri 1 Messy Data Example
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More informationLocal Whittle Likelihood Estimators and Tests for non-gaussian Linear Processes
Local Whittle Likelihood Estimators and Tests for non-gaussian Linear Processes By Tomohito NAITO, Kohei ASAI and Masanobu TANIGUCHI Department of Mathematical Sciences, School of Science and Engineering,
More informationNumerical Methods with Lévy Processes
Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:
More informationModeling Data with Linear Combinations of Basis Functions. Read Chapter 3 in the text by Bishop
Modeling Data with Linear Combinations of Basis Functions Read Chapter 3 in the text by Bishop A Type of Supervised Learning Problem We want to model data (x 1, t 1 ),..., (x N, t N ), where x i is a vector
More informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
More informationExtremogram and Ex-Periodogram for heavy-tailed time series
Extremogram and Ex-Periodogram for heavy-tailed time series 1 Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia) and Yuwei Zhao (Ulm) 1 Jussieu, April 9, 2014 1 2 Extremal
More informationMa 3/103: Lecture 24 Linear Regression I: Estimation
Ma 3/103: Lecture 24 Linear Regression I: Estimation March 3, 2017 KC Border Linear Regression I March 3, 2017 1 / 32 Regression analysis Regression analysis Estimate and test E(Y X) = f (X). f is the
More informationInference For High Dimensional M-estimates. Fixed Design Results
: Fixed Design Results Lihua Lei Advisors: Peter J. Bickel, Michael I. Jordan joint work with Peter J. Bickel and Noureddine El Karoui Dec. 8, 2016 1/57 Table of Contents 1 Background 2 Main Results and
More informationSwitching Regime Estimation
Switching Regime Estimation Series de Tiempo BIrkbeck March 2013 Martin Sola (FE) Markov Switching models 01/13 1 / 52 The economy (the time series) often behaves very different in periods such as booms
More informationInference For High Dimensional M-estimates: Fixed Design Results
Inference For High Dimensional M-estimates: Fixed Design Results Lihua Lei, Peter Bickel and Noureddine El Karoui Department of Statistics, UC Berkeley Berkeley-Stanford Econometrics Jamboree, 2017 1/49
More informationSTAT 100C: Linear models
STAT 100C: Linear models Arash A. Amini June 9, 2018 1 / 56 Table of Contents Multiple linear regression Linear model setup Estimation of β Geometric interpretation Estimation of σ 2 Hat matrix Gram matrix
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationNonparametric estimation of extreme risks from heavy-tailed distributions
Nonparametric estimation of extreme risks from heavy-tailed distributions Laurent GARDES joint work with Jonathan EL METHNI & Stéphane GIRARD December 2013 1 Introduction to risk measures 2 Introduction
More informationScaling exponents for certain 1+1 dimensional directed polymers
Scaling exponents for certain 1+1 dimensional directed polymers Timo Seppäläinen Department of Mathematics University of Wisconsin-Madison 2010 Scaling for a polymer 1/29 1 Introduction 2 KPZ equation
More informationWeak quenched limiting distributions of a one-dimensional random walk in a random environment
Weak quenched limiting distributions of a one-dimensional random walk in a random environment Jonathon Peterson Cornell University Department of Mathematics Joint work with Gennady Samorodnitsky September
More informationTime Series. Anthony Davison. c
Series Anthony Davison c 2008 http://stat.epfl.ch Periodogram 76 Motivation............................................................ 77 Lutenizing hormone data..................................................
More informationBickel Rosenblatt test
University of Latvia 28.05.2011. A classical Let X 1,..., X n be i.i.d. random variables with a continuous probability density function f. Consider a simple hypothesis H 0 : f = f 0 with a significance
More informationA Modern Look at Classical Multivariate Techniques
A Modern Look at Classical Multivariate Techniques Yoonkyung Lee Department of Statistics The Ohio State University March 16-20, 2015 The 13th School of Probability and Statistics CIMAT, Guanajuato, Mexico
More informationRadial Basis Functions I
Radial Basis Functions I Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo November 14, 2008 Today Reformulation of natural cubic spline interpolation Scattered
More informationExtremogram and ex-periodogram for heavy-tailed time series
Extremogram and ex-periodogram for heavy-tailed time series 1 Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia) and Yuwei Zhao (Ulm) 1 Zagreb, June 6, 2014 1 2 Extremal
More 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 information7 Day 3: Time Varying Parameter Models
7 Day 3: Time Varying Parameter Models References: 1. Durbin, J. and S.-J. Koopman (2001). Time Series Analysis by State Space Methods. Oxford University Press, Oxford 2. Koopman, S.-J., N. Shephard, and
More informationGibbs Sampling in Linear Models #2
Gibbs Sampling in Linear Models #2 Econ 690 Purdue University Outline 1 Linear Regression Model with a Changepoint Example with Temperature Data 2 The Seemingly Unrelated Regressions Model 3 Gibbs sampling
More informationOrder book resilience, price manipulation, and the positive portfolio problem
Order book resilience, price manipulation, and the positive portfolio problem Alexander Schied Mannheim University Workshop on New Directions in Financial Mathematics Institute for Pure and Applied Mathematics,
More informationTESTING PROCEDURES BASED ON THE EMPIRICAL CHARACTERISTIC FUNCTIONS II: k-sample PROBLEM, CHANGE POINT PROBLEM. 1. Introduction
Tatra Mt. Math. Publ. 39 (2008), 235 243 t m Mathematical Publications TESTING PROCEDURES BASED ON THE EMPIRICAL CHARACTERISTIC FUNCTIONS II: k-sample PROBLEM, CHANGE POINT PROBLEM Marie Hušková Simos
More informationGaussian with mean ( µ ) and standard deviation ( σ)
Slide from Pieter Abbeel Gaussian with mean ( µ ) and standard deviation ( σ) 10/6/16 CSE-571: Robotics X ~ N( µ, σ ) Y ~ N( aµ + b, a σ ) Y = ax + b + + + + 1 1 1 1 1 1 1 1 1 1, ~ ) ( ) ( ), ( ~ ), (
More informationStochastic volatility models: tails and memory
: tails and memory Rafa l Kulik and Philippe Soulier Conference in honour of Prof. Murad Taqqu 19 April 2012 Rafa l Kulik and Philippe Soulier Plan Model assumptions; Limit theorems for partial sums and
More informationAsymptotic Statistics-III. Changliang Zou
Asymptotic Statistics-III Changliang Zou The multivariate central limit theorem Theorem (Multivariate CLT for iid case) Let X i be iid random p-vectors with mean µ and and covariance matrix Σ. Then n (
More information5: MULTIVARATE STATIONARY PROCESSES
5: MULTIVARATE STATIONARY PROCESSES 1 1 Some Preliminary Definitions and Concepts Random Vector: A vector X = (X 1,..., X n ) whose components are scalarvalued random variables on the same probability
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationRandom Toeplitz Matrices
Arnab Sen University of Minnesota Conference on Limits Theorems in Probability, IISc January 11, 2013 Joint work with Bálint Virág What are Toeplitz matrices? a0 a 1 a 2... a1 a0 a 1... a2 a1 a0... a (n
More informationSums of exponentials of random walks
Sums of exponentials of random walks Robert de Jong Ohio State University August 27, 2009 Abstract This paper shows that the sum of the exponential of an oscillating random walk converges in distribution,
More informationStatistics of stochastic processes
Introduction Statistics of stochastic processes Generally statistics is performed on observations y 1,..., y n assumed to be realizations of independent random variables Y 1,..., Y n. 14 settembre 2014
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 informationVector Auto-Regressive Models
Vector Auto-Regressive Models Laurent Ferrara 1 1 University of Paris Nanterre M2 Oct. 2018 Overview of the presentation 1. Vector Auto-Regressions Definition Estimation Testing 2. Impulse responses functions
More informationVast Volatility Matrix Estimation for High Frequency Data
Vast Volatility Matrix Estimation for High Frequency Data Yazhen Wang National Science Foundation Yale Workshop, May 14-17, 2009 Disclaimer: My opinion, not the views of NSF Y. Wang (at NSF) 1 / 36 Outline
More informationGARCH Models Estimation and Inference. Eduardo Rossi University of Pavia
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationLectures on Structural Change
Lectures on Structural Change Eric Zivot Department of Economics, University of Washington April5,2003 1 Overview of Testing for and Estimating Structural Change in Econometric Models 1. Day 1: Tests of
More informationVAR Models and Applications
VAR Models and Applications Laurent Ferrara 1 1 University of Paris West M2 EIPMC Oct. 2016 Overview of the presentation 1. Vector Auto-Regressions Definition Estimation Testing 2. Impulse responses functions
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationCentral Bank of Chile October 29-31, 2013 Bruce Hansen (University of Wisconsin) Structural Breaks October 29-31, / 91. Bruce E.
Forecasting Lecture 3 Structural Breaks Central Bank of Chile October 29-31, 2013 Bruce Hansen (University of Wisconsin) Structural Breaks October 29-31, 2013 1 / 91 Bruce E. Hansen Organization Detection
More informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationFFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations
FFTs in Graphics and Vision Homogenous Polynomials and Irreducible Representations 1 Outline The 2π Term in Assignment 1 Homogenous Polynomials Representations of Functions on the Unit-Circle Sub-Representations
More informationVector spaces. DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis.
Vector spaces DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Vector space Consists of: A set V A scalar
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 informationLECTURE 2 LINEAR REGRESSION MODEL AND OLS
SEPTEMBER 29, 2014 LECTURE 2 LINEAR REGRESSION MODEL AND OLS Definitions A common question in econometrics is to study the effect of one group of variables X i, usually called the regressors, on another
More informationGaussian Random Fields: Geometric Properties and Extremes
Gaussian Random Fields: Geometric Properties and Extremes Yimin Xiao Michigan State University Outline Lecture 1: Gaussian random fields and their regularity Lecture 2: Hausdorff dimension results and
More information2018 Fall 2210Q Section 013 Midterm Exam I Solution
8 Fall Q Section 3 Midterm Exam I Solution True or False questions ( points = points) () An example of a linear combination of vectors v, v is the vector v. True. We can write v as v + v. () If two matrices
More informationNonparametric Drift Estimation for Stochastic Differential Equations
Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,
More informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationNonparametric Function Estimation with Infinite-Order Kernels
Nonparametric Function Estimation with Infinite-Order Kernels Arthur Berg Department of Statistics, University of Florida March 15, 2008 Kernel Density Estimation (IID Case) Let X 1,..., X n iid density
More informationA Bootstrap View on Dickey-Fuller Control Charts for AR(1) Series
AUSTRIAN JOURNAL OF STATISTICS Volume 35 (26), Number 2&3, 339 346 A Bootstrap View on Dickey-Fuller Control Charts for AR(1) Series Ansgar Steland RWTH Aachen University, Germany Abstract: Dickey-Fuller
More informationDS-GA 1002 Lecture notes 10 November 23, Linear models
DS-GA 2 Lecture notes November 23, 2 Linear functions Linear models A linear model encodes the assumption that two quantities are linearly related. Mathematically, this is characterized using linear functions.
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 informationNotes on Time Series Modeling
Notes on Time Series Modeling Garey Ramey University of California, San Diego January 17 1 Stationary processes De nition A stochastic process is any set of random variables y t indexed by t T : fy t g
More informationThreshold Autoregressions and NonLinear Autoregressions
Threshold Autoregressions and NonLinear Autoregressions Original Presentation: Central Bank of Chile October 29-31, 2013 Bruce Hansen (University of Wisconsin) Threshold Regression 1 / 47 Threshold Models
More informationIn terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.
1. GAUSSIAN PROCESSES A Gaussian process on a set T is a collection of random variables X =(X t ) t T on a common probability space such that for any n 1 and any t 1,...,t n T, the vector (X(t 1 ),...,X(t
More informationHigher order weak approximations of stochastic differential equations with and without jumps
Higher order weak approximations of stochastic differential equations with and without jumps Hideyuki TANAKA Graduate School of Science and Engineering, Ritsumeikan University Rough Path Analysis and Related
More informationUCSD ECE250 Handout #27 Prof. Young-Han Kim Friday, June 8, Practice Final Examination (Winter 2017)
UCSD ECE250 Handout #27 Prof. Young-Han Kim Friday, June 8, 208 Practice Final Examination (Winter 207) There are 6 problems, each problem with multiple parts. Your answer should be as clear and readable
More informationA BOOTSTRAP VIEW ON DICKEY-FULLER CONTROL CHARTS FOR AR(1) SERIES
A BOOTSTRAP VIEW ON DICKEY-FULLER CONTROL CHARTS FOR AR(1) SERIES Ansgar Steland Ruhr-Universität Bochum, Germany Abstract: Dickey-Fuller control charts aim at monitoring a random walk until a given time
More informationRobust Testing and Variable Selection for High-Dimensional Time Series
Robust Testing and Variable Selection for High-Dimensional Time Series Ruey S. Tsay Booth School of Business, University of Chicago May, 2017 Ruey S. Tsay HTS 1 / 36 Outline 1 Focus on high-dimensional
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationBayesian non-parametric model to longitudinally predict churn
Bayesian non-parametric model to longitudinally predict churn Bruno Scarpa Università di Padova Conference of European Statistics Stakeholders Methodologists, Producers and Users of European Statistics
More informationEC3062 ECONOMETRICS. THE MULTIPLE REGRESSION MODEL Consider T realisations of the regression equation. (1) y = β 0 + β 1 x β k x k + ε,
THE MULTIPLE REGRESSION MODEL Consider T realisations of the regression equation (1) y = β 0 + β 1 x 1 + + β k x k + ε, which can be written in the following form: (2) y 1 y 2.. y T = 1 x 11... x 1k 1
More informationHETEROSCEDASTICITY AND AUTOCORRELATION ROBUST STRUCTURAL CHANGE DETECTION. By Zhou Zhou 1 University of Toronto February 1, 2013.
HETEROSCEDASTICITY AND AUTOCORRELATION ROBUST STRUCTURAL CHANGE DETECTION By Zhou Zhou 1 University of Toronto February 1, 2013 Abstract The assumption of (weak) stationarity is crucial for the validity
More informationRecent Advances in the Field of Trade Theory and Policy Analysis Using Micro-Level Data
Recent Advances in the Field of Trade Theory and Policy Analysis Using Micro-Level Data July 2012 Bangkok, Thailand Cosimo Beverelli (World Trade Organization) 1 Content a) Classical regression model b)
More informationClass 1: Stationary Time Series Analysis
Class 1: Stationary Time Series Analysis Macroeconometrics - Fall 2009 Jacek Suda, BdF and PSE February 28, 2011 Outline Outline: 1 Covariance-Stationary Processes 2 Wold Decomposition Theorem 3 ARMA Models
More informationMOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS
MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties
More informationPreliminary Exam: Probability 9:00am 2:00pm, Friday, January 6, 2012
Preliminary Exam: Probability 9:00am 2:00pm, Friday, January 6, 202 The exam lasts from 9:00am until 2:00pm, with a walking break every hour. Your goal on this exam should be to demonstrate mastery of
More informationBootstrapping the triangles
Bootstrapping the triangles Prečo, kedy a ako (NE)bootstrapovat v trojuholníkoch? Michal Pešta Charles University in Prague Faculty of Mathematics and Physics Actuarial Seminar, Prague 19 October 01 Overview
More informationIMC 2015, Blagoevgrad, Bulgaria
IMC 05, Blagoevgrad, Bulgaria Day, July 9, 05 Problem. For any integer n and two n n matrices with real entries, B that satisfy the equation + B ( + B prove that det( det(b. Does the same conclusion follow
More informationDetection of structural breaks in multivariate time series
Detection of structural breaks in multivariate time series Holger Dette, Ruhr-Universität Bochum Philip Preuß, Ruhr-Universität Bochum Ruprecht Puchstein, Ruhr-Universität Bochum January 14, 2014 Outline
More informationMinimax Rate of Convergence for an Estimator of the Functional Component in a Semiparametric Multivariate Partially Linear Model.
Minimax Rate of Convergence for an Estimator of the Functional Component in a Semiparametric Multivariate Partially Linear Model By Michael Levine Purdue University Technical Report #14-03 Department of
More informationWeek 9-10: Recurrence Relations and Generating Functions
Week 9-10: Recurrence Relations and Generating Functions April 3, 2017 1 Some number sequences An infinite sequence (or just a sequence for short is an ordered array a 0, a 1, a 2,..., a n,... of countably
More informationOn 1.9, you will need to use the facts that, for any x and y, sin(x+y) = sin(x) cos(y) + cos(x) sin(y). cos(x+y) = cos(x) cos(y) - sin(x) sin(y).
On 1.9, you will need to use the facts that, for any x and y, sin(x+y) = sin(x) cos(y) + cos(x) sin(y). cos(x+y) = cos(x) cos(y) - sin(x) sin(y). (sin(x)) 2 + (cos(x)) 2 = 1. 28 1 Characteristics of Time
More informationA New Approach to Robust Inference in Cointegration
A New Approach to Robust Inference in Cointegration Sainan Jin Guanghua School of Management, Peking University Peter C. B. Phillips Cowles Foundation, Yale University, University of Auckland & University
More informationMultivariate Normal-Laplace Distribution and Processes
CHAPTER 4 Multivariate Normal-Laplace Distribution and Processes The normal-laplace distribution, which results from the convolution of independent normal and Laplace random variables is introduced by
More informationSome Aspects of Universal Portfolio
1 Some Aspects of Universal Portfolio Tomoyuki Ichiba (UC Santa Barbara) joint work with Marcel Brod (ETH Zurich) Conference on Stochastic Asymptotics & Applications Sixth Western Conference on Mathematical
More informationAsymptotic results for empirical measures of weighted sums of independent random variables
Asymptotic results for empirical measures of weighted sums of independent random variables B. Bercu and W. Bryc University Bordeaux 1, France Workshop on Limit Theorems, University Paris 1 Paris, January
More informationEconometrics II - EXAM Answer each question in separate sheets in three hours
Econometrics II - EXAM Answer each question in separate sheets in three hours. Let u and u be jointly Gaussian and independent of z in all the equations. a Investigate the identification of the following
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationMonte Carlo simulation inspired by computational optimization. Colin Fox Al Parker, John Bardsley MCQMC Feb 2012, Sydney
Monte Carlo simulation inspired by computational optimization Colin Fox fox@physics.otago.ac.nz Al Parker, John Bardsley MCQMC Feb 2012, Sydney Sampling from π(x) Consider : x is high-dimensional (10 4
More information