data <- read_nhgis(
ipums_example("nhgis0972_csv.zip"),
verbose = FALSE
)
sf_data <- read_ipums_sf(ipums_example("nhgis0972_shape_small.zip"))
joined_data <- ipums_shape_inner_join(data, sf_data, by = "GISJOIN")
colnames(joined_data)
Run the code above in your browser using DataLab