momo的回答
首先想到的是空气污染这块。
基于“地区空气污染水平”数据，从各地年均值来看：
cd C:\Download
use 地理信息_2022, clear
keep P* C*
drop *ODE
duplicates drop
g q = substr(C,1,6)
save A0, replace
*
cd C:\Download
*
forv i=1/35{
import exc 二氧化硫年平均浓度, sheet("Sheet`i'") clear
keep B C E
drop if mi(E)
replace B = substr(B,1,4)
destring B E, force replace
ren (B C E) (year city SO2)
save `i', replace
}
clear
forv i=1/35{
append using `i'
}
drop if mi(y)
compress
save SO2, replace
*
forv i=1/35{
import exc 二氧化氮年平均浓度, sheet("Sheet`i'") clear
keep B C E
drop if mi(E)
replace B = substr(B,1,4)
destring B E, force replace
ren (B C E) (year city NO2)
save `i', replace
}
clear
forv i=1/35{
append using `i'
}
drop if mi(y)
compress
save NO2, replace
*
forv i=1/35{
import exc 可吸入颗粒物年平均浓度, sheet("Sheet`i'") clear
keep B C E
drop if mi(E)
replace B = substr(B,1,4)
destring B E, force replace
ren (B C E) (year city PM10)
save `i', replace
}
clear
forv i=1/35{
append using `i'
}
drop if mi(y)
compress
save PM10, replace
clear
forv i=1/35{
cap erase `i'.dta
}
*
use SO2, clear
merge m:m y c using NO2, nogen
merge m:m y c using PM10, nogen
g AirPollu = (1/3)*(S/0.06+N/0.04+P/0.07)/1000
keep y c A
duplicates drop y c, force
save 0, replace
*
use 0, clear
keep c
duplicates drop
forv i=2004/2022{
g y`i' = .
}
gather y*
drop val
g year = substr(v,2,.)
drop v
destring y, force replace
merge 1:1 c y using 0, nogen
encode c, gen(i)
xtset i y
foreach c of var A{
bys c: ipolate `c' y, g(`c'_n) epolate
drop `c'
ren `c'_n `c'
replace `c' = . if `c' <= 0
bys c: fillmissing `c'
}
drop i
la var A 地区空气污染水平
duplicates drop
drop if mi(A)
save 地区空气污染水平, replace
*
use 地区空气污染水平, clear
g q = substr(c,1,6)
merge m:m q using A0, nogen keep(3)
graph hbar A, asyvars over(P) legend(row(3) ring(6) title("省份")) scale(0.6) blabel(bar) ytitle("")
得到：
数据显示，河北省污染最重，海南省污染最轻。
污染严重的省份，大自然的味道估计也好不了哪去。
数据来源：中国环境统计年鉴。
参考文献：Li, T., Li, Y., An, D., Han, Y., Xu, S., Lu, Z., & Crittenden, J. (2019).Mining of the association rules between industrialization level and air quality to inform high-quality development in China. Journal of Environmental Management, 246, 564–574.
（完）