I want to plot the polygons of this shapefile without the colors (just white with black borders); I also do not want the legend. I have tried scale_fill_manual and a few other things but to no avail. I would also like to include the labels (i.e. 901, 902 ....943) within each polygon. Any assistance would be greatly appreciated.
statsf2134 %>%
ggplot()+
geom_sf(aes(fill=Statistica, border = "black"))

As already suggested by @jdobres in the comments set your desired fill color as a parameter outside of
aes()instead of mapping a variable on thefillaesthetic insideaes(). Also note that, there is noborderaesthetic or parameter. If you want to switch the border color use thecolor=parameter. Finally, you could add labels at the center of each polygon usinggeom_sf_text.Using some fake example data based on the default example from
?geom_sf: