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)