我是时间序列分析的新手。我想基于ARIMA模型预测一些值,并根据保持数据绘制它们。
这就是我所拥有的:
% Retrieve IBM closing prices
c = yahoo;
IMB_data = fetch(c,'IBM','Close','1-sept-2014','31-dec-2014');
IMB_date = IMB_data(:,1);
IMB_close = IMB_data(:,2); % Contains 85 data points
使用autocorr
和parcorr
,我推断出最合适的模型应该是:
mdl = arima(1,0,0);
使模型适应前60个观测值,并保留剩余的15个观测值以评估预测性能。
fit_model = estimate(mdl, IMB_close(1:60));
使用拟合模型预测15个周期范围:
[Y, MSE] = forecast(fit_model,15,'Y0',IMB_close(1:60));
直观地将预测与保留数据进行比较:
figure
plot(IMB_date, IMB_close,'Color',[.7,.7,.7])
datetick
hold on
plot(61:85,Y,'b');
hold off
问题是预测的值没有被绘制,我做错了什么?