ARIMA Procedure

The SCAN Method

The smallest canonical (SCAN) correlation method can tentatively identify the orders of a stationary or nonstationary ARMA process. Tsay and Tiao (1985) proposed the technique, and for useful descriptions of the algorithm, see Box, Jenkins, and Reinsel (1994); Choi (1992).

Given a stationary or nonstationary time series StartSet z Subscript t Baseline colon 1 less-than-or-equal-to t less-than-or-equal-to n EndSet with mean corrected form z overTilde Subscript t Baseline equals z Subscript t Baseline minus mu Subscript z with a true autoregressive order of p plus d and with a true moving-average order of q, you can use the SCAN method to analyze eigenvalues of the correlation matrix of the ARMA process. The following paragraphs provide a brief description of the algorithm.

For autoregressive test order m equals p Subscript m i n Baseline comma ellipsis comma p Subscript m a x Baseline and for moving-average test order j equals q Subscript m i n Baseline comma ellipsis comma q Subscript m a x Baseline, perform the following steps:

  1. Let upper Y Subscript m comma t Baseline equals left-parenthesis z overTilde Subscript t Baseline comma z overTilde Subscript t minus 1 Baseline comma ellipsis comma z overTilde Subscript t minus m Baseline right-parenthesis prime. Compute the following left-parenthesis m plus 1 right-parenthesis times left-parenthesis m plus 1 right-parenthesis matrix,

    StartLayout 1st Row 1st Column ModifyingAbove beta With caret left-parenthesis m comma j plus 1 right-parenthesis 2nd Column equals 3rd Column left-parenthesis sigma-summation Underscript t Endscripts upper Y Subscript m comma t minus j minus 1 Baseline upper Y prime Subscript m comma t minus j minus 1 right-parenthesis Superscript negative 1 Baseline left-parenthesis sigma-summation Underscript t Endscripts upper Y Subscript m comma t minus j minus 1 Baseline upper Y prime Subscript m comma t right-parenthesis 2nd Row 1st Column ModifyingAbove beta With caret Superscript asterisk Baseline left-parenthesis m comma j plus 1 right-parenthesis 2nd Column equals 3rd Column left-parenthesis sigma-summation Underscript t Endscripts upper Y Subscript m comma t Baseline upper Y prime Subscript m comma t right-parenthesis Superscript negative 1 Baseline left-parenthesis sigma-summation Underscript t Endscripts upper Y Subscript m comma t Baseline upper Y prime Subscript m comma t minus j minus 1 right-parenthesis 3rd Row 1st Column ModifyingAbove upper A With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis 2nd Column equals 3rd Column ModifyingAbove beta With caret Superscript asterisk Baseline left-parenthesis m comma j plus 1 right-parenthesis ModifyingAbove beta With caret left-parenthesis m comma j plus 1 right-parenthesis EndLayout

    where t ranges from j plus m plus 2 to n.

  2. Find the smallest eigenvalue, ModifyingAbove lamda With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis, of ModifyingAbove upper A With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis and its corresponding normalized eigenvector, normal upper Phi Subscript m comma j Baseline equals left-parenthesis 1 comma minus phi 1 Superscript left-parenthesis m comma j right-parenthesis Baseline comma minus phi 2 Superscript left-parenthesis m comma j right-parenthesis Baseline comma ellipsis comma minus phi Subscript m Superscript left-parenthesis m comma j right-parenthesis Baseline right-parenthesis. The squared canonical correlation estimate is ModifyingAbove lamda With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis.

  3. Using the normal upper Phi Subscript m comma j as AR(m) coefficients, obtain the residuals for t equals j plus m plus 1 to n, by following the formula: w Subscript t Superscript left-parenthesis m comma j right-parenthesis Baseline equals z overTilde Subscript t Baseline minus phi 1 Superscript left-parenthesis m comma j right-parenthesis Baseline z overTilde Subscript t minus 1 Baseline minus phi 2 Superscript left-parenthesis m comma j right-parenthesis Baseline z overTilde Subscript t minus 2 Baseline minus midline-horizontal-ellipsis minus phi Subscript m Superscript left-parenthesis m comma j right-parenthesis Baseline z overTilde Subscript t minus m.

  4. From the sample autocorrelations of the residuals, r Subscript k Baseline left-parenthesis w right-parenthesis, approximate the standard error of the squared canonical correlation estimate by

    v a r left-parenthesis ModifyingAbove lamda With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis Superscript 1 slash 2 Baseline right-parenthesis almost-equals d left-parenthesis m comma j right-parenthesis slash left-parenthesis n minus m minus j right-parenthesis

    where d left-parenthesis m comma j right-parenthesis equals left-parenthesis 1 plus 2 sigma-summation Underscript i equals 1 Overscript j minus 1 Endscripts r Subscript k Baseline left-parenthesis w Superscript left-parenthesis m comma j right-parenthesis Baseline right-parenthesis right-parenthesis.

The test statistic to be used as an identification criterion is

c left-parenthesis m comma j right-parenthesis equals minus left-parenthesis n minus m minus j right-parenthesis normal l normal n left-parenthesis 1 minus ModifyingAbove lamda With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis slash d left-parenthesis m comma j right-parenthesis right-parenthesis

which is asymptotically chi 1 squared if m equals p plus d and j greater-than-or-equal-to q or if m greater-than-or-equal-to p plus d and j equals q. For m greater-than p and j less-than q, there is more than one theoretical zero canonical correlation between upper Y Subscript m comma t and upper Y Subscript m comma t minus j minus 1. Since the ModifyingAbove lamda With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis are the smallest canonical correlations for each left-parenthesis m comma j right-parenthesis, the percentiles of c left-parenthesis m comma j right-parenthesis are less than those of a chi 1 squared; therefore, Tsay and Tiao (1985) state that it is safe to assume a chi 1 squared. For m less-than p and j less-than q, no conclusions about the distribution of c left-parenthesis m comma j right-parenthesis are made.

A SCAN table is then constructed using c left-parenthesis m comma j right-parenthesis to determine which of the ModifyingAbove lamda With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis are significantly different from zero (see Table 7). The ARMA orders are tentatively identified by finding a (maximal) rectangular pattern in which the ModifyingAbove lamda With caret Superscript asterisk Baseline left-parenthesis m comma j right-parenthesis are insignificant for all test orders m greater-than-or-equal-to p plus d and j greater-than-or-equal-to q. There might be more than one pair of values (p plus d comma q) that permit such a rectangular pattern. In this case, parsimony and the number of insignificant items in the rectangular pattern should help determine the model order. Table 8 depicts the theoretical pattern associated with an ARMA(2,2) series.

Table 7: SCAN Table

MA
AR 0 1 2 3 dot dot
0 c left-parenthesis 0 comma 0 right-parenthesis c left-parenthesis 0 comma 1 right-parenthesis c left-parenthesis 0 comma 2 right-parenthesis c left-parenthesis 0 comma 3 right-parenthesis dot dot
1 c left-parenthesis 1 comma 0 right-parenthesis c left-parenthesis 1 comma 1 right-parenthesis c left-parenthesis 1 comma 2 right-parenthesis c left-parenthesis 1 comma 3 right-parenthesis dot dot
2 c left-parenthesis 2 comma 0 right-parenthesis c left-parenthesis 2 comma 1 right-parenthesis c left-parenthesis 2 comma 2 right-parenthesis c left-parenthesis 2 comma 3 right-parenthesis dot dot
3 c left-parenthesis 3 comma 0 right-parenthesis c left-parenthesis 3 comma 1 right-parenthesis c left-parenthesis 3 comma 2 right-parenthesis c left-parenthesis 3 comma 3 right-parenthesis dot dot
dot dot dot dot dot dot dot
dot dot dot dot dot dot dot


Table 8: Theoretical SCAN Table for an ARMA(2,2) Series

MA
AR 0 1 2 3 4 5 6 7
0 * X X X X X X X
1 * X X X X X X X
2 * X 0 0 0 0 0 0
3 * X 0 0 0 0 0 0
4 * X 0 0 0 0 0 0
X = significant terms
0 = insignificant terms
* = no pattern


Last updated: June 19, 2025