R kde2d 二維核密度估計

說明

使用 axis-aligned 二元正態核進行二維核密度估計，在方形網格上進行評估。

用法

kde2d(x, y, h, n = 25, lims = c(range(x), range(y)))

參數

x	數據的 x 坐標
y	數據的 y 坐標
h	x 和 y 方向的帶寬向量。默認為正常參考帶寬(請參閱bandwidth.nrd)。將采用標量值應用於兩個方向。
n	每個方向上的網格點數。可以是標量或長度為 2 的整數向量。
lims	網格覆蓋的矩形的範圍為 c(xl, xu, yl, yu) 。

值

三個組件的列表。

x , y	網格點的 x 和 y 坐標，長度為 n 的向量。
z	估計密度的 n[1] by n[2] 矩陣：行對應於 x 的值，列對應於 y 的值。

例子

attach(geyser)
plot(duration, waiting, xlim = c(0.5,6), ylim = c(40,100))
f1 <- kde2d(duration, waiting, n = 50, lims = c(0.5, 6, 40, 100))
image(f1, zlim = c(0, 0.05))
f2 <- kde2d(duration, waiting, n = 50, lims = c(0.5, 6, 40, 100),
            h = c(width.SJ(duration), width.SJ(waiting)) )
image(f2, zlim = c(0, 0.05))
persp(f2, phi = 30, theta = 20, d = 5)

plot(duration[-272], duration[-1], xlim = c(0.5, 6),
     ylim = c(1, 6),xlab = "previous duration", ylab = "duration")
f1 <- kde2d(duration[-272], duration[-1],
            h = rep(1.5, 2), n = 50, lims = c(0.5, 6, 0.5, 6))
contour(f1, xlab = "previous duration",
        ylab = "duration", levels  =  c(0.05, 0.1, 0.2, 0.4) )
f1 <- kde2d(duration[-272], duration[-1],
            h = rep(0.6, 2), n = 50, lims = c(0.5, 6, 0.5, 6))
contour(f1, xlab = "previous duration",
        ylab = "duration", levels  =  c(0.05, 0.1, 0.2, 0.4) )
f1 <- kde2d(duration[-272], duration[-1],
            h = rep(0.4, 2), n = 50, lims = c(0.5, 6, 0.5, 6))
contour(f1, xlab = "previous duration",
        ylab = "duration", levels  =  c(0.05, 0.1, 0.2, 0.4) )
detach("geyser")

參考

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.