如何用apply()系列函数批量绘制多列数据的plot()图形

问题描述 #

需要复用相同的plot()结构，以数据框dat的第6到40列作为y值，dat$elevation作为x值，生成多个散点图（每列对应一个图）。尝试的代码如下：

multi.plot <- apply(x = dat[,6:40], plot(dat[,i] ~ dat$elevation, xlim = c(2700,3500), ylab=names(dat)[i]))

收到错误提示：

Error in match.fun(FUN) : argument "FUN" is missing, with no default

疑问：

1. 如何让apply()正确识别并调用plot()函数？
2. 应如何指定apply()要提取y值的列索引？

解答 #

1. 让apply正确调用plot的方法 #

报错核心是apply参数格式错误：apply必须明确指定FUN参数（即要执行的函数），且需要将plot逻辑包装为可接收列数据的函数（匿名函数或自定义函数），同时设置MARGIN=2指定按列处理数据。

方法一：用apply结合匿名函数

# 按列处理第6-40列，执行自定义plot逻辑
apply(dat[, 6:40], MARGIN = 2, FUN = function(y_col) {
  # 获取当前列名称，用于设置y轴标签
  col_name <- colnames(dat)[which(colnames(dat) == deparse(substitute(y_col)))]
  plot(y_col ~ dat$elevation, xlim = c(2700, 3500), ylab = col_name)
})

方法二：用lapply遍历列名（推荐，可读性更高）

# 先获取要绘制的目标列名列表
target_cols <- colnames(dat)[6:40]
# 遍历列名，逐个生成图形
lapply(target_cols, function(col) {
  plot(dat[[col]] ~ dat$elevation, xlim = c(2700, 3500), ylab = col)
})

2. 指定y值列索引的方法 #

• 直接用dat[,6:40]即可选中第6到40列作为y值来源，这部分你的写法是正确的；
• 如果需要按列名筛选，可改为dat[, c("col_name1", "col_name2")]或dat[, grepl("关键词", colnames(dat))]；
• 提前将目标列名存入变量（如上述target_cols），再遍历变量会让代码更易维护，方便后续调整列范围。

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