# NOT RUN {
data("loaloa")
plot(loaloa, main="loaloa villages")
# elevation
plot(elevationLoa, col=terrain.colors(100), main="elevation")
points(loaloa)
# vegetation index
plot(eviLoa, main="evi")
points(loaloa)
plot(tempLoa, main="temperature")
points(loaloa)
# land type, a categorical variable
plot(ltLoa)
if(requireNamespace("mapmisc")){
mapmisc::legendBreaks("bottomleft",ltLoa, bty='n')
}
points(loaloa)
# }
Run the code above in your browser using DataLab