目前我正在尝试做一些空间统计。目前,我正在玩Getis和Ord的当地空间统计数据来指示"热点"和" coldspots"我的数据集。 gem $ avginc是每个城市的平均收入,gem.listw.20km是我的空间权重矩阵。
a<-0
for(i in 1:2379) {a[i] <- abs[i]}
abs <- as.data.frame(a)
abs$id <- gem_codes$code
abs <- left_join(abs,coords,by=c("id"="id"))
ggplot(abs,aes(long,lat)) +
geom_polygon(aes(group=id,fill=a)) +
scale_fill_gradient(low = "white", high = "darkred")
第一个问题:结果我得到了来自班级&#34; localG&#34;的东西。任何人都可以向我解释如何将其转换为普通整数或data.frame?或者这门课的优势在哪里?
到目前为止,这么好。因为我不知道另一种方法,所以我将其转换为数据帧并用ggplot绘制。
awk
从情节中可以看出,有两个城市不在他们应该的位置。在暗红色的中间,如图中的那个,应该不可能达到如此低的(因此是白色的)值。看起来好像混淆了,但我不知道怎么做。我试着使用&#34; spChk&#34; localG函数中的参数,但我没有得到它如何使用它。
我希望任何人都可以提供帮助!