之前给大家分享过历年日度夜间灯光数据的爬取方法，以及基于日度夜间灯光栅格数据处理的省市区县面板数据：
- 名师讲堂｜使用 R 语言下载和处理日度夜间灯光栅格数据：https://rstata.duanshu.com/#/brief/course/270178a4d6bb4ee2ba8b66ae71de3b8f
- 2012年1月19日～2024年12月31日中国各省市区县日度夜间灯光亮度面板数据：https://rstata.duanshu.com/#/course/49a7c949a7b04ef2b2006b46841f9139
今天给大家分享一份基于这个数据计算的 2012～2020 年税调企业周边夜间灯光亮度及是否加班面板数据。
数据概览
经过计算得到了 2012年1月19日～2020年12月31日税调企业周边夜间灯光亮度及是否加班面板数据（分年份），以 2020 年为例，数据概览如下：
这个数据分别根据时间和空间维度标准判断了加班情况，变量名称的含义一目了然。
为了便于大家使用，我还提供了 2012～2020 年税调企业年加班天数面板数据，数据概览如下：
就是根据 2012年1月19日～2020年12月31日税调企业周边夜间灯光亮度及是否加班面板数据统计的每个公司每年加班的天数，如果想计算比例，除以该年的总天数即可。
注意：该夜间灯光亮度的单位是 0.1 nWatts/(cm^2 sr) 如果想转换成 nWatts/(cm^2 sr)，需要把数值乘以 0.1 。
下图展示了各年税调企业的平均加班天数：
除此之外，大家还可以计算一些拓展指标，例如亮度差值：
计算方法
该数据的计算方法参考了附件中的文献：「超时加班与劳动收入份额：基于卫星夜间灯光的经验证据」：
论文中提供了两个判断税调企业是否加班的标准：
- 时间维度标准：该税调企业处的夜间灯光亮度大于该年节假日夜间灯光亮度的中位数；
- 空间维度标准：该税调企业处的夜间灯光亮度大于周边8个网格的中位数；
- 时空双维度：同时满足上面的时间和空间标准的日子才被视为加班日。
法定节假日数据之前也有给大家分享过：
2007~2024 年法定节假日及星期数据：https://rstata.duanshu.com/#/brief/course/63147c1f06924ddb89b152d1df84b78c
栅格数据的本质是一个个的格子：
高分辨率的栅格数据，这个夜间灯光数据的分辨率是 500mx500m的，很难看出来这种本质，但是低分辨率的栅格数据大家一看就明白了，例如这个栅格数据：
来源于课程「使用 R 语言绘制历年中国省市区县地图（小地图版本+长版）」：https://rstata.duanshu.com/#/course/379d19770956478ebc4d919910fca74c
每个税调企业都会落在其中的某个格子里面，例如 2020 年 12 月 31 日这天的某税调企业：
代码
为方便大家学习，附件中还提供了数据处理和绘图代码以供参考：
数据引用格式
由于该数据包含较多 RStata 处理的内容，在研究中使用该数据请使用清晰的方式注明数据来源于 RStata 或者 RStata 数据中心，并使用如下格式引用：
RStata 数据中心: 2012年1月19日～2020年12月31日税调企业周边夜间灯光亮度及是否加班面板数据. 2025. https://tidyfriday.cn/rsdb2/
英文文献可以使用下面的格式引用：
RStata Data Center: Panel Data on Nighttime Light Brightness in the Vicinity of and Overtime Work Indicator for Tax Surveyed Enterprises, January 19, 2012 – December 31, 2020. 2025. https://tidyfriday.cn/rsdb2/
点击这里跳转到 RStata 短书平台获取附件：2012年1月19日～2020年12月31日税调企业周边夜间灯光亮度及是否加班面板数据
评论