Корректировка некоторых соседних пространственных точек
Пожалуйста, найдите коды R, прилагаемые к настоящему документу.
library(maptools)
library(spdep)
library(leaflet)
library(RColorBrewer)
#Read polygon shapefile from the directory
chi.poly <- readShapePoly('Immunization.shp')
class(chi.poly)
str(slot(chi.poly,"data"))
summary(chi.poly@data$bcg)
plot(chi.poly)
chi.ols<-lm(bcg~ib+pnc+full_anc, data=chi.poly@data)
summary(chi.ols)
list.queen<-poly2nb(chi.poly, queen=TRUE)
W<-nb2listw(list.queen, style="W", zero.policy=T)
W
moran.lm<-lm.morantest(chi.ols, W, alternative="two.sided")
print(moran.lm)
** Основная проблема заключается в том, что у районов есть несколько соседних пространственных точек. Таким образом, когда я пытаюсь вычислить индекс Морана I через сгенерированную пространственную весовую матрицу, он показывает ошибку соседних областей.
Я не понимаю, как настроить пространственную весовую матрицу в случае, если у некоторых пространственных точек нет соседей.
Пожалуйста, помогите мне. заранее спасибо