使用for循环结合ggdensity按HIV状态批量绘制密度图的问题

错误原因解析 #

1. attempt to apply non-function 报错
   你用 file$toString(names(file[i]))[HIV=="pos"] 的写法不符合R的语法：$ 运算符后只能跟列名的字面量，不能直接嵌套函数调用或变量。R会把 toString(names(file[i])) 当成一个函数去执行，自然报错。正确的动态取列方式是用 [[ ]] 运算符，比如 file[[var_name]][file$HIV == "pos"]。

2. undefined columns selected 报错
   通常是循环时引用了不存在的列，或是列索引/列名拼写错误。比如如果你的循环包含了HIV列本身，或是索引超出了数据集的列范围，都会触发这个错误。

正确实现代码 #

步骤1：加载包并准备数据 #

library(ggpubr)
# 假设你的数据集名为file，HIV状态列名为HIV，取值为pos/neg

步骤2：定义要绘制的变量（排除HIV列） #

# 获取除HIV外的所有变量名，避免循环到非数值变量或状态列
plot_vars <- setdiff(names(file), "HIV")

步骤3：循环绘制分组密度图 #

方式1：按HIV状态分开绘制单组密度图

for (var in plot_vars) {
  # 绘制HIV阳性样本的密度图
  p_pos <- ggdensity(
    data = file[file$HIV == "pos", ],  # 直接筛选阳性样本子集
    x = var,  # 直接传入变量名字符串，ggdensity支持这种写法
    main = paste("Density of", var, "(HIV Positive)"),
    xlab = var,
    fill = "#1f77b4",
    alpha = 0.6
  )
  print(p_pos)  # 循环中必须print才能显示图形
  
  # 绘制HIV阴性样本的密度图
  p_neg <- ggdensity(
    data = file[file$HIV == "neg", ],
    x = var,
    main = paste("Density of", var, "(HIV Negative)"),
    xlab = var,
    fill = "#2ca02c",
    alpha = 0.6
  )
  print(p_neg)
}

方式2：单图展示同一变量的两组密度对比（更直观）

for (var in plot_vars) {
  p <- ggdensity(
    data = file,
    x = var,
    color = "HIV",  # 按HIV状态分组设置线条颜色
    fill = "HIV",   # 按HIV状态分组设置填充色
    main = paste("Density of", var, "by HIV Status"),
    xlab = var,
    alpha = 0.6,
    palette = c("#1f77b4", "#2ca02c")  # 自定义分组颜色
  )
  print(p)
}

关键注意事项 #

• R中动态取列优先用 [[ ]] 而非 $，因为 $ 不支持变量替换，而 [[ ]] 可以直接传入变量形式的列名。
• 筛选样本时尽量用 file[file$HIV == "pos", ] 这种数据框子集的方式，比单独取向量更清晰，也避免环境变量冲突。
• 循环绘制图形时必须用 print() 输出，否则R不会自动显示循环中的ggplot类图形。

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