回归分析答案.docx
《回归分析答案.docx》由会员分享,可在线阅读,更多相关《回归分析答案.docx(53页珍藏版)》请在冰豆网上搜索。
回归分析答案
试验设计作业
1、下表为小麦栽培试验的产量结果(kg),随机区组设计,小区计产面积为12m2,试作分析。
在表示最后结果时需化为每亩产量(kg)。
假定该试验为一完全随机设计,试分析后将其试验误差与随机区组时的误差作一比较,看看划分区组的效果如何?
处理
区组
Ⅰ
Ⅱ
Ⅲ
Ⅳ
A
6.2
6.6
6.9
6.1
B
5.8
6.7
6.0
6.3
C
7.2
6.6
6.8
7.0
D
5.6
5.8
5.4
6.0
E
6.9
7.2
7.0
7.4
F
7.5
7.8
7.3
7.6
完全随机设计的程序如下:
datali_1;
doi=1to6;
doj=1to4;
inputx@@;
output;
end;
end;
cards;
6.26.66.96.1
5.86.766.3
7.26.66.87
5.65.85.46
6.97.277.4
7.57.87.37.6
;
procanova;
classi;
modelx=i;
meansi;
run;
SAS输出结果如下:
Sumof
SourceDFSquaresMeanSquareFValuePr>F
Model58.972083331.7944166720.87<.0001
Error181.547500000.08597222
CorrectedTotal2310.51958333
R-SquareCoeffVarRootMSExMean
0.8528934.4064150.2932106.654167
SourceDFAnovaSSMeanSquareFValuePr>F
i58.972083331.7944166720.87<.0001
随机区组设计的程序如下:
datali_3;
doi=1to6;
doj=1to4;
inputx@@;
output;
end;
end;
cards;
6.26.66.96.1
5.86.766.3
7.26.66.87
5.65.85.46
6.97.277.4
7.57.87.37.6
;
procanova;
classij;
modelx=ij;
run;
结果如下:
Sumof
SourceDFSquaresMeanSquareFValuePr>F
Model89.243333331.1554166713.58<.0001
Error151.276250000.08508333
CorrectedTotal2310.51958333
R-SquareCoeffVarRootMSExMean
0.8786794.3835760.2916906.654167
SourceDFAnovaSSMeanSquareFValuePr>F
i58.972083331.7944166721.09<.0001
j30.271250000.090416671.060.3943
结果分析:
随机区组设计的误差要小一些。
2、下表为水稻品种比较试验的产量结果(kg),5×5拉丁方设计,小区计产面积30m2,试分析。
B
25
E
23
A
27
C
28
D
20
D
22
A
28
E
20
B
28
C
26
E
18
B
25
C
28
D
24
A
25
A
26
C
26
D
22
E
19
B
24
C
23
D
23
B
26
A
33
E
20
拉丁方设计的程序如下:
datali_4;
doi=1to5;
doj=1to5;
inputfz$x@@;
output;
end;
end;
cards;
B25E23A27C28D20
D22A28E20B28C26
E18B25C28D24A25
A26C26D22E19B24
C23D23B26A33E20
;
procanova;
classijfz;
modelx=ijfz;
run;
结果如下:
TheANOVAProcedure
DependentVariable:
x
Sumof
SourceDFSquaresMeanSquareFValuePr>F
Model12255.280000021.27333337.400.0008
Error1234.48000002.8733333
CorrectedTotal24289.7600000
R-SquareCoeffVarRootMSExMean
0.8810056.9585011.69509124.36000
SourceDFAnovaSSMeanSquareFValuePr>F
i48.56000002.14000000.740.5798
j444.560000011.14000003.880.0302
fz4202.160000050.540000017.59<.0001
结果分析:
列间,品种间差异是显著的。
3、左下表为玉米播期试验结果,右下表为油菜品比试验结果,皆为随机区组设计,试试分析。
播期
区组
Ⅰ
Ⅱ
Ⅲ
A
20.3
22.1
20.7
B
19.8
19.0
18.6
C
18.4
16.8
17.4
D
16.0
16.6
18.1
E
15.2
14.9
15.3
F
14.9
15.9
14.0
G
14.0
18.8
15.3
品种
区
组
Ⅰ
Ⅱ
Ⅲ
Ⅳ
A
3.9
3.9
3.6
4.5
B
5.8
6.3
5.9
6.8
C
4.4
4.4
5.6
4.5
D
5.5
5.2
5.4
6.7
E
6.8
6.9
7.4
6.0
F
7.3
7.2
7.5
7.0
玉米播期试验的产量(kg)结果油菜品比试验的产量(kg)结果
玉米的随机区组设计程序如下:
datali_3;
doi=1to7;
doj=1to3;
inputx@@;
output;
end;
end;
cards;
20.322.120.7
19.81918.6
18.416.817.4
1616.618.1
15.214.915.3
14.915.914
1418.815.3
;
procanova;
classij;
modelx=ij;
run;
结果如下:
Sumof
SourceDFSquaresMeanSquareFValuePr>F
Model890.694285711.33678577.610.0011
Error1217.87714291.4897619
CorrectedTotal20108.5714286
R-SquareCoeffVarRootMSExMean
0.8353427.0786301.22055817.24286
SourceDFAnovaSSMeanSquareFValuePr>F
i688.1714285714.695238109.860.0005
j22.522857141.261428570.850.4529
油菜随机区组设计程序如下:
datali_3;
doi=1to6;
doj=1to4;
inputx@@;
output;
end;
end;
cards;
3.93.93.64.5
5.86.35.96.8
4.44.45.64.5
5.55.25.46.7
6.86.97.46
7.37.27.57
;
procanova;
classij;
modelx=ij;
run;
结果如下:
Sumof
SourceDFSquaresMeanSquareFValuePr>F
Model831.275000003.9093750014.18<.0001
Error154.134583330.27563889
CorrectedTotal2335.40958333
R-SquareCoeffVarRootMSExMean
0.8832359.0977020.5250135.770833
SourceDFAnovaSSMeanSquareFValuePr>F
i530.817083336.1634166722.36<.0001
j30.457916670.152638890.550.6535
结果分析:
两播种期试验均表明,不同播种期对试验结果的影响是显著的。
4、下表为水稻栽培试验的小区产量(kg)结果,5×5拉丁方设计,,试分析。
。
B14
E15
C25
A12
D16
E18
D12
B15
C14
A11
C21
A13
D13
B13
E19
A10
C24
E18
D14
B12
D12
B15
A11
E20
C26
拉丁方设计程序如下:
datali_4;
doi=1to5;
doj=1to5;
inputfz$x@@;
output;
end;
end;
cards;
B14E15C25A12D16
E18D12B15C14A11
C21A13D13B13E19
A10C24E18D14B12
D12B15A11E20C26
;
procanova;
classijfz;
modelx=ijfz;
run;
结果如下:
Sumof
SourceDFSquaresMeanSquareFValuePr>F
Model12401.920000033.49333334.410.0078
Error1291.12000007.5933333
CorrectedTotal24493.0400000
R-SquareCoeffVarRootMSExMean
0.81518717.529262.75560015.72000
SourceDFAnovaSSMeanSquareFValuePr>F
i423.04000005.76000000.760.5717
j417.04000004.26000000.560.6954
fz4361.840000090.460000011.910.0004
结果分析:
品种间差异对结果的影响是显著的。
5、调查某队元麦及元麦和蚕豆混种、间种田块的产量(混、间种者为麦、豆产量合计),得结果(kg/0.1亩)于表:
(1)设以元麦单种为对照,试以LSD法作多重比较;
(2)设预定要作的比较是单种对混、间种,混种对间种,2麦1豆间种对3麦2豆间种,试作单一自由度的独立比较。
元麦单种
麦豆混种
2麦1豆间种
3麦2豆间种
20
24
30
30
24
23
28
33
22
28
34
31
18
21
32
36
21
24
31
35
先作方差分析:
datali_1;
doi=1to4;
doj=1to5;
inputx@@;
output;
end;
end;
cards;
2024221821
2423282124
3028343231
3033313635
;
procanova;
classi;
modelx=i;
meansi;
run;
结果如下:
Sumof
SourceDFSquaresMeanSquareFValuePr>F
Model3483.7500000161.250000028.04<.0001
Error1692.00000005.7500000
CorrectedTotal19575.7500000
R-SquareCoeffVarRootMSExMean
0.8402088.7996912.39791627.25000
SourceDFAnovaSSMeanSquareFValuePr>F
i3483.7500000161.250000028.04<.0001
TheANOVAProcedure
Levelof--------------x--------------
iNMeanStdDev
1521.00000002.23606798
2524.00000002.54950976
3531.00000002.23606798
4533.00000002.54950976
处理
均数
差数
均数-21
均数-24
均数-31
3麦2豆间种
33
12(极显著)
9(极显著)
2
2麦1豆间种
31
10(极显著)
7(极显著)
麦豆混种
24
3
元麦单种
21
6.有一大豆试验,A因素为品种,有A1、A2、A3、A44个水平,B因素为播期,有B1、B2、B33个水平,随机区组设计,重复3次,小区计产面积25平方米,其田间排列和产量(kg)如下图,试作方差分析。
检验:
品种、播期,品种×播期的效应是否显著?
区组Ⅰ
A1B1
12
A2B2
13
A3B3
14
A4B2
15
A2B1
13
A4B3
16
A3B2
14
A1B3
13
A4B1
16
A1B2
12
A3B1
14
A2B3
14
区组Ⅱ
A4B2
16
A1B3
14
A2B1
14
A3B3
15
A1B2
12
A2B3
13
A4B1
16
A3B2
13
A2B2
13
A3B1
15
A1B1
13
A4B3
17
区组Ⅲ
A2B3
13
A3B1
15
A1B2
11
A2B1
14
A4B3
17
A3B2
14
A2B2
12
A4B1
15
A3B3
15
A1B3
13
A4B2
15
A1B1
13
DATACaP;
DOa=1to4;
DOb=1to3;
DOn=1to3;
inputy@@;
output;
end;
end;
end;
DROPn;
CARDS;
121613
131415
141411
151514
131217
161314
141612
131315
161315
121513
141315
141713
;
PROCANOVA;CLASSAB;MODELy=ABA*B;
RUN;
结果如下:
TheANOVAProcedure
DependentVariable:
y
Sumof
SourceDFSquaresMeanSquareFValuePr>F
Model119.333333330.848484850.310.9778
Error2466.666666672.77777778
CorrectedTotal3576.00000000
R-SquareCoeffVarRootMSEyMean
0.12280711.904761.66666714.00000
SourceDFAnovaSSMeanSquareFValuePr>F
a32.888888890.962962960.350.7918
b20.500000000.250000000.090.9142
a*b65.944444440.990740740.360.8989
结果分析:
品种、播期,品种×播期对试验的影响均不显著。
7.有一小麦裂区试验,主区因素A,分A1(深耕)、A2(浅)两水平,副区因素B,分B1(多肥)、B2(少肥)两水平,重复3次,小区计产面积15平方米,其田间排列和产量(假设数字)如下图,试作方差分析。
因素A,因素B及其交互效应是否显著?
A1
A2
A2
A1
A2
A1
B1
9
B1
7
B2
3
B1
11
B2
1
B2
4
B2
6
B2
2
B1
5
B2
4
B1
6
B1
12
区组Ⅰ
区组Ⅱ
区组Ⅲ
解:
程序如下:
title'裂区试验的统计分析';
datalq;
doa=1to2;
dob=1to2;
dog=1to3;
inputyield@;output;
end;
end;
end;
cards;
91112
644
756
231;
procanova;
classgab;
modelyield=gab;
meansab/duncan;
run;
procanova;
classgab;
modelyield=gaba*bg*ag*a*b;
testH=gaE=g*a;
testH=ba*bE=g*a*b;
run;
试验结果:
裂区试验的统计分析208:
07Wednesday,July3,2002
AnalysisofVarianceProcedure
DependentVariable:
YIELD
SourceDFSumofSquaresFValuePr>F
Model4115.5000000014.270.0018
Error714.16666667
CorrectedTotal11129.66666667
R-SquareC.V.YIELDMean
0.89074624.387545.83333333
SourceDFAnovaSSFValuePr>F
G20.166666670.040.9599
A140.3333333319.930.0029
B175.0000000037.060.0005
裂区试验的统计分析608:
07Wednesday,July3,2002
AnalysisofVarianceProcedure
DependentVariable:
YIELD
SourceDFSumofSquaresFValuePr>F
Model11129.66666667..
Error0.
CorrectedTotal11129.66666667
R-SquareC.V.YIELDMean
1.00000005.83333333
SourceDFAnovaSSFValuePr>F
G20.16666667..
A140.33333333..
B175.00000000..
A*B13.00000000..
G*A21.16666667..
G*A*B410.00000000..
TestsofHypothesesusingtheAnovaMSforG*Aasanerrorterm
SourceDFAnovaSSFValuePr>F
G20.166666670.140.8750
A140.3333333369.140.0142
TestsofHypothesesusingtheAnovaMSforG*A*Basanerrorterm
SourceDFAnovaSSFValuePr>F
B175.0000000030.000.0054
A*B13.000000001.200.3349
结果分析:
因素A、因素B的显著性水平为0.0029、0.0005,是显著的,而AXB的F值是1.2,顾AXB是不显著的。
B1与B2有显著的差别,A1与A2有显著的差别,最好的是A1B1的组合。
8.设若上题小麦耕深与施肥量试验为条区设计,田间排列和产量将相应如下图,试作分析,并与裂区设计结果相比较)。
A1
A2
A2
A1
A2
A1
B1
9
7
B1
5
11
B2
1
4
B2
6
2
B2
3
4
B1
6
12
解:
程序如下:
title'条区试验的统计分析';
datatq;
dog=1to3;
doa=1to2;
dob=1to2;
inputyield@;output;
end;
end;
end;
cards;
9672
11453
12461;
pro