使用 ets() 时，为什么 R 没有响应并崩溃？答案

【问题标题】：When ets() is used, why R is not responding and crashes?使用 ets() 时，为什么 R 没有响应并崩溃？
【发布时间】：2021-12-16 01:09:35
【问题描述】：

我正在努力寻找预测特定地区月平均降雨量的最佳模型。到目前为止，我使用了一种季节性的幼稚方法和 SARIMA。但是当尝试运行ets() 时，R 崩溃而没有产生输出。

【问题讨论】：

如果 R 没有响应，etc 调用仍在运行。没有您的数据，很难看出是否存在问题。
@phiver 我该如何在此处附加数据？
要么在你的 github 中加载数据，要么分享源链接。降雨应该来自气象提供者之一，如 noaa、dwd、knmi 或气象局
@phiver github.com/noob-mun/Rainfall-data/blob/main/Copy.xlsx

标签： r time-series forecasting smoothing state-space

【解决方案1】：

我倾向于使用 fable 和 fabletools。预测的跟进。使用包 fpp3 会加载处理 tsibbles、dplyr 和日期对象所需的所有包。

我对您的数据运行任何预测方法没有任何问题。我尝试了寓言和预测并得到相同的结果。请参阅下面的代码。

# load your data
df1 <- readxl::read_excel("datasets/Copy.xlsx")
colnames(df1) <- c("date", "rainfall")


library(fpp3)
fit <- df1 %>% 
  mutate(date = yearmonth(date)) %>%
  as_tsibble() %>% 
  model(ets = ETS(rainfall)) 
  
report(fit)

Series: rainfall 
Model: ETS(M,N,A) 
  Smoothing parameters:
    alpha = 0.002516949 
    gamma = 0.0001065384 

  Initial states:
    l[0]      s[0]     s[-1]     s[-2]    s[-3]    s[-4]    s[-5]    s[-6]     s[-7]     s[-8]     s[-9]    s[-10]
 86.7627 -77.53686 -57.90353 -18.72201 86.57944 150.0896 166.8125 60.45602 -39.25331 -55.94238 -68.85851 -70.52719
    s[-11]
 -75.19377

  sigma^2:  0.1109

     AIC     AICc      BIC 
2797.766 2799.800 2850.708

使用预测：

library(forecast)
fit <- forecast::ets(ts(df1[, 2], frequency = 12))

fit

ETS(M,N,A) 

Call:
 forecast::ets(y = ts(df1[, 2], frequency = 12)) 

  Smoothing parameters:
    alpha = 0.0025 
    gamma = 1e-04 

  Initial states:
    l = 86.7627 
    s = -77.5369 -57.9035 -18.722 86.5794 150.0896 166.8125
           60.456 -39.2533 -55.9424 -68.8585 -70.5272 -75.1938

  sigma:  0.333

     AIC     AICc      BIC 
2797.766 2799.800 2850.708

【讨论】：

同样的问题仍然存在。我尝试运行代码，然后 R 没有响应并崩溃。也许这与我使用的系统有关。运行ets()或hw()时R崩溃的原因是我的系统不支持吗？
@noob，除非你的记忆有问题，否则不应该有任何原因。您是否尝试过在 R 控制台中运行而不使用 Rstudio 或任何其他 ide？span>
我一直在 R 控制台而不是 Rstudio 中运行它。但后来我在 Jupyter Notebook 中试了一下，得到了输出。