在 R 中创建联合核密度函数

Question

我有两个向量，A 和 B。我可以在 R 中凭经验创建一个联合核密度函数，以便我有一些来自向量 A 和 B 的函数 f(x,y)？以下示例不起作用，因为它不是联合概率分布

z <- cbind(A,B)
approxfun(density(z))

Answer 1

MASS::kde2d 在给定两个向量（x 和 y 坐标）的情况下进行二维核密度估计。但是，它不是返回可以在任意位置 ({newx,newy}) 进行评估的函数，而是返回在方形网格上评估的函数。

一旦您完成了诸如选择带宽之类的繁琐工作，在单个点 x0,y0 处进行内核密度估计的实际计算并不难，我认为它会是这样的

sum(dnorm((x0-x)/h)*dnorm((y0-y)/h)

MASS::kde2d 用outer() 和tcrossprod() 巧妙地计算了从所有数据点到评估网格上所有点的距离以及所有总和，在少数顶层操作，但我认为我上面所说的是它的症结所在。