在R中为直方图叠加正态曲线却生成直线的技术问询

嘿，我太懂你遇到的这个坑了！你的代码里出现直线的核心原因是直方图和正态曲线的尺度不匹配：

• 默认情况下，hist(g)绘制的是频数直方图（y轴是每个区间的观测数量）
• 而dnorm()输出的是概率密度值（取值范围通常很小，远小于频数），直接画出来就会变成贴在x轴的直线。

下面给你两种靠谱的解决办法，按需选择：

方法1：让直方图使用密度刻度（最常用） #

这种方法让直方图和正态曲线的y轴尺度统一，代码更简洁：

g <- unesco$Infant.Deaths
# 绘制密度直方图，freq=FALSE指定用密度而非频数
hist(g, freq = FALSE, main = "婴儿死亡数分布 + 正态拟合曲线", xlab = "婴儿死亡数")
# 生成覆盖数据范围的x序列，步长越小曲线越平滑
x_seq <- seq(min(g, na.rm = TRUE), max(g, na.rm = TRUE), by = 0.5)
# 绘制正态拟合曲线，注意处理NA值
lines(x_seq, 
      dnorm(x_seq, mean = mean(g, na.rm = TRUE), sd = sd(g, na.rm = TRUE)), 
      col = "blue", 
      lwd = 2) # lwd设置线宽，让曲线更显眼

方法2：保留频数刻度，缩放正态曲线 #

如果你必须保留频数的y轴，可以把密度值转换成频数尺度：

g <- unesco$Infant.Deaths
# 先绘制直方图并保存结果，方便获取组距等信息
hist_obj <- hist(g, main = "婴儿死亡数分布 + 正态拟合曲线", xlab = "婴儿死亡数")
# 计算缩放因子：组距 × 总观测数（把密度转成频数）
scale_factor <- diff(hist_obj$breaks)[1] * length(g)
# 生成x序列
x_seq <- seq(min(g, na.rm = TRUE), max(g, na.rm = TRUE), by = 0.5)
# 缩放后绘制曲线
lines(x_seq, 
      dnorm(x_seq, mean = mean(g, na.rm = TRUE), sd = sd(g, na.rm = TRUE)) * scale_factor, 
      col = "blue", 
      lwd = 2)

额外注意点 #

• 一定要加上na.rm = TRUE：如果你的数据里有缺失值，mean()和sd()会返回NA，导致曲线画不出来
• 不要硬编码x序列范围（比如你之前的seq(0,200,by=5)），用数据的实际最值能让曲线覆盖更完整的分布

内容的提问来源于stack exchange，提问作者R.N