您的訂閱是我製作影片的動力
訂閱點這裡~
若內容有誤,還請留言指正,謝謝您的指教
#影片程式碼
#延續DAY25
fit <- arima(時序格式 ,order=c(1,1,2),include.mean=FALSE)
fit=arima(h,order=c(10,1,12),seasonal=list(order=c(2,1,0),period=12))
#有差分中間請放1且include.mean=FALSE
#模型檢查
tsdisplay(residuals(fit), lag.max=50, main='殘差')
shapiro.test(fit$residuals) # > 常態
Box.test(fit$residuals, lag=10, type="Ljung-Box") # > 獨立
#非季節性=10 季節性=period*2
檔案下載處:
https://drive.google.com/drive/u/1/folders/1QIRQIvqYF0C37gfGURHMoSElh_6BYUr6
資料來源:
https://www.twse.com.tw/zh/