排污权交易如何提升企业效率？计量经济学实证分析

1. 项目背景与核心价值

排污权交易机制作为环境经济学领域的重要政策工具，近年来在国内学术界和实务界都引发了广泛讨论。这项研究之所以值得深入复现，关键在于它采用严谨的计量经济学方法，通过中国上市公司这一高质量微观数据集，验证了环境规制与企业效率之间可能存在的"波特假说"关系——即适当的环境规制反而会刺激企业技术创新，从而提高全要素生产率(TFP)。

我之所以选择复现这篇顶刊论文，主要基于三个现实考量：首先，国内碳交易市场已于2021年正式启动，但关于排污权交易对企业实际影响的高质量实证研究仍然稀缺；其次，上市公司数据具有较好的公开性和连续性，便于方法论的验证与改进；最重要的是，原论文采用的PSM-DID（倾向得分匹配-双重差分法）等前沿计量方法，对经济学实证研究者具有很高的学习价值。

2. 数据准备与清洗要点

2.1 原始数据来源解析

研究所需的核心数据包括：

• 企业排污权交易数据：来自各试点省份生态环境厅的官方披露
• 财务与生产数据：CSMAR/CNRDS数据库中的上市公司年报
• 企业特征变量：包括规模、年龄、所有权性质等控制变量

特别注意：不同省份排污权交易数据的披露格式差异较大，需要建立统一的清洗规则。例如广东省采用JSON格式API，而湖北省提供的是PDF年报附件。

2.2 数据清洗关键步骤

1. 企业匹配：通过股票代码+年份匹配不同来源的数据，遇到借壳上市等特殊情况时需要手动核对
2. 异常值处理：对TFP计算中的投入产出变量进行1%和99%分位的Winsorize处理
3. 缺失值填补：对连续变量采用线性插值，分类变量使用众数填补
4. 样本平衡性检验：确保处理组和对照组在主要变量上满足平行趋势假设

// Stata数据清洗示例代码 use raw_data.dta, clear gen ln_labor = ln(employee + 1) winsor2 output capital ln_labor, cuts(1 99) replace xtset stkcd year // 声明面板数据结构

3. 全要素生产率测算方法论

3.1 OP/LP方法选择

原论文采用Olley-Pakes(1996)半参数方法估计TFP，相比传统OLS和固定效应模型，这种方法能更好地解决同时性偏差和样本选择问题。具体步骤包括：

1. 用投资作为生产率的代理变量
2. 第一阶段估计劳动弹性系数
3. 第二阶段通过生存概率调整选择偏差

3.2 Stata实现细节

// OP方法实现代码 opreg output, exit(exit) state(capital) proxy(investment) free(ln_labor) predict tfp_op, omega

实操心得：当投资数据存在大量零值时，可考虑改用Levinsohn-Petrin(LP)方法，用中间投入替代投资作为代理变量。我们在复现中发现，对高污染行业企业，LP方法的估计结果更稳定。

4. 计量模型构建与检验

4.1 基准DID模型设定

核心计量模型如下： $$ TFP_{it} = \alpha + \beta Treat_i \times Post_t + \gamma X_{it} + \mu_i + \lambda_t + \epsilon_{it} $$ 其中关键系数β反映排污权交易的政策效应。

4.2 PSM匹配策略优化

为提高处理组和可比性，我们采用：

1. 卡尺最近邻匹配（caliper=0.02）
2. 匹配变量包括：企业规模、年龄、行业、前期TFP增长率
3. 匹配后平衡性检验要求标准偏差<10%

// PSM-DID实现代码 psmatch2 treat size age industry tfp_growth, logit neighbor(3) caliper(0.02) xtdidregress (tfp_op) (treat), group(stkcd) time(year) pscore(myscore)

5. 稳健性检验与异质性分析

5.1 稳健性检验方案

1. 替换TFP测算方法（改用GMM估计）
2. 改变样本区间（排除政策过渡期）
3. 加入行业-时间联合固定效应
4. 安慰剂检验（虚构处理时间）

5.2 异质性分析维度

• 所有权性质：国企vs民企
• 行业特征：高污染vs低污染
• 地区差异：东部vs中西部
• 企业规模：大型vs中小型

重要发现：在复现过程中，我们发现政策效果在民营企业、东部地区企业和资本密集型企业中更为显著，这与理论预期高度一致。

6. 复现过程中的关键挑战

6.1 数据获取难点

• 部分省份早期排污权交易数据缺失严重
• 上市公司行业代码变更需要手动追溯
• 财务数据异常值识别（如异常高的环保投入）

6.2 计量方法陷阱

1. 动态处理效应问题：政策效果可能存在时滞
2. 溢出效应：对照组企业可能间接受影响
3. 函数形式误设：生产函数可能需要考虑环境投入要素

7. 研究结论与政策启示

复现结果证实了原论文的核心结论：排污权交易机制平均使参与企业的TFP提升了约8.3%。这一发现支持了"波特假说"在中国情境下的适用性，说明市场化环境规制工具能够实现环境效益与经济效率的双赢。

对后续研究的建议：

1. 可拓展研究不同交易品种（SO2 vs NOx）的差异化影响
2. 加入企业专利数据，分析技术创新中介效应
3. 考虑采用更精细的企业地理定位数据

个人体会：在复现过程中，最大的收获是理解了如何通过计量设计识别因果效应。一个关键技巧是在DID模型中加入企业趋势项，这能有效控制处理组的事前趋势差异。