, 尚一斐* 通讯作者(Corresponding author): 魏永胜, E-mail:wysh70@nwsuaf.edu.cn 第一作者联系方式: E-mail:zhenglifei@nwsuaf.edu.cn
收稿日期:2016-12-07 接受日期:2017-05-10网络出版日期:2017-06-05基金:
摘要
关键词:冬小麦; 产量; 农艺性状; 结构方程模型
Structural Equation Model for Analyzing Relationship between Yield and Agronomic Traits in Winter Wheat
ZHENG Li-Fei
, SHANG Yi-FeiFund:
Abstract
Keyword:Winter wheat; Yield; Agronomic traits; Structural equation model
Show Figures
Show Figures
产量是小麦生产的主要目标之一, 不仅受到遗传因素的影响, 也受到土壤肥力、气候条件以及农耕措施等因素的影响。产量构成三要素穗数、穗粒数和千粒重一直是高产育种栽培的研究重点, 其分析模型与方法也多种多样, 如多元线性回归[1, 2]、通径分析[3, 4, 5, 6, 7, 8]、主成分分析(principal component analysis, PCA)[9]、简单相关分析[2, 5, 8, 10]和典型相关分析[11], 这些分析方法的根本目的是揭示冬小麦主要农艺性状与产量的关系, 为高产育种提供理论依据。
目前, 有关产量构成三要素对籽粒产量的贡献已有大量报道, 因为采用品种、样本量及试验条件不同, 其结论也不一致, 在特定条件下千粒重[3, 4, 7]、穗粒数[5]或单位面积穗数[6]都可能成为产量的主要贡献者。另外, 某些品种株高则是产量的主要决定因素[12]。为进一步明确冬小麦性状与产量的关系, 我们利用结构方程模型(structural equation model, SEM, 有时covariance structure analysis, covariance structure modeling and analysis of covariance structures[13])作为一种多变量分析方法。最初该法主要用于社会科学领域, 后来应用于生态学中复杂问题的解析[14]。Lamb等[15]将SEM初步用于作物分析。相对于“ 第一代” 多变量统计法PCA和聚类分析(cluster analysis, CA)等, SEM可称为“ 第二代” 多变量统计法, 能够同时分析系统内多个变量间的因果关系, 并明确给出各关系的强度大小[16]。本研究以19个冬小麦品种在19个不同地点的379套产量与主要农艺性状数据, 利用SEM分析主要农艺性状对产量的直接作用和间接作用, 以及各性状之间的关系, 为冬小麦高产品种的选育提供理论依据。
1 数据来源与分析方法1.1 数据来源数据来自2010— 2011年度国家冬小麦品种试验中长江上游组区域试验, 共19个品种或稳定品系, 分别是15610、07-129、07间551-55、08RC2525、09J6白、09J76、BL228、BL379、CD09-2542、SW20812、川07217、川08品32(W)、川09023、川麦42 (2次对照)、川麦53、川重组104、绵06-374、太空7号和云麦59, 这些材料均为各育种单位选育的遗传性稳定的品系, 依据农业部颁发的《农作物品种(小麦)区域试验技术规程》新标准进行各试点的小区设计、田间管理及性状调查。选取19个栽培点的379条数据。栽培点为甘肃省成县, 贵州省毕节县、贵阳市、铜仁市、遵义市、湖北省襄阳市、陕西省安康市、勉县, 重庆市北碚区、万州区、永川区, 四川省成都市、绵阳市、内江市、平昌县、西昌市及云南省楚雄州、德宏州、曲靖市。调查性状共7个, 分别是生育期、基本苗、单位面积穗数、株高、千粒重、穗粒数及产量。
1.2 数据标准化为了消除数据单位对分析的影响, 原始数据进行Z-score标准化处理。
$Z(x)=\frac{x-\bar{x}}{S(x)}$
式中, Z(x)为新的数据, x为原始数据, $\bar{x}$为x均值, S(x)为原始数据的标准差。
利用IBM SPSS AMOS 21 (Amos Development Corporation, Chicago, IL, USA)建立结构方程模型, 并利用SPSS 16软件进行数据基本特征、正态性检验和相关性分析。
1.3 初始结构方程模型结构方程模型是一种先验模型[14], 需要先在实践经验与理论的基础上建立假设路径, 再通过结构方程模型分析、验证和优化各因素间的直接和间接关系, 最终确定因素间合理的因果关系[15]。本研究的初始模型如图1所示。
图1
Fig. 1
| Figure OptionViewDownloadNew Window | |
 | 图1 冬小麦主要农艺性状与产量关系的初始模型 图内共25个箭头, 表示需估计25个参数; e1~e7为相应的误差项。GY: 产量; GNP: 穗粒数; BS: 基本苗; SN: 单位面积穗数; GD: 生育期; TGW: 千粒重; PH: 株高。Fig. 1 Initial modeling for the relationship between yield and agronomic traits in winter wheat The initial model contains 25 parameters shown by arrows, e1 to e7 indicate errors. GY: grain yield; GNP: grain number per spike; BS: density of basic seedlings; SN: spike number per ha; GD: growth duration; TGW: thousand-grain weight; PH: plant height. |
根据文献[6, 12, 15]及多年育种经验, 初始模型中假设千粒重、株高、穗粒数、生育期、基本苗、单位面积穗数对产量都有直接影响; 株高、生育期、单位面积穗数和穗粒数对千粒重有直接影响; 株高、生育期和单位面积穗数对穗粒数有直接影响; 生育期和基本苗对株高和单位面积穗数有直接影响, 因为生育期包括冬前生长和次年生长, 因此假设生育期长可得到更多分蘖, 进而获得更多的穗数。设定初始模型后需要估计25个参数(图1), 少于可以建立的方程总数[7(7+1)/2=28], 因此该模型是可以识别的。
在结构方程模型分析之前, 作为比较, 采用向后逐步回归方法对数据进行多元回归分析, 根据P值大小逐一剔除变量, P > 0.05则剔除该变量。本研究的变量P值均小于0.05, 因此只进行一步回归分析。
1.4 模型检验采用χ 2检验整体评价模型[17], χ 2值越小, 模型的拟合效果越好。若P> 0.05, 则表明模型没有遗漏重要参数。本研究中共有379个数据点, 满足结构方程模型χ 2检验所需要的最小推荐样本数, 即模型参数的5倍以上[14]。用RMSEA (root mean square error of approximation)和AIC (akaike information criterion)辅助检测模型的拟合度[18]。RMSEA < 0.05, 模型拟合非常好; 0.05 ≤ RMSEA < 0.08, 模型拟合较好; 0.08 ≤ RMSEA < 0.10, 模型拟合一般; RMSEA ≥ 0.10, 模型拟合不好。AIC值越小, 模型的拟合度越高。
2 结果与分析2.1 冬小麦数据的基本统计量如表1所示, 最高产量为SW20812的9270.0 kg hm-2, 出现在云南德宏试点, 而在重庆永川点的云麦59产量最低, 为1861.5 kg hm-2, 极差达7408.5 kg hm-2, 变异系数是所有性状中最大的。千粒重变异范围为27.3~62.2 g, 最低为种植于四川成都的15610, 最高为种植于四川西昌的07间551-55。株高分布范围为47~118 cm, 最低为种植于甘肃成县的07-129, 最高为种植于重庆永川的云麦59。穗粒数最少的是种植于四川平昌的07间551-55, 为每穗20.6粒, 最多的是种植于云南德宏州的07-129, 为每穗72.9粒。生育期变幅为142~232 d, 是测定性状中变异最小的, 除遗传特性外, 与试点纬度等关系密切, 其中在甘肃成县生育期最长(230.0± 0.8 d), 在云南德宏生育期最短(154.0± 5.4 d)。基本苗最低的是种植于甘肃省成县的15610, 为每公顷102万株, 最高的是种植于四川省西昌市的川重组104, 为每公顷292.5万株。单位面积穗数最低的为种植于贵州省铜仁市的BL228, 为每公顷240万穗, 最高的为种植于湖北省襄阳市的川07217, 为每公顷630万穗。
表1
Table 1
表1(Table 1)
 表1 本研究中冬小麦各性状的基本统计量 Table 1 Basic statistics of winter wheat traits used in this study (n = 379)
| 表1 本研究中冬小麦各性状的基本统计量 Table 1 Basic statistics of winter wheat traits used in this study (n = 379) |
2.2 冬小麦主要农艺性状与产量的相关性分析除株高外的其余性状与产量的相关性均达到显著水平, 其中穗粒数、基本苗和单位面积穗数与产量的相关系数较高, 生育期与产量呈显著负相关; 除产量外的各农艺性状之间也存在显著相关(表2)。因此, 有必要对各农艺性状与产量的直接和间接关系作进一步分析。
表2
Table 2
表2(Table 2)
表2 冬小麦性状的相关系数 Table 2 Correlation coefficients among winter wheat traits
| 表2 冬小麦性状的相关系数 Table 2 Correlation coefficients among winter wheat traits |
2.3 冬小麦主要农艺性状与产量的回归分析冬小麦产量(Ŷ )与千粒重(x1)、株高(x2)、穗粒数(x3)、生育期(x4)、基本苗(x5)、单位面积穗数(x6)的多元回归方程为Ŷ = 109.5x1 + 15.1x2 + 99.8x3 + 12.8x4+ 12.7x5 + 11.4x6 -13720.1 (R2 = 0.719), 自变量的标准化系数分别为0.440、0.125、0.557、0.156、0.301和0.601。因此, 各性状对冬小麦产量的影响, 由大到小依次为单位面积穗数、穗粒数、千粒重、基本苗、生育期和株高。
2.4 冬小麦主要农艺性状与产量的结构方程模型分析初始模型拟合结果, χ 2为36.665, 自由度(df)为4, χ 2/df为9.164, P< 0.05, 而RMR和RMSEA均大于0.05 (表3), 因此模型需要修订。修订时保留通径系数显著的, 删除不显著的。结构方程模型作为一种先验模型, 允许根据作物的生长实际, 保留一些统计上不显著但有一定生物学意义的路径。最终在初始模型的基础上, 我们增加了单位面积穗数到株高的路径, 去除了单位面积穗数到穗粒数的路径。修订模型参数见表3, P值为0.199, RMR和RMSEA均小于0.05。
表3
Table 3
表3(Table 3)
表3 初始与最终结构方程模型的主要拟合指标 Table 3 Initial and final fit indices of structural equation model
| 表3 初始与最终结构方程模型的主要拟合指标 Table 3 Initial and final fit indices of structural equation model |
在最终建立的结构方程模型中, 各农艺性状对产量的直接效应均为正效应, 其中单位面积穗数、穗粒数、千粒重对产量的影响较大, 通径系数分别为0.60、0.56和0.44; 基本苗对产量的通径系数为0.30; 而株高和生育期则对产量的影响较小。除了直接效应以外, 各性状还通过彼此的间接效应对最终产量发挥作用(图2)。如单位面积穗数和穗粒数对千粒重的影响均为负向, 标准化通径系数分别为-0.45和-0.18; 生育期则对株高和基本苗有负效应, 标准化通径系数分别为-0.56和-0.50; 单位面积穗数对穗粒数的影响则很小。因此, 各性状对产量的影响应按结构方程模型所示, 需要综合直接效应与间接效应进行评价, 讨论总效应。单位面积穗数对产量的直接效应为0.60, 通过千粒重、穗粒数和株高的间接效应分别为-0.20、-0.02和0.01, 总效应为0.39(表4)。千粒重对产量的最终效应系数是0.44。穗粒数对产量的直接效应为0.56, 但穗粒数与千粒重呈负相关, 即穗粒数越多, 籽粒越小, 这与生产实践中的结果是一致的。由于穗分化在先、灌浆在后, 因此穗粒数通过千粒重对产量产生间接影响(-0.08), 穗粒数对产量最终的效应为0.48。株高对产量的直接效应系数为0.13 (图2), 但株高分别通过千粒重(-0.03)和穗粒数(-0.06)对产量有负的影响, 因此, 其对产量的最终效应系数仅为0.04。基本苗的直接效应系数为0.30, 通过单位面积穗数、千粒重和株高效应系数分别为0.12、0.08和-0.01, 最终的效应系数为0.49。生育期对产量的直接效应为0.16, 通过株高、穗粒数、单位面积穗数和基本苗的效应系数分别为-0.02、-0.15、0.11和-0.25, 最终作用为-0.15。从决定系数看, 各性状随机误差的决定系数均高(0.74~0.99), 而产量的决定系数达到0.72, 与回归分析结果一致, 产量随机误差决定系数为0.28 (图2)。这表明生产中冬小麦所处环境的随机变化对各农艺性状的影响较大, 而这些性状特征一旦确定, 最终的产量也就基本确定了。
表4
Table 4
表4(Table 4)
表4 冬小麦主要农艺性状对产量的直接及间接效应 Table 4 Direct and indirect effects of agronomic traits on yield of winter wheat
| 表4 冬小麦主要农艺性状对产量的直接及间接效应 Table 4 Direct and indirect effects of agronomic traits on yield of winter wheat | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
图2
Fig. 2
| Figure OptionViewDownloadNew Window | |
 | 图2 结构方程模型中冬小麦主要性状对产量的影响及各性状间的关系 实线表示显著的路径, 实线的粗细表示通径系数绝对值的大小, 虚线表示该路径不显著, 但从冬小麦生理过程讲需要保留的路径。线旁的数值为标准化通径系数, 矩形框及椭圆框旁的数值为决定系数。GY: 产量; GNP: 穗粒数; BS: 基本苗; SN: 单位面积穗数; GD: 生育期; TGW: 千粒重; PH: 株高。Fig. 2 Effects of agronomic traits on yield of winter wheat and relationships among agronomic traits as results of structural equation modeling The solid lines show the significant paths, whose thick and thin show the great or small of absolute path coefficient. Dashed lines show the path is not significant, but these paths need to be conserved based on winter wheat’ s physiological process. The values beside the lines are standard path coefficients, and the values beside the rectangular or oval are determination coefficients. GY: grain yield; GNP: grain number per spike; BS: density of basic seedlings; SN: spike number per ha; GD: growth duration; TGW: thousand-grain weight; PH: plant height. |
3 讨论本研究所涉及的冬小麦产量基本数据(表1)变异系数最大的为产量(0.246), 其次为穗粒数(0.208)和单位面积穗数(0.199), 最小的是生育期(0.086), 表明产量因品种和栽培地不同存在较大差异, 这种差异更多的可能是由穗粒数和单位面积穗数的变异造成的, 而生育期相对变化小, 对产量的影响可能也较小。这与相关性分析结果(表2)所得判断有相似之处, 即穗粒数和单位面积穗数影响较大, 相关系数分别0.48和0.45, 但相关分析结果显示基本苗也对产量影响较大, 相关系数为0.48, 而株高对产量影响最小。
通过结构方程模型分析的结果表明, 各性状对产量总的效应与相关分析和回归分析的结果是不同的。结构方程模型中各性状的直接通径系数与多元回归方程中的标准化系数相同, 评估后的总效应系数则发生了变化, 其中单位面积穗数、穗粒数、株高和生育期的效应下降, 尤其是生育期, 其效应由正转负, 系数由0.16转变-0.16。而基本苗的总效应则由原来0.30增加到0.50。总的效应系数与性状产量相关系数判断的正负方向是一致的, 但具体数值不同(表4)。效应排序为基本苗> 穗粒数> 千粒重> 单位面积穗数> 株高> 生育期。
基本苗对产量的直接效应低于单位面积穗数、穗粒数和千粒重, 但加上通过单位面积穗数和千粒重对产量的间接效应后, 总的效应上升为第1位。因此推论, 本研究中所涉及的地区、品种在栽培时播种量有进一步提高空间。穗粒数和千粒重同为穗部特征, 在假设光合产物供应一定的情况, 穗粒数的增加可能会导致千粒重下降(表2)。从发育过程看, 穗粒数的分化在先, 籽粒灌浆在后, 所以在结构方程模型分析过程中允许考虑穗粒数对千粒重的影响, 不需要考虑千粒重对穗粒数的影响, 这也是先验模型的优势所在。理论上千粒重主要决定于灌浆程度, 而本研究中缺少评价灌浆的相关指标, 因此其直接效应即为最终效应。单位面积穗数的效应因通过千粒重对产量的间接负效应而削弱, 最终低于千粒重的直接效应。单位面积穗数分别通过千粒重和穗粒数对产量的间接作用均为负值, 尤其是通过千粒重的负向影响较大, 与千粒重对产量的直接效应相同。生育期对产量的直接效应是正向的, 从生理中可解释为长的生育期有利于延长总光合时间, 进而增加光合产物, 但生育期通过穗粒数和基本苗对产量有负效应, 最终导致生育期对产量也是负向影响。无论是哪一种分析方法, 株高对产量的影响都是最小的, 表明这些育成品种株高可能已经接近最适状态。
4 结论通过结构方程模型分析, 长江上游冬小麦区域试验19个冬小麦品种在栽培中播种量有进一步提高空间, 穗部特征对产量的影响高于单位面积穗数、株高和生育期, 在选育工作中应重点关注穗部性状。结构方程模型作为一种先验模型可应用于作物研究, 分析多变量关系, 尤其是多响应变量, 探讨其内部作用过程与机制。
The authors have declared that no competing interests exist.
参考文献View Option
原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
| [1] | |
| [2] | |
| [3] | |
| [4] | |
| [5] | |
| [6] | |
| [7] | |
| [8] | |
| [9] | |
| [10] | |
| [11] | |
| [12] | |
| [13] | |
| [14] | |
| [15] | |
| [16] | |
| [17] |
