Risk Modelling and Survival Analysis · Core concepts of time series models
Autoregressive (AR) Models: Stationarity, Yule-Walker, ACF and PACF
Updated 11 October 2026 · Fact-checked
An AR(p) model says today's value is a linear combination of the previous p values plus white noise. It is stationary when all roots of the characteristic polynomial 1 − α₁z − … − αₚzᵖ = 0 lie outside the unit circle. Yule-Walker equations give the ACF from the coefficients. The ACF decays; the PACF cuts off after lag p.
Understand Autoregressive (AR) Models
An autoregressive model of order p, written AR(p), regresses a series on its own past. Each value X_t depends on the last p values and a fresh random shock e_t. The shock is white noise: mean 0, variance σ², uncorrelated over time.
The simplest case is AR(1): X_t = μ + α(X_{t-1} − μ) + e_t. The parameter α controls memory. If α is near 1, shocks fade slowly and the series wanders. If α is near 0, the series looks like noise. If α is negative, the series flips sign often.
For the series to be stationary, its mean, variance and autocovariances must not change with time. An AR(p) process is stationary only if all roots of its characteristic polynomial lie outside the unit circle (modulus greater than 1). For AR(1) this means |α| < 1. If α = 1, you have a random walk, which is not stationary.
The ACF (autocorrelation function) measures correlation between X_t and X_{t-k}. For an AR process it never cuts off. It decays geometrically (AR(1)) or as a mixture of decays and damped waves (higher orders). The PACF measures the correlation at lag k after removing the effect of lags 1 to k−1. For AR(p), the PACF is non-zero up to lag p and zero after it. This cut-off is how you identify the order p from data.
The Yule-Walker equations link the autocorrelations to the coefficients. Multiply the model by X_{t-k}, take expectations, and divide by the variance. You get a recursion for ρ_k. You can solve it either way: from the α values to the ACF, or from sample autocorrelations to estimates of the α values.
Key rules to remember
- AR(p) model
- X_t = μ + α₁(X_{t-1} − μ) + … + αₚ(X_{t-p} − μ) + e_t
- e_t is white noise with mean 0 and variance σ². Written with the mean μ, as in IAI notation.
- Stationarity condition
- All roots of 1 − α₁z − α₂z² − … − αₚzᵖ = 0 satisfy |z| > 1
- For AR(1): |α| < 1. Roots may be complex; use the modulus.
- AR(1) ACF
- ρ_k = αᵏ for k ≥ 0
- Decays geometrically. Alternates in sign if α < 0.
- AR(1) variance
- γ₀ = σ² ÷ (1 − α²)
- Needs |α| < 1.
- Yule-Walker recursion for ACF
- ρ_k = α₁ρ_{k-1} + α₂ρ_{k-2} + … + αₚρ_{k-p} for k ≥ 1
- Use ρ₀ = 1 and ρ_{-k} = ρ_k.
- AR(2) Yule-Walker, first two lags
- ρ₁ = α₁ + α₂ρ₁ and ρ₂ = α₁ρ₁ + α₂
- So ρ₁ = α₁ ÷ (1 − α₂).
- AR(2) stationarity triangle
- α₁ + α₂ < 1, α₂ − α₁ < 1, |α₂| < 1
- All three must hold. This is equivalent to the root condition.
- Variance from autocovariances
- γ₀ = σ² ÷ (1 − α₁ρ₁ − … − αₚρₚ)
- Gives the variance once the ρ values are known.
- PACF of AR(p)
- φ_kk ≠ 0 for k ≤ p, φ_kk = 0 for k > p; φ₁₁ = ρ₁ and φ_pp = α_p
- Cut-off at lag p identifies the order.
How to solve Autoregressive (AR) Models questions
Use this order for most exam questions on AR models. Write out each step so you earn method marks.
- 1Write the model in standard form with the mean μ and noise e_t. State that e_t is white noise with variance σ².
- 2Test stationarity. For AR(1) check |α| < 1. For AR(2) check the three triangle conditions, or find the roots of the characteristic polynomial and check their modulus.
- 3Write the Yule-Walker equations for k = 1, 2, … using ρ_k = Σ α_i ρ_{k-i}, with ρ₀ = 1 and ρ_{-k} = ρ_k.
- 4Solve the first p equations for the unknown ρ values (or for the α values if you are given sample ρ).
- 5Use the recursion to extend to higher lags.
- 6Find the variance from γ₀ = σ² ÷ (1 − Σ α_i ρ_i) if needed, then γ_k = ρ_k γ₀.
- 7For identification questions, describe the ACF (decays, never cuts off) and the PACF (cuts off after lag p), and state the order.
- 8State assumptions and give the final answer with units or a clear conclusion.
Quickest way: Fast route for AR(1) and AR(2)
When to use it: Use under time pressure when the question asks for autocorrelations, a stationarity check or a model order.
- AR(1): check |α| < 1, then write ρ_k = αᵏ and γ₀ = σ² ÷ (1 − α²). Stop.
- AR(2): compute ρ₁ = α₁ ÷ (1 − α₂) directly, then ρ₂ = α₁ρ₁ + α₂, then ρ₃ = α₁ρ₂ + α₂ρ₁.
- Check the triangle conditions instead of solving for roots.
- To estimate coefficients from sample ρ₁ and ρ₂: α₂ = (ρ₂ − ρ₁²) ÷ (1 − ρ₁²) and α₁ = ρ₁(1 − α₂).
- For identification: PACF cuts off at lag p means AR(p). Do not count the lag-0 value.
Common mistakes in Autoregressive (AR) Models
Using |α| > 1 or the wrong side of the root condition.
Students confuse the condition on the coefficients with the condition on the roots of the polynomial in z.
Fix: For AR(1), the root of 1 − αz = 0 is z = 1/α. It must lie outside the unit circle, so |α| < 1. Check which polynomial you are using.
Checking only α₁ + α₂ < 1 for AR(2).
Students remember one of the three triangle conditions and stop.
Fix: Check all three: α₁ + α₂ < 1, α₂ − α₁ < 1 and |α₂| < 1.
Saying the ACF of an AR process cuts off.
It is mixed up with the MA(q) process, where the ACF cuts off.
Fix: Remember: AR has a decaying ACF and a cut-off PACF. MA is the reverse.
Forgetting ρ_{-k} = ρ_k in the Yule-Walker equations.
For k = 1 the AR(2) equation involves ρ_{-1}, which looks unfamiliar.
Fix: Replace ρ_{-1} with ρ₁ and ρ₀ with 1 before solving.
Using σ² as the variance of X_t.
σ² is the variance of the noise e_t, and students mix the two.
Fix: Use γ₀ = σ² ÷ (1 − α²) for AR(1). The process variance is larger than σ² when α ≠ 0.
Dropping the mean μ when writing the model.
Many examples are centred at zero, so the mean is forgotten.
Fix: Write the model in terms of (X_t − μ). The mean of a stationary AR process is μ.
Worked examples
Example 1
An AR(1) process is X_t = 10 + 0.6(X_{t-1} − 10) + e_t, where e_t is white noise with variance 16. (a) Is it stationary? (b) Find the mean, variance and ρ₃.
Show the solution
- Here μ = 10 and α = 0.6.
- Stationarity: |α| = 0.6 < 1, so the process is stationary.
- Mean: E[X_t] = μ = 10.
- Variance: γ₀ = σ² ÷ (1 − α²) = 16 ÷ (1 − 0.36) = 16 ÷ 0.64 = 25.
- ACF: ρ_k = αᵏ, so ρ₃ = 0.6³ = 0.216.
Answer: Stationary; mean 10; variance 25; ρ₃ = 0.216.
Example 2
A stationary AR(2) process (X_t − μ) = 0.5(X_{t-1} − μ) + 0.3(X_{t-2} − μ) + e_t has white noise variance σ² = 9. Find ρ₁, ρ₂, ρ₃ and the variance of X_t.
Show the solution
- Check stationarity: α₁ + α₂ = 0.8 < 1; α₂ − α₁ = −0.2 < 1; |α₂| = 0.3 < 1. Stationary.
- Yule-Walker at k = 1: ρ₁ = α₁ + α₂ρ₁, so ρ₁(1 − 0.3) = 0.5 and ρ₁ = 0.5 ÷ 0.7 = 5/7 ≈ 0.7143.
- At k = 2: ρ₂ = α₁ρ₁ + α₂ = 0.5 × 5/7 + 0.3 = 0.35714 + 0.3 = 0.65714.
- At k = 3: ρ₃ = α₁ρ₂ + α₂ρ₁ = 0.5 × 0.65714 + 0.3 × 0.71429 = 0.32857 + 0.21429 = 0.54286.
- Variance: γ₀ = σ² ÷ (1 − α₁ρ₁ − α₂ρ₂) = 9 ÷ (1 − 0.35714 − 0.19714) = 9 ÷ 0.44571.
- 9 ÷ 0.44571 ≈ 20.19.
Answer: ρ₁ ≈ 0.714, ρ₂ ≈ 0.657, ρ₃ ≈ 0.543; Var(X_t) ≈ 20.19.
Exam tips
- Always state the stationarity check explicitly. Examiners award a mark for it even when the answer is obvious.
- In identification questions, describe both ACF and PACF. Saying only 'PACF cuts off at lag 2' may lose the mark for the ACF behaviour.
- In computer-based papers, show the R call (for example arima or ar with the order) and quote the fitted coefficients, then interpret them. Check whether the output reports the mean or the intercept, as these differ.
- For AR(2) parameter estimation from sample ρ values, write the Yule-Walker equations in full before solving, so method marks are safe.
- Keep four decimals in intermediate steps. Rounding ρ₁ early shifts ρ₂ and the variance.
Practice questions from Core concepts of time series models
- A stationary AR(1) process has X_t - mu = 0.5 (X_{t-1} - mu) + e_t with mu = 40. Given X_10 = 48, what is the best forecast of X_12 at time …
- A time series {e_t} is described as Gaussian white noise with mean zero and variance sigma squared. Which statement about this process is co…
- For a sample from an AR(2) process, which pattern of sample autocorrelation function (ACF) and partial autocorrelation function (PACF) is mo…
- Which feature of the sample partial autocorrelation function (PACF) would suggest that an AR(2) model is appropriate for a stationary series…
- X_t is a random walk with X_0 = 0 and shock variance sigma squared. For s < t, what is Corr(X_s, X_t) when s = 9 and t = 36?
Autoregressive (AR) Models in other exams
The same ground in other exams, if you are preparing for more than one or want another angle on it.
Autoregressive (AR) Models: frequently asked questions
What is the stationarity condition for an AR(1) process?
The process X_t = μ + α(X_{t-1} − μ) + e_t is stationary when |α| < 1. This is the same as saying the root of 1 − αz = 0 lies outside the unit circle. If α = 1 it becomes a random walk and is not stationary.
How do I solve the Yule-Walker equations for AR(2)?
Write ρ₁ = α₁ + α₂ρ₁ and ρ₂ = α₁ρ₁ + α₂. Solve the first for ρ₁ = α₁ ÷ (1 − α₂), then substitute to get ρ₂. Higher lags follow from ρ_k = α₁ρ_{k-1} + α₂ρ_{k-2}.
What do the ACF and PACF of an AR(p) process look like?
The ACF decays towards zero, either geometrically or with damped oscillations, and does not cut off. The PACF is non-zero up to lag p and is zero afterwards. The lag where the PACF cuts off gives the order.
How is an AR model different from an MA model?
An AR model uses past values of the series. An MA model uses past noise terms. Their ACF and PACF patterns are reversed: AR has a cut-off PACF, MA has a cut-off ACF.