if (requireNamespace("splancs")) {
data("imdepi")
imdepiB <- subset(imdepi, type == "B")
mainpoly <- coordinates(imdepiB$W@polygons[[1]]@Polygons[[5]])
SGRID <- c(10, 25, 50, 100, 150)
TGRID <- c(1, 7, 14, 21)
B <- 19 # limited here for speed
if (!interactive()) B <- 4
imdBstKtest <- stKtest(imdepiB, eps.s = SGRID, eps.t = TGRID, B = B,
cores = 2, seed = 1, poly = list(mainpoly))
print(imdBstKtest)
plot(imdBstKtest)
}
Run the code above in your browser using DataLab