在 R 中，如何從 auto.arima 結果中捕獲 ARIMA 元素

Question

我有一堆系列要使用 forecast::auto.arima 函數進行預測。 我喜歡保存 auto.arima 適合的模型類型。 如果您運行以下代碼：

library(forecast)

set.seed(123)

y <- sin(seq(-pi,pi,0.05))+(rnorm(length(seq(-pi,pi,0.05)))/4)

arima.model <- auto.arima(y)

華宇模型

最后一行執行的結果顯示

Series: y 

**ARIMA(1,1,2)** 

Coefficients:
         ar1      ma1     ma2

      0.9594  -1.7285  0.7740

s.e.  0.0380   0.0745  0.0658

sigma^2 estimated as 0.06534:  log likelihood=-6.1

AIC=20.2   AICc=20.53   BIC=31.51

如何捕獲 ARIMA(1,1,2) 並保存結果？ 我希望做類似arima.model$的事情並捕獲我需要的東西，但我無法弄清楚。

Answer 1

您可以嘗試summary(arima.model) 、 arima.model$coef 、 arima.model$aic 、 arima.model$bic 。

如果你想要一個整潔的格式，你可以像這樣使用 broom 包：

library(broom) 
tidy(arima.model) #ar/ma terms
glance(arima.model) #information criteria

tidy(arima.model)
# A tibble: 3 x 3
  term  estimate std.error
  <fct>    <dbl>     <dbl>
1 ar1      0.959    0.0380
2 ma1     -1.73     0.0745
3 ma2      0.774    0.0658

glance(arima.model)
# A tibble: 1 x 4
  sigma logLik   AIC   BIC
  <dbl>  <dbl> <dbl> <dbl>
1 0.256  -6.10  20.2  31.5