Using sf and tidyverse to get weighted mean of neighbors for all columns

gistfile1.txt

library(tidyverse)
library(sf)

st_queen <- function(a, b = a){
  st_relate(a, b, pattern = "F***T****")
}

new_data <- original_data %>%
  mutate(neighbors = st_queen(.)) %>% 
  st_drop_geometry() %>%
  select(unit_id, neighbors) %>%
  tidyr::unnest(neighbors) %>%
  mutate(neighbor_id = original_data$unit_id[neighbors]) %>%
  left_join(original_data, by = c("neighbor_id"="unit_id")) %>%
  select(-neighbors, -neighbor_id) %>%
  group_by(unit_id) %>% 
  summarize_all(~mean(.)) %>%
  left_join(original_data, by = "unit_id")