CHARLS 是一项具备中国大陆 45 岁及以上人群代表性的追踪调查,旨在建设一个高质量的公共微观数据库,采集的信息涵盖社会经济状况和健康状况等多维度的信息,以满足老龄科学研究的需要。
为利用国际上最佳的数据采集方式,并确保研究结果的国际可比性CHARLS 参照包括美国的健康与退休研究(HRS)在内的系列国际老龄调查研究开展调查设计。其全国基线调查于 2011-12 年进行,于 2013 年、2015 年、2018 年和 2020 年分别开展了 4 轮常规问卷的追踪调查,并于 2014 年完成了中国中老年人生命历程调查。为确保样本的代表性,CHARLS 基线调查覆盖了全国 150 个国家/地区、450 个村庄/城市社区,涉及 10,257户家庭的 17,708 人,反映了中国中老年人群的总体情况。2019 年底到 2020 年初,新冠疫情在中国爆发,为及时记录新冠疫情对中国中老年人生活和健康的影响,在 2020 年的第 5 轮调查中增加采集了疫情相关的信息。
关于Charls数据的文章中,我们看到很多文章有中国地图的展示,便于展示一些疾病的分布人员分布等。
我目前专门编写了一个scimp包,目前只有一个charlsmap函数,顾名思义,这个是一个专门用于绘制charls数据地图的函数,但并不是说它只能绘制charls地图,基本中国省市级地图都能绘制,肯定没有其他的专业地图绘制包专业和功能详尽,但是可以快速出图,目的可以帮助新手小白快速出图,节省你的时间。
目前scimp包目前升级到了2.2版,为了方便粉丝发sci地图,增加了省名拼音表示,下面我来演示一下
先导入数据和R包
library(scimap)
library(haven)
data2021<-read_dta('E:/r/yanshi/2.dta') #data2011
names(data2021)
这是我自己提取的一个charls数据,我想研究吸烟者在全国的分布,先画一个常规的
charlsmap(data = data2021,username=username,token=token,sfig = T)
画个省级的
charlsmap(data = data2021,x="smoking",username=username,token=token,type = "province",pro.name = T,
sfig = T)
在参数e.pro.name = T就可以把省名改成拼音
charlsmap(data = data2021,x="smoking",username=username,token=token,type = "province",pro.name = T,e.pro.name = T,
sfig = T)
还可以把显示比例改成百分比
charlsmap(data = data2021,x="smoking",username=username,token=token,type = "province",pro.name = T,e.pro.name = T,
sfig = T,Legend.percentage=T)
更改图例标题
charlsmap(data = data2021,x="smoking",username=username,token=token,type = "province",pro.name = T,e.pro.name = T,
sfig = T,Legend.percentage=T,title = "吸烟比例")
支持三颜色,但是要是没调好比例配色,看起来有点像个非主流
charlsmap(data = data2021,x="smoking",username=username,token=token,type = "province",pro.name = T,e.pro.name = T,
sfig = T,Legend.percentage=T,lowcol="yellow",highcol="red",midcol = "blue")