世界经济 2025 年 1 月刊论文「田利辉 等 - 企业上云对其创新的影响：数据要素化的视角」中提到了使用雷电频率作为企业上云的工具变量。类似的，该数据也可以用作其他数字经济基础设施相关研究的工具变量。
今天给大家分享的就是这个工具变量数据：1995～2023年各省市区县雷击频率、雷击频率与长途光缆线路增长率交乘面板数据。
雷击频率的原始数据来源于 https://cmr.earthdata.nasa.gov/search/concepts/C1995863244-GHRC_DAAC.html。下载的到的是 netCDF 格式的 nc 文件，里面比较有用的是 HRFC_COM_FR、HRFC_LIS_FR 和 HRFC_OTD_FR 三个变量，分别表示：
- HRFC_LIS_FR：热带降雨测量任务（TRMM）卫星上的闪电成像传感器（LIS）监测得到的全球闪电总闪光频率；
- HRFC_OTD_FR：MicroLab-1卫星搭载的光学瞬变探测器（OTD）监测得到的全球闪电总闪光频率；
- HRFC_COM_FR：根据上面两个来源的数据合成的更优质的闪电频率结果。
这些指标记录的是 1995-05-04 到 2014-12-31 间的平均结果，是个截面数据。使用各省市区县矢量数据对这个栅格数据（nc 文件在 R 语言中可以使用 terra 包读取为栅格数据）进行分区域平均就可以得到各省市区县闪电频率的截面数据了。直接这些截面数据纳入模型回归会与企业固定效应产生共线性。
基于这部分的计算得到了下面三个文件：
- 各城市闪电频率数据.dta
- 各区县闪电频率数据.dta
- 各省份闪电频率数据.dta
以城市数据为例（数据里面有小数点，预览图是程序生成的，为了统一性只保留了整数部分）：
里面的变量含义就是上面介绍的含义。
按照参考文献的介绍，可以将雷击频率与全国长途光缆线路增长率（g）交乘，得到 “雷击频率×g”这一具有时变性质的交互项，并将其作为企业上云的工具变量进行估计。
全国长途光缆线路增长率数据可以使用「中国统计年鉴2024年」里面的数据计算得到的。附件中也提供了这个数据：「中国统计年鉴2024-电信主要通信能力.csv」
基于此，又计算了三个数据文件：
- 1995～2023年各省份雷击频率与长途光缆线路增长率交乘面板数据.dta
- 1995～2023年各区县雷击频率与长途光缆线路增长率交乘面板数据.dta
- 1995～2023年各城市雷击频率与长途光缆线路增长率交乘面板数据.dta
这些数据文件里面又包含了三个新变量：
- HFC_COM_FR_g = HRFC_COM_FR * 全国长途光缆线路增长率
- HRFC_LIS_FR_g = HRFC_LIS_FR * 全国长途光缆线路增长率
- HRFC_OTD_FR_g = HRFC_OTD_FR * 全国长途光缆线路增长率
这样得到的就是面板数据了。
（数据里面有小数点，预览图是程序生成的，为了统一性只保留了整数部分）：
下图展示了各省市区县的平均雷电频率及交互项的变化趋势：
数据引用格式
在论文中使用该数据请声明数据来源于RStata 数据中心，或者 RStata Data Center，并使用下述的格式引用：
RStata 数据中心: 1995～2023年各省市区县雷击频率、雷击频率与长途光缆线路增长率交乘面板数据. 2025. https://tidyfriday.cn/rsdb2/
英文引用格式：
RStata Data Center: Panel Data on Lightning Strike Frequency by Province, City, and District/County (1995–2023) and Its Interaction with Long-Distance Optical Cable Line Growth Rates. 2025. https://tidyfriday.cn/rsdb2/
另外也建议在说明数据来源的同时，根据本文的介绍在论文中描述数据的详细处理方法。
点击这里跳转到 RStata 短书平台获取附件：1995～2023 年各省市区县雷击频率、雷击频率与长途光缆线路增长率交乘面板数据
评论