如何计算数据框中所有列对的滚动相关性并生成结果表格？

以下是用R语言实现需求的完整方案，基于你提到的rollapply函数（来自zoo包）：

步骤1：加载依赖包与准备数据 #

首先加载所需工具包，然后构造你的示例数据：

library(zoo)
library(dplyr)

# 构造示例数据框
df <- data.frame(
  A = c(1, 7),
  B = c(8, 5),
  C = c(3, 5),
  D = c(4, 7)
)

步骤2：生成所有列对组合 #

用combn函数自动生成所有2列组合，并拼接成目标列名（如AB、AC）：

# 获取所有两列组合的列名
col_pairs <- combn(names(df), 2, FUN = paste0, collapse = "")

步骤3：定义滚动相关性计算函数 #

这个函数会接收窗口内的子数据框，计算指定两列的相关系数：

roll_cor <- function(window_data, col1, col2) {
  # use = "complete.obs" 用于处理窗口内的缺失值
  cor(window_data[, col1], window_data[, col2], use = "complete.obs")
}

步骤4：批量计算滚动相关性 #

遍历所有列对，用rollapply计算滚动相关，注意设置by.column = FALSE以传递整窗口数据：

# 计算每对列的滚动相关性，这里窗口宽度设为2（可根据你的需求调整）
roll_results <- lapply(col_pairs, function(pair) {
  # 拆分列对为单个列名
  cols <- strsplit(pair, "")[[1]]
  rollapply(df, width = 2, FUN = roll_cor, col1 = cols[1], col2 = cols[2], by.column = FALSE)
})

步骤5：合并结果为目标数据框 #

将计算结果合并成数据框，并设置对应的列名：

result_df <- do.call(cbind, roll_results) %>% 
  setNames(col_pairs)

# 查看最终结果
print(result_df)

示例输出 #

基于你提供的两行数据，窗口宽度设为2时，输出结果如下：

注意事项 #

• 调整width参数可以改变滚动窗口的大小，如果你的实际数据行数多于窗口大小，会生成对应长度的滚动相关向量
• 如果数据存在缺失值，可以调整cor函数的use参数来适配（比如use = "pairwise.complete.obs"）
• 若你的数据量较大，可考虑用data.table包的frollapply函数提升效率

内容的提问来源于stack exchange，提问作者tpscp