# NOT RUN {
data("Andropogon_shape", package = "bayesLopod")
ld_Shape = shapeLopodData(Shapefile = Andropogon_shape, fieldN = "sampEffort",
fieldY = "detections", Adjacency = TRUE, keepFields = FALSE)
mLopodShape = modelLopod(LopodData = ld_Shape, varP = TRUE, q = NULL,
pmin = 0, CAR = TRUE, nChains = 4,warmup = 500,sampling = 100,nCores =4)
psiShape_95 = lopodShape(mLopodShape, "psi_i", extrapolate = FALSE, quant = 0.95)
psiShape_05 = lopodShape(mLopodShape, "psi_i", extrapolate = TRUE, quant = 0.05)
#Visualize results
sp::spplot(psiShape_05, zcol = "psi_i")
sp::spplot(psiShape_95, zcol = "psi_i")
# }
Run the code above in your browser using DataLab