I am creating a map with different provinces and have 4 point layers over this map each from their own source. When I add the layers to create the map, I only get a legend for the provinces of the map, not for the rest of the layers.
Province %>% group_by(ProvinceNa) %>%
summarize() %>%
ggplot() +
geom_sf(aes(fill = ProvinceNa)) +
scale_fill_grey() +
geom_sf(data = It, col = "yellow") +
geom_sf(data = Io, col = "blue") +
geom_sf(data = Ib, col = "brown") +
geom_sf(data = Ie, col = "green")

