if (interactive() &&
bru_safe_inla() &&
bru_safe_sp() &&
require("sp") &&
require(ggplot2, quietly = TRUE) &&
requireNamespace("terra")) {
# plot all the nests, mesh and boundary
ggplot() +
gg(gorillas_sf$mesh) +
geom_sf(
data = gorillas_sf$boundary,
alpha = 0.1, fill = "blue"
) +
geom_sf(data = gorillas_sf$nests)
# Plot the elevation covariate
gorillas_sf$gcov <- terra::rast(
system.file(gorillas_sf$gcov_file, package = "inlabru")
)
plot(gorillas_sf$gcov$elevation)
# Plot the plot sample
ggplot() +
geom_sf(data = gorillas_sf$plotsample$plots) +
geom_sf(data = gorillas_sf$plotsample$nests)
}
if (FALSE) {
gorillas_sf$gcov <- gorillas_sf_gcov()
}
Run the code above in your browser using DataLab