Chapter 7

AR-, MA-, and ARMA-models

[15/0%]
Show
LaTeX
Problem 7:1

Suppose that an AR(1)-process, {Xt,t=0,±1,…}\left\{ X_{t}, t = 0, \pm 1, \ldots \right\}, has the covariance function r(k)=10⋅0.5∣k∣r(k) = 10 \cdot 0.5^{\left|k\right|}. Which AR(1)-process is it? What is the variance of the innovations?

?
Problem 7:2

Find the spectral density, R(f)R(f), and the covariance function, r(k)r(k), at time lag zero, one and two, of the AR(2)-process

Xt−Xt−1+0.5Xt−2=et, X_{t} - X_{t-1} + 0.5 X_{t-2} = e_{t},

where {et}\left\{ e_{t}\right\} is a white innovation process with variance Var⁡[et]=4\operatorname {Var}\left[e_{t}\right] = 4.

?
Problem 7:3

The covariance function of a signal, {Yt}\left\{ Y_{t}\right\}, in discrete time, has been well estimated at a few time lags:

τ:01234rY(τ):2.41.6−0.4−0.6 \begin{array}{c|ccccc} \tau : & 0 & 1 & 2 & 3 & 4 \\ \hline r_{Y}(\tau ): & 2.4 & 1.6 & & -0.4 & -0.6 \end{array}

We decide to use the following AR(2)-process to model the signal,

Yt−Yt−1+a2Yt−2=et. Y_{t} - Y_{t-1} + a_{2} Y_{t-2} = e_{t}.

Find the innovation variance, the second AR-parameter, and the covariance function at time lag two.

?
Problem 7:4

The MA(q)-process

Yt=et+b1et−1+⋯+bqet−q, Y_{t} = e_{t} + b_{1} e_{t-1} + \cdots + b_{q} e_{t-q},

where {et}\left\{ e_{t}\right\} is a white innovation process with unit variance, has covariance function r(0)=54r(0) = \frac{5}{4}, r(1)=0r(1) = 0, r(2)=12r(2) = \frac{1}{2}, and r(k)=0r(k) = 0, ∣k∣≥3\left|k\right| \geq 3. Determine the order of the process, the MA-coefficients, and the spectral density.

?
Problem 7:5

Let et,t∈Ze_{t}, t \in \mathbb {Z}, be independent random variables with expectation mm and variance σ2\sigma^{2}. Define, for a constant MM,

Xt=et+2et−1−et−3,Zt=Xt−M. X_{t} = e_{t} + 2 e_{t-1} - e_{t-3}, \quad Z_{t} = X_{t} - M.
?
(a)

Choose MM, such that {Zt}\left\{ Z_{t}\right\} is an MA(3)-process. Determine the expectation, covariance function, and spectral density of {Zt}\left\{ Z_{t}\right\}.

(b)

Find expectation, covariance function, and spectral density of {Xt}\left\{ X_{t}\right\}.

Problem 7:6

Let {Xt,t=0,±1,…}\left\{ X_{t}, t = 0, \pm 1, \ldots \right\} be a stationary AR(1)-process,

Xt+0.8Xt−1=et, X_{t} + 0.8 X_{t-1} = e_{t},

and let {Yt,t=0,±1,…}\left\{ Y_{t}, t = 0, \pm 1, \ldots \right\} be an MA(2)-process,

Yt=ut+0.8ut−1+0.4ut−2, Y_{t} = u_{t} + 0.8 u_{t-1} + 0.4 u_{t-2},

where {et}\left\{ e_{t}\right\} and {ut}\left\{ u_{t}\right\} are independent white Gaussian processes with unit variance. Realizations and estimated covariance functions of {Xt}\left\{ X_{t}\right\} and {Yt}\left\{ Y_{t}\right\} are shown in Figure 7.13, but which subplot shows what?

Figure 7.13 Realizations and estimated covariance functions of one AR(1)- and one MA(2)-process in Exercise 7:6.Figure 7.13 Realizations and estimated covariance functions of one AR(1)- and one MA(2)-process in Exercise 7:6.

?
Problem 7:7
?
(a)

Combine, in Figure 7.14, realizations, covariance functions, and spectral densities for two stationary processes.

Figure 7.14 Realizations, covariance functions, and spectral densities to Exercise 7:7a.Figure 7.14 Realizations, covariance functions, and spectral densities to Exercise 7:7a.

(b)

Combine the processes with the correct pole-zero plot in Figure 7.15, where 'x' depicts poles and 'o' zeros. The processes are of the type AR or MA. State the type and order of the processes.

Figure 7.15 Pole-zero plots to Exercise 7:7b.Figure 7.15 Pole-zero plots to Exercise 7:7b.

Problem 7:8

Figure 7.16 shows pole-zero plots, spectral densities, and covariance functions for three AR(2)-processes. Determine which figures that describe the same process.

Figure 7.16 Pole-zero plots, covariance functions, and spectral densities of two processes in Exercise 7:8.Figure 7.16 Pole-zero plots, covariance functions, and spectral densities of two processes in Exercise 7:8.

?
Problem 7:9

(This is a classical example of time series analysis in hydrology.) The Wabash river water level, WtW_{t} (meter), is measured once every month. The following AR-type model has been suggested:

Wt−1.559Wt−1+0.81Wt−2=3.5+et, W_{t} - 1.559 W_{t-1} + 0.81 W_{t-2} = 3.5 + e_{t},

where {et}\left\{ e_{t}\right\} is a white noise process with variance 4. What is the expected water level according to this model? What is the standard deviation?

?
Problem 7:10

Let {Xk,k∈Z}\left\{ X_{k}, k \in \mathbb {Z}\right\} be an AR(1)-process, Xk+a1Xk−1=ekX_{k} + a_{1} X_{k-1} = e_{k}, with Var⁡[ek]=σ2\operatorname {Var}\left[e_{k}\right] = \sigma^{2}. Construct a new process, {Ym,m∈Z}\left\{ Y_{m}, m \in \mathbb {Z}\right\}, by using every second element of {Xk}\left\{ X_{k}\right\}, that is: Ym=X2mY_{m} = X_{2m}. Prove that {Ym}\left\{ Y_{m}\right\} is an AR(1)-process and that it thus can be written in the following form:

Ym+aYm−1=um. Y_{m} + a Y_{m-1} = u_{m}.

Also, compute aa and Var⁡[um]\operatorname {Var}\left[u_{m}\right].

?
Problem 7:11

Compute the cross-covariance, Cov⁡[Xt+k,et]\operatorname {Cov}\left[ X_{t+k}, e_{t} \right], k=…,−1,0,1,…k = \ldots , -1, 0, 1, \ldots, where {Xt}\left\{ X_{t}\right\} is an AR(1)-process: Xt=0.4Xt−1+etX_{t} = 0.4 X_{t-1} + e_{t}, and Var⁡[et]=2\operatorname {Var}\left[e_{t}\right] = 2.

?
Problem 7:12

A random process, {Xt}\left\{ X_{t}\right\}, is generated by filtering white noise:

Xt+a1Xt−1=et,E[et]=0,Var⁡[et]=σ2, X_{t} + a_{1} X_{t-1} = e_{t}, \quad \mathbb {E}\left[e_{t}\right] = 0, \quad \operatorname {Var}\left[e_{t}\right] = \sigma ^{2},

where a1a_{1} and σ2\sigma^{2} are unknown. A simulation gives the following values:

x1,…,x8−0.15−0.971.350.22−0.79−0.52−0.57−0.92x9,…,x15−0.201.35−1.32−0.270.020.59−0.83 \begin{array}{c|cccccccc} x_{1}, \ldots , x_{8} & -0.15 & -0.97 & 1.35 & 0.22 & -0.79 & -0.52 & -0.57 & -0.92 \\ \hline x_{9}, \ldots , x_{15} & -0.20 & 1.35 & -1.32 & -0.27 & 0.02 & 0.59 & -0.83 & \end{array}
?
(a)

Estimate rX(0)r_{X}(0) and rX(1)r_{X}(1), using the given observations.

(b)

Use your estimates in (a) to estimate a1a_{1} and σ2\sigma^{2}.

Problem 7:13

Let {Xt}\left\{ X_{t}\right\} be stationary with mean value mm and assume

(Xt−m)+0.25(Xt−1−m)=et, (X_{t} - m) + 0.25 (X_{t-1} - m) = e_{t},

where {et}\left\{ e_{t}\right\} is Gaussian white noise with mean value zero and unit variance. The following values of XtX_{t} are observed:

x1,…,x50.5310.9441.693−0.3900.171x6,…,x101.126−0.4970.4312.1451.193 \begin{array}{c|ccccc} x_{1}, \ldots , x_{5} & 0.531 & 0.944 & 1.693 & -0.390 & 0.171 \\ \hline x_{6}, \ldots , x_{10} & 1.126 & -0.497 & 0.431 & 2.145 & 1.193 \end{array}

Estimate mm and compute a 95 confidence interval for your estimate.

?
Problem 7:14

Find the best one-step ahead predictor of the ARMA(1,1)-process defined by Xt−0.5Xt−1=et+0.5et+1X_{t} - 0.5 X_{t-1} = e_{t} + 0.5 e_{t+1}.

?
Problem 7:15

Figure 7.17 shows the mean wave height, XtX_{t}, outside the Pakistan coast, as measured by a satellite passing the Indian Ocean on May 21, 1993. The process is well described by an AR(3)-model:

Xt+0.2Xt−1+0.27Xt−2+0.0824Xt−3=σet, X_{t} + 0.2 X_{t-1} + 0.27 X_{t-2} + 0.0824 X_{t-3} = \sigma e_{t},

where ete_{t} is white noise with unit variance. Compute rX(4)r_{X}(4) and rX(5)r_{X}(5) from the first four covariances:

τ:0123rX(τ):0.3326−0.0474−0.07640.0007 \begin{array}{c|cccc} \tau : & 0 & 1 & 2 & 3 \\ \hline r_{X}(\tau ): & 0.3326 & -0.0474 & -0.0764 & 0.0007 \end{array}

Figure 7.17 Wave height, spectrum, and covariance function estimated using data, and spectrum and covariance function of the model in Exercise 7:15.Figure 7.17 Wave height, spectrum, and covariance function estimated using data, and spectrum and covariance function of the model in Exercise 7:15.

?