使用ggplot2绘制直方图上的密度

Question

我有一个“长”格式的数据框，它包含两列：第一列值，第二列性别[Male-1/Female-2]。我编写了一些代码来制作整个数据集的直方图（下面的代码）。

ggplot(kz6, aes(x = values)) + 
  geom_histogram()

但是，我还想在直方图上添加一个密度，以强调性别之间的差异，即我想组合3个图：整个数据集的直方图，以及每个性别的2个密度图。我尝试使用一些示例（一、二、三、四），但仍然不起作用。“密度”代码仅起作用，而“历史密度”的组合不起作用。

density <- ggplot(kz6, aes(x = x, fill = factor(sex))) + 
  geom_density()

both <- ggplot(kz6, aes(x = values)) + 
  geom_histogram() +
  geom_density()

both_2 <- ggplot(kz6, aes(x = values)) + 
  geom_histogram() +
  geom_density(aes(x = kz6[kz6$sex == 1,]))

P. S.一些例子包含y=...密度...什么意思？对此如何解读？

Answer 1

要绘制直方图并叠加由分类变量定义的两个密度，请在调用geom_density时使用适当的美学，如group或颜色。

ggplot(kz6, aes(x = values)) +
  geom_histogram(aes(y = ..density..), bins = 20) +
  geom_density(aes(group = sex, colour = sex), adjust = 2)

数据创建代码。

我将从内置数据集iris创建一个测试数据集。

kz6 <- iris[iris$Species != "virginica", 4:5]
kz6$sex <- "M"
kz6$sex[kz6$Species == "versicolor"] <- "F"
kz6$Species <- NULL
names(kz6)[1] <- "values"
head(kz6)