R语言三重嵌套for循环打印列表与变量名不符合预期问题排查

问题背景 #

尝试用三重嵌套for循环完成以下操作：

1. 遍历数据框列表
2. 对每个数据框，遍历treatment列的唯一值
3. 对每个唯一treatment值，对4个变量（QY_max、FvFo、NPQ_Lss、QY_Lss）执行shapiro.test正态性检验

但运行自定义函数后遇到三个问题：

• 所有参数名连续打印，未按「处理组→参数→检验结果」的顺序输出
• 检验结果中数据标识为df_treat[[parameter]]，而非具体数据框+变量名
• 用deparse(substitute(df))打印数据框名称时结果混乱

问题分析 #

1. 循环代码块范围错误：最内层for循环未加花括号，导致只有print(parameter)在循环内，print(shapiro.test(...))会在参数循环结束后一次性执行，打乱输出顺序。
2. 数据标识缺失：shapiro.test仅接收向量输入，无法自动关联数据框和变量名，需要手动添加上下文信息。
3. 数据框名称获取错误：直接遍历列表元素时，df是临时变量，deparse(substitute(df))无法获取列表中数据框的原始名称，需遍历列表的名称索引来获取对应名称。

修正后的代码 #

treatment_normality <- function(df_list) {
  # 遍历数据框列表的名称和对应元素
  for (df_name in names(df_list)) {
    df <- df_list[[df_name]]
    cat("=== 数据框:", df_name, " ===\n")
    
    for (treat in unique(df$treatment)) {
      cat("\n--- 处理组:", treat, " ---\n")
      df_treat <- dplyr::filter(df, treatment == treat)
      
      # 给内层循环加上花括号，确保代码块都在循环内
      for (parameter in c("QY_max", "FvFo", "NPQ_Lss", "QY_Lss")) {
        cat("\n参数:", parameter, "\n")
        # 执行检验并输出，同时添加自定义标识
        test_result <- shapiro.test(df_treat[[parameter]])
        # 修改检验结果的data.name，替换为具体的名称
        test_result$data.name <- paste0(df_name, "$", parameter, " (treatment: ", treat, ")")
        print(test_result)
      }
    }
  }
}

关键修正点说明 #

• 遍历列表名称：用names(df_list)获取数据框的原始名称，再通过df_list[[df_name]]提取对应数据框，解决名称打印混乱的问题。
• 添加循环花括号：给最内层for循环加上{}，确保参数打印和检验结果输出都在循环迭代中执行，保证输出顺序正确。
• 自定义检验结果标识：修改shapiro.test返回结果的data.name字段，替换为「数据框名$变量名 (处理组: X)」的格式，让检验结果的标识更清晰。

测试运行 #

使用提供的可复现示例数据调用函数：

treatment_normality(list_of_dataframes)

输出会按「数据框→处理组→参数→检验结果」的顺序展示，且每个检验结果的data.name会显示具体的上下文信息，比如：

Shapiro-Wilk normality test

data:  df1$QY_max (treatment: A)
W = 0.97942, p-value = 0.9409

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