全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
10970 3
2016-04-24
请问一下ARIMA预测中auto.arima函数返回的结果中的(p,d,q)参数如何提取啊?因为要用循环做大量的预测,所以手工提取实在不合适。再者,我也尝试了很多方法均无法提取出来这个参数,attr查看变量内部结构和内容竟然也没有发现有这个变量。求高手指教下,谢谢

  • 问题补充 : 3小时前

    >m<-auto.arima(apply.weekly(ts_300,mean))

    >m

    Series: apply.weekly(ts_300, mean)

    ARIMA(1,1,1)                    

    Coefficients:

             ar1      ma1

          0.4111  -0.8993

    s.e.  0.1490   0.0750

    sigma^2 estimated as 38.07:  log likelihood=-203.48

    AIC=412.96   AICc=413.37   BIC=419.39

    就是想提取标红的结果.

    我查看str(m)的结果中没有这个参数

    > str(m)

    List of 17

    $ coef     : Named num [1:2] 0.411 -0.899

      ..- attr(*, "names")= chr [1:2] "ar1" "ma1"

    $ sigma2   : num 38.1

    $ var.coef : num [1:2, 1:2] 0.02219 -0.0069 -0.0069 0.00562

      ..- attr(*, "dimnames")=List of 2

      .. ..$ : chr [1:2] "ar1" "ma1"

      .. ..$ : chr [1:2] "ar1" "ma1"

    $ mask     : logi [1:2] TRUE TRUE

    $ loglik   : num -203

    $ aic      : num 413

    $ arma     : int [1:7] 1 1 0 0 1 1 0

    $ residuals: Time-Series [1:64] from 1 to 64: 0.00367 -0.58768 4.29421 -3.38195 2.15109 ...

    $ call     : language auto.arima(x = structure(list(x = structure(c(3.66666666666667, 3, 7.85714285714286,  2.71428571428571, 5.71428571428571, 44.5714285714286, 17.2857142857143,  ...

    $ series   : chr "ts_300_weekly"

    $ code     : int 0

    $ n.cond   : int 0

    $ nobs     : int 63

    $ model    :List of 10

      ..$ phi  : num 0.411

      ..$ theta: num -0.899

      ..$ Delta: num 1

      ..$ Z    : num [1:3] 1 0 1

      ..$ a    : num [1:3] -1.77 1.64 2.57

      ..$ P    : num [1:3, 1:3] 0.00 0.00 7.10e-29 0.00 1.79e-07 ...

      ..$ T    : num [1:3, 1:3] 0.411 0 1 1 0 ...

      ..$ V    : num [1:3, 1:3] 1 -0.899 0 -0.899 0.809 ...

      ..$ h    : num 0

      ..$ Pn   : num [1:3, 1:3] 1.00 -8.99e-01 2.66e-23 -8.99e-01 8.09e-01 ...

    $ bic      : num 419

    $ aicc     : num 413

    $ x        : Time-Series [1:64] from 1 to 64: 3.67 3 7.86 2.71 5.71 ...

    - attr(*, "class")= chr [1:2] "ARIMA" "Arima"




二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2016-4-25 10:56:35
没人愿意解答啊,好头疼啊
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-11-30 14:45:12
可以用arimaorder(m)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2020-3-16 09:11:49
用names()查看结果中参数名
用m$...提取
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群