D 57 MTW 11 Solwb
D 57 MTW 11 Solwb
D 57 MTW 11 Solwb
Please do NOT open this booklet until you are asked to do so.
There are 9 pages including this page and t tables. Please check to see if you have all the pages.
Aids allowed: You are allowed to use the textbook (Time Series Analysis and Its Applications
with R Examples by R. H. Shumway and D. S. Stoffer), the class notes and a non-
communicating calculator. No other material will be allowed during the test.
All your work must be presented clearly in order to get credit. Just an answer with no other work
shown will only qualify for zero credit. Show your work and answer in the space provided, in
ink. Pencil may be used, but then any re-grading will NOT be allowed.
Last name:____________________
First name:____________________
Other names:___________________
Page 1 of 9
1) [ 6 points] Let { X t } be a weakly stationary time series and let’s assume that we have made
three observations from this series i,e, n = 3. Let = 5, = 3 and = 4. Use this data to
estimate γ (1). (i.e. calculate γˆX (1) )
Sol
1 n−h
γˆX (h) = ∑ ( xt −x )( xt + h − x )
n t =1
1 n −1
γˆX (1) = ∑ ( xt −x )( xt + h − x )
n t =1
1 2
= ∑ ( xt −x )( xt + h − x )
3 t =1
1
= ( ( 5 − 4 ) (3 − 4) + (3 − 4)((4 − 4) )
3
1
=−
3
Page 2 of 9
2) Consider the AR(2) process given by X t = X t −1 − 0.5 X t − 2 + wt .
a) [5 points] Is this process weakly stationary? Give reasons for your answer.
Sol
The characteristic equation is
The roots of this equation (regarding B as a variable) are complex, namely, 1± i. As the rnodulus
of both roots exceeds one, the roots are both outside the unit circle and so the process is
stationary.
Sol:
Since this process is stationary Corr ( X 1 , X 2 ) = ρ (1) and Corr ( X 5 , X 7 ) = ρ (2)
p
Yule-Walker equations γ (h) − ∑ φiγ (h − i ) = σ w2 I (h = 0) .
i =1
h=1
2
γ (1) − ∑ φiγ (1 − i ) = 0
i =1
h= 2
2
γ (2) − ∑ φiγ (2 − i ) = 0
i =1
Page 3 of 9
3) [7 points] Consider the stationary ARMA(1, 1) process given by
X t = 0.5 X t −1 + wt − 0.3wt −1 .
Since the process is (weakly) stationary, it has the Wold representation. Calculate the first four
coefficients of the Wold representation (i.e. ψ 0 , ψ 1 , ψ 2 , and ψ 4 , )
Sol
Page 4 of 9
4) [ 6 points] A stationary process { X t } has ρ (1) = 0.97 and ρ (2) = 0.90. Calculate φ2,2 (i.e. the
partial autocorrelation at lag 2).
Sol
−1
1 0.97 0.97 1.64129
φ2,2 is the second element of the vector =
0.97 1 0.90 -0.69205
Data Display
Matrix M2
16.9205 -16.4129
-16.4129 16.9205
.
MTB > Read 2 1 m3.
DATA> 0.97
DATA> 0.90
2 rows read.
Data Display
Matrix M3
0.97
0.90
Data Display
Matrix M5
1.64129
-0.69205
MTB >
∑α
i =1
*
ih γ ( j − i ) = γ ( j ), j = 1,… , h can use ρ’s in place of γ’s.
Page 5 of 9
In this question h = 2 and so
j =1
α12* ρ (0) + α 22
*
ρ (1) = ρ (1)
α12* + 0.97α 22
*
= 0.97
j= 2
α12* ρ (1) + α 22
*
ρ (0) = ρ (2)
0.97α12* + α 22
*
= 0.97
*
and φ2,2 = α 22
Page 6 of 9
5) [ 6 points] A time series { X t } is given by X t = 2 wt + wt −1 + vt where wt and vt are
uncorrelated (i.e. Cov( ws , vt ) = 0 for all s, t.) white noise processes, each with mean 0 and
variance 3.
Calculate the cross-correlation function between {wt } and { X t } .
Sol
V ( X t ) = 4V ( wt ) + V ( wt −1 ) + V ( vt ) = 4 × 3 + 3 + 3 = 18
γ w, x (h) = Cov( wt , X t + h )
= Cov( wt , 2 wt + h + wt + h−1 + vt + h )
= 2Cov( wt , wt + h ) + Cov( wt , wt + h−1 )
6 if h = 0
= 3 if h =1
0 otherwise
Corr ( wt , X t + h )
ρ w , x ( h) =
VwVX t t+h
6 / 3 ×18 if h = 0
= 3 / 3 ×18 if h =1
0 otherwise
Page 7 of 9
6) [ 6 points] A time series { X t } is given by X t = U cos(ct ) + V sin(ct ) where U and V are iid
N(0, σ ), σ > 0 and c > 0 is a constant. Calculate Cov( , ) for h = 0, ±1, ±2, ⋯. Is { X t }
a weakly stationary process?
( Hint: cos( A − B ) = cos( A) cos( B ) + sin( A) sin( B ) )
Sol
Page 8 of 9
7) Consider the time series model given by X t = 0.5 + 0.75 X t −1 + wt where wt are i.i.d.
(0, σ ). You may assume that this process is weakly stationary. [Note: This is not exactly the
AR(1) process we discussed in class, but the necessary formulas can be developed in a similar
manner.]
a) [4 points] If =1.5, find a one-step ahead forecast (i.e. )
Sol
( )
X tt+l = E X t + l | X t
= E ( 0.5 + 0.75 X t + l −1 + wt + l | X t )
0.5 + 0.75 X t if l = 1
= t
0.5 + 0.75 X t +l −1 if l ≥ 2
And so = 0.5+0.75*1.5 = 1.625
b) [ 6 points] If =1.5, find a 10-step ahead forecast (i.e. )
Sol
X tt+ 2 = 0.5 + 0.75 X tt+1
= 0.5 + 0.75 ( 0.5 + 0.75 X t )
= 0.5 + 0.5 × 0.75 + 0.752 X t
( )
X tt+l = 0.5 + 0.5 × 0.75 + 0.5 × 0.752 + + 0.5 × 0.75l −1 + 0.75l X t
0.5(1 − 0.75l )
= + 0.75l X t
1 − 0.75
= 2(1 − 0.75l ) + 0.75l X t
Sol
E ( X t ) = µ ∀t since the process is stationaty .
µ = 0.5/(1-0.75) = 2
Page 9 of 9