On Fri, 14 Mar 2014, Sven Schreiber wrote:

> Am 14.03.2014 09:50, schrieb michel.pouchain(a)univ-paris13.fr:
>> Dear gretl users,
>> I have an example for a
>> little dynamic model like
>>   :
>> y(t) =a0+a1*y(t-1)+a2*x(t)
>> When I use, via the scroll menu, "prévisions" I have
>> the choice beetween static oy dynamic forecats.
>> When I use "dyn", the forecast for T+1 is not equal to the
>> forecats(static) for T+1. Why ?
>
> Do you mean the point forecasts or the standard errors / intervals? If
> you really mean the point forecasts, my guess would be that there is
> some small difference in the specification (sample) that you haven't
> noticed, otherwise it would indeed be strange.
>
> Showing the output would be helpful.
>
>> And I like to know which
>> formulas used for the standards error?
>> I have use with (R) program the formulas given by Pagan Nicholls
>> (1984). Like microfit5
>>
>
> The 'fcast' help says: "For static linear models standard errors are
> computed using the method outlined by Davidson and MacKinnon (2004);
> they incorporate both uncertainty due to the error process and parameter
> uncertainty (summarized in the covariance matrix of the parameter
> estimates). For dynamic models, forecast standard errors are computed
> only in the case of a dynamic forecast, and they do not incorporate
> parameter uncertainty. For nonlinear models, forecast standard errors
> are not presently available."

Here's an example:

<hansl>
open data9-7
smpl ; -2
ols QNC const QNC(-1) INCOME
smpl --full
# static forecast
fcast 1990:3 1990:4 --static --quiet
series fc1 = $fcast
series fe1 = $fcerr
# default forecast (dynamic out of sample)
fcast 1990:3 1990:4 --quiet
series fc2 = $fcast
series fe2 = $fcerr
smpl 1990:3 1990:4
print fc1 fe1 fc2 fe2 -o
</hansl>

<partial-output>
                 fc1          fe1          fc2          fe2

1990:3     2561.975     281.4759     2561.975     272.0654
1990:4     2461.011     280.6808     2572.629     310.1355
</partial-output>

As regards point values, the two forecasts for 1990:3 (the first 
out-of-sample observation) are equal. But the forecasts for 1990:4 are not 
equal, because the static one uses observed QNC on the right-hand side 
while the dynamic one uses forecast QNC.

The forecast standard errors are not equal even for 1990:3. This is 
because of the point noted by Sven: our standard errors for dynamic 
forecasts do not incorporate parameter uncertainty. Therefore, for 1990:3 
the fe2 value is slightly smaller than fe1. For 1990:4, however, the fe2 
is larger than fe1 -- naturally, since the uncertainty due to the error 
process is compounded; the latter uncertainty comes to dominate as the 
forecast horizon increases.

Allin Cottrell

Reply via email to