互相关帧

Question

我有三个数据帧Tushka，ARQE和ARQW

我使用互相关的代码为每一个数据帧分别这样

# cross correlations with US column 

#Tushka 

head(Tushka) 

# US PW1 PW2 PW3 PW4 

#1 173.62 153.01 144.65 152.53 137.05 

#2 173.57 152.97 144.64 152.52 137.10 

#3 173.52 152.95 144.64 152.52 137.11 

lag = 200 

t3 = data.frame(lag = c(-lag:lag)) 

vars = names(Tushka[, -1]) 

for (j in vars) { 

t4 = ccf(Tushka[, 1], Tushka[, j], lag = lag, na.action = na.pass) 

t3[, j] = data.frame(t4$acf)[, 1] 

} 

但是，如果我创建一个三个数据帧的列表，我如何使用交叉相关的数据帧的列表使用for-loop或任何其他功能

Answer 1

这里是一些可重现的样本数据：

n <- 100 
data_list <- list(
    d1 = data.frame(x = runif(n)), 
    d2 = data.frame(x = runif(n)), 
    d3 = data.frame(x = runif(n)) 
) 

诀窍是使用expand.grid来获得数据帧对。

pairs <- expand.grid(first = names(data_list), second = names(data_list)) 

现在你可以遍历的pairs行，运行互相关的每对行中命名的数据帧。

lapply(
    seq_len(nrow(pairs)), 
    function(i) 
    { 
    first_data <- data_list[[pairs[i, "first"]]] 
    second_data <- data_list[[pairs[i, "second"]]] 
    ccf(first_data$x, second_data$x, plot = FALSE) 
    } 
)