我想将季节性成分重新组合为由 STL 分解的时间序列的季节性调整成分。例如:
library(fpp)
fit <- stl(elecequip, s.window=7)
plot(fit)
sa_fit <- seasadj(fit)
如何将季节性分量合并到 sa_fit 中，以便获取原始季节性数据？
提前致谢。
最佳答案
为了获取 R 中(对您而言)不熟悉的类的属性，请使用 attributes。
快速浏览一下就会发现您需要 time.series 属性，如下所示
unadjusted <- fit$time.series[, 1] + sa_fit
然后，您可以使用 all.equal 检查未调整是否等于原始数据
all.equal(unadjusted, elecequip)
## [1] TRUE
关于r - 如何在 R 中重新组合季节性分解的 STL 组件，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/27468366/