我正在使用一个Shapefile(包含311个元素的大型SpatialPoygonsDataFrame)在R中工作,这是1970年印度的地区 map .我可以通过数据列b_all(特定于地区的回归系数)的值在splot中绘制一张 map ,使用:
brk1 = c(-130, -5, 0, 5, 311)
spplot(distVDS_betas,zcol="b_all",at=brk1)
我可以通过另一个数据列t_all(每个回归系数的t值)绘制类似的图:
brk2 = c(0, 1.645, 1.96, 53)
spplot(distVDS_betas,zcol="t_all",at=brk2)
我想要做的是(A)用divergent的 colored颜色 模式绘制beta系数(比如橙红色上升fr/零,蓝色紫色下降fr/零),然后(B)用t值使 colored颜色 饱和.有人能帮我做这件事吗?
我知道ssplot2已经过时了,我应该转到ggplot2,我很想学习如何使用它.已try 按照代码here进行操作.但是,当我try 运行此代码时:
ggplot(data = distVDS_betas) +
geom_sf()
我得到的错误是:stat_sf() requires the following missing aesthetics: geometry.
所以我甚至不确定如何从ggplot2+sf开始.