如何在R中按分组列循环计算多列间的cor.test相关性？

针对你的需求——按Result.Test.Mnemonic分组，计算Run-in、Control、TRE三列两两之间的相关性，以下是两种可靠的实现方式，可解决你之前遇到的字符串/非数值输入问题：

方法一：使用Tidyverse工具链（推荐） #

这种方法代码更整洁，结果输出为结构化数据框，方便后续分析：

library(dplyr)
library(tidyr)
library(purrr)
library(broom)

# 定义需要计算相关性的列
cor_target_cols <- c("Run-in", "Control", "TRE")
# 生成所有两两列组合
column_pairs <- combn(cor_target_cols, 2, simplify = FALSE) %>%
  map(set_names, c("x_col", "y_col"))

# 分组计算并整理结果
correlation_results <- change %>%
  group_by(Result.Test.Mnemonic) %>%
  group_modify(function(group_data, group_info) {
    # 对每一对列计算相关性
    map_dfr(column_pairs, function(pair) {
      # 通过列名字符串提取数值向量
      x_values <- group_data[[pair$x_col]]
      y_values <- group_data[[pair$y_col]]
      
      # 执行相关性检验，自动过滤缺失值
      test_result <- cor.test(x_values, y_values, use = "complete.obs")
      
      # 将检验结果转为结构化数据框，并添加列对信息
      tidy(test_result) %>%
        mutate(x_column = pair$x_col, y_column = pair$y_col, .before = 1)
    })
  }) %>%
  ungroup()

代码说明： #

• 用group_modify对每个分组的数据单独处理，避免全局循环的混乱
• 通过[[列名字符串]]直接提取数据框中的数值向量，这是解决"非数值输入"错误的核心
• broom::tidy()将cor.test()的结果转为标准化数据框，便于查看和后续操作

方法二：基础R循环实现 #

如果你更习惯基础R语法，可采用以下循环方式：

# 定义目标列和两两组合
cor_cols <- c("Run-in", "Control", "TRE")
col_pairs <- combn(cor_cols, 2)
# 获取所有分组名称
all_groups <- unique(change$Result.Test.Mnemonic)
# 初始化结果存储列表
result_list <- list()

# 遍历每个分组
for (group_name in all_groups) {
  # 提取当前分组的数据
  current_data <- change[change$Result.Test.Mnemonic == group_name, ]
  
  # 遍历每一对列
  for (pair_idx in 1:ncol(col_pairs)) {
    x_col <- col_pairs[1, pair_idx]
    y_col <- col_pairs[2, pair_idx]
    
    # 提取数值向量
    x <- current_data[[x_col]]
    y <- current_data[[y_col]]
    
    # 计算相关性
    cor_test <- cor.test(x, y, use = "complete.obs")
    
    # 整理结果并添加到列表
    result_list[[paste(group_name, x_col, y_col, sep = "_")]] <- data.frame(
      Group = group_name,
      X_Column = x_col,
      Y_Column = y_col,
      Correlation_Coeff = cor_test$estimate,
      P_Value = cor_test$p.value,
      CI_Lower = cor_test$conf.int[1],
      CI_Upper = cor_test$conf.int[2]
    )
  }
}

# 将列表合并为最终数据框
final_results <- do.call(rbind, result_list)

你之前出错的原因 #

• 直接在cor.test()中传入paste0生成的字符串：cor.test需要的是数值向量，而非列名字符串，因此会报错"非数值类型"
• 用paste0生成执行字符串但未执行：即使生成了类似cor.test(change$Run-in, change$Control)的字符串，也需要用eval(parse(text=...))来执行，但这种方式易出错且代码可读性差，不推荐使用

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