data(grave)
out = spscan.test(grave, case = 2, alpha = 0.1, nsim = 49)
plot(out, chars = c(1, 20), main = "most likely cluster",
border = "orange", ccol = NA)
# change color, lty, lwd of circles
set.seed(2)
out2 = spscan.test(grave, case = 2, alpha = 0.8, nsim = 49)
plot(out2, chars = c(1, 20), border = "blue")
plot(out2, chars = c(1, 20), border = c("blue", "orange"),
clwd = c(3, 2), clty = c(2, 3))
Run the code above in your browser using DataLab