if (FALSE) {
plot_mapbox(res_mn)
plot_mapbox(res_mn, color = ~INDRESNAME)
map_data("world", "canada") %>%
  group_by(group) %>%
  plot_mapbox(x = ~long, y = ~lat) %>%
  add_polygons() %>%
  layout(
    mapbox = list(
      center = list(lat = ~median(lat), lon = ~median(long))
    )
  )
}
Run the code above in your browser using DataLab