time-series - 尝试对非平稳和季节性时间序列进行预白化

Question

我正在尝试在两个系列都是非静止的地方进行预白化。在这种情况下，我应该这样写吗？

x <- diff(rtc.ts)
y<- diff(precipitation.ts- fitted(Arima(diff(precipitation.ts), model = mod_rtc))
ccf(as.numeric(x), as.numeric(y), 12)

如果我再次处理两个非平稳时间序列，如果我使用 TSA 的 prewhiten 函数，如何编写代码？