如何在 R 中重新组合季节性分解的 stl 组件

Question

我想将季节性分量重新组合到由 strl 分解的时间序列的季节性调整分量。例如：

    library(fpp)
    fit <- stl(elecequip, s.window=7)
    plot(fit)
    sa_fit <- seasadj(fit)

如何将季节性成分合并到sa_fit以便我获得原始季节性数据？

提前谢谢。

Answer 1

为了在 R 中获取不熟悉的（对您来说）类的属性，请使用 attributes .

快速浏览将显示您需要 time.series 属性，如

unadjusted <- fit$time.series[, 1] + sa_fit

然后，您可以使用all.equal检查unadjusted是否等于原始数据

all.equal(unadjusted, elecequip)
## [1] TRUE