基于python语言的冬小麦品种稳定性分析

doi:10.11924/j.issn.1000-6850.casb19010094

摘要/Abstract

摘要：

旨在挑选出适宜的冬小麦品种稳定性分析参数及其计算方法。利用冬小麦实验数据,首次通过python语言计算了5个稳定性参数,并采用秩相关性和主成分分析,研究了参数的相关性。结果表明,python语言可以方便、灵活的计算各类参数,计算相关性和进行主成分分析。相关分析表明,PCOA和产量显著正相关,CV和Bi和产量极显著负相关。研究表明,PCOA模型能够同时评价品种产量及其稳定性。多元统计方法是单变量参数模型的有益补充。

关键词: 冬小麦, 品种, 稳定性, python语言, 程序, 秩相关性

Abstract:

The study aims at selecting suitable stability analysis parameters and calculation methods for winter wheat varieties. Using the experimental data of winter wheat, five stability parameters were calculated by Python for the first time. The correlations of parameters were studied by rank correlation and principal component analysis. The results showed that the Python language could easily and flexibly calculate various parameters and rank correlation and analyze principal components. The correlation analysis showed that PCOA was significantly and positively correlated with yield; CV and Bi were significantly and negatively correlated with yield. The results indicate that the PCOA model can simultaneously evaluate yield and stability of varieties, and the multivariate method can be a useful complement to the univariate parameter model.

Key words: winter wheat, variety, stability, python language, program, rank correlation

中图分类号:

S11+4

李晓航, 王映红. 基于python语言的冬小麦品种稳定性分析[J]. 中国农学通报, 2020, 36(14): 156-164.

Li Xiaohang, Wang Yinghong. Stability Analysis of Winter Wheat Varieties Based on Python Language[J]. Chinese Agricultural Science Bulletin, 2020, 36(14): 156-164.

图/表 10

参考文献 27

[1]	Mohammadi R, Amri A . Comparison of parametric and non-parametric methods for selecting stable and adapted durum wheat genotypes in variable environments[J]. Euphytica, 2008,159(3):419-432. doi: 10.1007/s10681-007-9600-6 URL
[2]	Liu Y J, Chen Q M . Modelling the impacts of climate change and crop management on phenological trends of spring and winter wheat in China[J]. Agricultural and Forest Meteorology, 2018,246(2):518-526.
[3]	Farshadfar E . Incorporation of AMMI Stability Value and Grain Yield in a Single Non-Parametric Index (GSI) in Bread Wheat[J]. Pakistan Journal of Biological Sciences, 2008(8):1790-1796.
[4]	Farshadfar E, Mohammadi R, Aghaee M , et al. GGE biplot analysis of genotype x environment interaction in wheat-barley disomic addition lines[J]. Australian Journal of Crop Science, 2012,6(6):1074-1079.
[5]	Lempereur I, RouauJ X, Abecassis J . Genetic and Agronomic Variation in Arabinoxylan and Ferulic Acid Contents of Durum Wheat (Triticum durum L.) Grain and Its Milling Fractions[J]. Journal of Cereal Science, 1997,25(2):103-110.
[6]	Akçura M, Kaya Y, Taner S . Evaluation of durum wheat genotypes using parametric and non parametric stability statistics[J]. Turkish J. of Field Crops, 2009,14(2):111-122.
[7]	Mohammadi R, Amri A . Comparison of parametric and non-parametric methods for selecting stable and adapted durum wheat genotypes in variable environments[J]. Euphytica, 2008,159:419-432.
[8]	Lucila O, Susanna-Assunta S . Finding useful data across multiple biomedical data repositories using DataMed[J]. Nature genetics , 2017,49(6), 816-820.
[9]	Leonardo M, Vanessa B, Juliana F , et al. noWorkflow: capturing and analyzing provenance of scripts. In International Provenance and Annotation Workshop[J]. Springer, 2007(14):71-83.
[10]	Jake V . python数据科学手册[M]. 北京: 人民邮电出版, 2018(1):3-4.
[11]	金建猛, 赵鹏飞, 赵国建 . 开麦20农艺特性分析及小麦育种方向思考[J]. 农业与技术. , 2018,2018(03).
[12]	张嵩午, 刘党校 . 冷型小麦品质稳定性的研究[J]. 自然科学进展, 2007,17(1):29-34.
[13]	Shukla G K . Some Statistical Aspects of Partitioning Genotype-Environmental Components of Variability[J]. Heredity, 1972(29):237-245.
[14]	Francis T R, Kannenberg L W . Yield Stability Studies in Short Season Maize: I. Descriptive Method for Grouping Genotypes[J]. Plant Sci, 1978(58):1029-1034.
[15]	Hasan K . Assessment of Parametric and Non-parametric Methods for Selecting Stable and Adapted Durum Wheat Genotypes in Multi-Environments[J]. Journal of Animal & Plant Sciences, 2012.
[16]	Akçura M, Kaya Y . Non-parametric stability methods for interpreting genotype by environment interaction of bread wheat genotypes (Triticum aestivum L.) Genetics and Molecular Biology, 2008,31(4):906-913.
[17]	吴元奇, 潘光堂 . 作物稳定性研究进展[J]. 四川农业大学学报, 2005,12(4):482-489.
[18]	Alagarswamy G, Chandra S . Pattern analysis of international sorghum multi-environment trials for grain-yield adaptation[J]. Theoretical and Applied Genetics, 1998(96):397-405.
[19]	朱文辉 . 离散隐马尔科夫模型及其齿轮箱动态故障诊断[D]. 太原:中北大学, 2018: 1-73.
[20]	Chandra S, Sohoo M S, Singh K P . Genotype environment interaction for yield in ram[J]. Journal of Research, 1974(8):165-168.
[21]	Roy D . Plant Breeding, Analysis and Exploitation of Variation[J]. Alpha Science International Ltd. Pang Bourne, India, 2000(11):163-195.
[22]	Mohammadi R, Amri A . Analysis of Genotype × Environment interaction for grain yield in durum wheat[J]. Crop Science, 2008(49):1177-1186.
[23]	Eberhart S A, Russell W A . Stability Parameters for Comparing Varieties[J]. Crop Science, 1966(6):36-40.
[24]	Fernando P, Brian E . Project Jupyter: Computational narratives as the engine of collaborative data science[EB/OL]. http://archive.ipython.org/JupyterGrantNarrative-2015.pdf.
[25]	Oinn T, Greenwood M, Goble A , et al. lessons in creating a workflow environmentfor the life sciences[J]. 2006(10):1067-1100.
[26]	Pinthus M J . Estimates of Genotypic Value: A Proposed Method[J]. Euphytica, 1973(22):345-51.
[27]	Flores F, Moreno M T, Cubero J I . A comparison of univariate and multivariate methods to analyze G × E interaction[J]. Field Crops Research, 1998,56(3):271-286.

品种	地点
品种	D1	D2	D3	D4	D5	D6	D7	D8	D9	D10	D11	D12	D13	D14
s156	580.5	502.7	366.3	537.9	380.3	401.4	444.8	405.3	411.7	372.6	512.3	458.5	420.2	464.7
l146	5416.6	539.3	394.2	537.9	424.7	352.2	439.1	436.6	405.2	386.6	532.4	503.9	393.5	482
w182	621.6	523.8	443.2	526.4	458.8	390.5	435.7	457.3	417.8	367.1	529.3	423.6	431.2	498.7
gy13028	531.1	544.2	406.7	491.7	416.3	414.9	449.8	457	401.5	342.1	467.6	438.8	408.1	443.3
z875	606.2	557.4	450.7	552.6	421.4	346.4	473.5	392.4	401.7	398.8	523.1	357.7	424.5	393.1
w803	586	548.2	335	529.4	365.5	389.5	482.5	426.4	414.5	358.3	574.1	398.4	425.9	457.9
l899	582.6	568.7	413	551.1	422.4	370.3	395.3	445.1	423.7	361.5	513.9	418.8	435.5	467
a1303	554.4	540.7	398.7	527.5	380.8	419.9	445.2	453.4	392.5	369.9	472.2	429.2	442.6	496.3
z162	517.1	473.9	417.2	500	380.8	389.2	420.8	391.1	440.2	365.8	501.5	427．1	416.4	447.9
z18ck	520.6	484.4	398.2	513	410.9	367.6	404.2	399.4	414.1	372.5	470.7	403.8	367.4	441.3

品种	地点
品种	D1	D2	D3	D4	D5	D6	D7	D8	D9	D10	D11	D12	D13	D14
s156	580.5	502.7	366.3	537.9	380.3	401.4	444.8	405.3	411.7	372.6	512.3	458.5	420.2	464.7
l146	5416.6	539.3	394.2	537.9	424.7	352.2	439.1	436.6	405.2	386.6	532.4	503.9	393.5	482
w182	621.6	523.8	443.2	526.4	458.8	390.5	435.7	457.3	417.8	367.1	529.3	423.6	431.2	498.7
gy13028	531.1	544.2	406.7	491.7	416.3	414.9	449.8	457	401.5	342.1	467.6	438.8	408.1	443.3
z875	606.2	557.4	450.7	552.6	421.4	346.4	473.5	392.4	401.7	398.8	523.1	357.7	424.5	393.1
w803	586	548.2	335	529.4	365.5	389.5	482.5	426.4	414.5	358.3	574.1	398.4	425.9	457.9
l899	582.6	568.7	413	551.1	422.4	370.3	395.3	445.1	423.7	361.5	513.9	418.8	435.5	467
a1303	554.4	540.7	398.7	527.5	380.8	419.9	445.2	453.4	392.5	369.9	472.2	429.2	442.6	496.3
z162	517.1	473.9	417.2	500	380.8	389.2	420.8	391.1	440.2	365.8	501.5	427．1	416.4	447.9
z18ck	520.6	484.4	398.2	513	410.9	367.6	404.2	399.4	414.1	372.5	470.7	403.8	367.4	441.3

品种	YIELD	CV	Bi	pi
s156	446.4	0.1	1	1200
l146	454.6	0.1	1	876.2
w182	466.1	0.1	1	568
gy13028	443.8	0.1	1	1496
z875	450.0	0.2	1	1728.4
w803	449.4	0.2	1	1484.4
1899	454.9	0.2	1	1001.6
a1303	451.7	0.1	1	1257.3
z162	434.9	0.1	1	1947.8
z18ck	426.3	0.1	1	2399.3

品种	YIELD	CV	Bi	pi
s156	446.4	0.1	1	1200
l146	454.6	0.1	1	876.2
w182	466.1	0.1	1	568
gy13028	443.8	0.1	1	1496
z875	450.0	0.2	1	1728.4
w803	449.4	0.2	1	1484.4
1899	454.9	0.2	1	1001.6
a1303	451.7	0.1	1	1257.3
z162	434.9	0.1	1	1947.8
z18ck	426.3	0.1	1	2399.3

品种	YIELD	CV	Bi	pi	PCOA
s156	4	7	6	4	4
1146	5	6	5	2	1
w186	6	5	7	1	6
gy13028	9	2	2	7	9
z875	2	9	9	8	3
w803	1	10	10	6	5
1899	3	8	8	3	2
a1303	7	4	4	5	7
z162	10	1	1	9	10
zl18ck	8	3	3	10	8