# NOT RUN {
library(raster)
r <- raster(nrows=824, ncols=767, xmn=2451905, xmx=3218905,
ymn=-2744771, ymx=-1920771, resolution = 5000)
r[] <- runif(ncell(r), 0, 1)
# Modified z-score
z <- raster.Zscore(r)
# P-value
p <- raster.Zscore(r, p.value = TRUE)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab