我有2个时间序列,正在使用ccf
查找它们之间的互相关。ccf(ts1, ts2)
列出所有时滞的互相关。如何在不手动查看数据的情况下找到导致最大相关性的滞后?
最佳答案
发布答案http://r.789695.n4.nabble.com/ccf-function-td2288257.html
Find_Max_CCF<- function(a,b)
{
d <- ccf(a, b, plot = FALSE)
cor = d$acf[,,1]
lag = d$lag[,,1]
res = data.frame(cor,lag)
res_max = res[which.max(res$cor),]
return(res_max)
}