Tutorials › Exam SRM - Statistics for Risk Modeling › Using ACF and PACF to Screen Candidate Orders

Exam SRM · Time-Series Model Selection · Tutorial 512 of 926

Using ACF and PACF to Screen Candidate Orders

Use ACF and PACF patterns, seasonal lags, and preliminary residuals to shortlist plausible time-series orders without mistaking sample cutoffs for certainty.

Advanced 13 min read

What You'll Learn

  • Define sample ACF and PACF and distinguish what each measures.
  • Recognize the ideal ACF and PACF patterns for stationary AR and invertible MA models.
  • Use approximate white-noise reference bands while accounting for their limitations.
  • Read seasonal dependence at seasonal lags and distinguish it from ordinary short-lag structure.
  • Use residual patterns to shortlist candidate orders without treating them as proof of a model.
  • Avoid common errors caused by trends, multiple lag comparisons, and over-interpreting sample cutoffs.

Why ACF and PACF Patterns Help Screen Orders

When a time series has been made suitable for stationary modeling, its dependence across time can help narrow the plausible model orders. The autocorrelation function (ACF) summarizes correlation between observations separated by each lag. The partial autocorrelation function (PACF) measures the correlation at a lag after accounting for the intervening lags. Their patterns can suggest AR or MA structure, but sample plots are screening tools: finite samples fluctuate, and real series need not display ideal theoretical cutoffs.

The screening should be done on a series whose trend, level changes, and seasonality have been considered. A persistent trend can make many ACF values large even when the short-memory order is small. If a preliminary trend or seasonal component has been removed, the ACF and PACF of the resulting residuals can help expose remaining dependence. They do not, by themselves, establish that the preliminary treatment or a proposed model is adequate.

Definition: For a stationary series \(X_t\) with mean \(\mu\), the theoretical autocorrelation at lag \(k\) is \(\rho_k=\operatorname{Corr}(X_t,X_{t-k})\). The partial autocorrelation at lag \(k\), denoted \(\phi_{kk}\), is the correlation between \(X_t\) and \(X_{t-k}\) after each has been linearly adjusted for \(X_{t-1},\ldots,X_{t-k+1}\). The sample ACF \(r_k\) and sample PACF \(\hat{\phi}_{kk}\) estimate these quantities from observed data.

A common sample ACF definition, using observations \(x_1,\ldots,x_n\) and their sample mean \(\bar{x}\), is

$$ r_k= \frac{\displaystyle\sum_{t=k+1}^{n}(x_t-\bar{x})(x_{t-k}-\bar{x})} {\displaystyle\sum_{t=1}^{n}(x_t-\bar{x})^2}, \qquad k=0,1,\ldots $$

In particular, \(r_0=1\). Software can use slightly different finite-sample conventions, so exact plotted values may differ a little. The PACF at lag \(k\) is commonly estimated by fitting the linear prediction of \(X_t\) from its preceding \(k\) values and reading the coefficient on \(X_{t-k}\). At lag 1, the PACF equals the ACF.

Theoretical Patterns and Their Limits

For a stationary autoregressive model of order \(p\), written AR(\(p\)), the theoretical PACF is zero beyond lag \(p\), while the ACF generally tails off. For an invertible moving-average model of order \(q\), written MA(\(q\)), the ACF is zero beyond lag \(q\), while the PACF generally tails off. A stationary, invertible ARMA model usually has both functions tail off rather than showing a clean cutoff. “Tails off” may mean gradual decay, possibly with alternating signs or damped oscillation; it does not require a smooth monotone decline.

These are population patterns under the stated model conditions, not rules that every sample plot must obey. A finite-sample ACF may have a nonzero value past a theoretical cutoff, and a true nonzero correlation may look small in a short series. If the process is nonstationary, these stationary-model patterns are not a sound basis for selecting orders until the nonstationarity has been addressed.

Screening guide: ACF tails off and PACF appears to cut off after \(p\): consider AR(\(p\)). ACF appears to cut off after \(q\) and PACF tails off: consider MA(\(q\)). Both appear to tail off: consider ARMA candidates. Treat each pattern as evidence for a shortlist, not as a definitive identification.

One reason for the MA cutoff follows directly from the model. Suppose \(X_t=\varepsilon_t+\theta\varepsilon_{t-1}\), where the innovations \(\varepsilon_t\) are uncorrelated, have mean zero, and have variance \(\sigma^2\). Then

$$ \gamma_0=\operatorname{Var}(X_t)=(1+\theta^2)\sigma^2, \qquad \gamma_1=\operatorname{Cov}(X_t,X_{t-1})=\theta\sigma^2. $$

For \(k\geq 2\), \(X_t\) and \(X_{t-k}\) contain no common innovation, so their covariance is zero. Dividing each covariance by \(\gamma_0\) gives

$$ \rho_1=\frac{\theta}{1+\theta^2}, \qquad \rho_k=0 \quad (k\geq 2). $$

That is the theoretical ACF cutoff for an MA(1). In observed data, the estimated \(r_k\) values beyond lag 1 will generally not be exactly zero. By contrast, for a stationary AR(1), the ACF follows \(\rho_k=\phi^k\), where \(|\phi|<1\), so it decays rather than cutting off; its PACF cuts off after lag 1.

Reference Bands: A Useful Approximation, Not a Verdict

A standard visual reference for an ACF under a white-noise null is the approximate pointwise band

$$ -\frac{1.96}{\sqrt{n}} \ \text{to}\ \frac{1.96}{\sqrt{n}}. $$

Under white noise and suitable large-sample conditions, an individual sample autocorrelation at a fixed nonzero lag has approximate standard error \(1/\sqrt{n}\). The factor 1.96 gives an approximate 95% pointwise reference interval. The same rough band is often used when scanning PACFs, although more careful uncertainty calculations can depend on the setting and estimation method.

These bands are not simultaneous confidence bands for all lags. When many lags are inspected, some estimates may fall outside by chance. Also, if the series has meaningful autocorrelation, the white-noise approximation need not describe the uncertainty at later lags. Use the band to judge whether an isolated value is striking relative to sampling noise, not as a mechanical keep-or-reject test.

Worked Example: A Decaying ACF and a PACF Cutoff

A stationary, seasonally adjusted series has \(n=64\) observations. Its sample ACF at lags 1 through 6 is \(0.62, 0.38, 0.23, 0.14, 0.08, 0.04\). Its sample PACF at those lags is \(0.62, 0.02, -0.05, 0.03, -0.01, 0.04\). Use the reference band to shortlist orders.

Method. Calculate the approximate white-noise band, compare each displayed value with it, and then assess the overall ACF/PACF shapes. The method avoids treating a single plotted bar as an exact order test.

Step 1: Calculate the band. Here \(\sqrt{64}=8\), so \(1.96/\sqrt{64}=1.96/8=0.245\). The approximate band is \(-0.245\) to \(0.245\).

Step 2: Read the ACF. The lag-1 and lag-2 values, \(0.62\) and \(0.38\), exceed \(0.245\). The lag-3 value, \(0.23\), is inside the band, as are lags 4 through 6. The ACF declines across the lags rather than stopping sharply at lag 1.

Step 3: Read the PACF. The lag-1 value \(0.62\) exceeds the band. Each displayed PACF value at lags 2 through 6 lies within it. This is consistent with a PACF that cuts off after lag 1.

Conclusion and check. A stationary AR(1) is a reasonable first candidate because its ACF tails off and its PACF appears to stop after lag 1. The values do not prove the process is AR(1); for example, a mixed model or sampling variation could create a similar plot. The calculation is internally consistent: the stated threshold is \(0.245\), and only the stated ACF and PACF values exceed it.

Seasonal Lags and Preliminary Residuals

For data observed at a regular seasonal interval \(s\), inspect not just the first few lags but also \(s,2s,3s,\ldots\). For monthly data with annual seasonality, these are lags 12, 24, 36, and so on. For quarterly data with annual seasonality, they are lags 4, 8, 12, and so on. Repeated dependence at those lags can indicate that seasonal structure remains. A seasonal AR pattern may show a PACF feature at a seasonal lag and a declining ACF at seasonal multiples; a seasonal MA pattern may instead show an ACF feature at a seasonal lag. Mixed or short samples can blur these ideal patterns.

A residual ACF is useful only with context: identify which trend, seasonal, or other preliminary structure was removed before calculating it. For example, residual autocorrelation at lags 12 and 24 after removing a monthly series’ trend may signal that seasonal dependence remains. Conversely, residual values near zero at those lags are not proof of adequate seasonal modeling. Later model comparisons and residual diagnostics are needed to assess a fitted candidate.

Worked Example: Screening a Seasonal Pattern in Residuals

For a monthly series with 144 observations, residuals from a preliminary trend and nonseasonal adjustment have ACF values \(r_{12}=0.58\), \(r_{24}=0.35\), \(r_{36}=0.20\), and \(r_{48}=0.04\). Their seasonal-lag PACF values are \(\hat{\phi}_{12,12}=0.60\), \(\hat{\phi}_{24,24}=-0.08\), and \(\hat{\phi}_{36,36}=0.04\). The ordinary-lag residual ACF values at lags 1 through 4 are \(0.12,-0.08,0.05,0.10\). Shortlist seasonal structure.

Method. Compare seasonal-lag values with the approximate band, then compare the pattern in the ACF and PACF at seasonal multiples. The method avoids confusing a remaining seasonal signal with ordinary short-lag dependence.

Step 1: Find the band. Since \(\sqrt{144}=12\), the approximate band is \(\pm 1.96/12=\pm 0.1633\), rounded to four decimals.

Step 2: Check seasonal ACF values. At lags 12, 24, and 36, \(0.58\), \(0.35\), and \(0.20\) all exceed \(0.1633\). At lag 48, \(0.04\) is inside the band. The seasonal ACF appears to decay over successive annual lags.

Step 3: Check seasonal PACF and ordinary lags. The seasonal PACF at lag 12, \(0.60\), exceeds the band; the values at lags 24 and 36, \(-0.08\) and \(0.04\), are inside it. Each stated ordinary-lag ACF value is also inside the band.

Conclusion and check. A seasonal AR(1) component with period 12 is a reasonable candidate to include in a shortlist: the seasonal ACF declines across multiples of 12 while the seasonal PACF has its main feature at 12. This is not a unique identification. The band calculation is \(1.96/12=0.1633\), so the comparisons above use the same threshold consistently.

Common Trap: Calling a Sample Cutoff Exact

Suppose an analyst sees a PACF bar outside a reference band at lag 1 and no later bars outside it. The incorrect conclusion is, “The series is definitely AR(1), so every other order is ruled out.” The correction is, “The pattern supports AR(1) as a candidate, subject to model comparison and further checking.” A sample PACF can hide smaller true effects, and an ARMA process can resemble a pure AR process over a limited range of lags.

A related calculation error is to change the denominator of the common white-noise band at every lag. For \(n=64\), an incorrect lag-6 band would use \(1.96/\sqrt{64-6}=1.96/\sqrt{58}\approx0.257\). The usual rough white-noise pointwise band is instead \(1.96/\sqrt{64}=0.245\) at each fixed lag. More importantly, neither version makes a visual cutoff conclusive; the band is approximate and pointwise.

Other practical cautions are to avoid reading a trending raw series as a short-memory AR model, to inspect seasonal multiples as well as small lags, and to avoid selecting a high order merely because one late-lag bar crosses a band. ACF and PACF plots are most useful when they narrow a small, defensible set of candidates.

1
Establish the series being screened.
Check whether trend, changing level, and seasonal structure require attention before applying stationary-model patterns.
2
Read ordinary and seasonal lags.
Compare ACF and PACF shapes at small lags and at seasonal multiples such as \(s,2s,3s\).
3
Form a shortlist, not a verdict.
Use apparent cutoffs and tails to propose a few AR, MA, ARMA, or seasonal candidates; retain plausible alternatives when the pattern is ambiguous.
4
Carry candidates forward for comparison.
Visual screening does not estimate parameters or establish adequacy. Candidate models must be assessed using later model-comparison and diagnostic steps.
Key takeaway: For stationary series, a PACF cutoff with a tailing ACF suggests an AR order, while an ACF cutoff with a tailing PACF suggests an MA order. Seasonal multiples and residual patterns add useful clues. All are provisional evidence, not exact sample rules.

Check Your Understanding

For each problem, use the approximate band \( \pm 1.96/\sqrt{n} \) when a sample size is given. Treat the result as a screening conclusion.

  1. A stationary series has \(n=49\), ACF values \(0.31,0.10,-0.07\) at lags 1–3, and PACF values \(0.31,0.02,-0.04\). What is the band and a reasonable first candidate? (Answer: The band is approximately \(\pm0.2800\); AR(1) is a reasonable first candidate.)
  2. For \(n=100\), the ACF at lags 1–4 is \(0.48,0.05,-0.02,0.03\), while the PACF is \(0.48,-0.25,0.14,-0.08\). What pattern is suggested? (Answer: The band is \(\pm0.1960\); the ACF appears to cut off after lag 1 while the PACF tails off, suggesting MA(1) as a candidate.)
  3. A monthly residual series has \(n=121\), with ACF values \(0.52,0.30,0.18\) at lags 12, 24, and 36, and PACF values \(0.54,0.03,-0.02\) at those lags. What seasonal candidate is suggested? (Answer: The band is approximately \(\pm0.1782\); a seasonal AR(1) candidate with period 12 is suggested.)
  4. A raw series has a slowly declining ACF beginning at \(0.91\), but the series visibly rises over time. Should that pattern alone determine a stationary AR order? (Answer: No; address the trend or other nonstationarity first, then screen the appropriate transformed or residual series.)
  5. For a quarterly residual series, the ACF has notable values at lags 1 and 4, with another notable value at lag 8. What should be included in the shortlist? (Answer: Consider both ordinary-lag and seasonal structure with period 4; the repeated lag-4 and lag-8 signal may support a seasonal candidate, but the plots do not uniquely identify its order.)
Clinical Biostats