R：在经度/纬度图上叠加状态图

Question

我在经度和纬度的网格上绘制了邮政编码区域。但是，我想在此图中添加宾夕法尼亚州的轮廓。我不确定该怎么做。这是我到目前为止的代码：

library(tidyverse)
library(zipcode)
data(zipcode)


toplot <- merge(zipcode, nhmzipuse, sort = FALSE, by.x = "zip", by.y = "zip")

ggplot(toplot, aes(longitude, latitude)) + geom_point(aes(color=count))

这是我当前情节的样子：

我知道我需要以某种方式获取宾夕法尼亚州边界的数据并将其与我的数据合并，但是我不确定该怎么做。

任何建议都将不胜感激！