Econ0022 Problem Set 3
Econ0022 Problem Set 3
Econ0022 Problem Set 3
Autumn 2018
1 Practice Problems
1. Consider iterative forecasts over the horizon T + 1; :::; T + h using an AR(1)
model, Yt = + Yt 1 + "t , given data up to time T . We ignore parameter
uncertainty and treat the parameters and as known such that the iterative
forecast takes the form
(a) Show that the …nal h-step ahead forecast, YT +hjT , can be expressed as
n o
YT +hjT = h YT + 1 + + 2 + ::: + h 1 :
If j j < 1 : YT +hjT = Y;
1
If = 1 : YT +hjT ! 1;
eT +h = YT +h YT +hjT :
1
(d) Show that
n o
2 2 2(h 1)
E [eT +h ] = 0; Var (eT +h ) = " 1+ + :::: :
What is the h-step MSFE for the iterative procedure using an AR(1)
model? What happens to the MSFE as h grows large?
(e) Explain how you could use the results in question (d) to construct a 95%
forecast interval for your h-period ahead point forecast.
2. Suppose we have obtained the following estimates of an ADL model for annual
US unemployment (unemp) using in‡ation (inf l) as additional preditor,
d t = 1:304 + :647unempt
unemp 1 + :184inf lt 1; ^ " = :883: (1)
(:490) (:084) ((:041)
(a) Suppose that unemp2012 = 7:8 and inf l2012 = 1:3. What is your forecast
of unemployment in 2013?
(b) Suppose the errors in the ADL model are normally distributed. Compute
a 95% and 67% forecast interval.
(c) Can you use the above ADL model to forecast employment for 2014?
Explain.
(d) You decide to build a model for in‡ation in addition to the above model
for unemployment. You …nd that an AR(1) model provides an adequate
description of in‡ation,
c lt = 1:277 + :665inf lt
inf 1; ^ " = :653: (2)
(:558) (:107)
(a) Let W be a random variable with mean W and variance 2 . Show that
W
for any constant c:
h i
E (W c)2 = 2W + ( W c)2 :
2
denote the unconditional mean squared forecast error associated with
fT +1jT . Show that M SF E is minimized w.r.t. fT +1jT at fT +1jT = E [YT +1 ].
Interpret the result.
(c) Consider instead the conditional mean squared forecast error associated
with fT +1jT ,
h i
2
M SF ET = E YT +1 fT +1jT jYT ; YT 1 ; ::: :
(a) As a …rst way of choosing the number of lags for the AR model for Yt ,
compute the AIC and BIC and LR criteria for the given sample [this can
be done with the STATA command varsoc dlgdp - note that STATA call
the BIC criterion SBIC]. Which model(s) does each of the three criteria
recommend?
(b) As a second way of choosing the appropriate number of lags, we can use
the pseudo-forecast evaluation criterion.
To this end, we wish to compute rolling pseudo-forecasts as described in
the slides. To perform this in STATA, you need to use the STATA com-
mands that you …nd in the …le ’rolling forecast US-dgp.txt’ posted
on Moodle/Problem Sets which performs the forecast evaluation using an
AR(3) model.
The …rst part of the code says:
//Pseudo-forecast of dlgdp using AR(3) model for T=150,...,252.
gen f_dlgdp3=.
local start 150
local end 252
forv T=‘start’/‘end’ {
quietly regress dlgdp L(1/3).dlgdp in 1/‘T’
quietly replace f_dlgdp3 = _b[_cons] + _b[L1.dlgdp]*dlgdp[‘T’]
+ _b[L2.dlgdp]*dlgdp[‘T’-1] + _b[L3.dlgdp]*dlgdp[‘T’-2] in ‘T+1’/‘T+1’
The above code runs a loop over T = 150; 151; :::; 252, where for each value
of T STATA does the following:
3
i. Estimate an AR(3), Y^t = ^ T + ^ 1;=t Yt 1 + ^ 2;T Yt 2 + ^ 3;T Yt 3 , using
only observations fY1 ; :::; YT g.
ii. Compute Y^T +1jT = ^ T + ^ 1;T YT + ^ 2;T YT 1 + ^ 3;T YT 2 and save this
value in f_dlgdp3.
To perform steps (i)-(ii), copy the above part of the code from ’rolling
forecast US-dgp.txt’and paste it into the command prompt in STATA
(and press enter...).
If you open the Data Editor, you’ll see that the variable f_dlgdp3 has
been created containing the rolling pseudo-forecasts, Y^T +1jT .
(c) To visualize the performance of the AR(3) as forecast model, plot the
actual observed value of YT +1 together with its pseudo-forecasts. This can
be done with the command ’tsline dlgdp f_dlgdp3 if tin(1984q2,)’.
Does the AR(3) model appear to do a good job?
(d) Next, compute pseudo-RMSFE in STATA by de…ning a new variable which
is the di¤erence between by Y and Y^ , and compute its sample variance
[gen e3 = dlgdp - f_dlgdp3; summarize e3]
(e) Next, re-do questions (b)-(c), this time using an AR(1) model. This can
be done by copy-pasting the second part of the code in the txt …le rolling
forecast US-dgp.txt into the command line in STATA. This generates a
new variable f_dlgdp1 that contains the rolling forecasts using an AR(1)
model. Plot this new forecast vs. actual outcome, and compute the cor-
responding pseudo-RMSFE.
(f) We now have two pseudo-RMSFE’s - one using an AR(1) model, and one
using an AR(3) model. In terms of pseudo-RMSFE, which of the two
models would you recommend? How does this recommendation compare
to the model(s) chosen by AIC, BIC and LR, respectively, in question (a)?