changing ggmap projection

Hello,

I am trying to create a study area map using ggmap and geom_sf. I originally had some issues but the ideas in [this](https://stackoverflow.com/a/50844502/3653140) post helped line things up. But now I am struggling with the geographic projection of ggmap. The map I am trying to produce is large (>1,000,000 km2) and northern (Canada), so the distortion from Mercator is apparent. I was hoping to transform the projection to a projected coordinate system, such as Albers Equal Area. But can not figure out how to transform the ggmap object. Is this possible?

I tried transforming the ggplot object, but coord_map doesn't work with geom_sf, and coord_sf(crs=st_crs(3005)) didn't seem to do what I wanted (nothing happened).

Thanks