Как получить список значимых сигналов / лагов на графике ccf, используя пакет прогноза?

Question

Я сделал кросс-корреляционный график с использованием пакета прогноза.

Существует довольно много значимых лагов / сигналов, и я хочу выяснить, с каким номером лага они находятся, учитывая, что мой диапазонот -120 до 120. Есть ли какая-то функция в прогнозе для создания вывода, которое сделает список примерно таким:

В цериях есть эквивалентная функцияпакет называется 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")

Answer 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