免费文献传递   相关文献

Correlation and regression analysis of trait and yield of soybean

大豆植株性状相关性与产量回归分析



全 文 :中国生态农业学报 2008 年 11 月 第 16 卷 第 6 期
Chinese Journal of Eco唱Agriculture,Nov.2008,16(6):1429唱1433   
DOI: 10.3724/SP.J.1011.2008.01429
 
大豆植株性状相关性与产量回归分析倡 倡
韩秉进 潘相文 金 剑 王光华 刘长江 刘晓冰
(中国科学院东北地理与农业生态研究所 黑土生态重点实验室 哈尔滨 150081)
摘 要 本文利用 3 个大豆品种、4 种施肥模式条件下,植株产生的 17 个形态指标及产量差异,进行相关性分
析与产量回归分析。 结果表明,株高 x1、分枝数 x5、1 粒荚数 x9、2 粒荚数 x10、3 粒荚数 x11、4 粒荚数 x12、百粒重
x15、经济系数 x16等 8 个大豆植株性状因子对产量有显著性影响;由标准化回归系数得知,对产量作用的大小
顺序为:x11 >x12 >x9 >x10 >x16 >x15 >x1 >x5。 说明提高大豆产量首先要考虑增加有效荚数,其次是提高经济
系数,而增加百粒重和株高对于增产的效应较小。 育种上需加强生物量大、3 粒荚和 4 粒荚数多、结荚节位低
的性状选择;栽培上需注重促进生物量生长、降低结荚节位的促控措施。
关键词 大豆 植株性状 相关性 产量回归
中图分类号:S513  文献标识码:A  文章编号:1671-3990(2008)06-1429-05
Correlation and regression analysis of trait and yield of soybean
HAN Bing唱Jin, PAN Xiang唱Wen, JIN Jian, WANG Guang唱Hua, LIU Chang唱Jiang, LIU Xiao唱Bing
(Key Laboratory for Black Soil Ecology, Northeast Institute of
Geography and Agro唱ecology, Chinese Academy of Sciences, Harbin 150081, China)
Abstract Correlation and regression analyses were performed based on variations in 17 morphological traits and yield obtained
from field experiments on three soybean varieties with four fertilization modes.Results show that eight morphological traits inclu唱
ding plant height (x1 ), branch number ( x5 ), uni唱grain pod ( x9 ), di唱grain pod ( x10 ), tri唱grain pod ( x11 ), tetra唱grain pod
( x12), 100 -grain weight (x15) and economic coefficient (x16) have significant impact on yield with relative effects ( from stand唱
ard regression coefficient) in the following order: x11 >x12 >x9 >x10 >x16 >x15 >x1 >x5.Enhancing effective pod number is
an important first step for improving soybean yield, followed by improving economic coefficient.In comparison, increasing 100 -
grain weight and plant height are less significant for improving yield.Selection of traits with high biomass, more tri唱grain pod, tet唱
ra唱grain pod and low podding node site should be emphasized in breeding and developing regulation measures to enhance soybean
biomass.Descending podding node site should also be adopted in soybean cultivation.
Key words Soybean, Plant trait, Correlation, Yield regression
(Received Sept.19, 2007; accepted Dec.30, 2007)
  选育高产优质的作物品种、采用高效的栽培技
术,一直是农业科研工作者研究的主题[ 1 -7] 。 多年
以来,人们为了解作物植株相关性状与产量的关
系,以便做到更有效的性状选择和采取相应的栽培
措施,进行了一系列的研究[ 8 -19] 。 关于作物产量与
施肥、土壤条件、气候因子等的关系研究报道也较
多 [20 -26] 。 但在作物植株诸多性状与产量相关性分
析基础上,剔除自相关性强的自变量、建立较高精
度和稳定产量回归方程分析的报道较少[27 -31] 。 本
研究利用相关分析,剖析大豆植株性状间的相互关
系及其对产量的影响,并通过共线性诊断、逐步回
归的方法,建立产量回归模型,进行产量分析,对各
回归变量的作用大小做出评价,从中找出产量的主
导性状因子,为大豆高产育种中各农艺性状的选择
和高产栽培技术提供理论依据。
1 材料与方法
1.1 数据来源与处理
数据来源于黑龙江省重点攻关课题“大豆超高
产试验示范与配套技术研究”。 试验设 3 个大豆品
种、4 种施肥模式与用量,共 24 组考种数据,每组数
据是 15 株大豆植株的平均值。 设计包含有品种、
施肥上的差异,使包括产量在内的植株性状指标充
倡 黑龙江省重点科技攻关课题(GA06B101 -3)和中国科学院东北地理与农业生态研究所学科前沿领域项目(KZCX3 -SW -NA3 -22)资助
韩秉进(1957 ~) ,男,博士,研究员,主要从事作物高产栽培理论与技术研究。 E唱mail:hanbj@neigaehrb.ac.cn
收稿日期:2007-09-19 接受日期:2007-12-30
中国生态农业学报 2008 第 16卷
分表达。 将这 24 组考种测定的 17 个植株性状指
标,采用计算机 SAS9.0 统计软件处理,进行相关性
和多元回归分析。 多元回归首先利用共线性诊断
去掉了自变量中自相关性强、作用较小的变量,使
保留下来的所有自变量之间尽可能相互独立,将剩
余的自变量性状与因变量产量性状,采用逐步回归
法进行多元回归分析[ 32,33] 。
1.2 调查项目
调查项目包括株高、结荚部位、主茎节数、结荚
节位、分枝数、地上干重、空瘪荚、1 粒荚数、2 粒荚
数、3 粒荚数、4 粒荚数、有效荚数、每荚粒数、百粒
重、经济系数、单株粒数、单株粒重。
2 结果与分析
2.1 大豆植株性状间的相关性
由表 1 可以看出,单株产量与单株粒数、有效
荚数、地上干重、3 粒荚数、每荚粒数、4 粒荚数、经
济系数、百粒重都呈极显著的正相关,与分枝数呈
显著正相关,而与结荚节位呈极显著的负相关。 说
明要获得较高的产量,首先要考虑增加单株粒数、
增加百粒重、提高经济系数、降低结荚节位;通过增
加有效荚数、3 粒荚数和 4 粒荚数,从而增加每荚粒
数和单株粒数。 从表 1 还可看出,大豆的百粒重 x15
与分枝数、地上干重、3 粒荚数、有效荚数、经济系数
呈显著或极显著正相关,与结荚部位、结荚节位呈
显著或极显著负相关。 欲提高百粒重,栽培上需考
虑增加 3 粒荚数和地上干重、降低结荚部位性状的
措施,育种上则要注意相关性状的选择。 经济系数
与 3 粒荚数、有效荚数呈显著或极显著正相关,而
与株高、结荚部位、主茎节数、结荚节位成显著或极
显著负相关。 欲提高品种的经济系数,在育种性状
选择上要考虑 3 粒荚数多、低株高、低结荚部位性
状的选择。 3 粒荚数和 4 粒荚数都与地上干重极显
著正相关,与结荚节位显著或极显著负相关,因此
欲增加 3 粒荚数和 4 粒荚数,育种上也需加强生物
量大、结荚节位低的性状选择,栽培上需注重相应的
促进生物量生长、降低结荚节位的促控措施。
2.2 共线性诊断
当回归模型中自变量之间存在高度线性关系
时,会造成模型估计量的稳定性差,必须通过共线
性诊断,排除自相关性强的自变量。 本研究以产量
为因变量,其他 16 个性状为自变量,输入 SAS 程序
进行共线性诊断。 结果是条件指数 K =64 784,远
远大于 30,而在变量 1 粒荚数、2 粒荚数、3 粒荚数、
4 粒荚数、有效荚数等 5 个性状上的方差比率都在
0.99 以上,远超过 50%,说明这 5 个自变量性状之
间有严重的多重共线性[33] 。 去掉 1 个综合性性状
有效荚数,再次诊断。 这样反复进行,每次去掉 1
个综合性强或作用较小的变量。 共进行了 6 次诊
断,先后剔除了有效荚数、单株粒数、每荚粒数、地
上部干重、结荚部位等 5 个变量。 诊断后保留下来
的 11 个自变量性状之间基本相互独立,然后利用
这 11 个自变量性状因子为自变量,产量为因变量
进行逐步回归筛选变量,进行多元回归分析。
2.3 变量筛选及回归方程建立
经过共线性诊断保留下来的 11 个自变量,对
因变量的影响不一定都有显著性意义,需要根据各
自变量对因变量的贡献大小进行变量筛选,从而求
出精炼、稳定的回归方程。 本研究筛选变量的方法
是采用最常用的逐步回归法。 模型中的变量从无
到有,根据各变量的 F 统计量及 P 值决定是否进入
模型。 规定选变量进入方程的显著性水平为 SLE =
0.3,当 P <0.3,则该变量入选,否则不能入选;规定
从方程中剔除变量的显著性水平为 SLS =0.1,当 P
<0.1,则将此变量保留在方程中,否则将变量从方
程中剔除[20,21,32] 。 并要求求出标准化回归参数的
估计值(STB)。
经过逐步回归法筛选,11 个自变量中有 8 项指
标对产量的贡献达到显著水平,其余指标的作用均
不显著。 表 2 列出了逐步回归法选出的自变量与
系数。 产量回归方程为:
y =-11.559 5 +0.024 9x1 -1.260 1x5 +
0.383 7x9 +0.152 6x10 +0.439 8x11 +0.336 9x12 +
0.406 9x15 +7.767 7x16  (n =24,R =0.992 9,
F =131.14,Pr ﹤ 0.000 1) (1)
2.4 模型检验及解释
为直观检验产量回归方程的拟合程度,绘制了
残差图(图 1)。 由因变量预测值与残差的标准差的
关系(残差图)可以看出,方程的拟合较好,实际观
测值与预测值的拟合检验时,尽管标准化残差中尚有
1 个大于 2 的点,但模型的决定系数 R2 =0.985 9,检
验统计量 F =131.14,Pr﹤ 0.000 1,达到极显著水
平,剩余标准差 S 也仅为 0.058 6 。 说明模型的精
确度较高,98.59%大豆产量变异是株高、分枝数、
有效荚数、百粒重和经济系数影响的结果 [ 34,35 ] 。
因此,可以实际应用模型进行因变量与自变量关
系的解析,对各回归变量的作用大小做出评价,利
用回归方程对因变量进行预测、对自变量进行控
制等。
由表 2 和回归方程可以看出,株高 x1 、分枝数
x5 、1 粒荚数 x9、2 粒荚数 x10 、3 粒荚数 x11 、4 粒荚数
x12 、百粒重 x15 、经济系数 x16等 8 个大豆植株性状因
子对产量有显著性影响。 方程中的回归系数表明,
0341

1 












Tab
.1 
Cor
rela
tion
coef
ficie
nta
mon
gso
ybe
anp
lant
trait


目 Item

高 Pla
nt
hei
ght x 1




Low
est
pod
din

hei
ght x 2




Nod
es
num
ber
of
mai
ns
tem x 3




Nod
esi
te
ofp
odd
ing
x 4



Bra
nch
ing
num
ber x 5




Sho
otd
ry
wei
ght x 6




Nil
唱gra
in
pod num
ber x 8
1粒


Uni
唱gra
in
pod num
ber x 9
2粒


Di唱
gra
in
pod num
ber x 10
3粒


Tri唱
grai

pod num
ber x 11
4粒


Tet
ra唱g
rain
pod num
ber x 12




Effe
ctiv

pod num
ber x 13




Gra
in
num
ber
per
pod x 14



100
唱gra
in
wei
ght x 15




Eco
nom
ical
coe
ffic
ien

x 16




Gra
in
num
ber
per
pla
nt
x 17




Gra
in
wei
ght
per
pla
nt
x 18
x 1
1.0
00

x 2
0.8
87
7倡倡
1.0
00

x 3
0.4
55
7倡
0.4
48
5倡
1.0
00

x 4
0.2
48

0.4
39
2倡
0.8
03
1倡倡
1.0
00

x 5
-0
.02
99
-0
.15
89
-0
.25
83
-0
.37
85
1.0
00

x 6
0.4
36
7倡
0.0
88

-0
.00
90
-0
.36
06
0.3
84

1.0
00

x 8
-0
.02
04
-0
.00
79
-0
.06
19
-0
.11
60
0.4
09
5倡
0.0
21

1.0
00

x 9
-0
.29
36
-0
.25
35
0.1
55

0.2
75

0.0
43

-0
.12
48
0.0
37

1.0
00

x 10
-0
.10
34
-0
.18
65
0.1
24

0.1
15

0.2
35

0.1
67

0.5
93
8倡倡
0.5
19
4倡倡
1.0
00

x 11
0.2
17

-0
.11
53
-0
.26
81
-0
.59
77
倡倡
0.5
14
8倡
0.8
87
9倡倡
0.0
13

-0
.31
44
0.0
06

1.0
00

x 12
0.2
79

0.0
68

-0
.15
10
-0
.43
31

0.2
01

0.7
43
2倡倡
-0
.27
38
-0
.34
56
-0
.38
06
0.7
21
2倡倡
1.0
00

x 13
0.1
10

-0
.20
51
-0
.11
90
-0
.40
09
0.5
25
3倡倡
0.8
74
1倡倡
0.2
29

0.1
88

0.5
12
0倡
0.8
08
0倡倡
0.5
23
4倡倡
1.0
00

x 14
0.2
86

0.0
39

-0
.25
49
-0
.55
33
倡倡
0.2
93

0.7
51
2倡倡
-0
.17
88
-0
.59
88倡

-0
.41
17

0.8
61
6倡倡
0.8
98
7倡倡
0.4
99
4倡
1.0
00

x 15
-0
.37
33
-0
.58
90
倡倡
-0
.31
75
-0
.41
42

0.5
16
4倡倡
0.4
08
4倡
0.1
01

0.0
11

0.2
50

0.5
49
2倡倡
0.2
60

0.5
67
9倡倡
0.3
49

1.0
00

x 16
-0
.51
46
倡 -
0.6
48
0倡倡
-0
.52
60
倡倡
-0
.56
34
倡倡
0.3
73

0.3
21

-0
.02
38
0.2
92

0.2
06

0.4
96
6倡
0.3
00

0.5
93
1倡倡
0.2
78

0.7
19
8倡倡
1.0
00

x 17
0.2
03

-0
.12
96
-0
.18
48
-0
.51
50

0.5
05
5倡
0.9
46
3倡倡
0.0
98

-0
.09
64
0.2
16

0.9
36
8倡倡
0.7
51
2倡倡
0.9
39
4倡倡
0.7
64
4倡倡
0.5
52
6倡倡
0.5
43
4倡倡
1.0
00

x 18
0.1
33

-0
.20
53
-0
.21
60
-0
.51
76
倡倡
0.4
61
6倡
0.9
21
3倡倡
0.0
08

0.0
18

0.2
24

0.9
08
2倡倡
0.6
97
2倡倡
0.9
34
1倡倡
0.7
00
5倡倡
0.6
26
2倡倡
0.6
59
8倡倡
0.9
68
7倡倡
1.0
00

  
倡表


5%




;倡倡



1%




。倡
mea
ns
sign
ific
ant
diff
ere
nce
at0
.05
;倡倡
mea
ns
sign
ific
ant
diff
ere
nce
at0
.01

中国生态农业学报 2008 第 16卷
表 2 逐步回归法选出的自变量与系数
Tab.2 Parameters and coefficients selected from the method of stepwise regress analysis
变量 自由度 系数 标准误 T 值 显著水平 Significant level 标准化回归系数
 Variable DF Coefficient Standard error T value Pr > |t| Standardized regress coefficient
 截距 Intercept 1 -11.559 5 2.412 7 -4.79 0.000 2 0.000 0
 x1 1 0.024 9 0.009 7 2.56 0.021 6 0.130 8
 x5 1 -1.260 1 0.472 2 -2.67 0.017 5 -0.106 7
 x9 1 0.383 7 0.089 0 4.31 0.000 6 0.213 9
 x10 1 0.152 6 0.038 2 4.00 0.001 2 0.184 6
 x11 1 0.439 8 0.051 3 8.58 <0.000 1 0.644 1
 x12 1 0.336 9 0.068 8 4.90 0.000 2 0.276 7
 x15 1 0.406 9 0.158 9 2.56 0.021 8 0.137 6
 x16 1 7.767 7 3.257 7 2.38 0.030 7 0.164 2
8 个大豆植株性状对产量增产效应大小顺序为:
x16 >x11 >x15 >x9 >x12 >x10 >x1;其中作用最大的是
经济系数 x16 ,作用最小的是株高 x1 ;而分枝数 x5 对
产量有负向效应。 由于方程中的偏回归系数受单
位的干扰,要比较各自变量对因变量产量 y 作用的
大小,需要用标准化回归系数[ 29,32,33] 。 由表 2 可以
看出,标准化回归系数大小顺序为:x11 >x12 >x9 >
x10 >x16 >x15 >x1 。 说明增加大豆有效荚数、经济系
数、百粒重和株高对于增产的效应依次减小。 而 x5
分枝对产量起到了负向作用,说明该地区适合栽培
主茎型大豆品种而不是分枝型大豆品种,主茎型大
豆品种的合理密植是获取产量的重要因素。
图 1 产量回归分析残差图
Fig.1 Standard residual plot of regress analysis on yield
3 讨论
大豆的单株产量与单株粒数、有效荚数、地上
干重、3 粒荚数、每荚粒数、4 粒荚数、经济系数、百
粒重和分枝数都显著或极显著正相关。 说明要获
得较高的产量,首先要考虑增加单株粒数、百粒重,
提高经济系数;通过增加有效荚数、3 粒荚数、4 粒
荚数,从而增加每荚粒数和单株粒数。 大豆的百粒
重与分枝数、地上干重、3 粒荚数、有效荚数、经济系
数呈显著或极显著正相关,与结荚部位、结荚节位
呈显著或极显著负相关,因此欲提高百粒重,栽培
上需考虑增加 3 粒荚数和地上干重、降低结荚部位
性状的措施,育种上则要注意相关性状的选择。
从大豆的单株产量与诸多性状逐步回归选元
结果分析,影响大豆产量的主要因子是经济系数、
分枝数、有效荚数(2 粒荚数、1 粒荚数、3 粒荚数、4
粒荚数)、百粒重、株高;从多元回归的标准化回归
系数分析,欲提高大豆产量,首先要考虑增加大豆
有效荚数(尤其是增加 3 粒荚数和 4 粒荚数),其次
是提高经济系数,而增加百粒重和株高对增产的效
应较小。
由于大豆的 3 粒荚数和 4 粒荚数都与地上干重
极显著正相关,与结荚节位显著或极显著负相关,
因此欲增加 3 粒荚数和 4 粒荚数,育种上也需加强
生物量大、结荚节位低的性状选择;栽培上需注重
采取相应的促进生物量生长、降低结荚节位的促控
措施。 由于经济系数与 3 粒荚数、有效荚数呈显著
或极显著正相关,而与株高、结荚部位、主茎节数、
结荚节位呈显著或极显著负相关关系,因此欲提高
品种的经济系数,在育种性状选择上要考虑增加 3
粒荚数和降低株高、降低结荚部位性状的选择。 研
究表明,大豆分枝对产量起到了负向作用,说明该
地区更适合栽培主茎型大豆品种,配合合理密植是
获取较高产量的有效措施。
参考文献
[1] 赵团结, 盖钧镒, 李海旺, 等.超高产大豆育种研究的
进展与讨论[J].中国农业科学,2006,39(1):29 -37
[2] 满为群,杜维广,张桂茹, 等.高光效大豆几项光合生
理指标的研究[J].作物学报, 2003,29(5):697 -700
2341
第 6期 韩秉进等:大豆植株性状相关性与产量回归分析
[3] 王铁军, 张海泉.大豆育种工作的现状与展望[ J].沈
阳农业大学学报,2000,31(4):375 -379
[4] Rosenzweig V.E., Goloenko D.V., Davydenko O.G.,
et al.Breeding strategies for early soybeans in Belarus
[J].Plant Breeding,2003, 122: 456 -458
[5] 李新民.大豆群体结构的研究[ J].大豆科学,1990,
9(3):185 -190
[6] 刘忠堂.大豆窄行密植高产栽培技术的研究[J].大豆
科学,2002,21(2):117 -122
[7] Cui S.Y., Yu D.Y.Estimates of relative contribution of
biomass, harvest index and yield components to soybean
yield improvements in China[ J].Plant Breeding, 2005,
124 (5): 473 -476
[8] 战勇,罗赓彤,刘胜利,等.春大豆性状的多元遗传分析
及应用分析[J].新疆农业科学, 2000 (2):55 -59
[9] 张玉先,张瑞朋.主成分分析在大豆产量与相关性状中
的应用[J].河北农业科学,2004,8(3):107 -109
[10] 宁海龙,李文霞,潘相文,等.大豆氨基酸组分影响蛋
白质含量的通径分析[ J].大豆科学,2002,21 (4):
259 -262  
[11] 宁海龙,李文霞,王继安,等.黑龙江省大豆蛋白质油
分及蛋白质组分类型[ J].作物学报,2003,29 (4):
551 -556 
[12] 孙宪印,吴科,钱兆国,等.黄淮冬麦区北片水地组供
试小麦品种(系)主要品质性状的主成分分析和聚类
分析[J].山东农业科学, 2006 (1):24 -26
[13] 宋志伟,杨首乐.春性小麦品种主要农艺性状与产量
的相关及通径分析[J].中国农学通报, 2006,22(5):
174 -176
[14] 宋启建,吴天侠,邱家驯,等.夏大豆群体结构对不同
类型品种产量及农艺性状的影响[ J].大豆科学,
1995,14(1):40 -46
[15] 余泽高, 方燕妮.小麦经济系数与其它若干性状关系
的研究[J].湖北农学院学报, 2002, 22(6):486 -493
[16] 季彪俊,孙依斌,郑喻.旱作对籼粳稻产量及相关性状
的影响[J].中国农学通报, 2006,22(7):178 -182
[17] 庄萍萍,李伟,魏育明,等.波斯小麦农艺性状相关性
及主成分分析 [ J].麦类作物学报,2006,26 (4):
11 -14 
[18] Board J.E., Kang M.S., Harville B.G.Path analyses
of the yield formation process for late唱planted soybean
[J].Agronomy Journal, 1999, 91:128 -135
[19] Machikowa T., Waranyuwat A., Laosuwan P.Relation唱
ships between seed yield and other characters of different
maturity types of soybean grown in different environments
and levels of fertilizer [ J].Science Asia, 2005, 31:
37 -41 
[20] 季彪俊.影响水稻产量因子的研究[ J].西南农业大
学学报:自然科学版,2005,27(5):579 -583.
[21] 付凤玲,周树峰,潘光堂,等.玉米耐旱系数的多元回
归分析[J].作物学报,2003,29(3):468 -472
[22] 李卫东,卢为国,梁慧珍,等.大豆蛋白质含量与生态
因子关系的研究 [ J].作物学报, 2004,30 ( 10 ):
1605 -1068 
[23] Lobell D.B., Ortiz唱Monasterio J.I., Asner G.P., et
al.Analysis of wheat yield and climatic trends in Mexico
[J].Field Crops Research, 2005, 94: 250 -256
[24] Kalivas D.P., Kollias V.J.Effects of soil, climate and
cultivation techniques on cotton yield in Central Greece,
using different statistical methods[J].Agronomie, 2001,
21:73 -89
[25] Cox M.S., Gerard P.D., Wardlaw M.C., et al.Vari唱
ability of selected soil properties and their relationships
with soybean yield[ J].Soil Science Society of America
Journal, 2003,67: 1296 -1302
[26] 韩秉进.松嫩平原黑土区玉米生产氮磷配合肥效优化
模型的研究[J].土壤学报, 1998,35(3):392 -397
[27] 谭禾平,王桂跃,胡贤女,等.影响玉米产量效应因子
的多元回归与通径分析[ J].浙江农业学报, 2006,
18(4):238 -240
[28] 黄大辉,彭懿紫,黄天进,等.杂交水稻主要性状的多重
逐步回归和通径分析[ J].广西农业生物科学,2004,
23 (2) :100 -103
[29] 刘廷辉.食用大麦八个农艺性状与其产量的多元回归
分析[J].大麦科学,2004 (1):7 -9
[30] Board J.E.A regression model to predict soybean culti唱
var yield performance at late planting dates[ J].Agrono唱
my Journal, 2002, 94(3): 483 -492
[31] Natarajan S., Kumaresan G.R., Sasikala S.Correlation
and regression analysis among some agronomic characters
in soybean[J].Res.on Crops, 2006, 7 (1): 363 -365
[32] 胡良平.现代统计学与 SAS应用[M].北京:军事医
学科学出版社,2002:248 -263
[33] 岳朝龙,黄永兴,严忠.SAS 系统与经济统计分析
[M].合肥:中国科学技术大学出版社, 2003:
265 -310 
[34] 马育华.试验统计 [M].北京:农业出版社,
1998:272 -318
[35] 陶勤南.回归分析与回归设计[ J].北京农业科学,
1984 (专集):1 -76
3341