<< Chapter < Page Chapter >> Page >

Estimation of arma models

The estimation of ARMA models are relatively easy in Stata . The basic command to estimate an ARMA model is: .arima depvar [varlist], ar( numlist ) ma( numlist ) . ARIMA means A uto R egressive I ntegrated M oving A verage. See Enders (2005: 67) for a discussion of what integrated means. We can ignore it given our limited purposes. The first thing to notice in the command that this command can apply to either to a single variable or to an equation. If [varlist] is omitted, Stata will produce an estimate of the ARMA model for that variable; if the list is included, it will estimate the model with the disturbances allowed to have the ARMA structure specified in the command. Figure 14 reports the estimation of an ARMA model for real investment levels. Notice that we write AR(1/2) so that Stata knows to include both the first and second autoregressive term. A command of AR(2) would include only the second autoregressive term. In Figure 15 we report the ARMA (2, 1) estimation of (1).

Output from the estimation of an ARMA(2, 2) model of real investment.
Estimation of an ARMA(2, 2) model of real investment.

Output from the estimation of Equation (1) using an ARMA(2, 1) model.
Estimation of Equation (1) using an ARMA(2, 1) model.

Estimation of various arma models of real investment.
ARMA(1, 1) ARMA(2, 1) AR(1) AR(2) MA(1)
Intercept 185.307 185.6556 184.8208 185.2092 189.373
(10.06) (10.83) (9.27) (10.25) (18.09)
AR (L1) 0.70936 1.76342 0.80307 0.95257
(3.12) (5.27) (5.51) (4.47)
AR (L2) -0.81715 -0.18963
(-3.21) (-0.91)
MA (L1) 0.26236 -0.99998 0.87262
(0.90) (-0.00) (2.97)
Log likelihood -86.1791 -85.8702 -86.47780 -86.21224 -88.48713
Wald χ2 26.96 422.60 30.36 31.65 8.81
Probability>χ2 0.0000 0.0000 0.0000 0.0000 0.0000
Sample size 19 19 19 19 19
(14,1) 1964-1982 1964-1982 1964-1982 1964-1982 1964-1982

The interpretation of these results is not obvious. We check the sensitivity of these results by estimation some other models. The results of these estimations are reported in Table 2 and Table 3. Based purely on ML tests, it would appear that AR(1) model in Table 2 is as good as any of the models describing the ARMA structure of real investments. On the other hand, the results reported in Table 3 suggests that the ARMA(2, 1) appears to be the best model to assume for the disturbance term in the estimates of Equation (1).

Various arma estimates of equation (1).
AR(1) ARMA(1, 1) ARMA(2, 1)
Intercept -14.49489 -13.37455 -16.89182
(-0.26) (-0.23) (-1.68)
Real GNP 0.17006 0.16912 0.17253
(3.96) (3.78) (20.18)
Real interest rate -0.82517 -0.92007 -0.33692
(-0.46) (-0.33) (-0.25)
AR (L1) 0.27953 -0.02028 0.85619
(0.60) (-0.02) (1.46)
AR (L2) -0.70702
(-2.64)
MA (L1) 0.41151 -1.00000
(0.42) (-2.98)
Log likelihood -78.7868 -78.4279 -72.94569
Wald χ 2 26.30 31.86 980.18
Probability>χ 2 0.0000 0.0000 0.0000
Sample size 19 19 19
Sample period 1964-1982 1964-1982 1964-1982

Other time-series concepts

There are a large number of additional time-series methods and issues that are not discussed in this module. These topics include, among others, ARCH and GARCH estimators, unit roots, the Dickey-Fuller test, and vector autoregression (VAR) models. There is no way to do justice to these topics in notes as short as these are. Moreover, it is necessary to discuss difference equations (the discrete version of differential equations) if one wants to understand many of these topic at anything more than an intuitive level. Those interested in these topics should enroll in the forecasting course (Economics 422) or, if they cannot, plan to read several textbooks on whatever econometric tool they need to understand.

Exercise

This exercise is designed to be sure you know how to use Stata in analyzing time-series data sets; there is no economic content in the exercise. The MS Excel file Rabun County Temperature Data reports the morning temperature (MornTemp) observed in Rabun County, Georgia for every day between March 15, 2005 to November 2, 2008. The data set includes a variable “edate” that is the daily date in Stata notation. The data set also includes dummy variables for the season, the month, and the year of each observation (with the Winter, the December, and the 2008 dummy variables omitted).

a. Create a graph of (a) the data set morntemp, (2) the autocorrelations of morntemp, and (3) the partial autocorrelations of morntemp (you will have to set the matrix size to some number greater than 43 using the command .set matsize # ).

b. Estimate the following models:

  1. ARMA(2,2) for morntemp.
  2. ARMA(2,2) for morntemp as a function of the season dummy variables.
  3. ARMA(2,2) for morntemp as a function of the monthly dummy variables.
  4. ARMA(2,2) for morntemp as a function of the monthly dummy variables and the annual dummy variables.
  5. ARMA(1,2) for morntemp as a function of the monthly dummy variables and the annual dummy variables.
  6. ARMA(1,1) for morntemp as a function of the monthly dummy variables and the annual dummy variables.

c. Arrange the parameter estimates in a table and comment on them. Include the results of estimating (6) using OLS; what is the DW-statistic for this regression?

References

Cochran, D. and G. Orcutt (1949). Application of Least Squares Regression to Relationships Containing Autocorrelated Error Terms. Journal of the American Statistical Association 44 : 32-61.

Enders, Walter (1995). Applied Econometric Time Series (New York: John Wiley&Sons, Inc.).

Greene, William H. (1990). Econometric Analysis (New York: Macmillan Publishing Company).

StataCorp (2003). Stata Statistical Software: Release 8.0: Stata Time-Series Reference Manual (College Station, TX: Stat Corporation).

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Econometrics for honors students. OpenStax CNX. Jul 20, 2010 Download for free at http://cnx.org/content/col11208/1.2
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Econometrics for honors students' conversation and receive update notifications?

Ask