r - 在R中使用 "density"函数查找新数据点的概率密度

Question

我正在尝试使用R中的“密度”函数查找具有未知分布的连续数据的最佳PDF。现在，给定一个新的数据点，我想基于核密度找到该数据点的概率密度从“密度”函数结果得到的估计量。
我怎样才能做到这一点？

Answer 1

如果您的新点在density生成的值的范围内，则很容易做到-我建议使用approx(或approxfun，如果需要将其作为函数)来处理网格值之间的插值。

这是一个例子:

set.seed(2937107)
x <- rnorm(10,30,3)
dx <- density(x)
xnew <- 32.137
approx(dx$x,dx$y,xout=xnew)

如果我们绘制密度和新点，我们可以看到它正在执行您需要的操作:

如果需要外推新值，则将返回 NA。如果要处理外推，我建议针对该点直接计算KDE(使用您拥有的KDE的带宽)。