一、数据概述
本数据集提供了1998年至2023年中国各地级市的地方生产总值及政府竞争程度指标面板数据。
该数据基于《财政研究》2019年第4期论文《经济分权、地方政府竞争与城市全要素生产率》
中提出的政府竞争程度衡量方法，采用地区经济增长率与空间权重矩阵的乘积（Wgdpg）
来衡量地方政府竞争程度。
二、理论背景
1. 指标来源
   - 文献出处：《财政研究》2019年第4期
   - 论文标题：《经济分权、地方政府竞争与城市全要素生产率》
   - 核心指标：Wgdpg（地方政府竞争程度指标）
2. 计算方法
   - 政府竞争程度 = 空间权重矩阵(W) × 地区经济增长率(gdpg)
   - 空间权重矩阵基于城市间的地理邻接关系构建
   - 地区经济增长率使用地区生产总值的年度增长率衡量
3. 指标含义
   - Wgdpg反映了周边地区经济增长对本地区的竞争压力
   - 该指标越高，表明周边地区经济增长越快，本地区面临的竞争压力越大
   - 可用于研究地方政府间的"标尺竞争"、"晋升锦标赛"等问题
三、数据文件说明
本数据集包含以下两个核心数据文件：
【文件1】1998~2023年各城市政府竞争程度指标面板数据.dta
  - 文件格式：Stata数据文件（.dta）
  - 观测数量：7201条
  - 时间范围：1998年～2023年
  - 变量列表：

变量名          变量说明                           
年份             数据年份（1998-2023）            
省代码          省级行政区划代码                  
省                省级行政区名称                     
市代码          地级市行政区划代码                 
市                地级市名称                        
Wgdpg         政府竞争程度指标                    
【文件2】1998～2023年各地级市地区生产总值面板数据(填补后).dta
  - 文件格式：Stata数据文件（.dta）
  - 观测数量：7346条
  - 时间范围：1998年～2023年
  - 变量列表：
变量名                                        变量说明                       
省                                               省级行政区名称                 
省代码                                         省级行政区划代码               
市                                               地级市名称                     
市代码                                         地级市行政区划代码            
年份                                            数据年份（1998-2023）         
地区生产总值_当年价格_亿元_全市   地区生产总值（当年价格，亿元）
地区生产总值增长率_百分比_全市     地区生产总值增长率（%）        

四、算法逻辑说明
Wgdpg指标的计算分为以下五个步骤：
【步骤1】数据预处理 - 填补缺失的GDP增长率
  (1) 读取中国城市统计年鉴地级市面板数据
  (2) 保留关键变量：省、省代码、市、市代码、年份、地区生产总值、GDP增长率
  (3) 删除2017年数据（该年数据存在问题），用单独处理的2017年数据替换
  (4) 创建城市唯一标识符（格式：省_市）
  (5) 设定面板数据结构，对缺失的GDP增长率进行填补：
填补公式：增长率 = (当年GDP - 上年GDP) / 当年GDP × 100
【步骤2】处理城市接壤矩阵
  (1) 读取城市接壤矩阵数据（2019年各市是否接壤矩阵）
  (2) 矩阵结构：n×n 的方阵，n为城市数量
  (3) 矩阵元素：
      - W[i,j] = 1 表示城市i与城市j地理上接壤
      - W[i,j] = 0 表示城市i与城市j不接壤
      - 对角线元素均为0（城市不与自身接壤）
【步骤3】单年份计算过程（以2021年为例）
  (1) 筛选当年数据，删除GDP增长率缺失的观测
  (2) 将GDP数据与城市接壤矩阵进行匹配合并
  (3) 仅保留两个数据集中都存在的城市（确保矩阵维度一致）
  (4) 重新排列矩阵，确保行列顺序一致
  (5) 构建计算所需的矩阵：
      - m1：GDP增长率向量（n×1维列向量）
      - m2：城市接壤矩阵（n×n维方阵）
  (6) 执行矩阵乘法计算Wgdpg：   
      Wgdpg = W × gdpg
      即：Wgdpg = Σ(W[i,j] × gdpg[j])，对所有j求和
      结果含义：城市i的Wgdpg等于其所有接壤城市GDP增长率之和
【步骤4】循环计算所有年份
  (1) 创建结果数据存储文件夹
  (2) 循环处理1998-2023年每一年的数据
  (3) 对每一年执行步骤3的计算过程
  (4) 将每年的计算结果保存为单独的数据文件
【步骤5】合并生成最终面板数据
  (1) 读取1998年数据作为基础
  (2) 循环追加1999-2023年数据
  (3) 添加年份变量，调整变量顺序
  (4) 保存最终的面板数据集
【计算示例】
假设有3个城市A、B、C，其接壤关系和GDP增长率如下：
接壤矩阵W：          GDP增长率向量gdpg：
     A  B  C              A: 5%
A [  0  1  0  ]           B: 8%
B [  1  0  1  ]           C: 6%
C [  0  1  0  ]
计算结果：
Wgdpg_A = 0×5 + 1×8 + 0×6 = 8  （A只与B接壤，所以等于B的增长率）
Wgdpg_B = 1×5 + 0×8 + 1×6 = 11 （B与A、C都接壤，等于A和C增长率之和）
Wgdpg_C = 0×5 + 1×8 + 0×6 = 8  （C只与B接壤，所以等于B的增长率）

五、相关代码
本数据集附带了计算政府竞争程度指标的相关代码，位于"相关代码"文件夹中：
1. main.do    - Stata代码，用于计算Wgdpg指标
2. main.py    - Python代码，提供替代计算方案
辅助数据文件：
- 2019年各市是否接壤矩阵（对角线为0）.dta - 空间权重矩阵数据
- 其他辅助数据文件
六、绑图展示
"绑图"文件夹中包含了数据可视化的相关代码和地图数据：
1. main.do / main.py - 绑图代码
2. chinacity2019mini_*.dta - 中国地级市地图数据
3. chinaprov2019mini_*.dta - 中国省级地图数据
4. nicecut.ado / nicecut.py - 分组切分辅助程序
可绑制的图形包括：
- 各城市地区生产总值增长率分布图
- 各城市政府竞争程度指标分布图