我已經為大型數據集(幾年的每小時數據)開發了一個 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 模型的用例和強大功能。
任何幫助是極大的讚賞!