# NOT RUN {
library(sp)
s <- readRDS(system.file("nigeria/SR_Naija.rds", package = "phyloregion"))
k <- 10
COLOUR <- hcl.colors(k, "RdYlBu")
y <- choropleth(s$SR, k)
plot(s$SR, y)
## To plot and color according to some metric:
plot(s, col = COLOUR[y])
# }
Run the code above in your browser using DataLab