3414肥料试验结果统计分析方法研究.docxVIP

3414肥料试验结果统计分析方法研究 自2005年以来，国家在全国范围内开展了广泛的土壤测试和施肥试验，将实验与生产相结合，指导科学施肥，促进粮食高产和农民增收，降低肥料成本，提高经济效益和土壤肥力，促进环境效应。“3414”肥料试验设计吸收了回归最优设计处理少、效率高的优点, 是测土配方施肥工作主要推荐的肥料施肥方案。通过“3414”试验可以建立包括肥料效应函数法、养分平衡法、土壤养分丰缺指标法等在内的测土配方体系。 肥料效应函数法是测土配方施肥技术中最常用的方法。通过统计分析, 不仅可以作为一个完整的三因素试验用于建立三元二次肥料效应函数, 而且还可以建立三个两因素的二元二次肥料效应函数和三个单因素的一元二次肥料效应函数。即使某一个或几个处理出问题, 仍可以获得一些用于肥料决策的价值信息, 提高了试验效率。 近期, 作者查阅了大量关于“3414”肥料试验及相关内容的文章。由于不同的工作者对试验分析的侧重点不同, 肥料试验结果的分析方法存在很大的差异, 甚至个别文章分析方法存在误区, 结果是错误的。肥料效应函数既关系到肥料效应, 又属于数学函数范畴, 问题较为复杂。因此, 把肥料效应函数分析存在的误区及本人的理解归纳总结如下, 与大家共同探讨。 1 对f检验的检验 利用SPSS、DPS、SAS、Excel等数据统计软件, 对试验结果进行回归分析, 可以获得三元二次肥料方程: ?Y=b0+b1X1+b2X21+b3X2+b4X22+b5X3+b6X23+b7X1X2+b8X1X3+b9X2X3(1)Y?=b0+b1X1+b2X21+b3X2+b4X22+b5X3+b6X23+b7X1X2+b8X1X3+b9X2X3(1) 同时获得对方程的F检验值 (相关系数检验与F值检验存在一致性) 。 根据肥料报酬递减律和二次数学函数曲线的特点, 一般认为方程F检验达到显著水平 (P0.05) 以上且符合肥料报酬递减律 (b1、b3和b5均为正值, b2、b4和b6均为负值) 时, 认为该三元二次方程拟合成功。这是最为理想的结果, 一般称为典型肥料效应函数。 实际上, 田间试验影响因素较多, 模拟的肥料效应方程会出现多种情况。可能不符合报酬递减律, 即方程中二次项系数可能为正值, 一次项系数可能为负值。但是只要方程F检验达显著水平 (P0.05) 以上, 也可认为拟合成功, 一般称这类函数成为非典型肥料效应函数。 如果拟合三元二次方程不成功, 方程F检验不显著, 也不要放弃对数据的分析和有价值信息的挖掘, 可以拟合单因素效应、两因素效应方程。 2 计算的最大功率、最佳功率和相应的含肥量 拟合三元二次方程成功后, 需要进行最高产量、最佳产量及对应施肥量的计算。 2.1 最佳产量值的确定及验证 根据肥料报酬递减律, 当边际产量 (dy/dx) 等于零时, 作物产量达到最高点, 此时的施肥量为最大施肥量或称最高产量施肥量。对一元二次函数采用微分求导的方法, 求解到的最大施肥量是正确的。因此, 许多工作者依此类推, 对于三元二次函数也采用了微分求导的方法。即对方程 (1) 求各因素的偏导数并使其等于零, 获得三元一次方程组: {?y?x1=b1+2b2X1+b7X2+b8X3=0?y?x2=b3+2b4X2+b7X1+b9X3=0?y?x3=b5+2b6X3+b8X1+b9X2=0????????????y?x1=b1+2b2X1+b7X2+b8X3=0?y?x2=b3+2b4X2+b7X1+b9X3=0?y?x3=b5+2b6X3+b8X1+b9X2=0 求解, 获得最高产量施肥量X1、X2、X3值, 代入方程 (1) 即获得最高产量值。 当边际效益等于边际成本 (dyPy=dxPx, 即dy/dx=Px/Py时) , 肥料的边际利润为零, 获得最佳产量及对应施肥量。也采用类似一元二次函数微分求导的方法, 对方程 (1) 求各因素的偏导数, 并使其等于Px/Py, 获得三元一次方程组: {?y?x1=b1+2b2X1+b7X2+b8X3=Ρx1Ρy?y?x2=b3+2b4X2+b7X1+b9X3=Ρx2Ρy?y?x3=b5+2b6X3+b8X1+b9X2=Ρx3Ρy????????????y?x1=b1+2b2X1+b7X2+b8X3=Px1Py?y?x2=b3+2b4X2+b7X1+b9X3=Px2Py?y?x3=b5+2b6X3+b8X1+b9X2=Px3Py 求解, 获得最佳产量施肥量X1、X2、X3值, 代入方程 (1) 即获得最佳产量值。Px1、Px2、Px3、Py分别表示肥料X1、X2、X3和农产品y的单价, 一般肥料和农产品的价格取前三年价格的均值。 2.2 最高产量函数 模拟获得的三元二次肥料函数即使符合报酬递减律, 即二次项为