今天数据中心将要介绍一个独特但正如火如荼使用的地市级面板数据——中国省/地级市夜间灯光数据（1992-2013的面板数据）。

    当前应用最为广泛的全球夜间灯光数据由美国国防气象卫星计划的多个卫星观测所得。Henderson et al在2012年将夜间灯光数据引入经济学的分析框架，夜间灯光亮度的变化率被作为GDP增长率的代理变量，将经济统计数据和夜间灯光数据结合起来更准确地度量经济增长绩效。近年来，夜间灯光数据被中国社科学者运用于经济活动总量测算等领域，改进了经济增长等统计指标的准确性，产生了大量的研究成果。

   以下这个是中国各个省份和各个地级市的夜间灯光数据，可以构造一个面板数据用于实证研究。数据完整期间是1992-2013，这里只展示了1998-2013的两个地级市。

根据年份对以下代表性省份做的时间趋势图中，可以看出上海的夜间灯光确实最亮，而处于中国西部的四川在夜晚就显得暗淡一些。只展示了1998-2013年的数据，完整的为1992-2013总计22年数据。

把画这个图的程序也列出来，大家可以根据数据进行操作演练，对于面板数据这是很实用的技巧。

use lightdata.dta, clear
preserve
collapse (mean) lightmean, by(year)
sort year
save lightmean,replace
restore
preserve
collapse (mean) lightmean if province==11, by(year)
rename lightmean lightmean1a
sort year
save lightmean1a,replace
restore
preserve
collapse (mean) lightmean if province==31, by(year)
rename lightmean lightmean2a
sort year
save lightmean2a,replace
restore
preserve
collapse (mean) lightmean if province==44, by(year)
rename lightmean lightmean3a
sort year
save lightmean3a,replace
restore
preserve
collapse (mean) lightmean if province==33, by(year)
rename lightmean lightmean4a
sort year
save lightmean4a,replace
restore
preserve
collapse (mean) lightmean if province==51, by(year)
rename lightmean lightmean5a
sort year
save lightmean5a,replace
restore

use lightmean,clear
merge year using lightmean1a
drop _merge
sort year
merge year using lightmean2a
drop _merge
sort year
merge year using lightmean3a
drop _merge
sort year
merge year using lightmean4a
drop _merge
sort year
merge year using lightmean5a
drop _merge

twoway(line lightmean year)(line lightmean1a year,lpattern(dash)) ///
(line lightmean2a year,lpattern(dot)lwidth(0.4)color(green)) ///
(line lightmean3a year,lpattern(shortdash) color(blue))(line lightmean4a year,lpattern(solid)color(pink)) ///
(line lightmean5a year,lpattern(dash_dot)color(brown)) ,graphregion(fcolor(white) ///
lcolor(white))plotregion(fcolor(white)) scheme(s2mono) legend(position(1) ////
ring(0)c(1) symxsize(6)textwidth(6)lab(1 "全国")lab(2 "北京")lab(3 "上海") ////
lab(4 "广东")lab(5 "浙江")lab(6 "四川")) xtitle("year")ytitle("灯光数据的时间趋势变化图")xlabel(1998 (2) 2013) ///
ylabel(,nogrid) xscale(r(1998 2013)) name(graph1b, replace)其他所有地级市1992-2013的完整夜间灯光数据欢迎学院师生到本中心查看拷贝。

其他所有地级市1992-2013的完整夜间灯光数据欢迎学院师生到本中心查看拷贝。

本文部分图文来源于微信公众号”计量经济圈“，版权属于原作者或网站。如有版权问题，请与管理员邮箱联系，将立刻进行相应处理。