试验数据处理及配方施肥简介（上）

1、“3414”田间试验数据处理 问题与讨论,河北省农林科学院农业资源环境研究所,报告内容,报告内容,“3414”试验简介 “3414”试验数据处理与分析 土壤养分丰缺指标的建立 实际操作举例,“3414”是指氮、磷、钾3个因素、4个水平、14个处理。 4个水平的含义： 0水平指不施肥 2水平指当地最佳施肥量 1水平=2水平0.5 3水平=2水平1.5（该水平为过量施肥水平）。 “3414”方案设计吸收了回归最优设计处理少、效率高的优点。,3414简介,“3414”施肥方案设计,测土配方施肥工作的目的是 通过多年、多点的3414试验，,1、建立/更新土壤肥力丰缺指标 2、确定推荐施肥量,土壤养分丰

2、缺指标法是国际上最通用的测土配方施肥方法！,为什么要建立/更新测土配方施肥技术指标体系？,1）由于土壤肥力、作物产量等因素发生了很大变化，二十世纪八十年代建立的测土配方施肥指标体系已不能适用于目前的生产条件，而全国很多地区还没有建立新的指标体系。,通过田间试验建立测土配方施肥技术指标体系,全国土壤养分丰缺指标协作组：,2）近十多年来，新的作物生产体系不断出现，蔬菜、果树、花卉等经济作物的种植面积和重要性越来越突出，但就全国而言，尚未建立这些作物测土施肥的技术体系。,为什么要建立/更新测土配方施肥技术指标体系？,通过田间试验建立测土配方施肥技术指标体系,河北省“3414”试验各因素施肥水平（kg

3、/mu）,报告内容,“3414”试验简介 “3414”试验数据处理与分析 土壤养分丰缺指标的建立 实际操作举例,试验处理（三元二次效应函数）,y = b0 + b1x1 + b2x2 + b3x3 + b4x12 + b5x22 + b6x32 +b7x1x2 + b8x1x3 + b9x2x3,系数均为正值,系数均为负值,最高施肥量 自变量xj (j = 1,2.3.) 的边际效应等于0,b1 + 2b4x1 +b7x2 + b8x3 = 0,b2 + 2b5x2 + b7x1 + b9x2 = 0,b3 + 2b6x32 +b8x1 + b9x2 = 0,求解上述方程组，获取最高施肥量 x

4、1,2,3,最佳施肥量 自变量xj (j = 1,2.3.) 的边际效应等于投入价格与产出价格比,b1 + 2b4x1 +b7x2 + b8x3 = Px1/Py,b2 + 2b5x2 + b7x1 + b9x2 = Px2/Py,b3 + 2b6x32 +b8x1 + b9x2 = Px3/Py,求解上述方程组，获取最高施肥量 x1,2,3,3413三元二次效应函数处理软件：,肥料效应函数的拟合和推荐施肥量计算,1 、利用网络工具进行运算,3414 三元二次效应函数处理软件： ,中国肥料信息网,试验数据的分析,第一步： 数据准备,第二步：进行网络运算,第三步：分析结果的判读,2、利用EXCE

5、L运算,以Microsoft Office2000为例,第一步： 数据准备,第二步：数据计算,第三步：打开工具下拉菜单数据分析回归，确定,第四步：在“y值输入区域”输入产量，“x值输入区域”输入 x1, x2, x3等 ，若将表头也选中，在标志框处打“”，并选择任意位置为输出区域，确定。,第五步：数据运算结果判读,b1 + 2b2x1 + b7x2 + b8x3 = Px1/Py b3 + 2b4x2 + b7x1 + b9x3 = Px2/Py b5 + 2b6x3 + b8x1 + b9x2 = Px3/Py,上述方程拟合成功，根据边际收益等于边际成本原则计算经济最佳施肥量,第六步：三元二

6、次方程,注意： EXCEL软件必须安装数据分析功能才能进行上述计算。,1、放入OFFICE安装盘；2、按下操作进行,注意：,第一步：打开office 2000安装盘，在Microsoft Office文件夹下有MSOWC.DLL文件，复制，粘贴到C:Program FilesMicrosoft OfficeOffice文件中。 第二步：点击开始菜单，点运行选项，在运行选项中，输入命令: RegSvr32 “C:Program FilesMicrosoft OfficeOfficeMSOWC.DLL”；其中 C:Program FilesMicrosoft Office 为 Office 的路径

7、，要根据实际情况进行修改，确定，会出现注册成功字样。 第三步： 如果在你的计算机中安装了Microsoft Office并且正确注册了 OWC (Office Web Component) 组件，仍然打不开土肥工具，请在网页工具中，IE 中的安全级别设置为“中”、“中低”或“低”；具体步骤：工具 Internet 选项 安全 调整 Internet安全级别,因为许多版本的office 2003中不含有MSOWC.DLL文件，所以在在打开网页后，表格不能输入数据，从而无法计算提供最佳施肥量。,3 、一元二次方程,第一步：输入处理,第二步：打开工具下拉菜单数据分析回归，确定 第三步：在“y值输入区

8、域”输入产量，“x值输入区域” 输入 x1, x2, x3等 ，若将表头也选中，在标志框处打“”，并选择任意位置为输出区域，确定。 第四步：数据运算结果判读 第五步：写出一元二次公式 第六步：计算最佳氮肥用量,N， y = ax2 +bx +c P， y = ax2 +bx +c K， y = ax2 +bx +c,分别计算氮磷钾的一元二次方程,最高施肥量 自变量x 的边际效应等于0,N，b + 2ax = 0,求解上述方程组，获取最高施肥量 x,最佳施肥量 自变量xj (j = 1,2.3.) 的边际效应等于投入价格与产出价格比,求解上述方程组，获取最高施肥量 x,P，b + 2ax = 0

9、,K，b + 2ax = 0,N，b + 2ax = Px/Py,P，b + 2ax = Px/Py,K，b + 2ax = Px/Py,4、 线性加平台模型拟合,线性加平台模型,第一步：数据准备、保存,第二步，打开SAS程序,第三步：修改参数,第三步：修改参数,（1）修改文件保存位置和文件名（infile），例中文件位置为d:book2.prn。 注意：SAS无法识别中文系统，故要求大家将待处理的数据均以book文件名保存在根目录下，便于SAS程序寻找待处理问题。 （2）修改第七行中作物产量范围和施肥量变化范围及其步长，A值代表作物产量范围和步长；K值代表施肥量变化范围和步长，B值无需变动。

10、例中产量范围是60009200 kg/ha，运算的补充为200 kg/ha，施肥量变化范围为0400 kg/ha，运算的步长为20 kg/ha。 注意：程序中的作物产量和施肥量范围一定要大于试验处理结果的变化范围。,运行完结果表示： PLATEAU=* 即为平台产量（最高产量） INTERCEPT=* 表示直线的起点（最低产量） SLOPE=* 表示直线的斜率 JOINT=* 表示直线与平台交会的X轴坐标（优化肥料用量） RSQ=* 为R2 F=* 为显著性检验的F值 PROBF=* 为显著性检验结果是否显著,第四步：结果运算；第五步 结果判读与记录,第六步：添加散点,第七步：添加趋势线,第八

11、步：修饰图形,不同产量增加趋势需要不同的模型、任何如何数据都不应随意修改或舍弃,能用三元二次模型拟合的，可以直接利用三元二次模型提供的公式计算施肥量。 不能用三元二次模型拟合的，要考虑用一元二 次、线性加平台等模型进行拟合，计算最佳施肥量。 不能用任何模型拟合的，要根据各自的情况进行判断，给出相应的推荐施肥量,推荐施肥量的确定,b1 + 2b2x1 + b7x2 + b8x3 = Px1/Py b3 + 2b4x2 + b7x1 + b9x3 = Px2/Py b5 + 2b6x3 + b8x1 + b9x2 = Px3/Py,对用三元二次模型拟合成功的，根据边际收益等于边际成本原则计算经济最

12、佳施肥量,既适合线性加平台模型又适合一元二次模型,当在前三个施肥水平作物产量随施肥量增加而增加，第四个施肥水平与第三个施肥水平相差不大时 ，采用任何模型计算最佳施氮量都是可以的。,仅适合线性加平台模型、仅适合一元二次方程,当第二个施肥水平已接近或达到最高产量时，应用一元二次方法拟合最佳产量施肥量明显高于应用线性加平台模型计算施肥量，故应该选用线性加平台模型。,当在前三个施肥水平产量随施肥量增加而增加，第四个施肥水平与第三个施肥水平相比，有明显降低时，应该选用一元二次方程拟合最佳产量施肥量,异常反应曲线1,作物产量没有随施肥量的增加而变化，或作物产量随施肥量的增加而出现减产趋势。,最大产量施肥量(kg/ha) 优化产量施肥量(kg/ha),194 -294,2811 1288,优化肥料用量：0,作物产量随施肥量的增加而呈跳跃式变化。,异常反应曲线2,检查是否是产量计算或数据输入错误； 若没有错误，对于高肥力田块可定义优化施肥量为0，对于无法判断肥 力水平的田块应舍弃,优化肥料