试验设计及其统计建模幻灯片.ppt

《试验设计及其统计建模幻灯片.ppt》由会员分享，可在线阅读，更多相关《试验设计及其统计建模幻灯片.ppt（30页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、试验设计及其统计建模第1页，共30页，编辑于2022年，星期二3414试验设计第2页，共30页，编辑于2022年，星期二3414试验简介“3414”试验方案类似二次回归D-最优设计。该实验设计由氮、磷、钾3个因素、4个水平、14个处理组成。一般认为，“3414”试验设计方案设计吸收了回归最优设计处理少、效率高的优点。4个水平的含义指：0水平为不施肥，2水平为当地最佳施肥量的近似值，1水平=2水平0.5，3水平=2水平1.5(该水平为过量施肥水平)。第3页，共30页，编辑于2022年，星期二试验设计基本原则均匀分散；整齐可比（如正交设计示例）第4页，共30页，编辑于2022年，星期二和3414试

2、验近似的几个试验设计4因素D饱和设计（D饱和，取前3列，15个处理）Box-Behnken设计（BB，不完全的3水平试验设计，15个处理）。5因素4水平的正交试验设计(正交，取前3列，16次试验)。3因素(1/2实施)复合中心设计(CCD，中心点4个，14次试验)。第5页，共30页，编辑于2022年，星期二几个试验设计的散点图比较第6页，共30页，编辑于2022年，星期二试验设计两个优良指标比较试验E-最优(条件数)行列式值34142.04080.0244D饱和1.28000.0237B-B1.01060.0241正交1.00000.0201CCD1.06590.0292第7页，共30页，编辑

3、于2022年，星期二3因子(1/2)CCD设计示意图第8页，共30页，编辑于2022年，星期二几个3因素CCD设计的水平值表正交正交表表-101+试验试验点点中中心心点点条件数条件数行列式值行列式值0(正交正交)0.36 1.50 2.6431441.06591.06590.02920.02920(旋转旋转)0.44 1.50 2.5631441.18750.0292全全0(正交正交)0.27 1.50 2.7331511.04350.0241全全0(旋转旋转)0.61 1.50 2.3931511.14290.1162全全0(正交正交)0.33 1.50 2.6631621.05710.17

4、36全全0(旋转旋转)0.61 1.50 2.3931621.18750.0836第9页，共30页，编辑于2022年，星期二从统计学意义上讲，建议的试验设计：3因子CCD设计（中心点为2的正交回归组合设计）。No.No.c1 c1 c2 c2 c3 c31 11 11 11 12 21 11 1-1-13 31 1-1-11 14 41 1-1-1-1-15 5-1-11 11 16 6-1-11 1-1-17 7-1-1-1-11 18 8-1-1-1-1-1-19 9-1.2872-1.2872 0 0 0 0 10101.2872 1.2872 0 0 0 0 11110 0-1.287

5、2-1.2872 0 0 12120 0 1.2872 1.2872 0 0 13130 0 0 0-1.2872-1.2872 14140 0 0 0 1.2872 1.2872 15150 0 0 0 0 0 16160 0 0 0 0 0 第10页，共30页，编辑于2022年，星期二04水平试验设计表x1x2x33.55383.55383.55383.5538 3.5538 0.44623.55380.44623.55383.55380.44620.44620.44623.55383.55380.44623.55380.44620.4462 0.4462 3.55380.44620.44

6、620.4462022422202242220224222222第11页，共30页，编辑于2022年，星期二统计建模（反应面分析）：拟合二次多项式回归模型第12页，共30页，编辑于2022年，星期二3414试验，一般建模过程的数据格式第13页，共30页，编辑于2022年，星期二数据格式（有重复的情形）第14页，共30页，编辑于2022年，星期二3414试验，一般建模过程的用户界面第15页，共30页，编辑于2022年，星期二三元二次多项式回归方程F值=16.4746Df=(9,4)p值=0.0080剩余标准差=0.9058，复相关R=0.9868，决定系数R2=0.9737回归方程:y=115.

7、2177+8.4457N+11.6216P-4.3366K-0.2140N2-0.4208P2-0.5060K2-0.2395NP+0.4582NK+0.2991PK优化结果 产量最高 效益最大 目标y 257.3996196.5864x1 24.505212.1522x2 10.34188.6085x3 9.86671.7851第16页，共30页，编辑于2022年，星期二选取部分重要肥料因子建模当三元二次回归方程不能建立时，根据3414试验设计思想，可用14个处理还可分别进行氮、磷、钾中任意二元二次或一元二次肥料效应函数的拟合。第17页，共30页，编辑于2022年，星期二模型拟合三元二次方程

8、一元二次方程PK成功()9.037.332.8不成功()91.062.767.2山东省冬小麦66个3414试验方程拟合成功率比较引自孙义祥等应用3414试验建立冬小麦测土配方施肥指标体系，植物营养和肥料学报，2009，15（1）：197203第18页，共30页，编辑于2022年，星期二当前建模方法的不足没有将区组（小区控制）因素所产生的效应从试验误差中去掉，可能是造成拟合三元二次回归方程拟合率的一个重要原因。进行氮、磷、钾中任意二元二次或一元二次肥料效应函数的拟合，没有考虑肥料因子间互作，因为3414试验各个因子间并不是相互正交。这种方法看似合理，实际上没有充分利益实验结果的所有信息。因此，用

9、统计学的方法挑选“重要”因子建模更合理。第19页，共30页，编辑于2022年，星期二3414试验数据统计分析方法改进区组设计的区组效应分解，可提高模型拟合精度。逐步回归分析，选取部分二次项、交互项作用显著的变量（肥料因子），建立回归模型。对试验结果进行典范分析，确定优化方向，在下一年试验中，以现有试验为基础进行进一步试验。第20页，共30页，编辑于2022年，星期二一个区组设计例子No.NPK10003.43.33.820224.04.14.431224.74.55.442024.64.65.352124.84.85.362224.84.44.772324.55.05.082203.54.64

10、.392214.14.84.9102233.94.63.7113225.35.46.2121123.94.04.3131214.64.54.8142114.95.04.9第21页，共30页，编辑于2022年，星期二3414试验结果统计建模（目前方法）回归方程y=3.5019+0.3090N+0.4675P+0.4029K+0.1153N2-0.0093P2-0.3696K2-0.3649NP+0.2353NK+0.1408PK复相关系数R=0.8329，决定系数R2=0.6937F值=8.0535，Df=(9,32)，p值0.0001第22页，共30页，编辑于2022年，星期二试验区组效应的处

11、理3414试验区组设计，以往试验优化分析工具，都没有提供可处理区组设计功能。当3414试验有重复、且是区组设计时，须充分利用区组的小区控制效应，进行建模分析，一提高数据分析的精度。统计建模时，我们可增加了处理含有若干区组的二次正交回归组合(中心复合)设计试验数据分析建模功能。DPS中，该功能模块在“试验统计”“试验优化分析”“区组设计二次多项式回归”里面。第23页，共30页，编辑于2022年，星期二考虑区组效应的3414试验建模Y=3.31+0.1649b1+0.4159b2+0.3089x1+0.4674x2+0.4029x3+0.1152x1*x1-0.00934*x2*x2-0.3696

12、x3*x3-0.3649x1*x2+0.2353x1*x3+0.1408*x2*x3方差分析表变异来源 平方和df均方F值p值回归11.01 11 1.009.72 0.0001模型9.7891.0910.55 0.0001区组1.2320.615.960.0066残差3.09300.10总变异14.10 41相关系数R=0.8837决定系数R2=0.7809第24页，共30页，编辑于2022年，星期二点间、年间效应类似处理将一年多点数据或一点多年数据中的试验点或年份作为分类变量，类似区组效应分解的方式进行处理，可改进建模的精度，提高建模的效果。第25页，共30页，编辑于2022年，星期二挑选

13、重要因子的逐步回归第26页，共30页，编辑于2022年，星期二第27页，共30页，编辑于2022年，星期二分析结果Y=3.295+0.1649*b1+0.4159*b2+0.0513*x1+0.5864*x2+0.5376*x3+0.1276*x1*x1-0.3650*x3*x3-0.3057*x1*x2+0.2929*x1*x3方差分析表变异来源平方和df均方 F值 p值回归10.9491.22 12.300.0001模型9.7171.3914.040.0001区组1.23 20.61 6.220.0066残差3.16 320.10总变异14.10 41第28页，共30页，编辑于2022年，星期二典范分析的应用当典范分析得到的稳定点和数值优化分析得到的稳定点一致时，最优方案方可应用。如果当前的优化点不是最优点，典范分析可提示其优化方向。根据优化方向，可作为制定来年试验设计计划的依据。第29页，共30页，编辑于2022年，星期二试验的序贯设计第30页，共30页，编辑于2022年，星期二