r - 在 `forecast()` 中包装工作 `tryCatch()` 语句时没有结果

Question

我有大量时间序列需要生成预测。为了自动生成最佳预测，我想应用一些模型，如 auto.arima、ets、(s)naive、神经网络等。不幸的是，当它循环遍历时间序列时，一些模型会失败，从而停止R脚本的执行。为了使它更健壮，我开始使用 tryCatch();我这样做的主要目标是我希望脚本继续运行，而不必捕获错误。执行代码时，tryCatch() 中的 forecast() 无法生成正确的预测。

请在下面找到我遇到的错误的可重现示例。

历史时间序列:

ts <- structure(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 9, 10, 10, 16, 7, 13, 0, 9, 1, 11, 2, 11, 3, 
11, 4, 1, 20, 13, 13, 13, 9, 14, 16, 16, 18, 17, 20, 18, 19, 
16, 16, 16, 15, 14, 27, 24, 35, 8, 18, 21, 20, 19, 22, 18, 21, 
19, 24, 33, 23, 18, 26, 18, 17, 19, 19, 22, 19, 24, 29, 29, 18
), .Tsp = c(2025.25, 2032.16666666667, 12), class = "ts")

以下代码行产生正确的预测:

fcast_arima <- forecast(auto.arima(ts),h=4)
fcast_arima
         Point Forecast    Lo 80    Hi 80    Lo 95    Hi 95
Apr 2032       24.69032 18.57094 30.80971 15.33153 34.04911
May 2032       25.00539 18.84018 31.17061 15.57651 34.43428
Jun 2032       25.32046 19.10975 31.53118 15.82200 34.81893
Jul 2032       25.63554 19.37966 31.89141 16.06800 35.20307

当我将相同的代码包装到 tryCatch 行时，它无法生成预测，请参见下面的示例:

fcast_arima <- tryCatch({ fcast_arima <- forecast(auto.arima(ts),h=4)}, warning = function(warningcondition) {message(warningcondition)}, error = function(errorcondition) {message(errorcondition)}, finally={})
p-value smaller than printed p-value
fcast_arima
NULL

谁能解释为什么在没有 tryCatch() 的情况下工作正常的代码在 tryCatch() 中却不能工作？

Answer 1

我在想有什么东西通过设置 options(warn=-1) 来抑制警告，并生成警告。我可以模拟您在这里看到的行为。

创建一个生成警告并返回 99 但警告被选项抑制的函数:

> foo=function(){o=options()$warn; options(warn=-1);warning("Yikes!");options(warn=o);return(99)}

这似乎运行得很好:

> foo()
[1] 99

我可以在一个简单的 tryCatch 中运行它:

> tryCatch(foo())
[1] 99

但是如果我添加一个 warning 子句，警告就会被捕获:

> tryCatch(foo(),warning=function(w){message(w)})
Yikes!> 

虽然它似乎仍然在没有警告的情况下运行:

> foo()
[1] 99

我不确定解决方案是什么......抱歉......也许让编写试图抑制警告的代码的人重写它以使用 suppressWarnings 代替:

> foo=function(){suppressWarnings({warning("Yikes!")}) ; return(99)}
> foo()
[1] 99
> tryCatch(foo(),warning=function(w){message(w)})
[1] 99