scimp包1.7版发布,用于绘制charls加权地图
CHARLS 是一项具备中国大陆 45 岁及以上人群代表性的追踪调查,旨在建设一个高质量的公共微观数据库,采集的信息涵盖社会经济状况和健康状况等多维度的信息,以满足老龄科学研究的需要。
在文章《scimp包发布,轻松秒绘制charls数据文章中的地图》中,咱们已经介绍了scimp包绘制chalrs地图,但是由粉丝提出,地图要加入权重,不然人数太少了,不符合实际情况,我想了一下,确实如此,因此在1.7版加入了权重参数,下面我来演示一下
导入数据和R包
library(scimap)
library(haven)
data2021<-read_dta('E:/r/yanshi/3.dta') #data2011
names(data2021)
数据是我自己提取的charls2011年的基线表数据,我们可以看到有城市,和省的变量分布
如果咱们想了解吸烟患者的分布,不加权的话
charlsmap(data = data2021,x="smoking",username=username,token=token,type = "province")
上图咱们可以看到,虽然按比例显示了吸烟患者的分布,但是一个省的吸烟人数才几百例,显然不符合实际情况,因此咱们可以给他加上一个权重
charlsmap(data = data2021,x="smoking",username=username,token=token,weights = "weight",type = "province")
这样就相对合理多了。
有不少粉丝私信我提醒我南海九段线的问题,非常感谢,我还不至于犯这么明显的错误,我们只要调下坐标,南海九段线就会出来。默认的坐标是xlim = c(-3000000,2700000),ylim = c(-2787082,1654989),
我们只要把Y轴的坐标稍微往下调一点
charlsmap(data = data2021,weights = "weight",type = "province",
xlim = c(-3000000,2700000),ylim = c(-3887082,1654989))
这样完整的南海九段线就出来了。