如何用R将数据框列表合并为完整的相关系数矩阵？

基础R实现方式 #

我们可以通过以下步骤构建目标相关系数矩阵：

1. 提取所有变量并初始化全0矩阵
   先从所有数据框中提取唯一变量名称并排序，以此作为矩阵的行名和列名，同时创建一个全0的初始矩阵：

# 提取所有唯一变量并排序
all_vars <- sort(unique(unlist(lapply(myList, function(x) x$name))))
# 初始化全0相关矩阵
corMatrix <- matrix(0, 
                    nrow = length(all_vars), 
                    ncol = length(all_vars), 
                    dimnames = list(all_vars, all_vars))

2. 确定每个数据框对应的主变量
   每个数据框的第一个name值就是该数据框对应的主变量（比如A1对应a1），提取这些主变量：

row_vars <- sapply(myList, function(x) x$name[1])

3. 填充相关系数到矩阵中
   循环遍历每个数据框，找到主变量在矩阵中的行位置、关联变量的列位置，填充对应系数，同时利用相关矩阵的对称性填充对称位置：

for (i in seq_along(myList)) {
  current_df <- myList[[i]]
  # 找到当前主变量的行索引
  row_idx <- match(row_vars[i], all_vars)
  # 找到当前数据框中变量的列索引
  col_idx <- match(current_df$name, all_vars)
  # 填充行方向系数
  corMatrix[row_idx, col_idx] <- current_df$cor
  # 填充对称的列方向系数
  corMatrix[col_idx, row_idx] <- current_df$cor
}

# 保留两位小数，与示例格式一致
corMatrix <- round(corMatrix, 2)

运行上述代码后，corMatrix就会生成你需要的结果。

tidyverse 实现方式 #

如果你习惯使用tidyverse工具链，可以用更简洁的代码完成：

library(tidyverse)

# 给每个数据框添加主变量列
df_list <- map2(myList, sapply(myList, function(x) x$name[1]), 
                ~ mutate(.x, row_var = .y))

# 合并所有数据框并转换为宽格式，缺失值用0填充
wide_df <- bind_rows(df_list) %>%
  pivot_wider(names_from = name, values_from = cor, values_fill = 0) %>%
  column_to_rownames("row_var")

# 转换为矩阵并确保对称性（避免重复填充对角线）
corMatrix <- as.matrix(wide_df)
corMatrix <- corMatrix + t(corMatrix) - diag(diag(corMatrix))
corMatrix <- round(corMatrix, 2)

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