的R - 獲取從二維核密度估計

Question

聯合概率我有兩個向量S和V，並使用功能kde2d，我得到他們的聯合密度的以下情節：

利用這些數據，是否有可能以P（S [i]，V [j]）的形式獲得聯合概率的經驗估計？

在問題How to find/estimate probability density function from density function in R中，建議我們使用approxfun來獲取1D KDE圖中值的高度。有沒有辦法將這個想法擴展到2維？

Answer 1

一種方法是使用bilinear interpolation由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)