Am 23.11.2024 um 18:02 schrieb Brian Revell:
Is there any reason why recursive forecasts cannot be included in the ARIMA option for univariate modelling together with their confidence intervals when the model is truly univariate with no exogenous varisbles included in the specification  Clearly post sample data the MA terms would drop out of the forecasts that would effectively only require the AR parameters and any differencing.

I'm wondering whether what you're trying to achieve is the same thing as what's called "recursive" in gretl. Maybe it is, but the terminology can be complex and not always universal.

So: If you have a base sample from t=T1 to t=T2, and you estimate your Arima model on that sample, and then you want to create forecasts for the out-of-sample range T2+1 through T2+h, for a certain positive integer h, you can do that, but we wouldn't call it recursive. Example:

<hansl>
open AWM18.gdt

h = 3
smpl +0 -h    # leave some obs for the forecasting range, T2 is now 2017:1
arima 1 1 1; log(YED)

fcast --dynamic --out-of-sample # prints out h=3 forecast values up to 2017:4

</hansl>

This is a forward-iterated forecast. Again, maybe this is _not_ what you're actually trying to do, I just want to make sure there are no misunderstandings, because sometimes people call this thing recursive.

In contrast, what gretl calls "recursive" --and I hope I'm getting this right-- entails updating/re-estimating the model coefficients for every new value, starting from very early in the original base sample. So T2 is not fixed anymore (and this could become computationally expensive for non-OLS estimators). Example:

<hansl>
open AWM18.gdt

h = 3
smpl --full
ols ldiff(YED) const ldiff(YED(-1))

fcast h --recursive # prints out many forecast values, each h-step

</hansl>

I'm noticing that the printout here has no values for the final two observations - to be honest, I don't know why that is so, all the necessary ingredients should be there. Can anybody explain this, please?

thanks

sven