Use ggplot2::facet_wrap().
library(sf)
#> Linking to GEOS 3.9.1, GDAL 3.2.1, PROJ 7.2.1
library(ggplot2)
load(url("https://github.com/mgimond/Spatial/raw/main/Data/Sample1.RData"))
p=ggplot(data = s.sf, aes(fill = Income)) +
geom_sf(col = "red")
p+ facet_wrap(.~ NAME, ncol = 4)

Created on 2021-12-28 by the reprex package (v2.0.1)