内核密度估计-将图例比例更改为每m²的密度

Question

我已经为R中的一个物种计算了内核密度估计（KDE），为简化起见，我在这里使用spatstat的ants数据集（蚂蚁巢），我会希望有一个带有每个测量单位（例如m²）的蚁巢数量的图例标度。如果我正确理解KDE图例显示的概率密度在0到1之间，该如何以“实际”密度（例如，每m²的点（蚁巢））转换该概率？

enter image description here

这里是示例：

require("spatstat")
data(ants)
dat <- ants

# estimate bandwith
h_cox <- bw.diggle(dat)

# calculate KDE
kd_cox <- density(dat, h_cox, diggle=TRUE, se=TRUE, eps=diff(dat$window$xrange)/500)

# Plot KDE, contours and points
plot(kd_cox$estimate, main="KDE ants bw.diggle")
contour(kd_cox$estimate, labels="", add=TRUE, col=gray(.5)) 
points(dat)

Answer 1

感谢您对示例数据有一个明确的问题。