current code works with ggplot 3.3.2, but not with 3.3.3 a fix by aaron is to add the fixed length of land types: + scale_shape_manual(values = c(1:18))