I am following the example from here: https://www.r-spatial.org/r/2018/10/25/ggplot2-sf.html
How can I overlay a shape file on top of this that has a boundary box?