R localG get';身份证有问题吗

R localG get';身份证有问题吗,r,geospatial,R,Geospatial,目前我正在尝试做一些空间统计。目前,我正在使用Getis和Ord的本地空间统计数据来表示数据集的“热点”和“冷点”。gem$avginc是每个城市的平均收入,gem.listw.20km是我的空间权重矩阵 abs <- localG(gem$avginc,gem.listw.20km,zero.policy=F) abs您可以尝试使用class函数强制更改类:class(abs)不幸的是,这不起作用…R返回的消息是什么?R将其转换为data.frame,但其中没有观察值。localG对象

目前我正在尝试做一些空间统计。目前,我正在使用Getis和Ord的本地空间统计数据来表示数据集的“热点”和“冷点”。gem$avginc是每个城市的平均收入,gem.listw.20km是我的空间权重矩阵

abs <- localG(gem$avginc,gem.listw.20km,zero.policy=F)

abs您可以尝试使用class函数强制更改类:
class(abs)不幸的是,这不起作用…R返回的消息是什么?R将其转换为data.frame,但其中没有观察值。localG对象可能是一个列表,比data.frame更通用。使用str()函数查看对象中包含的内容。可能有更多的信息,而不仅仅是数据。但是,列表中的一个项目将包含数据。您可以使用$符号提取此项。您可以尝试使用class函数强制更改类:
class(abs)。不幸的是,这不起作用……R返回的消息是什么?R将其转换为data.frame,但其中没有观察值。localG对象可能是一个列表,比data.frame更通用。使用str()函数查看对象中包含的内容。可能有更多的信息,而不仅仅是数据。但是,列表中的一个项目将包含数据。您可以使用$符号提取此项目。
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")