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 

を推定私は{stats}でdensity機能を使用していて、それが出力に特定の値での濃度を確認することが可能です場合、私は思ったんだけど？たとえば、現時点では、噴出値が[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