<< 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).

Questions & Answers

if three forces F1.f2 .f3 act at a point on a Cartesian plane in the daigram .....so if the question says write down the x and y components ..... I really don't understand
Syamthanda Reply
hey , can you please explain oxidation reaction & redox ?
Boitumelo Reply
hey , can you please explain oxidation reaction and redox ?
Boitumelo
for grade 12 or grade 11?
Sibulele
the value of V1 and V2
Tumelo Reply
advantages of electrons in a circuit
Rethabile Reply
we're do you find electromagnetism past papers
Ntombifuthi
what a normal force
Tholulwazi Reply
it is the force or component of the force that the surface exert on an object incontact with it and which acts perpendicular to the surface
Sihle
what is physics?
Petrus Reply
what is the half reaction of Potassium and chlorine
Anna Reply
how to calculate coefficient of static friction
Lisa Reply
how to calculate static friction
Lisa
How to calculate a current
Tumelo
how to calculate the magnitude of horizontal component of the applied force
Mogano
How to calculate force
Monambi
a structure of a thermocouple used to measure inner temperature
Anna Reply
a fixed gas of a mass is held at standard pressure temperature of 15 degrees Celsius .Calculate the temperature of the gas in Celsius if the pressure is changed to 2×10 to the power 4
Amahle Reply
How is energy being used in bonding?
Raymond Reply
what is acceleration
Syamthanda Reply
a rate of change in velocity of an object whith respect to time
Khuthadzo
how can we find the moment of torque of a circular object
Kidist
Acceleration is a rate of change in velocity.
Justice
t =r×f
Khuthadzo
how to calculate tension by substitution
Precious Reply
hi
Shongi
hi
Leago
use fnet method. how many obects are being calculated ?
Khuthadzo
khuthadzo hii
Hulisani
how to calculate acceleration and tension force
Lungile Reply
you use Fnet equals ma , newtoms second law formula
Masego
please help me with vectors in two dimensions
Mulaudzi Reply
how to calculate normal force
Mulaudzi
Got questions? Join the online conversation and get instant answers!
Jobilize.com Reply

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