问题补充 : 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"