# NOT RUN {
library(raster)
# create a raster object with NA
data <- runif(100,0,11)
data <- matrix(data, ncol=10)
data <- raster::raster(data)
# replace NA with zero in raster object
data <- calc(data,intFun.theta)
plot(data); text(data)
# }
Run the code above in your browser using DataLab