使用R中的核密度估计的热点地图
答
spatstat
包使这容易。由于我不知道你的数据是什么样的,我将首先创建一些演示数据。
x <- rnorm(327)
y <- runif(327)
现在数据使用ppp()
功能从spatstat
转换为地理空间格式。
library(spatstat)
dta <- ppp(x, y, window = owin(c(-5, 5), c(-5, 5)))
最后,计算密度,然后绘制结果。
# Compute the density function
dta <- density(dta)
# Plot the density
plot(dta, main = "Density plot of sample dataset")
+0
非常感谢。我试图使用一个不同的数据,它有几个热点,而不仅仅是你的例子中的一个,但我总是得到1个热点。我尝试了MASS库中有3个热点的“间歇泉”数据。为什么它不起作用? –
如果你希望有人真正帮助您完成此,请[提供一个可重复的例子(https://*.com/questions/5963269/how-to-make-a-great -r重现-例子)。例如用你用MASS的geyser数据试过的代码更新你的问题,所以我们可以看到错误是什么。 –