# NOT RUN {
file <- system.file('external/dem_example.grd',package='elsa')
r <- raster(file)
plot(r,main='a continuous raster map')
mo <- lisa(r,d2=2000,statistic='i') # local moran's I (Z.Ii value)
plot(mo, main="local Moran's I (Z.Ii)")
mo <- lisa(r,d2=2000,statistic='i',mi='I') # local moran's I (Ii value (non-standardized))
plot(mo, main="local Moran's I (Ii))")
gc <- lisa(r,d2=2000,statistic='c') # local Geary's c
plot(gc, main="local Geary's c")
g <- lisa(r,d2=2000,statistic='g') # local G
plot(g, main="local G")
# }
Run the code above in your browser using DataLab