Я использую модель GEE, чтобы увидеть пространственную автокорреляцию.Но когда я использую функцию correlog() в пакете ncf в R, он говорит: Ответ многомерный: корреляционная матрица будет центрирована на нуле.Используйте correlog.nc() для нецентрированной коррелограммы
Затем я переключаюсь на correlog.nc(), и он говорит:
Ошибка в if (moran [1] <0) {:пропущенное значение, где требуется TRUE / FALSE </p>
Я проверил все переменные в моей модели GEE, и пропущенное значение отсутствует.Может ли кто-нибудь любезно сказать мне, где я ошибся?
Ранее я использовал обычный glm() с распределением Пуассона, и такой проблемы не было.И я мог бы также построить графики пространственных коррелограмм.
Большое спасибо за любую помощь!
Мой код такой:
mgee = geeglm(death0814 ~ pop_dst + 5 SES variables + ndvi2011 + FAD + offset(logpop) + factor(age) + factor(sex), family=poisson, data=data_corr, id=fTPU, corstr="ar1")
rp<-resid(mgee,type="pearson")
correlog1 <- correlog.nc(data_corr$X, data_corr$Y, rp, na.rm=T, increment=1, resamp=0)