st_coordinates()
sf
inla()
ggplot() + geom_sf(data = nc_subset) + geom_sf(data = nc_points) + theme_bw()