One step ahead forecast from an estimated model - error term

Question

na ja 2016-9-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/303411-one-step-ahead-forecast-from-an-estimated-model-error-term

评论： Lorenzo 2023-1-21

I have estimated the model for my series and it is arima(1,1,1). Instead of available k-step ahead option, I need to do one-step ahead forecast. The model is:

 ARIMA(1,1,1) Model:
    --------------------
    Conditional Probability Distribution: Gaussian
                                  Standard          t     
     Parameter       Value          Error       Statistic 
    -----------   -----------   ------------   -----------
     Constant    8.96034e-05    0.00121444      0.0737815
        AR{1}       0.531086     0.0802215        6.62025
        MA{1}      -0.917878     0.0394706       -23.2548
     Variance      0.0378884    0.00419511        9.03158
so, the equation will be

y(k) = .000089 + (0.531086*y(k-1)) + e(k) + (-0.917878*e(k-1))

What will be the value of 'e' term? what should be 'e(k)'? Correct me if I have interpreted anything wrongly.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Brendan Hamm 2016-9-20

在 MATLAB Online 中打开

Yes the model w ould be:

y(k) = .000089 +  y(k-1) + [0.531086*(y(k-1) - y(k-2))] + e(k) + (-0.917878*e(k-1))

e(k) in this example is simply a Normal random variable with mean zero and variance 0.038 such e(k) is independent of e(k-t)) for all t.

Look at the forecast method, you should provide the pre-sample response 'Y0' and innovations 'E0'.

doc arima\forecast

Forecast provides MLE (i.e. all e sampled will be zero), but you can also use the simulate method to sample from that distribution:

doc arima\simulate

na ja 2016-9-22

在 MATLAB Online 中打开

Thanks Brendan! Will do the same. I have one more doubt regarding statistical concept of ARIMA. As I have discussed above, this model is ARIMA(1,1,1). For one step ahead forecast I am using the equation

 [y(k) - y(k-1)] = .000089 + [0.531086*(y(k-1) - y(k-2))] + e(k) + (-0.917878*e(k-1))

What will be the value of e(k-1) term to get 1st forecast (i.e. for k =1)? Is it the last value of residual matrix of training data or something else? I am bothered because according to the equation, i am predicting differences, not the actual values. The residual table gives the residuals for (actual values- predicted values).

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Brendan Hamm 2016-9-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/303411-one-step-ahead-forecast-from-an-estimated-model-error-term#answer_235744

编辑：Brendan Hamm 2016-9-24

在 MATLAB Online 中打开

1. The arima\estimate method will return to you the data on the original scale. That is if you do:

Mdl = arima(1,1,1);
Mdl = estimate(Mdl,data);
res = infer(Mdl,data); % Retrieve inferred residuals (innovations)
foreValues = forecast(Mdl,1,'Y0',data','E0',res) % forecast

Your forecasted data will be on the same scale as data and not the differenced data. The difference operator is just applied and you have the model you wrote above, but are simply returned y(k),y(k+1),...

2. The residual e(k-1) is he last value in the variable res above. That is, if you pass 'E0' to the forecast method it will use the residuals which were infered from the model and data.

5 个评论
显示 3更早的评论隐藏 3更早的评论

na ja 2016-9-23

编辑：na ja 2016-9-23

That's right!

Actually I predicted the model using the same approach.

What I want to implement is one step ahead forecast. For this, the equation is as discussed above. As I am predicting the difference and not the actual value, my error term should be the residuals obtained from differences and not the actual value. 'infer' fits model in terms of original series, which is not appropriate in my case. I have found the previous residual terms using one step ahead approach, as discussed below.

The training data 'y' is of 200 values. So, total differenced terms 'dy' will be 199. To forecast 200th differenced value, e(k-1) = e(199) is needed..

ARIMA Equation is [y(k) - y(k-1)] = .000089 + [0.531086 * (y(k-1) - y(k-2))] + e(k) + [-0.917878*e(k-1)]

forecast = 0; // As no past value is available for forecast, initial value is 0 // _

for

k = 1:199 // total 199 differenced terms//

error(k) = dy(k) - forecast; // find residual//

next_forecast(k) = .000089 + (0.531086 * dy(k) ) +(-0.917878) *error(k); //next forecast

forecast = next_forecast(k); // update forecast value

end

Now last value of error(k) will be the first residual value e(k-1) for next forecast , which is implemented for 50 points of test data, one-step ahead .

Is this right??

na ja 2016-9-26

Thanks Brendan! That made it all clear. :)

Lorenzo 2023-1-21

Thank you for this answer! Can you then plot the forecasted values combined with the observed values so that we get a graph were we see the continuity between observed and forecasted values?

请先登录，再进行评论。

One step ahead forecast from an estimated model - error term

3 个评论
显示 1更早的评论隐藏 1更早的评论

采纳的回答

5 个评论
显示 3更早的评论隐藏 3更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

One step ahead forecast from an estimated model - error term

3 个评论 显示 1更早的评论隐藏 1更早的评论

采纳的回答

5 个评论 显示 3更早的评论隐藏 3更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

3 个评论
显示 1更早的评论隐藏 1更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论