我使用预测软件包制作了一个互相关图。
有很多滞后/信号很重要,考虑到我的范围是-120到120,我想找出它们的滞后数。在预测中是否有某种函数可以创建输出这会使列表有点像这样:
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")
答案 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