我有两个向量 S 和 V,使用函数 kde2d,我得到了它们联合密度的下面这张图: 使用这些数据,是否可以获得联合概率的经验估计,形式为P(S[i],V[j])?在问题如何从R中的密度函数找到/估计概率密度函数中建议我们使用approxfun来获取1D KDE图中值的高度。是否有一种方法将此想法扩展到2维?
一种方法是使用由kde2d返回的网格的双线性插值:library(fields) points <- data.frame(x=0:2, y=c(0, 5, 5)) interp.surface(k, points) # [1] 0.066104795 0.040191482 0.001943069 数据: library(MASS) set.seed(144) x <- rnorm(1000) y <- 5*x + rnorm(1000) k <- kde2d(x, y)