如何从R语言生成的相关矩阵中提取指定相关系数？

以下提供两种适配不同需求的实现方案，可直接适配300个变量的大规模相关矩阵场景：

方法1：单个提取指定变量的相关系数 #

如果你只需要查询少量指定变量对的相关系数，直接按行列名索引即可：

# 构建数据集并计算相关矩阵
df <- data.frame(a = c(2, 3, 3, 5, 6, 9, 14, 15, 19, 21, 22, 23),
                 b = c(23, 24, 24, 23, 17, 28, 38, 34, 35, 39, 41, 43),
                 c = c(13, 14, 14, 14, 15, 17, 18, 19, 22, 20, 24, 26),
                 d = c(6, 6, 7, 8, 8, 8, 7, 6, 5, 3, 3, 2))
cor_mat <- cor(df[, c('a', 'b', 'c')])

# 提取a和b的相关系数
cor_ab <- cor_mat["a", "b"]
# 按要求格式输出
cat("a和b的相关系数为：", cor_ab, "\n")

方法2：批量提取所有非重复两两相关系数 #

如果需要输出全部变量对的相关系数，可通过上三角矩阵索引自动去重，无需手动列配对：

# 获取上三角矩阵索引（自动排除对角线的自相关系数、避免正反两次重复配对）
upper_ind <- which(upper.tri(cor_mat), arr.ind = TRUE)

# 拼接配对信息和对应相关系数
result <- data.frame(
  var1 = rownames(cor_mat)[upper_ind[, "row"]],
  var2 = colnames(cor_mat)[upper_ind[, "col"]],
  cor_value = cor_mat[upper_ind]
)

# 批量输出要求格式的结果
for (i in 1:nrow(result)) {
  cat(i, ". ", result$var1[i], "和", result$var2[i], "的相关系数为：", result$cor_value[i], "\n", sep = "")
}

运行以上代码将直接输出你需要的格式：

1. a和b的相关系数为：0.9279869
2. a和c的相关系数为：0.9604329
3. b和c的相关系数为：0.8942139

补充：批量提取指定部分配对的相关系数 #

如果仅需要从大规模矩阵中提取部分指定配对，可先构造目标配对表再批量索引：

# 构造你需要提取的变量配对表
target_pairs <- data.frame(var1 = c("a", "b"), var2 = c("c", "c"))
# 批量提取对应相关系数
target_pairs$cor_value <- mapply(function(x, y) cor_mat[x, y], target_pairs$var1, target_pairs$var2)

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