r - 从密度对象中提取概率

Question

我目前正在为我正在进行的生态项目构建密度对象。 density函数在 R非常适合将密度函数拟合到我的数据中。

我如何使用 density 的示例功能如下图所示:

dens.iris <- density(iris$Sepal.Length, bw = "bcv")

这非常有效。然而， predict功能似乎不适用于 density对象。有谁知道提取特定点的密度值的方法(例如，在 iris 数据集中，提取 Sepal.Length 为 6.432)？重要的是我使用 biased cross validation为此的技术。

Answer 1

您可以使用 approxfun函数在密度结果给出的点之间线性插值。因此你可以使用

diris <- with(dens.iris, approxfun(x, y, rule=1))
diris(6.432)
# [1] 0.349344

还

curve(diris(x), from=4.0, to=7.9)

当然，您必须记住，密度曲线的值与概率不同。与任何连续分布一样，萼片长度恰好为 6.432 的概率为 0。