R 密度图：如何添加从 x 轴到密度曲线顶部的实线

Question

我有一个密度图，使用以下方法绘制：

plot(density(x))

我感兴趣的是创建一条从 x 轴到曲线上相应点的线，例如 x = 5。

像这样：

Answer 1

为此，您可以首先将密度值存储在对象中，然后从该对象中检索 x 和 y 元素。在以下示例中，我使用 findInterval 检索给定 x 值的 y 值：

x <- rnorm(1000) # Sample data
y <- density(x)  # Calculate and store density

x0 <- 2 # Desired position on x-axis
y0 <- y$y[findInterval(x0, y$x)] # Corresponding y value

plot(density(x))
segments(x0, 0, x0, y0)

Answer 2

如果它确实是 z 分数，那么只需实际绘制密度函数 dnorm() 即可。您似乎还希望实际的 x 轴位于 0 处。

zmax <- 4
curve(dnorm, -zmax, zmax, xaxt = 'n', bty = 'n')
axis(1, -zmax:zmax, pos = 0)

要绘制线条，您可以再次使用 dnorm 函数。

zscore <- 1.65
segments(zscore, 0, zscore, dnorm(zscore))

您甚至还可以很好地遮蔽它...:)

x <- seq(zscore, zmax, 0.01)
y <- c(0, dnorm(x)[1:(length(x)-2)],0)
polygon(x,y, density = 20)

您还可以使用segments和text命令将标签放在您的图形上，显示阴影和未阴影的内容地区的意思。