data:
![在这里插入图片描述](https://img-blog.csdnimg.cn/20210308164300599.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzIzOTI0Njkx,size_16,color_FFFFFF,t_70)
ChinaRD2<-readRDS("gadm36_CHN_1_sp.rds")
names(ChinaRD2)[4]<-c("province_name")
popDataSP<-merge(ChinaRD2, data_sp, by=c("province_name"),all.x=TRUE)
popDataSP<-popDataSP[!is.na(popDataSP$zao_va),]
summary(popDataSP)
coordinates(popDataSP)->coord
ID<-popDataSP@data$province_name
library(maptools)
library(spdep)
sh1<-knn2nb(knearneigh(coord,k=2),row.names=ID)
sh1_nb<-nb2listw(sh1)
soco_nbq <- poly2nb(popDataSP)
soco_nbq_w <- nb2listw(soco_nbq)
wm <- nb2mat(soco_nbq, style='B')
(moran(popDataSP$zao_va,sh1_nb,length(sh1),Szero(sh1_nb)))
moran.test(popDataSP$zao_va,sh1_nb,randomisation=F)
geary.test(popDataSP$zao_va, sh1_nb, randomisation=TRUE, zero.policy=NULL,
alternative="greater", spChk=NULL, adjust.n=TRUE)
oid<-order(ID)
resl<-localmoran(popDataSP$zao_va,sh1_nb,p.adjust.method="bonferroni")
printCoefmat(data.frame(resl[oid,],row.names =popDataSP$province_name),cheek.names=FALSE)
moran.plot (popDataSP$zao_va,sh1_nb,labels =(popDataSP$province_name),pch=19)
moran.plot (popDataSP$zao_va,soco_nbq_w,labels =(popDataSP$province_name),pch=19,xlab = "First Trimester",ylab = "Vitamin A")
moran.mc(popDataSP$zao_va, soco_nbq_w, nsim=99)
![在这里插入图片描述](https://img-blog.csdnimg.cn/20210308164920386.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzIzOTI0Njkx,size_16,color_FFFFFF,t_70#pic_center)
可以画出左上角的图
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)