Example 7.4 An Intervention Model for Ozone Data

This example fits an intervention model to ozone data as suggested by Box and Tiao (1975). Notice that the response variable, OZONE, and the innovation, X1, are seasonally differenced. The final model for the differenced data is a multiple regression model with a moving-average structure assumed for the residuals.

The model is fit by maximum likelihood. The seasonal moving-average parameter and its standard error are fairly sensitive to which method is chosen to fit the model, in agreement with the observations of Davidson (1981) and Ansley and Newbold (1980); thus, fitting the model by the unconditional or conditional least squares method produces somewhat different estimates for these parameters.

Some missing values are appended to the end of the input data to generate additional values for the independent variables. Since the independent variables are not modeled, values for them must be available for any times at which predicted values are desired. In this case, predicted values are requested for 12 periods beyond the end of the data. Thus, values for X1, WINTER, and SUMMER must be given for 12 periods ahead.

The following statements read in the data and compute dummy variables for use as intervention inputs:

   title1 'Intervention Data for Ozone Concentration';
   title2 '(Box and Tiao, JASA 1975 P.70)';
   data air;
      input ozone @@;
      label ozone  = 'Ozone Concentration'
            x1     = 'Intervention for post 1960 period'
            summer = 'Summer Months Intervention'
            winter = 'Winter Months Intervention';
      date = intnx( 'month', '31dec1954'd, _n_ );
      format date monyy.;
      month = month( date );
      year = year( date );
      x1 = year >= 1960;
      summer = ( 5 < month < 11 ) * ( year > 1965 );
      winter = ( year > 1965 ) - summer;
   datalines;
   2.7  2.0  3.6  5.0  6.5  6.1  5.9  5.0  6.4  7.4  8.2  3.9
   4.1  4.5  5.5  3.8  4.8  5.6  6.3  5.9  8.7  5.3  5.7  5.7
   3.0  3.4  4.9  4.5  4.0  5.7  6.3  7.1  8.0  5.2  5.0  4.7
   
   ... more lines ...   

The following statements produce Output 7.4.1 through Output 7.4.3:

   proc arima data=air;
   
      /* Identify and seasonally difference ozone series */
      identify var=ozone(12)
               crosscorr=( x1(12) summer winter ) noprint;
   
      /* Fit a multiple regression with a seasonal MA model */
      /*     by the maximum likelihood method               */
      estimate q=(1)(12) input=( x1 summer winter )
               noconstant method=ml;
   
      /* Forecast */
      forecast  lead=12 id=date interval=month;
   
   run;

The ESTIMATE statement results are shown in Output 7.4.1 and Output 7.4.2.