| 研究问题 | 农业机械化跨区协作(包括收割机跨省调度、移动式烘干服务、精量播种技术推广)是否显著提升了小麦主产区的粮食全要素生产率(TFP)?不同技术类型(收割 vs. 烘干 vs. 播种)的生产率效应是否存在异质性? |
| 理论框架 | 诱致性技术变迁理论(Hayami & Ruttan, 1971);农业发展阶段理论(Timmer, 1988);技术扩散理论(Griliches, 1957) |
| 研究方法 | 随机前沿分析(SFA)测算粮食TFP + 双向固定效应面板模型;以各省农机跨区调度协议签订年份作为准自然实验,构建渐进DID识别跨区协作的因果效应;样本期 2010-2025 年省级面板数据 |
| 数据来源 | 农业农村部农机跨区作业调度数据、国家统计局省级粮食产量与投入数据(劳动力、化肥、机械动力、播种面积)、Wind农业经济数据库、各省农事服务组织登记数据 |
| 创新点 | 首次将农业机械化"跨区协作"作为一个独立的技术扩散变量纳入TFP分析框架,突破了现有文献仅关注单一技术采用(如收割机械化率)的局限;利用跨省协议签订时点的外生冲击识别因果效应,区分收割、烘干、播种三类技术协同的生产率贡献 |
诱致性技术变迁理论(Hayami & Ruttan, 1971)指出,农业技术进步的方向由要素相对稀缺性内生决定——在劳动力成本持续上升的背景下,机械化替代劳动是要素价格信号引导的必然结果。然而,单一地区的机械化水平受制于本地装备存量和作业窗口期,跨区协作通过空间再配置打破了这一约束,使得机械这一"流动资本"在更广地理范围内实现最优配置。技术扩散理论(Griliches, 1957)进一步预测,跨区协作加速了先进机械从高采纳率地区向低采纳率地区的扩散速度,产生技术溢出效应。
据此提出以下假说:
- H1:农业机械化跨区协作程度每提高1个标准差,粮食TFP提升0.3-0.5个百分点,且该效应在收割期窗口较短的地区(如华北平原麦收期多雨区)更为显著。
- H2:移动式烘干服务对TFP的边际贡献大于收割机跨区调度——因为烘干环节直接减少产后损失,属于"止损型"技术,而收割替代仅改变生产方式。
- H3:跨区收割协作与精量播种技术推广存在协同效应——收割效率提升释放的劳动力窗口期使农户更有余力采纳精密播种技术。
核心内生性挑战在于:跨区协作程度本身可能与地区农业发展水平正相关(反向因果)——生产率高的地区更可能主动寻求跨区合作。本文以各省与邻省签订农机互助协议的时间点作为外生冲击,构建渐进DID(staggered DID)识别策略。识别假定为:(1)协议签订的时机主要由当年气候条件(抢收窗口期长短)决定,与长期TFP趋势不直接相关;(2)协议签订前,处理组与对照组在TFP增长率上满足平行趋势。采用Callaway & Sant'Anna(2021)的处理效应估计量,解决传统TWFE在渐进DID中的负权重问题。同时,以"协议签订省份距离邻省粮仓的距离"作为工具变量进行2SLS估计,利用地理外生性增强因果推断。
基准模型设定为渐进DID:
TFP_it = α + β TreatPost_it + γ X_it + μ_i + λ_t + ε_it
其中 TFP_it 为i省在t年的粮食全要素生产率(SFA方法测算),TreatPost_it 为处理指示变量(签订跨区协议后=1),β 为核心待估参数(跨区协作的平均处理效应ATT)。X_it 包括:有效灌溉面积占比、农业劳动力人均受教育年限、化肥投入强度、农业R&D支出占比。μ_i 为省份固定效应(吸收不变的地区特征),λ_t 为年份固定效应(吸收共同时间趋势)。使用Callaway & Sant'Anna(2021)方法估计组别-时间异质性处理效应。
被解释变量粮食TFP采用超越对数随机前沿生产函数(Translog-SFA)测算:以粮食产量为产出,投入要素包括劳动力(农业从业人数 × 粮食播种占比)、机械动力(农机总动力 × 粮食播种占比)、化肥施用量(折纯量 × 粮食播种占比)、播种面积。核心解释变量TreatPost:以农业农村部农机跨区调度协议数据库记录的各省首次与邻省签订互助协议的年份构造处理组虚拟变量。控制变量:有效灌溉面积占比(国家统计局)、农业劳动力平均受教育年限(CGSS微观数据聚合)、化肥折纯施用量/亩(国家统计局)、农业R&D支出占农业GDP比(科技统计年鉴)。样本:13个小麦主产省(黑、吉、辽、冀、鲁、豫、皖、苏、鄂、川、陕、甘、新),2010-2025年面板。
- 替换TFP测算方法:以DEA-Malmquist指数替代SFA重新测算TFP,检验结论是否对生产前沿估计方法敏感。
- 安慰剂检验:将各省协议签订时间分别前移3年和后移3年作为虚假处理时点,检验处理效应是否仅在真实时点存在——若虚假时点的估计系数不显著,则加强了因果推断的可信度。
- Bacon-Decomp分解:将TWFE估计量分解为各组别-时间的2×2 DID加权平均,检验是否存在"已处理组作为对照组"的负权重污染(Goodman-Bacon, 2021)。
- 排除极端天气干扰:在回归中额外控制各省当年"连续降雨天数"和"高温天数"(中国气象局站点数据),检验跨区协作效应是否在控制天气冲击后仍然稳健。