回归分析答案.doc

试验设计作业1、下表为小麦栽培试验的产量结果(kg)，随机区组设计，小区计产面积为12m2，试作分析。在表示最后结果时需化为每亩产量(kg)。假定该试验为一完全随机设计，试分析后将其试验误差与随机区组时的误差作一比较，看看划分区组的效果如何？处 理区 组A6.26.66.96.1B5.86.76.06.3C7.26.66.87.0D5.65.85.46.0E6.97.27.07.4F7.57.87.37.6完全随机设计的程序如下：data li_1; do i=1 to 6; do j=1 to 4; input x; output; end; end; cards;6.2 6.6 6.9 6.15.8 6.7 6 6.37.2 6.6 6.8 75.6 5.8 5.4 66.9 7.2 7 7.47.5 7.8 7.3 7.6;proc anova;class i;model x=i;means i;run;SAS输出结果如下： Sum of Source DF Squares Mean Square F Value Pr F Model 5 8.97208333 1.79441667 20.87 F i 5 8.97208333 1.79441667 20.87 F Model 8 9.24333333 1.15541667 13.58 F i 5 8.97208333 1.79441667 21.09 F Model 12 255.2800000 21.2733333 7.40 0.0008 Error 12 34.4800000 2.8733333 Corrected Total 24 289.7600000 R-Square Coeff Var Root MSE x Mean 0.881005 6.958501 1.695091 24.36000 Source DF Anova SS Mean Square F Value Pr F i 4 8.5600000 2.1400000 0.74 0.5798 j 4 44.5600000 11.1400000 3.88 0.0302 fz 4 202.1600000 50.5400000 17.59 F Model 8 90.6942857 11.3367857 7.61 0.0011 Error 12 17.8771429 1.4897619 Corrected Total 20 108.5714286 R-Square Coeff Var Root MSE x Mean 0.835342 7.078630 1.220558 17.24286 Source DF Anova SS Mean Square F Value Pr F i 6 88.17142857 14.69523810 9.86 0.0005 j 2 2.52285714 1.26142857 0.85 0.4529油菜随机区组设计程序如下：data li_3; do i=1 to 6; do j=1 to 4; input x; output; end; end;cards;3.9 3.9 3.6 4.55.8 6.3 5.9 6.84.4 4.4 5.6 4.55.5 5.2 5.4 6.76.8 6.9 7.4 67.3 7.2 7.5 7;proc anova; class i j; model x=i j;run;结果如下： Sum of Source DF Squares Mean Square F Value Pr F Model 8 31.27500000 3.90937500 14.18 F i 5 30.81708333 6.16341667 22.36 F Model 12 401.9200000 33.4933333 4.41 0.0078 Error 12 91.1200000 7.5933333 Corrected Total 24 493.0400000 R-Square Coeff Var Root MSE x Mean 0.815187 17.52926 2.755600 15.72000 Source DF Anova SS Mean Square F Value Pr F i 4 23.0400000 5.7600000 0.76 0.5717 j 4 17.0400000 4.2600000 0.56 0.6954 fz 4 361.8400000 90.4600000 11.91 0.0004结果分析：品种间差异对结果的影响是显著的。5、 调查某队元麦及元麦和蚕豆混种、间种田块的产量(混、间种者为麦、豆产量合计)，得结果(kg/0.1亩)于表：(1)设以元麦单种为对照，试以LSD法作多重比较；(2)设预定要作的比较是单种对混、间种，混种对间种，2麦1豆间种对3麦2豆间种，试作单一自由度的独立比较。元麦单种麦豆混种2麦1豆间种3麦2豆间种2024303024232833222834311821323621243135先作方差分析：data li_1; do i=1 to 4; do j=1 to 5; input x; output; end; end; cards;20 24 22 18 2124 23 28 21 2430 28 34 32 3130 33 31 36 35;proc anova;class i;model x=i;means i;run;结果如下： Sum of Source DF Squares Mean Square F Value Pr F Model 3 483.7500000 161.2500000 28.04 F i 3 483.7500000 161.2500000 28.04 F Model 11 9.33333333 0.84848485 0.31 0.9778 Error 24 66.66666667 2.77777778 Corrected Total 35 76.00000000 R-Square Coeff Var Root MSE y Mean 0.122807 11.90476 1.666667 14.00000 Source DF Anova SS Mean Square F Value Pr F a 3 2.88888889 0.96296296 0.35 0.7918 b 2 0.50000000 0.25000000 0.09 0.9142 a*b 6 5.94444444 0.99074074 0.36 0.8989结果分析：品种、播期，品种播期对试验的影响均不显著。7. 有一小麦裂区试验，主区因素A，分A1(深耕)、A2(浅)两水平，副区因素B，分B1(多肥)、B2(少肥)两水平，重复3次，小区计产面积15平方米，其田间排列和产量(假设数字)如下图，试作方差分析。因素A, 因素B及其交互效应是否显著?A1A2A2A1A2A1B19B17B23B111B21B24B26B22B15B24B16B112区组区组区组解：程序如下：title 裂区试验的统计分析;data lq;do a=1 to 2; do b=1 to 2; do g=1 to 3; input yield ; output; end; end;end;cards;9 11 126 4 47 5 62 3 1;proc anova;class g a b;model yield= g a b ;means a b/duncan;run;proc anova;class g a b;model yield= g a b a*b g*a g*a*b;test H=g a E=g*a;test H=b a*b E=g*a*b;run;试验结果：裂区试验的统计分析208:07 Wednesday, July 3, 2002 Analysis of Variance ProcedureDependent Variable: YIELDSource DF Sum of Squares F Value Pr FModel 4 115.50000000 14.27 0.0018Error 7 14.16666667Corrected Total 11 129.66666667 R-Square C.V. YIELD Mean 0.890746 24.38754 5.83333333Source DF Anova SS F Value Pr FG 2 0.16666667 0.04 0.9599A 1 40.33333333 19.93 0.0029B 1 75.00000000 37.06 0.0005裂区试验的统计分析608:07 Wednesday, July 3, 2002 Analysis of Variance ProcedureDependent Variable: YIELDSource DF Sum of Squares F Value Pr FModel 11 129.66666667 . .Error 0 .Corrected Total 11 129.66666667 R-Square C.V. YIELD Mean 1.000000 0 5.83333333Source DF Anova SS F Value Pr FG 2 0.16666667 . .A 1 40.33333333 . .B 1 75.00000000 . .A*B 1 3.00000000 . .G*A 2 1.16666667 . .G*A*B 4 10.00000000 . .Tests of Hypotheses using the Anova MS for G*A as an error termSource DF Anova SS F Value Pr FG 2 0.16666667 0.14 0.8750A 1 40.33333333 69.14 0.0142Tests of Hypotheses using the Anova MS for G*A*B as an error termSource DF Anova SS F Value Pr FB 1 75.00000000 30.00 0.0054A*B 1 3.00000000 1.20 0.3349结果分析：因素A、因素B的显著性水平为0.0029、0.0005，是显著的，而AXB的F值是1.2，顾AXB是不显著的。B1与B2有显著的差别，A1与A2有显著的差别，最好的是A1B1的组合。8. 设若上题小麦耕深与施肥量试验为条区设计，田间排列和产量将相应如下图，试作分析，并与裂区设计结果相比较)。A1A2A2A1A2A1B197B1511B214B262B234B1612解：程序如下：title 条区试验的统计分析;data tq;do g=1 to 3; do a=1 to 2; do b=1 to 2; input yield ; output; end; end;end;cards;9 6 7 211 4 5 312 4 6 1;proc anova;class g a b;model yield= g a b ;means a b/duncan;run;proc anova;class g a b;model yield= g a b a*b g*a b*g g*a*b;test H=g a E=g*a;test H=b a*b E=g*a*b;run; 试验结果：条区试验的统计分析208:42 Wednesday, July 3, 2002 Analysis of Variance ProcedureDependent Variable: YIELDSource DF Sum of Squares F Value Pr FModel 4 115.50000000 14.27 0.0018Error 7 14.16666667Corrected Total 11 129.66666667 R-Square C.V. YIELD Mean 0.890746 24.38754 5.83333333Source DF Anova SS F Value Pr FG 2 0.16666667 0.04 0.9599A 1 40.33333333 19.93 0.0029B 1 75.00000000 37.06 0.0005条区试验的统计分析608:42 Wednesday, July 3, 2002 Analysis of Variance ProcedureDependent Variable: YIELDSource DF Sum of Squares F Value Pr FModel 11 129.66666667 . .Error 0 .Corrected Total 11 129.66666667 R-Square C.V. YIELD Mean 1.000000 0 5.83333333Source DF Anova SS F Value Pr FG 2 0.16666667 . .A 1 40.33333333 . .B 1 75.00000000 . .A*B 1 3.00000000 . .G*A 2 1.16666667 . .G*B 2 3.50000000 . .G*A*B 2 6.50000000 . .Tests of Hypotheses using the Anova MS for G*A as an error termSource DF Anova SS F Value Pr FG 2 0.16666667 0.14 0.8750A 1 40.33333333 69.14 0.0142Tests of Hypotheses using the Anova MS for G*A*B as an error termSource DF Anova SS F Value Pr FB 1 75.00000000 23.08 0.0407A*B 1 3.00000000 0.92 0.4380结果分析：因素A、因素B的显著性水平为0.8750，0.0407，A*B的显著性水平为0.4380，故均属于不显著。9. 在药物处理大豆种子试验中，使用了大中小三种类型种子，分别用五种浓度、两种处理时间进行试验处理，播种后45天对每种各取两个样本，每个样本取10株测定其干物重，求其平均数，结果如下表。试进行方差分析。处理时间A种子类型C浓度BB1(010-6)B2(1010-6)B3(2010-6)B4(3010-6)B5(4010-6)A1(12小时)C1(小粒)7.012.822.021.324.46.511.421.820.323.2C2(中粒)13.513.220.419.024.613.814.221.419.623.8C3(大粒)10.712.422.621.324.510.313.221.822.424.2A2(24小时)C1(小粒)3.610.74.712.413.61.58.83.410.513.7C2(中粒)4.79.82.712.414.04.910.54.213.214.2C3(大粒)8.79.63.413.014.83.59.74.212.712.6解：程序如下：title 三因素随机区组试验的统计分析;data dsy;do a=1 to 2; do b=1 to 5; do c=1 to 3; do rep=1 to 2; input y ; output; end; end; end;end;cards;7.0 6.513.5 13.810.7 10.312.8 11.413.2 14.212.4 13.222.0 21.820.4 21.422.6 21.821.3 20.319.0 19.621.3 22.424.4 23.224.6 23.824.5 24.23.6 1.54.7 4.98.7 3.510.7 8.89.8 10.59.6 9.74.7 3.42.7 4.23.4 4.212.4 10.512.4 13.213.0 12.713.6 13.714.0 14.214.8 12.6;proc anova;cla