如何基于R中的相邻多边形填充缺少的值？

假设以下数据集并添加一些缺少的值（仅用于说明）：

这给了我一个每行相邻县的指数列表。现在我想用相邻多边形的平均面积来填充缺少的面积值。我将如何使用这些索引获取对应行的平均值

index <- st_touches(nc, nc)

output <- nc %>% 
  mutate(AREA = ifelse(is.na(AREA),
                             apply(index, 1, function(i){mean(.$AREA[i])}),
                             AREA))

output$AREA[c(30, 45)]

[1] 0.1510 0.1335

手动查找区域

mean(output$AREA[index[[30]]])

[1] 0.151

mean(output$AREA[index[[45]]])

[1] 0.1335

---

nc看起来像什么？

index <- st_touches(nc, nc)

output <- nc %>% 
  mutate(AREA = ifelse(is.na(AREA),
                             apply(index, 1, function(i){mean(.$AREA[i])}),
                             AREA))

output$AREA[c(30, 45)]

[1] 0.1510 0.1335

index[c(30, 45)]

[[1]]
[1] 13 14 29 37 48

[[2]]
[1] 44 87

mean(output$AREA[index[[30]]])

[1] 0.151

mean(output$AREA[index[[45]]])

[1] 0.1335