昨天群里有个小伙伴在找各省的互联网普及率，这个数据可以从中经网上找到，不过数据不是很完整，另外互联网普及率应该小于 1 ，而中经网的数据上福建省部分年份的数据超过了 100%，显然是错误的。于是我就把这个数据简单整理了下。
| 全国 2002～2019 年互联网普及率 |
 |
这幅图可以使用下面的 Stata 代码绘制：
use "各省互联网普及率(2002~2019).dta", clear
tw area 互联网普及率 年份 if 地区 == "全国",
color("gs60") ||
conn 互联网普及率 年份 if 地区 == "全国",
m(o) lp(solid) lc(gs80) mc(gs20) leg(off)
ti("2002～2019 年中国互联网普及率", size(*1.2))
subti("数据来源：中经网", size(*0.95))
caption("绘制：微信公众号 RStata", size(*0.8))
xla(2002(2)2018 2019)
yla(0 "0%" 20 "20%" 40 "40%" 60 "60%")
gr export "2002～2019年中国互联网普及率.png", width(1200) replace
gr export "2002～2019年中国互联网普及率.pdf", replace
|
各个省份的互联网普及率数据缺失比较严重：
| 各省互联网普及率(2002~2019) |
 |
drop if 地区 == "全国"
encode 地区, gen(region)
xtset region 年份
xtline 互联网普及率, byopts(xr yr ///
note("绘制：微信公众号 RStata", size(*0.8)) ///
ti("2002～2019 年全国及各省互联网普及率", size(*1.2)) ///
subti("数据来源：中经网", size(*0.95))) xti("")
gr export "各省互联网普及率(2002~2019).png", width(1200) replace
gr export "各省互联网普及率(2002~2019).pdf", replace
|
最后一个数据齐全的年份是 2016 年：
| 2016年中国各省互联网普及率 |
 |
其中福建省数据缺失正是因为概念的数据超过了 100%，被我删除了。
该地图的绘制代码为：
use china_prov_db.dta, clear
gen prov = substr(省, 1, 6)
save china_prov_db2.dta, replace
use "各省互联网普及率(2002~2019).dta", clear
drop if missing(互联网普及率) | 地区 == "全国" | 年份 != 2016
hist 互联网普及率
gen prov = substr(地区, 1, 6)
merge m:1 prov using china_prov_db2.dta
replace 互联网普及率 = -1 if missing(年份)
grmap 互联网普及率 using china_prov_coord, id(ID)
clmethod(custom) clbreaks(-1 0 40 50 60 80 120)
fcolor("gray" "237 248 251" "178 226 226" "102 194 164" "44 162 95" "0 109 44")
ocolor("gray" ...)
ti("2016 年各省互联网普及率", size(*1.1))
subtitle("数据来源：中经网")
graphr(margin(medium))
osize(vvthin ...)
legend(size(*1.1)
order(2 "无数据" 3 "< 40%"
4 "40%～50%" 5 "50%～60%"
6 "70%～80%" 7 "> 80%")
ti(互联网普及率, size(*0.5)
pos(11)))
xsize(18) ysize(12.5)
caption("绘制：微信公众号 RStata", size(*0.8))
line(data(china_city_line_coord.dta) size(*0.5 ...) color(black) select(keep if _ID <= 6))
label(data(china_prov_label.dta) x(X) y(Y) l(name) size(*0.8))
xsize(13) ysize(10)
polygon(data(polygon) by(value) fcolor(black)
osize(vvthin))
gr export "2016年中国各省互联网普及率.pdf", replace
gr export "2016年中国各省互联网普及率.png", replace width(1200)
|
点击这里跳转到 RStata 短书平台获取附件：各省互联网普及率(2002~2019)
评论