st_combine por variáveis ​​r

states %>%
  group_by(group_var) %>% 
  summarize(geometry = st_union(geometry))
marcofaga