R：如何获得特定的密度估计值

Question

# The Old Faithful geyser data
d <- density(faithful$eruptions, bw = "sj")
> head(d$x)
[1] 1.179869 1.188363 1.196857 1.205350 1.213844 1.222338

我在density中使用{stats}函数，我想知道是否可以在输出中看到特定值的密度？例如，目前，我在[1] 1.179869 1.188363 ...的喷发值处有密度估算值，但如果我想知道喷发值1 2 5 10 ...的密度估算值，该怎么办？有没有办法提取这些density对象，d？

Answer 1

如果我理解正确，你想要x值等于某个数字的概率（在我的解决方案中是3或4）？

d <- density(faithful$eruptions, bw = "sj")
densityDF <- data.frame(xVals = d$x, prob = d$y)
densityDF$xVals <- round(densityDF$xVals)

densitySearch <- densityDF[densityDF$xVals %in% c(3,4),]

结果：

   xVals       prob
157     3 0.11229482
158     3 0.10721410
159     3 0.10230912
160     3 0.09765156
161     3 0.09318662
162     3 0.08891621