具有用于季节性分解的海函数的自动绘图产生错误下标越界
autoplot with seas function for seasonal decomposition produces error subscript out of bounds
我正在尝试使用 seasonal::seas() 函数对我的时间序列数据进行季节性分解并绘制结果。将 autoplot() 函数包装在 seas() 周围通常效果很好，但对于下面的数据，我收到了一个错误。我看到一些关于季节性分解的下标越界错误的帖子，但 none 已经解决了这个问题。 seas() 函数有警告 Model used in SEATS is different: (1 1 2) 但我不知道为什么会导致错误。如果您能帮助解决此问题，我们将不胜感激。
代码：
library(forecast)
library(seasonal)
data <- c(11666, 18526, 14955, 13791, 13508, 15682, 23039, 24190, 23628, 23411, 25068, 20204, 16084, 15726, 13677, 14668, 17915, 25673, 30529, 31990, 31334, 30129, 24323, 19584, 20827, 19766, 21886, 25119, 31832, 31032, 26647, 28725, 26160, 22611, 20698, 21711, 20014, 21698, 27068, 27181, 23442, 25894, 22482, 21894, 20365, 20078, 20995, 22959, 25745, 33074, 29875, 27740, 29499, 27302, 19854, 19658, 20850, 20130, 21415, 25808, 35240, 31326, 31128, 29594, 25660, 26359, 25723, 26995, 25283, 26823)
time_series <- ts(data, freq=12, start=c(2013, 10))
autoplot(seas(time_series))
产生错误：
Model used in SEATS is different: (1 1 2)
Error in `[.default`(object$data, , c("trend", "seasonal", "irregular")) :
subscript out of bounds
发生这种情况是因为分解没有季节性成分，autoplot() 假设它有。
我现在已经更新了包的开发版本 (https://github.com/robjhyndman/forecast) 以允许这种情况。
我正在尝试使用 seasonal::seas() 函数对我的时间序列数据进行季节性分解并绘制结果。将 autoplot() 函数包装在 seas() 周围通常效果很好，但对于下面的数据，我收到了一个错误。我看到一些关于季节性分解的下标越界错误的帖子，但 none 已经解决了这个问题。 seas() 函数有警告 Model used in SEATS is different: (1 1 2) 但我不知道为什么会导致错误。如果您能帮助解决此问题，我们将不胜感激。
代码：
library(forecast)
library(seasonal)
data <- c(11666, 18526, 14955, 13791, 13508, 15682, 23039, 24190, 23628, 23411, 25068, 20204, 16084, 15726, 13677, 14668, 17915, 25673, 30529, 31990, 31334, 30129, 24323, 19584, 20827, 19766, 21886, 25119, 31832, 31032, 26647, 28725, 26160, 22611, 20698, 21711, 20014, 21698, 27068, 27181, 23442, 25894, 22482, 21894, 20365, 20078, 20995, 22959, 25745, 33074, 29875, 27740, 29499, 27302, 19854, 19658, 20850, 20130, 21415, 25808, 35240, 31326, 31128, 29594, 25660, 26359, 25723, 26995, 25283, 26823)
time_series <- ts(data, freq=12, start=c(2013, 10))
autoplot(seas(time_series))
产生错误：
Model used in SEATS is different: (1 1 2)
Error in `[.default`(object$data, , c("trend", "seasonal", "irregular")) :
subscript out of bounds
发生这种情况是因为分解没有季节性成分，autoplot() 假设它有。
我现在已经更新了包的开发版本 (https://github.com/robjhyndman/forecast) 以允许这种情况。