-
化肥过量施用导致农业面源污染与环境恶化是中国生态治理面临的严重问题之一[1-3]。据统计,2017年中国农用化肥施用总量为5 859.4万t,施用强度为352 kg·hm−2,超过国际警戒强度。近年来,中国政府先后出台了一系列政策文件,推动化肥减量化,测土配方施肥技术是促进化肥减量化的重要途经[4-5]。家庭农场是中国未来现代农业发展进程中最为合意和最具生命力的经营主体[6],家庭农场的生产决策行为直接关乎化肥减量化这一目标能否实现。中国目前测土配方施肥技术采用率不足1/3[7]。因此,以中国家庭农场等新型经营主体为对象,探究影响家庭农场采纳环境友好型技术的因素及影响强度,对促进中国农业绿色转型发展具有重要的现实意义。已有研究表明:创新意识、经营规模、外出务工等是影响家庭农场采纳测土配方施肥技术的重要因素[8-10],但从契约农业参与的视角来理论和实证分析其采纳测土配方技术的内在机制与效应还未见报道;同时,已有研究大多基于区域调查数据展开,鲜有基于全国大样本调查数据的实证检验,其结论的普适性有待进一步验证。鉴于此,本研究以2018年全国家庭农场监测数据为基础,分析作用机制,构建Probit-倾向得分匹配(PSM)模型,估计契约农业对家庭农场采纳测土配方施肥技术的实际影响,以期为政府制定相关政策提供依据。
HTML
-
家庭农场测土配方施肥技术采纳与否,所涉及到的被解释变量为二分类变量,即采纳测土配方施肥技术为1,反之为0。因此可采用Probit模型进行实证分析。同时,家庭农场选择是否参与契约农业还要结合自身需求和资源禀赋,即自选择会导致估计偏误。因此本研究利用基础模型分析,采用PSM构建“反事实推断模型”。
Probit模型如下:
式(1)中:Yi表示家庭农场i是否采纳测土配方施肥技术,采纳取值为1,未采纳取值为0;X1表示是否参与契约农业,参与取值为1,不参与取值为0;Xi表示影响家庭农场测土配方施肥技术采纳行为的特征变量,β1和βi分别为其估计系数;β0表示常数项,ε表示随机扰动项。
PSM匹配模型如下:
式(2)中:P(Zi)为倾向匹配得分值,表示家庭农场参与契约农业的倾向匹配得分值;Di为二值虚拟变量,即家庭农场i参与契约农业时取值为1,不参与取值为0;Zi为一系列的匹配变量,包括农场特征、农场主特征、经营特征等;β表示相应匹配变量的确定性系数。
基于匹配结果测算家庭农场参与契约农业对测土配方施肥技术采纳行为的平均处理效应(average treatment effect of the treated, EATT),计算公式如下:
式(3)中:y1i为参与契约农业的家庭农场采纳测土配方施肥技术的概率,y0i为匹配后得到的假如处理组未参与契约农业时家庭农场采纳测土配方施肥技术的概率。X1=1表示家庭农场参与契约农业。
-
基于上述机制分析,给出模型变量及测度。①被解释变量。被解释变量为家庭农场测土配方施肥技术采纳行为,是指家庭农场在全年的生产经营过程中是否采纳测土配方施肥技术。②关键解释变量。关键解释变量为家庭农场是否参与契约农业,是指家庭农场是否和合作社或农业企业签订了农产品销售合同。③其他控制变量。参考文献[9-10]和[15-16],选取其他可能影响家庭农场采纳测土配方施肥技术的变量,如农场主个人特征(性别、年龄、受教育程度、是否参加农业技术培训、从事农业规模经营年限、是否为外村户籍和从业经历)、农场特征(是否在工商部门登记注册、是否为农业部门认定的示范农场、农场土地经营总面积、农场经营的全部土地共有几块和家庭成员中在农场工作的人数)、经营特征(是否有日常收支记录、农场未来扩张意愿、是否有生产经营性借款、是否购买过农业保险)等。变量说明及描述性统计结果见表1。
类别 变量名 具体含义与赋值 均值 标准差 最小值 最大值 因变量 是否测土配方施肥技术 否为0;是为1 0.61 0.49 0 1 自变量 是否参与契约农业 否为0;是为1 0.32 0.46 0 1 农场主特征 性别 女为0;男为1 0.89 0.32 0 1 年龄 岁 46.14 8.52 18 73 受教育程度 未上过学为1;小学为2;初中为3;高中为4;中专为5;
职高为6;大专为7;本科为8;研究生及以上为93.89 1.38 1 9 是否接受过土肥培育技术培训 否为0;是为1 0.54 0.50 0 1 从事农业规模经营年限 a 7.75 4.81 0 36 是否为外村户籍 否为0;是为1 0.82 0.39 0 1 主要从业经历 非普通农民为0;普通农民为1 0.37 0.48 0 1 农场特征 是否在工商部门登记注册 否为0;是为1 0.80 0.40 0 1 是否为农业部门认定的示范农场 否为0;是为1 0.55 0.50 0 1 土地经营总面积 hm2 23.88 32.63 2.00 475.30 全部经营土地块数 块 13.30 28.32 1 385 家庭成员中在农场工作的人数 人 2.88 1.17 1 15 日常收支记录 否为0;是为1 0.76 0.43 0 1 土地未来扩张意愿 否为0;是为1 0.41 0.49 0 1 是否有生产经营性借款 否为0;是为1 0.60 0.49 0 1 是否购买农业保险 否为0;是为1 0.62 0.48 0 1 Table 1. variable assignment and description
2.1. 模型设定
2.2. 变量选取
-
由表2可知:模型整体拟合程度较好,参数估计符合预期。从关键解释变量来看,家庭农场参与契约农业对其采纳测土配方施肥技术的影响系数为正,且P<0.01。表明相比未参与契约农业的家庭农场,参与契约农业的家庭农场采纳测土配方施肥技术的可能性更大。
类别 变量 测土配方施肥技术采纳行为 系数 标准误 是否参与契约农业 0.621*** 0.083 农场主特征 性别 −0.013 0.108 年龄 0.006 0.004 受教育程度 0.005 0.027 是否参加土肥培育技术培训 0.516*** 0.069 从事农业规模经营年限 0.018** 0.007 是否为外村户籍 −0.144 0.091 主要从业经历 −0.425*** 0.076 农场特征 是否在工商部门登记注册 −0.091 0.093 是否为农业部门认定的示范农场 −0.123 0.075 土地经营总面积 0.000 0.000 全部经营土地块数 0.000 0.001 家庭成员中在农场工作的人数 0.007 0.029 是否有日常收支记录 0.613*** 0.087 土地未来扩张意愿 0.205*** 0.072 是否有生产经营性借款 −0.031 0.073 是否购买农业保险 0.196*** 0.074 常数项 −0.880*** 0.302 样本量 1 706 拟合优度 0.184 χ2 420.10 说明:
**、***分别代表在5%、1%的统计水平下显著Table 2. Probit model estimation results
就农场主特征来看,农场主参加土肥培育技术培训对家庭农场采纳测土配方施肥技术具有正向影响,影响系数为0.516,且P<0.01。说明参与土肥培育技术培训的家庭农场更愿意采纳测土配方施肥技术。农场主从事农业规模年限对家庭农场采纳测土配方施肥技术具有正向影响,影响系数为0.018,且P<0.05。说明农场主农业规模年限越长,越倾向于采纳测土配方施肥技术。另外,农场主的从业经历对家庭农场采纳测土配方施肥技术也有显著影响。
就农场特征来看,农场具有完整日常收支记录对家庭农场采纳测土配方施肥技术具有正向影响,影响系数为0.613,且P<0.01。说明具有完整日常收支记录的家庭农场更愿意采纳测土配方施肥技术;可能的原因是,家庭农场有完整日常收支记录会极大地提高其成为示范典型的可能性[17],从而提升农场主采纳测土配方施肥技术的意愿和动机。农场土地未来扩张意愿对家庭农场采纳测土配方施肥技术具有正向影响,影响系数为0.205,且P<0.01。说明具有土地未来扩张意愿的家庭农场采纳测土配方施肥技术的可能性更大;可能的原因是家庭农场具有未来继续扩张土地的意愿一定程度能反映出家庭农场土地经营的稳定性,较高的地权稳定性是促进家庭农场可持续发展的重要指标[18],对家庭农场采纳测土配方施肥技术意愿有利。
-
家庭农场是否参与契约是自主选择的结果,Probit模型回归结果可能存在选择性偏误。因此需要采用倾向得分匹配(PSM)来处理可能存在的自选择问题。PSM模型需要满足2个基本假定,即共同支撑假设和匹配变量的平衡性假定。
共同支撑假设要求对照组与实验组的倾向得分值重叠区间要足够大,否则将导致样本缺失。结果显示:1 706个家庭农场样本数据中,共有1 647个样本数据满足共同支撑假设,可进行匹配。由图1和图2可以看出:样本匹配前,对照组和实验组概率密度曲线重叠较少,吻合度不高,说明存在显著性差异;样本匹配后,概率密度曲线重叠增多,吻合度较高,说明2组样本无显著差异,各个维度特征基本趋于相似,满足共同支撑假设。
PSM的平衡性检验主要是考察协变量在控制组与实验组之间是否存在显著性差异。最近邻匹配下的平衡性检验假设结果显示:匹配后各个协变量的标准偏差均大幅度减少,除农场主年龄(8.5%)和农场家庭劳动力(22.9%)的标准偏差下降幅较小,其他协变量标准偏差下降幅度均超过60%,并且所有协变量的标准偏差均小于10%。2组样本均值非常接近,无显著性差异,即通过了平衡性检验。
-
本研究PSM模型分别采用最近邻匹配(K=1)、最近邻匹配(K=4)、半径匹配(R=0.01)和核匹配4种方法对样本进行匹配。由表3可知:4种匹配结果均通过了1%的显著性检验,且处理效应系数差异不大。一方面表明PSM模型估计结果较为稳健,另一方面说明家庭农场参与契约农业对其采纳测土配方施肥技术有显著的促进作用。
匹配方法 处理组/对照组(样本数) EATT 匹配方法 处理组/对照组(样本数) EATT 最近邻匹配(K=1) 538/1 168 0.195***(0.040) 半径匹配(R=0.01) 529/1 091 0.196***(0.029) 最近邻匹配(K=4) 533/1 114 0.183***(0.034) 核匹配 533/1 114 0.190***(0.026) 说明:括号内为标准误;***表示在1%的统计水平下差异显著 Table 3. Average processing effect of contract farming participation