1 / 43

Econ 240C

Econ 240C. Lecture Five. Part I: Time Series Components Model. The conceptual framework for inertial (mechanical) time series models: Time series (t) = trend + cycle + seasonal + residual We are familiar with trend models, e.g. Time series = a + b*t + e(t) , i.e.

yaholo
Télécharger la présentation

Econ 240C

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Econ 240C Lecture Five

  2. Part I: Time Series Components Model • The conceptual framework for inertial (mechanical) time series models: • Time series (t) = trend + cycle + seasonal + residual • We are familiar with trend models, e.g. • Time series = a + b*t + e(t) , i.e. • time series = trend + residual where e(t) is i.i.d. N(0, s2)

  3. Time Series Components Model • We also know how to deal with seasonality. For example, using quarterly data we could add a dummy zero-one variable, D1 that takes on the value of one if the observation is for the first quarter and zero otherwise. Similarly, we could add dummy variables for second quarter observations and for third quarter observations: Time series = a + b*t + c1*D1 + c2*D2 + c3*D3 + e(t)

  4. Time Series Components Models • So we have: time series = trend + seasonal + residual • But how do we model cycles? • Since macroeconomic variables are likely to be affected by economic conditions and the business cycle, this is an important question. • The answer lies in Box-Jenkins or ARIMA models.

  5. Time Series Components Models • ARIMA models are about modeling the residual • The simplest time series model is: time series(t) = white noise(t) • some other time series models are of the form: time series = A(Z)*white noise(t), where A(Z) is a polynomial in Z, a dynamic multiplier for white noise. • For the white noise model, A(Z) = Z0 =1

  6. Time Series Components Models • For the random walk, RW(t), RW(t) = A(Z)*WN(t) where A(Z) = (1+Z+Z2 +Z3 + …) • For the autoregressive process of the first order, ARONE(t), ARONE(t) = A(Z) * WN(t) where A(Z) = (1+b*Z+b2*Z2 +b3*Z3 +…) and -1<b<1 , i.e. b is on the real number line and is less than one in absolute value

  7. Time Series Components Models • So ARIMA models have a residual, white noise, as an input, and transform it with the polynomial in lag, A(Z), to model time series behavior. • One can think of ARIMA models in terms of the time series components model, where the time series, y(t), for quarterly data is modeled as: y(t)=a+b*t+c1*D1+c2*D2 +c3*D3+A(Z)WN(t)

  8. Time Series Components Models • But y(t), with a trend component, and a seasonal component, is evolutionary, i.e. time dependent, on two counts. So first we difference the time series, y(t), to remove trend, and seasonally difference it to remove the seasonal component, making it stationary. Then we can model it as an ARMA model, i.e. an autoregressive- moving average time series.

  9. Time Series Components Models • Symbolically, we difference, D, y(t) to remove trend, obtaining Dy(t) • Then we seasonally difference, DS, Dy(t) to remove the seasonality, obtaining D DS y(t). • Now we can model this stationary time series, D DS y(t) as ARMA, e.g. D DS y(t) = A(Z)*WN(t)

  10. Time Series Components Models • After modeling D DS y(t) as an ARMA process, we can recover the model for the original time series, y(t), by undoing the differencing and seasonal differencing. • This is accomplished by summation, i.e. integration, the inverse of differencing. Hence the name autoregressive integrated moving average, or ARIMA, for the model of y(t).

  11. Part II. Behavior of Autoregressive Processes of the First Order • From PowerPoint Lecture Three

  12. Model Three: Autoregressive Time Series of Order One • An analogy to our model of trend plus shock for the logarithm of the Standard Poors is inertia plus shock for an economic time series such as the ratio of inventory to sales for total business • Source: FRED http://research.stlouisfed.org/fred/

  13. Trace of Inventory to Sales, Total Business

  14. Behavior of ARONE Processes • So we have a typical trace of an ARONE • How about the histogram? • How about the correlogram?

  15. Histogram: Ratio of Inventory to Sales, Total Business

  16. Ratio of Inventory to Sales Sample: 1992:01 2003:01 Included observations: 133 Autocorrelation Partial Correlation AC PAC Q-Stat Prob .|*******| .|*******| 1 0.928 0.928 117.23 0.000 .|*******| .|* | 2 0.886 0.175 224.83 0.000 .|*******| .|* | 3 0.851 0.074 324.86 0.000 .|****** | *|. | 4 0.803 -0.087 414.54 0.000 .|****** | .|. | 5 0.764 0.019 496.51 0.000 .|****** | *|. | 6 0.715 -0.092 568.80 0.000 .|***** | .|. | 7 0.665 -0.048 631.90 0.000 .|***** | *|. | 8 0.611 -0.086 685.61 0.000 .|**** | .|. | 9 0.563 0.001 731.48 0.000 .|**** | .|. | 10 0.513 -0.038 769.87 0.000 .|**** | .|. | 11 0.462 -0.028 801.26 0.000 .|*** | .|. | 12 0.416 -0.006 826.93 0.000

  17. Part III. Characterizing Autoregressive Processes of the First Order • ARONE(t) = b*ARONE(t-1) + WN(t) • Lag by one • ARONE(t-1) = b*ARONE(t-2) + WN(t-1) • Substitute for ARONE(t-1) • ARONE(t) = b*[b*ARONE(t-2) + WN(t-1] + WN(t) • ARONE(t) = WN(t) + b*WN(t-1) +b2*ARONE(t-2)

  18. Characterize AR of 1st Order • Keep lagging and substituting to obtain • ARONE(t) = WN(t) +b*WN(t-1) + b2*WN(t-2) + ….. • ARONE(t) = [1+b*Z+b2Z2+…] WN(t) • ARONE(t) = A(Z)*WN(t)

  19. Characterize AR of 1st Order • ARONE(t) = WN(t) +b*WN(t-1) + b2*WN(t-2) + ….. • Note that the mean function of an ARONE process is zero • m(t) = E ARONE(t) = E{WN(t) + b*WN(t-1) + b2*WN(t-2) + …..} where E WN(t) =0, and EWN(t-1) =0 etc. • m(t) = 0

  20. WN(t) and WN(t-1)

  21. Autocovariance of ARONE • E{[ARONE(t) - EARONE(t)]*[ARONE(t-u)- EARONE(t-u)]}=E{ARONE(t)*ARONE(t-u)] since EARONE(t) = 0 = EARONE(t-u) • So gAR,AR(u) = E{ARONE(t)*ARONE(t-u)} • For u=1, i.e. lag one, gAR,AR(1) = E{ARONE(t)*ARONE(t-1)}, and • use ARONE(t) = b*ARONE(t-1) + WN(t) • and multiply byARONE(t-1)

  22. Autocovariance of ARONE • ARONE(t)*ARONE(t-1) = b*[ARONE(t-1)]2 +ARONE(t-1)*WN(t) • and take expectations, E • E{ARONE(t)*ARONE(t-1) = b*[ARONE(t-1)]2 +ARONE(t-1)*WN(t)} • where the LHS E{ARONE(t)*ARONE(t-1) is gAR,AR(1) by definition and • b*E *[ARONE(t-1)]2 is b*gAR,AR(0) , i.e. b* the variance by definition but how about

  23. Autocovariance of an ARONE • E{ARONE(t-1)*WN(t)} = ? • Note that ARONE(t) = WN(t) +b*WN(t-1) + b2*WN(t-2) + ….. • And lagging by one, ARONE(t-1) = WN(t-1) +b*WN(t-2) + b2*WN(t-3) + ….. • So ARONE(t-1) depends on WN(t-1) and earlier shocks, so that E{ARONE(t-1)*WN(t)} = 0, i.e. ARONE(t) is independent of WN(t).

  24. Autocovariance of ARONE • In sum, gAR,AR(1) = b*gAR,AR(0) • or in general for an ARONE, gAR,AR(u) = b*gAR,AR(u-1) which can be confirmed by taking the formula : • ARONE(t) = b*ARONE(t-1) + WN(t), multiplying by ARONE(t-u) and taking expectations. • Note rAR,AR (u) = gAR,AR(u) / gAR,AR(0)

  25. Autocorrelation of ARONE(t) • So dividing gAR,AR(u) = b*gAR,AR(u-1) by gAR,AR(0) results in • rAR,AR (u) = b* rAR,AR (u-1), u>0 • rAR,AR (1) = b* rAR,AR (0) = b • rAR,AR (2) = b* rAR,AR (1) = b*b = b2 • etc.

  26. Autocorrelation Function of an Autoregressive Process of the First Order

  27. Part IV. Forecasting-Information Sources • The Conference Board publishes a monthly, Business Cycle Indicators • A monthly series followed in the popular press is the Index of Leading Indicators

  28. Leading Index Components • Average Weekly Hours, anufacturing • Initial claims For Unemployment Insurance • Manufacturerers’ New Orders • Consumer goods • Vendor Performance • Building Permits • new private housing • Manufacturerers’ New Orders • nondefense capital goods

  29. Leading Index (Cont.0 • Stock Prices • 500 common stocks • Money Supply M2 • Interest Rate Spread • 10 Treasury bonds - Federal Funds Rate • Index of Consumer Expectations

  30. Article About Leading Indices • http://www.tcb-indicators.org/GeneralInfo/bci4.pdf • BCI Web page: http://www.tcb-indicators.org/

  31. Revised Versus Old Leading Index

  32. Part V: Autoregressive of 1st Order • The model for an autoregressive process of the first order, ARONE(t) is: • ARONE(t) = b*ARONE(t-1) + WN(t) • or, using the lag operator, ARONE(t) = b*ZARONE(t) + WN(t), i.e. ARONE(t) - b*ZARONE(t) = WN(t), and factoring out ARONE(t): [1 -b*Z]*ARONE(t) = WN(t)

  33. Autoregressive of 1st Order • Dividing by [1 - b*Z], • ARONE(t) = {1/[1 - b*Z]}WN(t) where the reciprocal of [1 - bZ] is: {1/[1 - b*Z]} = (1+b*Z+b2*Z2 +b3*Z3 +…) which can be verified by multiplying [1 - b*Z] by (1+b*Z+b2*Z2 +b3*Z3 +…) to obtain 1.

  34. Autoregressive of the First Order • Note that we can write ARONE(t) as 1/[1-b*Z]*WN(t), i.e. ARONE(t) ={1/B(Z)}* WN(t), where B(Z) = [1 - b*Z] is a first order polynomial in Z, • Or, we can write ARONE(t) as ARONE(t) = A(Z) * WN(t) where A(Z) = (1+b*Z+b2*Z2 +b3*Z3 +…) is a polynomial in Z of infinite order.

  35. Autoregreesive of the First Order • So 1/B(Z) = A(Z). A first order polynomial in the denominator can approximate an infinite order polynomial in the numerator. • Box and Jenkins achieved parsimony, i.e. the use of only a few parameters which you need to estimate by modeling time series using the ratio of low order polynomials in the numerator and denominator: • ARMA(t) = {A(Z)/B(Z)}* WN(t)

  36. Autoregressive of the First Order • For Now we will concentrate on the denominator: ARONE(t) = {1/B(Z)}*WN(t), where the polynomial in the denominator, B(Z) = [1 - b*Z], captures autoregressive behavior of the first order. • Later, we will turn our attention to the numerator, where A(Z) captures moving average behavior. • Then we will combine A(Z) and B(Z).

  37. Puzzles • Annual data on output per hour; all persons, manufacturing • measure of productivity

  38. Fractional Changes: Productivity

  39. Histogram

  40. Date: 04/15/03 Time: 16:59 Sample: 1949 2000 Included observations: 51 Autocorrelation Partial Correlation AC PAC Q-Stat Prob **| . | **| . | 1 -0.237 -0.237 3.0243 0.082 . | . | .*| . | 2 -0.036 -0.097 3.0953 0.213 .*| . | .*| . | 3 -0.059 -0.098 3.2906 0.349 ***| . | ****| . | 4 -0.402 -0.481 12.605 0.013 . |**** | . |**** | 5 0.588 0.461 32.933 0.000 .*| . | .*| . | 6 -0.187 -0.132 35.033 0.000 . | . | .*| . | 7 -0.041 -0.139 35.136 0.000 . | . | .*| . | 8 -0.017 -0.125 35.155 0.000 **| . | . |** | 9 -0.221 0.205 38.310 0.000 . |*** | .*| . | 10 0.351 -0.107 46.419 0.000 . | . | . |*. | 11 -0.030 0.114 46.480 0.000

More Related