Tutorials › Exam SRM - Statistics for Risk Modeling › Solving ARMA(1,1) Parameter and ACF Problems

Exam SRM · ARMA, ARIMA, and Seasonality · Tutorial 457 of 926

Solving ARMA(1,1) Parameter and ACF Problems

Learn to recover ARMA(1,1) parameters from autocorrelations, check the model conditions, and calculate autocovariances and selected ACF values.

Advanced 12 min read

What You'll Learn

  • State the ARMA(1,1) model and distinguish its mean, parameters, and innovation variance.
  • Check stationarity and invertibility using the model’s lag-polynomial roots.
  • Derive the lag-zero and positive-lag autocovariances from the causal representation.
  • Calculate the ARMA(1,1) ACF and use its lag pattern to recover parameters.
  • Resolve multiple parameter solutions by applying the invertibility condition.
  • Avoid common errors involving parameter signs, ACF ratios, and invalid model conditions.

Solving ARMA(1,1) Problems

An ARMA(1,1) process combines dependence on the previous observation with dependence on the current and previous innovations. Problems often ask for more than substitution into a formula: a candidate may need to recover parameters from autocorrelations, determine whether the model is stationary and invertible, or identify which parameter solution is admissible. The order of those steps matters. In particular, do not use a stationary ACF formula before checking the stationarity condition.

Throughout this tutorial, use the following sign convention. Let \(X_t\) be a process with mean \(\mu\), and let \(Z_t\) be white noise with mean zero and variance \(\sigma_Z^2\). The model is

$$ X_t-\mu=\phi(X_{t-1}-\mu)+Z_t+\theta Z_{t-1}. $$

White noise here means that the \(Z_t\) values have constant variance and are uncorrelated at distinct times. They need not be normally distributed for the second-moment results below. The sign convention is important: some sources write the moving-average term with a minus sign and therefore report a different signed MA parameter.

Definition: In the stated ARMA(1,1) convention, \(\phi\) is the autoregressive parameter, \(\theta\) is the moving-average parameter, and \(\sigma_Z^2\) is the innovation variance. A stationary process has constant mean and autocovariances that depend only on the lag. The model is causal when \(|\phi|<1\), and it is invertible when \(|\theta|<1\). These are distinct conditions.

Conditions Before Calculating the ACF

The autoregressive lag polynomial is \(1-\phi B\), where \(B X_t=X_{t-1}\). Its root is \(1/\phi\) when \(\phi\ne 0\). Causality requires that this root lie outside the unit circle, which gives \(|\phi|<1\). The moving-average polynomial is \(1+\theta B\), with root \(-1/\theta\) when \(\theta\ne 0\). Invertibility requires that root to lie outside the unit circle, giving \(|\theta|<1\). When a parameter is zero, its polynomial has no finite root to test and its corresponding condition is satisfied.

For the usual stationary ARMA(1,1) autocovariance calculation, first establish \(|\phi|<1\). Invertibility is not needed merely to calculate autocovariances, but it is important when identifying a unique, standard MA representation from the model. In parameter-recovery problems, use it to reject a noninvertible alternative when the question specifies an invertible model.

Worked Example: Recovering the MA Parameter From a Lag-One ACF

A stationary, invertible ARMA(1,1) process uses the convention above, has \(\phi=0.5\), and has theoretical lag-one autocorrelation \(\rho_1=9/14\). Find \(\theta\). State the method and the common mistake it avoids.

Method. Substitute the given \(\phi\) and \(\rho_1\) into the ARMA(1,1) formula for \(\rho_1\), solve the resulting quadratic, and apply invertibility to select the admissible solution. The innovations are assumed to be white noise with some variance \(\sigma_Z^2>0\); its value cancels from the ACF.

The formula derived below is \(\rho_1=(\phi+\theta)(1+\phi\theta)/(1+\theta^2+2\phi\theta)\). Substitution gives

$$ \frac{9}{14} = \frac{(1/2+\theta)(1+\theta/2)} {1+\theta^2+\theta}. $$

Multiplying by the denominator and expanding both sides,

$$ \frac{9}{14}+\frac{9}{14}\theta+\frac{9}{14}\theta^2 = \frac{1}{2}+\frac{5}{4}\theta+\frac{1}{2}\theta^2. $$

Moving the right-hand side to the left gives \(\frac{1}{7}-\frac{17}{28}\theta+\frac{1}{7}\theta^2=0\). Multiplying by 28 yields

$$ 4\theta^2-17\theta+4=0 \quad\Longrightarrow\quad \theta=\frac{17\pm\sqrt{17^2-4(4)(4)}}{8} =\frac{17\pm15}{8}. $$

The two solutions are \(\theta=4\) and \(\theta=1/4\). Only \(\theta=1/4\) satisfies \(|\theta|<1\), so the invertible solution is \(\theta=0.25\). The common mistake is to report both roots as equally acceptable, or to select a root without checking invertibility.

Check. With \(\theta=1/4\), the numerator is \((3/4)(9/8)=27/32\), and the denominator is \(1+1/16+1/4=21/16\). Their ratio is \((27/32)/(21/16)=9/14\), as required.

Deriving the ARMA(1,1) Autocovariances

When \(|\phi|<1\), the model has the causal representation

$$ X_t-\mu=Z_t+(\phi+\theta)Z_{t-1} +(\phi+\theta)\phi Z_{t-2} +(\phi+\theta)\phi^2 Z_{t-3}+\cdots. $$

Thus the coefficient of \(Z_t\) is \(\psi_0=1\), and for \(j\geq1\), the coefficient of \(Z_{t-j}\) is \(\psi_j=(\phi+\theta)\phi^{j-1}\). Since distinct white-noise terms are uncorrelated, the variance is \(\sigma_Z^2\sum_{j=0}^{\infty}\psi_j^2\). For lag \(k\geq1\), the covariance is \(\sigma_Z^2\sum_{j=0}^{\infty}\psi_j\psi_{j+k}\). These sums converge when \(|\phi|<1\).

1
Calculate the variance.
Use \(\psi_0=1\) and the geometric series \(\sum_{j=0}^{\infty}\phi^{2j}=1/(1-\phi^2)\). This gives \(\gamma_0=\sigma_Z^2[1+(\phi+\theta)^2/(1-\phi^2)]\), which simplifies to the expression below.
2
Calculate a positive-lag covariance.
For \(k\geq1\), separate the \(j=0\) term from the remaining geometric series. Factor the result to obtain \(\gamma_k=\sigma_Z^2(\phi+\theta)\phi^{k-1}(1+\phi\theta)/(1-\phi^2)\).
3
Divide by the variance.
By definition, \(\rho_k=\gamma_k/\gamma_0\). The common factor \(\sigma_Z^2\) cancels, leaving the ACF formulas below.

For completeness, the positive-lag simplification follows by summing the covariance terms:

$$ \begin{aligned} \gamma_k &=\sigma_Z^2\left[ (\phi+\theta)\phi^{k-1} +(\phi+\theta)^2\sum_{j=1}^{\infty}\phi^{2j+k-2} \right]\\ &=\sigma_Z^2(\phi+\theta)\phi^{k-1} \left[1+\frac{(\phi+\theta)\phi}{1-\phi^2}\right]\\ &=\sigma_Z^2\frac{(\phi+\theta)\phi^{k-1}(1+\phi\theta)}{1-\phi^2}. \end{aligned} $$

In the last equality, \(1-\phi^2+\phi(\phi+\theta)=1+\phi\theta\). The variance simplifies because \(1-\phi^2+(\phi+\theta)^2=1+\theta^2+2\phi\theta\). Therefore,

$$ \gamma_0=\sigma_Z^2\frac{1+\theta^2+2\phi\theta}{1-\phi^2}, \qquad \gamma_k=\sigma_Z^2\frac{(\phi+\theta)(1+\phi\theta)\phi^{k-1}}{1-\phi^2} \quad (k\geq1), $$

and

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

In particular, \(\rho_k=\phi\rho_{k-1}\) for \(k\geq2\). The ACF therefore tails off geometrically after lag one; it does not generally cut off after lag one. The denominator in \(\rho_1\) is positive under stationarity because \(1+\theta^2+2\phi\theta=(\theta+\phi)^2+1-\phi^2>0\).

Worked Example: Calculating Selected Autocorrelations and Covariances

Let \(\phi=-0.4\), \(\theta=0.6\), and \(\sigma_Z^2=4\). Calculate \(\gamma_0\), \(\rho_1\), \(\rho_2\), and \(\rho_3\), checking the model conditions. The method is to check stationarity, calculate the variance and lag-one ACF, then use the geometric ACF relation. This avoids treating an ARMA process as a pure MA(1) whose ACF cuts off.

The process is stationary because \(|-0.4|<1\), and it is invertible because \(|0.6|<1\). First calculate the denominator term: \(1+\theta^2+2\phi\theta=1+0.36+2(-0.4)(0.6)=0.88\). Also, \(1-\phi^2=1-0.16=0.84\). Hence

$$ \gamma_0=4\frac{0.88}{0.84} =4\frac{22}{21} =\frac{88}{21}\approx4.1905. $$

For the lag-one ACF, \(\phi+\theta=0.2\) and \(1+\phi\theta=1-0.24=0.76\). Thus

$$ \rho_1=\frac{(0.2)(0.76)}{0.88} =\frac{0.152}{0.88} =\frac{19}{110}\approx0.1727. $$

Now apply \(\rho_k=\phi^{k-1}\rho_1\):

$$ \rho_2=(-0.4)\frac{19}{110} =-\frac{19}{275}\approx-0.0691, \qquad \rho_3=(-0.4)\left(-\frac{19}{275}\right) =\frac{38}{1375}\approx0.0276. $$

The signs alternate after lag one because \(\phi<0\). As a check, \(\gamma_1=\rho_1\gamma_0=(19/110)(88/21)=76/105\), and \(\gamma_2=\phi\gamma_1=-152/525\). Their signs and ratios agree with the calculated autocorrelations.

Recovering Parameters From Several ACF Values

If \(\rho_1\ne0\), the lag relation gives a direct estimate of the autoregressive parameter from theoretical autocorrelations:

$$ \phi=\frac{\rho_2}{\rho_1}. $$

After finding \(\phi\), substitute it and \(\rho_1\) into the lag-one formula to solve for \(\theta\). The resulting equation is quadratic:

$$ (\rho_1-\phi)\theta^2+ \bigl(2\phi\rho_1-(1+\phi^2)\bigr)\theta+ (\rho_1-\phi)=0. $$

This equation can have two solutions. Check the original model conditions for each candidate; in an invertible model retain only a solution with \(|\theta|<1\). If \(\rho_1=0\), dividing \(\rho_2\) by \(\rho_1\) is invalid. Under stationarity, \(|\phi\theta|<1\), so \(1+\phi\theta\ne0\); the lag-one formula then implies \(\phi+\theta=0\). This is a cancellation case, not a reason to use the ratio formula.

Worked Example: Recovering Both Parameters and the Innovation Variance

A stationary, invertible ARMA(1,1) process has theoretical autocorrelations \(\rho_1=9/14\) and \(\rho_2=9/28\), and variance \(\gamma_0=7\). Recover \(\phi\), \(\theta\), and \(\sigma_Z^2\). Assume the stated model convention and white-noise innovations. The method is to recover \(\phi\) from consecutive ACF values, solve the lag-one equation, apply invertibility, and then use the variance formula. This avoids mistaking an algebraic root for an admissible model parameter.

Since \(\rho_1\ne0\),

$$ \phi=\frac{\rho_2}{\rho_1} =\frac{9/28}{9/14} =\frac{1}{2}. $$

With \(\phi=1/2\) and \(\rho_1=9/14\), the quadratic from the lag-one equation is \(4\theta^2-17\theta+4=0\), as calculated in the first worked example. Its roots are \(4\) and \(1/4\). Invertibility requires \(|\theta|<1\), so \(\theta=1/4\).

Rearrange the variance formula to solve for the innovation variance:

$$ \sigma_Z^2 =\gamma_0\frac{1-\phi^2}{1+\theta^2+2\phi\theta}. $$

Here \(1-\phi^2=3/4\), and \(1+\theta^2+2\phi\theta=1+1/16+1/4=21/16\). Therefore,

$$ \sigma_Z^2 =7\frac{3/4}{21/16} =7\left(\frac{3}{4}\right)\left(\frac{16}{21}\right) =4. $$

Check. Substitution gives \(\gamma_0=4(21/16)/(3/4)=7\). Also, \(\rho_1=9/14\) and \(\rho_2=(1/2)(9/14)=9/28\), matching both given autocorrelations.

Common Exam Traps

  • Using the wrong MA sign convention. The formulas here assume \(+\,\theta Z_{t-1}\). If the equation instead has \(-\,\theta Z_{t-1}\), replace \(\theta\) by \(-\theta\) before applying these formulas.
  • Using \(\rho_1/\rho_2\) for \(\phi\). The correct relation is \(\rho_2=\phi\rho_1\), so when \(\rho_1\ne0\), \(\phi=\rho_2/\rho_1\).
  • Confusing stationarity with invertibility. Check \(|\phi|<1\) for causality and \(|\theta|<1\) for invertibility. One condition does not imply the other.
  • Ignoring multiple parameter roots. The equation for \(\theta\) can produce more than one algebraic solution. Apply the stated invertibility condition and verify the selected root in the original formula.
  • Assuming the ACF cuts off at lag one. That cutoff is characteristic of an MA(1), not a general ARMA(1,1). For this model, \(\rho_k=\phi^{k-1}\rho_1\) for positive lags.

Parameter recovery from exact theoretical autocorrelations is an algebra problem. In applied work, estimated sample autocorrelations are subject to sampling variation and may not obey the theoretical relations exactly. The formulas remain useful for understanding the model, checking candidate parameters, and recognizing whether a proposed solution is consistent with stationarity and invertibility. These techniques are aligned with the SOA Exam SRM (Statistics for Risk Modeling) syllabus.

Check Your Understanding

For each problem, use the sign convention \(X_t-\mu=\phi(X_{t-1}-\mu)+Z_t+\theta Z_{t-1}\). Answers are given to four decimal places where needed.

  1. A stationary, invertible model has \(\phi=0.3\) and \(\theta=0.2\). Calculate \(\rho_1\) and \(\rho_2\). (Answer: \(\rho_1=0.4569,\ \rho_2=0.1371\))
  2. A stationary model has \(\phi=-0.5\), \(\theta=0.25\), and \(\sigma_Z^2=2\). Calculate \(\gamma_0\), \(\rho_1\), and \(\rho_2\). (Answer: \(\gamma_0=13/6,\ \rho_1=-7/26,\ \rho_2=7/52\))
  3. A stationary, invertible model has \(\rho_1=9/14\) and \(\rho_2=9/28\). Find \(\phi\) and \(\theta\), and explain why the other quadratic root is rejected. (Answer: \(\phi=1/2,\ \theta=1/4\); reject \(\theta=4\) because it is not invertible)
  4. A proposed model has \(\phi=1.1\) and \(\theta=0.2\). Is it stationary and is it invertible? Should the stationary ARMA(1,1) ACF formulas derived here be used? (Answer: not stationary because \(|\phi|>1\); invertible because \(|\theta|<1\); do not use the stationary ACF formulas)
  5. For a stationary model, \(\rho_1=0\). Is \(\rho_2/\rho_1\) a valid way to recover \(\phi\)? What parameter relationship follows from the lag-one formula? (Answer: no; division by zero is invalid, and \(\theta=-\phi\))
Clinical Biostats