library(DRHotNet)
library(spatstat.geom)
library(spatstat.linnet)
library(spdep)
library(raster)
rel_assault <- relpnet(X = chicago,
lixel_length = 50, h = 50, mark = "marks", category_mark = "assault")
hotspots_assault <- drhot(X = chicago, rel_probs = rel_assault,
k = 0.5, n = 4)
plothot(X = chicago, hotspots = hotspots_assault)
Run the code above in your browser using DataLab