ARMALIK Call

Chapter Contents

Language Reference

ARMALIK Call

computes the log likelihood and residuals for an ARMA model

CALL ARMALIK( lnl, resid, std, x, phi, theta);

The inputs to the ARMALIK subroutine are as follows:

x: is an n ×1 or 1 ×n matrix containing values of the time series (assuming mean zero).
phi: is a 1 ×(p+1) matrix containing the autoregressive parameter values. The first element is assumed to have the value 1.
theta: is a 1 ×(q+1) matrix containing the moving-average parameter values. The first element is assumed to have the value 1.

The ARMALIK subroutine returns the following values:

lnl: specifies a 3 ×1 matrix containing the log likelihood concentrated with respect to the innovation variance; the estimate of the innovation variance (the unconditional sum of squares divided by n); and the log of the determinant of the variance matrix, which is standardized to unit variance for the innovations.
resid: specifies an n ×1 matrix containing the standardized residuals. These values are uncorrelated with a constant variance if the specified ARMA model is the correct one.
std: specifies an n ×1 matrix containing the scale factors used to standardize the residuals. The actual residuals from the one-step-ahead predictions using the past values can be computed as std#resid.

The ARMALIK subroutine computes the concentrated log likelihood function for an ARMA model. The unconditional sum of squares is readily available, as are the one-step-ahead prediction residuals. Factors that can be used to generate confidence limits associated with prediction from a finite past sample are also returned.

The notational conventions for the ARMALIK subroutine are the same as those used by the ARMACOV subroutine. See the description of the ARMACOV call for the model employed. In addition, the condition $\sum_{i=0}^q \theta_{iz}^i \ne 0$ for |z|<1 should be satisfied to guard against floating-point overflow.

If the column vector x contains n values of a time series and the variance matrix is denoted $\Sigma = \sigma^2 V$ , where $\sigma^2$ is the variance of the innovations, then, up to additive constants, the log likelihood, concentrated with respect to $\sigma^2$ , is

$-\frac{n}2 \log ( x^' V^{-1} x ) - \frac{1}2 \log |{V}| .$

The matrix V is a function of the specified ARMA model parameters. If L is the lower Cholesky root of V (that is, V = LL'), then the standardized residuals are computed as resid= L^-1x. The elements of std are the diagonal elements of L. The variance estimate is x'V^-1x/n, and the log determinant is $\log|{V}|$ .See Ansley (1979) for further details. To compute the log-likelihood for the model

y_t - y_t-1 + 0.25y_t-2 = e_t + 0.5e_t-1

use the following IML code:

   proc iml;
      phi={ 1 -1 0.25} ;
      theta={ 1 0.5} ;
      x={ 1 2 3 4 5} ;
      call armalik(lnl,resid,std,x,phi,theta);
      print lnl resid std;

The printed output is

                           LNL     RESID       STD
                     -0.822608 0.4057513 2.4645637
                     0.8721154 0.9198158 1.2330147
                     2.3293833 0.8417343 1.0419028
                               1.0854175 1.0098042
                               1.2096421 1.0024125

Chapter Contents
Previous
Next
Top