for循环生成按Species分面的数值列带密度线直方图报错如何解决？

问题背景 #

需求：为Iris数据集的所有数值列（如Sepal.Length）生成带有密度线的直方图，并按Species列分面。

执行两段循环代码后均触发报错：

Error in `f()`:
StatBin requires a continuous x variable: the x variable is discrete. Perhaps you want stat="count"?

尝试代码1：

for (i in colnames(subset(iris, select = -`Species`))) {
plot=  ggplot(iris, aes(x= i))+
    geom_histogram()+
    geom_density(colour = "blue", size = 1) +
    facet_wrap(~ Species, scales = "free")

 print(plot)
}

尝试代码2：

for (i in colnames(subset(iris, select = -`Species`))) {
plot=  ggplot(subset(iris, select = -`Species`), aes(x= i))+
    geom_histogram()+
    geom_density(colour = "blue", size = 1) +
    facet_wrap(~ iris$Species, scales = "free")

 print(plot)
}

疑问：是否需要在geom_histogram()命令中添加参数解决该问题？

问题分析与解决 #

不是geom_histogram()参数的问题，核心错误是直接在aes(x = i)中传入字符串列名，ggplot会将其识别为离散变量，而我们需要的是将数值列作为连续变量传入。

修正后代码 #

library(ggplot2)

for (i in colnames(subset(iris, select = -Species))) {
  plot <- ggplot(iris, aes(x = .data[[i]], y = ..density..)) +
    geom_histogram(fill = "lightgray", color = "black") +
    geom_density(colour = "blue", size = 1) +
    facet_wrap(~ Species, scales = "free") +
    labs(x = i, title = paste("直方图+密度线：", i))
  
  print(plot)
}

关键调整说明 #

• 使用.data[[i]]引用循环中的列名，这是ggplot官方推荐的动态列引用方式，能正确识别列的数值类型，避免离散变量错误
• 给geom_histogram添加y = ..density..参数，将直方图的高度转换为密度值，和密度线的y轴尺度统一，解决两者比例失调的问题
• 新增填充色、边框和标题，优化可视化效果

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