从 R 中的 3D 内核密度估计和提取点值

Question

从 R 文档link 中的以下代码中，我想弄清楚如何为数据帧的每个原始点提取单个值（“密度”） 这样，我就可以将新列用于我想到的 3D 空间中的“点”密度的进一步应用。

我查看了这个 SO 问题：link，但该代码给出了一个 51 层的矩阵。我不知道如何获得 150 个点值。

在我的情况下，我担心使用的三元示例：

fhat <- kde(x=iris[,1:3])

这是代码，它给出了一个很好的情节，但我无法为 iris 数据集中的 150 行/点中的每一个获取一个值

library(ks)
library(MASS)
data(iris)

## univariate example
fhat <- kde(x=iris[,2])
plot(fhat, cont=50, col.cont="blue", cont.lwd=2, xlab="Sepal length")

## bivariate example
fhat <- kde(x=iris[,2:3])
plot(fhat, display="filled.contour2", cont=seq(10,90,by=10))
plot(fhat, display="persp", thin=3, border=1, col="white")

## trivariate example
fhat <- kde(x=iris[,2:4])
plot(fhat, drawpoints=TRUE)

Answer 1

我假设您指的是 ks 包中的 kde()？

如果是这样，请使用参数 eval.points 使用估计值计算点向量/矩阵的密度：

## univariate example

library(ks)

fhat <- kde(x = iris[,2])
plot(fhat, cont=50, col.cont="blue", cont.lwd=2, xlab="Sepal length")

estimate <- kde(iris[,2], eval.points = iris[,2])$estimate
points(iris[,2], estimate)

类似地适用于更高的输入维度。