我正在运行GEE模型以查看空间自相关。
但是,当我在R的ncf包中使用correlog()
函数时,它说:
响应是多元的:相关矩阵将以零为中心。将correlog.nc()
用于非中心相关图
然后我切换到correlog.nc()
,它说:
if(moran [1] <0){中的错误:需要TRUE / FALSE的缺失值
我检查了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)