How to combine two shapefiles using a common attribute in R?

Question

I have two shapefiles sf1 and sf2 for same area. sf1 is population data for year 2015 and has 6 variables.

sf2 is popul