R 时空中的随机点

我有一套枣子堆肥，在几年的死亡道路上。我的数据是坐标x和y。我需要创建一个随机分布的数据作为对照，以便能够将我的数据与随机对照进行比较，看看是否存在死亡热点。我知道可以用R创建一组随机点，但我不知道如何使两者在时间上具有可比性。我如何在几年内将这些积分装箱

再次感谢你

---

Aurora

为了识别热点，最好使用它。请参见下面的示例，了解道路交通事故死亡率。为了识别热点，使用了

等高线

绘图和

kde2d

功能。请参阅下面的代码：

library(lubridate)
mort <- read.csv("https://data.wa.gov/api/views/mcp7-tcwf/rows.csv?accessType=DOWNLOAD", stringsAsFactors = FALSE)

ll <- t(sapply(mort$Location, function(x) na.omit(as.numeric(unlist(strsplit(x, "\\(|\\,| |\\)"))))))
rownames(ll) <- NULL
colnames(ll) <- c("lat", "lon")
mort2 <- cbind(mort, ll)
mort2$Salvage.Date.Time2 <- mdy_hms(mort2$Salvage.Date.Time)
mort2$month <- month(mort2$Salvage.Date.Time2)
mort2$year <- year(mort2$Salvage.Date.Time2)
mort2 <- mort2[mort2$year> 2016, ]

mort3 <- mort2[with(mort2, lat > 45.5 & lat < 49.& lon > -125 & lon < -116), ]
f1 <- with(mort3, kde2d(lat, lon, n = 100))

plot(mort3$lat, mort3$lon, pch = 18, col = "lightblue")
contour(f1, levels  =  c(0.01, 0.05, 0.1, 0.2), add = TRUE, labcex = 1)

库（lubridate）
mort通常使用
runif
、
rexp
、
rnorm
等函数在R中构建随机数据。请查看这些函数。您可以通过在控制台中键入例如
？runif
来获得帮助。请访问“SO！”！在这个社区中，如果您可以发布一些数据作为示例，展示您的尝试和结果，那么您将更容易获得帮助。在这里，您应该使用
poisson
distribution来生成示例数据。非常感谢！！我会尝试所有这些，我会回到你身边，看看我是否做到了！谢谢