# NOT RUN {
data(nydf)
coords = with(nydf, cbind(x, y))
out = scan.test(coords = coords, cases = floor(nydf$cases),
pop = nydf$pop, nsim = 49,
longlat = TRUE, alpha = 0.12)
summary(out)
# }
Run the code above in your browser using DataLab