Question

我想将季节性成分重新组合到季节性调整后的成分中，这些时间序列由stl分解。例如：

    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

然后，您可以使用unadjusted

检查all.equal是否等于原始数据

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