我已经为大型数据集(几年的每小时数据)开发了一个 SARIMA 模型。 我一直在尝试获取模型的结果并手动计算(验证)一些每小时结果,但我一直无法这样做。 我只是希望能够获取生成值和实际值的表格,并根据模型的输入以及(已经估计的)参数和残差来计算预测一小时的时间。 我不想重新计算参数、残差、任何误差计算或拟合结果,因为我知道这些是该解决方案迭代收敛的结果,无法通过简单的一次性计算复制。
附加到这篇文章的是从 JMP 导出的结果以及我的手动计算(验证尝试)。
A 到 G 列是 SARIMA (1,1,1) (0,1,0) [24] 无截距模型的结果。该模型并不理想(我为自己的目的使用了更复杂的模型),但为了验证预测的计算方式,它应该就足够了。 (如果我能让这个工作,我可以让另一个工作)。
J 到 N 列是我根据模型参数信息和相关的先前实际值和残差值手动计算给定预测的尝试。 为此,我计算 J 列中的双差 y,然后计算 K 中预测的 AR 部分,L 部分中预测的 MA 部分,然后将它们放在一起并撤消 M 中的双差以获得计算预测。 N 列具有使用此手动预测和实际目标值计算的残差。
两个适用的参数在 P4 开始的表格中。 原始拟合被迫计算没有截距,所以它只是两个。
对于我的生活,我无法弄清楚他们为什么不结盟。 我是不是试图迭代地计算误差项,我只想看看给定的预测是如何从文件中的信息计算出来的,我不能。 当我强制它运行更简单的模型(例如 MA(1))时,我仍然无法解决问题。 我有一段时间认为这只是我的一个简单的计算错误,但我不再相信它是。 如果有人能告诉我我做错了什么计算,或者为什么这个任务不可能完成,我将不胜感激。 我怀疑这很简单,但我尝试在 MA 术语上颠倒符号并没有帮助,我相信我的差分和非差分计算是正确的。
全面披露:本作品既不是作业/考试/项目,也不是工作项目,这只是个人项目,数据为公开信息。 这样做的目的是能够向一位老同事展示 excel 工作表验证,以便他可以看到 SARIMA 模型的用例和强大功能。
任何帮助是极大的赞赏!