如何使用预测包获取ccf图中重要信号/滞后的列表?

时间:2018-12-20 18:54:57

标签: r forecasting cross-correlation

我使用预测软件包制作了一个互相关图。

有很多滞后/信号很重要,考虑到我的范围是-120到120,我想找出它们的滞后数。在预测中是否有某种函数可以创建输出这会使列表有点像这样: enter image description here

tseries软件包中有一个等效功能,称为ccfvalues,但我发现tseries不如预测灵活。

以下是该图的代码:

Ccf(diff(Site1V2OxygenWoNA$Oxygen), diff(Site1V2TempWoNA$Temp), type = "correlation", calc.ci = TRUE, level = 95, lag.max = 118, plot = TRUE, na.action = na.interp, main = "Site 1 OT")

enter image description here

1 个答案:

答案 0 :(得分:1)

是的,您想要的是acf元素:

out <- ccf(diff(Site1V2OxygenWoNA$Oxygen), diff(Site1V2TempWoNA$Temp), type = "correlation", calc.ci = TRUE, level = 95, lag.max = 118, plot = TRUE, na.action = na.interp, main = "Site 1 OT")$acf[, , 1]
names(out) <- -floor(length(out) / 2):floor(length(out) / 2) # Adding lag names