我可以阻止 sarima.for 进行绘图吗？答案

【问题标题】：Can I prevent sarima.for from plotting?我可以阻止 sarima.for 进行绘图吗？
【发布时间】：2019-04-26 16:02:47
【问题描述】：

我想使用 astsa 库中的函数 sarima.for 预测时间序列。然而，我想阻止情节并只保留数据。

        require("astsa")
        set.seed(1)
        x<-ts(rnorm(31),start=1980,end=2010)
        sarima.for(x, n.ahead = 10, p=0, d=1, q=0)

我在这里找不到可选性。有没有？或者有没有办法为函数定义一个合适的包装器？

【问题讨论】：

蛮力替代方案是：sarima.for(x, n.ahead = 10, p=0, d=1, q=0); dev.off()

标签： r plot arima

【解决方案1】：

您可以将输出发送到一个临时文件，然后删除该文件：

sarima.noplot = function(x, ...) {
  png(tf<-tempfile())
  out <- sarima.for(x, ...)
  dev.off()
  file.remove(tf)
  return(out)
}

sarima.noplot(x, n.ahead = 10, p=0, d=1, q=0)
# $pred
# Time Series:
#   Start = 2011 
# End = 2020 
# Frequency = 1 
# [1] 1.424851 1.491022 1.557193 1.623364 1.689535 1.755706 1.821877
# [8] 1.888048 1.954220 2.020391

【讨论】：

感谢您接受这个答案。虽然我个人觉得@JilberUrbina 的答案更好。它更干净，并且还返回s.e.

【解决方案2】：

您还可以使用 stats 包中的 arima 和 predict 函数来获得相同的结果而无需绘图。

你可以从帮助文件中看到sarima.for - is a wrapper for R's predict.Arima.

> n <- length(x)
> fit <- stats::arima(x, order = c(p=0, d=1, q=0), xreg=1:n)
> stats::predict(fit, n.ahead=10, newxreg = (n + 1):(n + 10))
$`pred`
Time Series:
Start = 2011 
End = 2020 
Frequency = 1 
 [1] 1.424851 1.491022 1.557193 1.623364 1.689535 1.755706 1.821877 1.888048 1.954220 2.020391

$se
Time Series:
Start = 2011 
End = 2020 
Frequency = 1 
 [1] 1.306001 1.846964 2.262060 2.612002 2.920307 3.199036 3.455353 3.693928 3.918002 4.129937

【讨论】：