AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

Similar documents
Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

Elementary Differential Equations and Boundary Value Problems

innovations shocks white noise

Institute of Actuaries of India

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

Wave Equation (2 Week)

4.1 The Uniform Distribution Def n: A c.r.v. X has a continuous uniform distribution on [a, b] when its pdf is = 1 a x b

Components Model. Remember that we said that it was useful to think about the components representation

EXERCISE - 01 CHECK YOUR GRASP

Copyright 2012 Pearson Education, Inc. Publishing as Prentice Hall.

Wisconsin Unemployment Rate Forecast Revisited

Control System Engineering (EE301T) Assignment: 2

The Science of Monetary Policy

GDP Advance Estimate, 2016Q4

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

CSE 245: Computer Aided Circuit Simulation and Verification

Chapter 17 Handout: Autocorrelation (Serial Correlation)

EXST Regression Techniques Page 1

Lecture 4: Laplace Transforms

On the Speed of Heat Wave. Mihály Makai

Microscopic Flow Characteristics Time Headway - Distribution

Chapter 12 Introduction To The Laplace Transform

Physics 160 Lecture 3. R. Johnson April 6, 2015

University of Kansas, Department of Economics Economics 911: Applied Macroeconomics. Problem Set 2: Multivariate Time Series Analysis

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

I) Title: Rational Expectations and Adaptive Learning. II) Contents: Introduction to Adaptive Learning

46. Let y = ln r. Then dy = dr, and so. = [ sin (ln r) cos (ln r)

Discussion 06 Solutions

Applied Statistics and Probability for Engineers, 6 th edition October 17, 2016

SOLUTIONS. 1. Consider two continuous random variables X and Y with joint p.d.f. f ( x, y ) = = = 15. Stepanov Dalpiaz

MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

Properties of Autocorrelated Processes Economics 30331

Midterm Examination (100 pts)

Estimation Uncertainty

Instructors Solution for Assignment 3 Chapter 3: Time Domain Analysis of LTIC Systems

Ratio-Product Type Exponential Estimator For Estimating Finite Population Mean Using Information On Auxiliary Attribute

Ali Karimpour Associate Professor Ferdowsi University of Mashhad. Reference: System Identification Theory For The User Lennart Ljung

H is equal to the surface current J S

Economics 302 (Sec. 001) Intermediate Macroeconomic Theory and Policy (Spring 2011) 4/25/2011. UW Madison

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

Logistic equation of Human population growth (generalization to the case of reactive environment).

5. An object moving along an x-coordinate axis with its scale measured in meters has a velocity of 6t

Real time estimation of traffic flow and travel time Based on time series analysis

3.4 Repeated Roots; Reduction of Order

Mundell-Fleming I: Setup

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

3(8 ) (8 x x ) 3x x (8 )

Economics 302 (Sec. 001) Intermediate Macroeconomic Theory and Policy (Spring 2011) 3/28/2012. UW Madison

Improving the Power of the Diebold- Mariano Test for Least Squares Predictions

Charging of capacitor through inductor and resistor

(10) (a) Derive and plot the spectrum of y. Discuss how the seasonality in the process is evident in spectrum.

Let s look again at the first order linear differential equation we are attempting to solve, in its standard form:

Lecture 2: Current in RC circuit D.K.Pandey

Chapter 6 Differential Equations and Mathematical Modeling

Transfer function and the Laplace transformation

EE 350 Signals and Systems Spring 2005 Sample Exam #2 - Solutions

ERROR ANALYSIS A.J. Pintar and D. Caspary Department of Chemical Engineering Michigan Technological University Houghton, MI September, 2012

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

Poisson process Markov process

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional

REPETITION before the exam PART 2, Transform Methods. Laplace transforms: τ dτ. L1. Derive the formulas : L2. Find the Laplace transform F(s) if.

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

Almost power law : Tempered power-law models (T-FADE)

14.02 Principles of Macroeconomics Fall 2005 Quiz 3 Solutions

Math 266, Practice Midterm Exam 2

( A) ( B) ( C) ( D) ( E)

LaPlace Transform in Circuit Analysis

CHAPTER CHAPTER14. Expectations: The Basic Tools. Prepared by: Fernando Quijano and Yvonn Quijano

fiziks Institute for NET/JRF, GATE, IIT JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES MATEMATICAL PHYSICS SOLUTIONS are

Problem 2. Describe the following signals in terms of elementary functions (δ, u,r, ) and compute. x(t+2) x(2-t) RT_1[x] -3-2 = 1 2 = 1

A Simple Procedure to Calculate the Control Limit of Z Chart

Themes. Flexible exchange rates with inflation targeting. Expectations formation under flexible exchange rates

Week 06 Discussion Suppose a discrete random variable X has the following probability distribution: f ( 0 ) = 8

Forecasting in functional regressive or autoregressive models

Review Lecture 5. The source-free R-C/R-L circuit Step response of an RC/RL circuit. The time constant = RC The final capacitor voltage v( )

Chapter 6 Test December 9, 2008 Name

symmetric/hermitian matrices, and similarity transformations

1973 AP Calculus BC: Section I

First Lecture of Machine Learning. Hung-yi Lee

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

Power communication network traffic prediction based on two dimensional prediction algorithms Lei Xiao Xue Yang Min Zhu Lipeng Zhu

Midterm exam 2, April 7, 2009 (solutions)

EE 434 Lecture 22. Bipolar Device Models

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Linear Systems

Digital Image Processing

PFC Predictive Functional Control

Lecture 1: Growth and decay of current in RL circuit. Growth of current in LR Circuit. D.K.Pandey

Continuous probability distributions

INFLATION EXPECTATIONS DERIVED FROM BUSINESS TENDENCY SURVEY OF THE CENTRAL BANK

a dt a dt a dt dt If 1, then the poles in the transfer function are complex conjugates. Let s look at f t H t f s / s. So, for a 2 nd order system:

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

Search sequence databases 3 10/25/2016

The transition:transversion rate ratio vs. the T-ratio.

Interpreting coefficients for transformed variables

Stationary Time Series

Transcription:

AR() Procss Th firs-ordr auorgrssiv procss, AR() is whr is WN(0, σ )

Condiional Man and Varianc of AR() Condiional man: Condiional varianc: ) ( ) ( Ω Ω E E ) var( ) ) ( var( ) var( σ Ω Ω Ω Ω E

Auocovarianc of AR() Tak h quaion And hn mulipl boh sids b -k Thn ak xpcaions. Sinc is whi nois, i is uncorrlad wih or k k k E( k ) E( k ) E( k ) γ ( k) γ ( k )

Auocorrlaion of AR() Dividing b h varianc, his implis W know Thn ρ( k) ρ( k ) ρ( 0) ρ() ρ() ρ( k) ρ(0) ρ() k

Auocorrlaion of AR() W hav drivd ρ ( k) k Th auocorrlaion of h saionar AR() is a simpl gomric dca ( < ) If is small, h auocorrlaions dca rapidl o zro wih k If is larg (clos o ) hn h auocorrlaions dca modral Th AR() paramr dscribs h prsisnc in h im sris

On-Sp-Ahad Forcas As w showd arlir Thus E E ( Ω Th opimal on-sp-ahad forcas is a linar funcion of h final obsrvd valu ) ( ΩT ) T T

-sp-ahad forcas B back-subsiuion Thus and ( ) ) ( ) ( Ω Ω E E T T T E ) ( Ω

-sp-ahad forcas This shows ha h opimal -sp-ahad forcas is also a linar funcion of h final obsrvd valu, bu wih h cofficin. E( Ω T ) T T

h-sp-ahad forcas Similarl So Opimal forcas: h h h h h h h h h h h E E Ω Ω ) ( ) ( T h T h T E Ω ) (

Invrsion of AR() B invring h lag opraor Which is h sam as found b back subsiuion ( ) ( ) i i i i i i L L L 0 0

Condiion for Invribili Th opraor (-L) is invribl whn < This is h sam as for h MA() modl is h invrs of h roo of h polnomial -L Th roo of a funcion is h valu whr i crosss h x-axis Th roo of -L is /, h invrs of h roo is Invribili rquirs ha h invrs of h roo b lss han on

AR() wih Inrcp An AR() wih inrcp is Taking xpcaions α Thus and E ( ) α E( ) E( µ α µ α µ )

Bs Linar Prdicor A linar prdicor of givn - is α Th forcas rror is α Th linar prdicor which minimizs h xpcd squard forcas rror solvs min α, E( α )

Las-Squars Th sima of h xpcd squard linar forcas rror is h sum of squard rrors Th las squars sima ˆ ˆ α minimizs h sum of squard rrors, so is h sima of h bs linar prdicor This is a linar rgrssion, raing - as a rgrssor. ˆ

Unmplomn Ra. rgrss ur L.ur Sourc SS df MS Numbr of obs 88 F(, 86) 7943. Modl 5337.574 5337.574 Prob > F 0.0000 Rsidual 45.6976 86.9744886 R-squard 0.9560 Adj R-squard 0.9559 Toal 558.85 87 6.750755 Roo MSE.54539 ur Cof. Sd. Err. P> [95% Conf. Inrval] ur L..976788.00796 33.95 0.000.96465.9904 _cons.30.0709754 3. 0.00.089969.360636

Fid AR() To plo a scar and fid rgrssion, wowa scar ur L.ur lfi ur L.ur, il( Unmplomn Ra, 0-4 ar olds ) Unmplomn Ra, 0-4 ar olds 0 5 0 5 0 0 5 0 5 0 ur, L ur Fid valus

GDP Growh Ras. rgrss gdp L.gdp Sourc SS df MS Numbr of obs 78 F(, 76) 44.3 Modl 585.7068 585.7068 Prob > F 0.0000 Rsidual 3643.9846 76 3.0835 R-squard 0.384 Adj R-squard 0.35 Toal 49.5453 77 5.677059 Roo MSE 3.6336 gdp Cof. Sd. Err. P> [95% Conf. Inrval] gdp L..37475.0557984 6.66 0.000.66307.48397 _cons.036509.86904 7.0 0.000.480006.59303

GDP Growh Ras -0 0 0 0-0 0 0 0 Ral Gross Domsic Produc, L Ral Gross Domsic Produc Fid valus

On-Sp-Ahad Forcas Th opimal forcas for T givn T is ˆ T T α T Th forcas using h simas is ˆ T T α ˆ ˆ T

Unmplomn Ra, 0-4 ar olds Th simas wr ur Cof. Sd. Err. P> [95% Conf. Inrval] ur L..976788.00796 33.95 0.000.96465.9904 _cons.30.0709754 3. 0.00.089969.360636 0. 0.977 Th valu for Jan 06 is 8.3%, so ˆ : Poin forcas is 8.3% (unchangd from currn) ˆ 0. 0.977 8.3 07 8.3

lis im ur L.ur p if im>m(06m) L. im ur ur p 87. 06m 8. 9.4 9.4056 88. 06m 8.6 8. 8.30487 89. 06m3 8.4 8.6 8.678 80. 06m4 8.8 8.4 8.45833 8. 06m5 8.3 8.8 8.8654 8. 06m6 8.6 8.3 8.3859 83. 06m7 8.9 8.6 8.678 84. 06m8 8. 8.9 8.9497 85. 06m9 8. 8. 8.383 86. 06m0 8.4 8. 8.30487 87. 06m 8. 8.4 8.45833 88. 06m 8. 8. 8.383 89. 07m 8.3 8. 8.30487 830. 07m. 8.3 8.3859

Exampl GDP Growh Th simas wr gdp Cof. Sd. Err. P> [95% Conf. Inrval] gdp L..37475.0557984 6.66 0.000.66307.48397 _cons.036509.86904 7.0 0.000.480006.59303.04 0.37 Th valu for 4 h quarr 06 is.9%, so ˆ ˆ07 :.04 0.37.9.7%

GDP Growh L. im gdp gdp p 7. 05q.3.89090 73. 05q.6.77946 74. 05q3.6 3.00345 75. 05q4.9.77946 76. 06q.8.9.370837 77. 06q.4.8.333689 78. 06q3 3.5.4.556575 79. 06q4.9 3.5 3.336673 80. 07q..9.743

On-Sp-Ahad Forcas Error Th forcas rror is ˆ T T T α T T ( α ) T T Th forcas varianc is var( ˆ T T T ) var( T ) σ

Forcas varianc simaion Avrag of squard rsiduals whr h las-squars rsiduals ar ˆ ˆ ˆ ˆ σ σ σ T T ˆ ˆ ˆ α

Unmplomn Ra, 0-4 ar olds Unmplomn Ra, 0-4 ar olds Rsiduals ur 0 5 0 5 0 Rsiduals -3 - - 0 950m 960m 970m 980m 990m 000m 00m 00m im 950m 960m 970m 980m 990m 000m 00m 00m im

GDP Growh Ral Gross Domsic Produc -0 0 0 0 GDP Growh 950q 960q 970q 980q 990q 000q 00q 00q im Rsiduals -0-5 0 5 0 5 Rsiduals 950q 960q 970q 980q 990q 000q 00q 00q im

On-Sp-Ahad Inrvals Normal Mhod Assum forcas rror is normall disribud Forcas inrval is poin sima, plus and minus h sandard dviaion of forcas muliplid b a normal prcnil For a 95% inrval: ˆ ˆ ˆ T T ± σ z.05 T For a 90% inrval T ± σ ˆ ˆ ˆ ± σ z.05 T T T T ± σ ˆ ˆ.96.645

Sandard Dv of Forcas Th forcas varianc is h varianc of h rror, plus h varianc of h sima Th sandard dviaion of h forcas is is squar roo In mos cass h major componn of h forcas varianc is h varianc of h rror Thus a simpl mhod o sima h sandard dviaion of h forcas is o us σ, roo man squard rror, from h rgrssion sima In STATA, for sdf ou can us h command prdic s, sdf Th diffrnc bwn h rms and sdf grows whn h numbr of rgrssors is larg rlaiv o h sampl siz.

Exampl: GDP Th Roo MSE is 3.63 Th sdf is 3.64, onl slighl highr!. rgrss gdp L.gdp Sourc SS df MS Numbr of obs 78 F(, 76) 44.3 Modl 585.7068 585.7068 Prob > F 0.0000 Rsidual 3643.9846 76 3.0835 R-squard 0.384 Adj R-squard 0.35 Toal 49.5453 77 5.677059 Roo MSE 3.6336 gdp Cof. Sd. Err. P> [95% Conf. Inrval] gdp L..37475.0557984 6.66 0.000.66307.48397 _cons.036509.86904 7.0 0.000.480006.59303

Forcas Inrval Consrucion sappnd, add() prdic p prdic s, sdf gn p p.645*s gn p p.645*s lis im p p p s if im>q(06q4) im p p p s 80. 07q -3.4689.743 8.7355 3.640853 Poin sima.7% Sdf 3.64 90% Inrval [-3.%, 8.7%]

Assignmns Rad Dibold hrough Chapr 7 Problm S # 5 Du Tusda (/) Rad Chapr 5 from Th Signal and h Nois Rading Rflcion Du Thursda (/3)