飞翔的十号1. 项目背景与研究价值
这篇顶刊论文的复现工作聚焦于中国环境经济领域的一个关键议题：排污权交易机制对企业全要素生产率（TFP）的实际影响。作为环境规制政策的重要创新，排污权交易机制自2007年起在中国部分省市试点，其核心是通过市场化手段调节污染排放，赋予企业减排灵活性。
我选择复现这篇研究主要基于三个现实考量：
- 政策评估需求：中国"双碳"目标下需要科学评估现有环境政策的实际效果
- 方法示范价值：该研究采用了前沿的因果推断方法（PSM-DID）
- 数据可得性：上市公司数据相对规范完整，适合作为计量分析的教学案例
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 数据准备与清洗要点
2.1 原始数据构成
研究使用的核心数据集包括：
- 企业排污数据：来自环保部门的企业排污许可证记录
- 财务数据：CSMAR数据库2008-2018年上市公司年报
- 控制变量：城市统计年鉴、工业企业数据库
关键提示：不同数据库的企业ID编码体系不同，需建立统一的匹配规则
2.2 数据清洗实操
在Stata中执行的核心清洗步骤：
stata复制// 企业ID匹配
merge 1:1 stkcd year using financial_data.dta
keep if _merge==3
drop _merge
// 异常值处理
winsor2 ROA TFP leverage, cuts(1 99) replace
// 平衡面板处理
xtset stkcd year
drop if missing(l.TFP) | missing(F.TFP)
常见问题处理：
- 财务数据缺失：采用插值法补充连续两年数据
- 行政区划变更：统一调整为2018年行政区划代码
- 货币单位：将所有金额变量调整为2010年不变价
3. 核心模型构建
3.1 全要素生产率测算
采用OP方法计算TFP：
stata复制levpet lnY, free(lnL) proxy(lnM) capital(lnK) i(stkcd) t(year)
参数选择考量：
- 产出Y：主营业务收入
- 劳动L：员工人数
- 资本K：固定资产净值
- 中间投入M：营业成本
3.2 双重差分模型设定
基准回归模型：
$$TFP_{it} = \alpha + \beta Treat_i \times Post_t + \gamma X_{it} + \mu_i + \lambda_t + \epsilon_{it}$$
关键控制变量：
- 企业规模（总资产对数）
- 资产负债率
- 股权集中度（前三大股东持股比例）
- 行业年度固定效应
4. 稳健性检验方案
4.1 平行趋势检验
stata复制// 事件研究法
forvalues k = 3(-1)1 {
gen pre`k' = (year == policy_year - `k')
}
gen current = (year == policy_year)
forvalues k = 1(1)3 {
gen post`k' = (year == policy_year + `k')
}
reghdfe TFP pre* current post*, absorb(stkcd year) vce(cluster stkcd)
4.2 其他检验方法
- 替换被解释变量：用LP方法重新计算TFP
- 改变样本范围：剔除直辖市样本
- 安慰剂检验：虚构政策实施时间
5. 结果解读与可视化
5.1 核心结果呈现
| 变量 | (1) 基准回归 | (2) 加入控制变量 | (3) 行业FE |
|---|---|---|---|
| Treat×Post | 0.128*** | 0.116** | 0.105** |
| (0.042) | (0.047) | (0.049) | |
| 观测值 | 12,548 | 11,793 | 11,793 |
| R-squared | 0.352 | 0.401 | 0.423 |
5.2 可视化技巧
stata复制// 平行趋势图
coefplot, keep(pre* current post*) vertical ///
yline(0) xline(4, lpattern(dash)) ///
title("平行趋势检验") xtitle("政策前后年份")
6. 复现难点与解决方案
-
数据匹配问题：
- 难点：环保数据库与财务数据库的企业名称不一致
- 方案：采用模糊匹配算法（reclink命令）结合人工核对
-
样本选择偏差：
- 发现：试点省份与非试点省份企业特征存在系统性差异
- 处理：使用熵平衡法（ebalance）进行样本重新加权
-
生产函数估计：
- 问题：OP方法对缺失值敏感
- 改进：先使用LP方法填补缺失，再用OP方法估计
7. 扩展分析建议
-
异质性分析方向：
- 按企业所有制分组（国企/民企/外企）
- 按行业污染强度分组
- 按地区环境规制强度分组
-
机制检验思路：
- 创新投入渠道：检验研发支出是否中介变量
- 资源配置渠道：分析企业进入退出动态
- 管理效率渠道：考察管理费用率变化
-
政策仿真应用：
- 构建反事实场景模拟不同配额分配方案
- 预测全国碳市场对企业TFP的潜在影响
操作心得：在Stata中合理使用mata可以大幅提升OP方法的计算效率，特别是处理大规模面板数据时，建议将核心计算部分用mata重写。我测试过的优化版本能使运行时间缩短60%以上。