ge_stats()
computes parametric and non-parametric stability statistics
given a data set with environment, genotype, and block factors.
Arguments
- .data
The dataset containing the columns related to Environments, Genotypes, replication/block and response variable(s).
- env
The name of the column that contains the levels of the environments.
- gen
The name of the column that contains the levels of the genotypes.
- rep
The name of the column that contains the levels of the replications/blocks.
- resp
The response variable(s). To analyze multiple variables in a single procedure use, for example,
resp = c(var1, var2, var3)
.- verbose
Logical argument. If
verbose = FALSE
the code will run silently.- prob
The probability error assumed.
Value
An object of class ge_stats
which is a list with one data
frame for each variable containing the computed indexes.
Details
The function computes the statistics and ranks for the following stability indexes.
"Y"
(Response variable),"CV"
(coefficient of variation)"ACV"
(adjusted coefficient of variation callingge_acv()
internally)POLAR
(Power Law Residuals, callingge_polar()
internally)"Var"
(Genotype's variance)"Shukla"
(Shukla's variance, callingShukla()
internally)"Wi_g", "Wi_f", "Wi_u"
(Annichiarrico's genotypic confidence index for all, favorable and unfavorable environments, respectively, callingAnnicchiarico()
internally )"Ecoval"
(Wricke's ecovalence,ecovalence()
internally)"Sij"
(Deviations from the joint-regression analysis) and"R2"
(R-squared from the joint-regression analysis, callingge_reg()
internally)"ASTAB"
(AMMI Based Stability Parameter),"ASI"
(AMMI Stability Index),"ASV"
(AMMI-stability value),"AVAMGE"
(Sum Across Environments of Absolute Value of GEI Modelled by AMMI ),"Da"
(Annicchiarico's D Parameter values),"Dz"
(Zhang's D Parameter),"EV"
(Sums of the Averages of the Squared Eigenvector Values),"FA"
(Stability Measure Based on Fitted AMMI Model),"MASV"
(Modified AMMI Stability Value),"SIPC"
(Sums of the Absolute Value of the IPC Scores),"Za"
(Absolute Value of the Relative Contribution of IPCs to the Interaction),"WAAS"
(Weighted average of absolute scores), callingammi_indexes()
internally"HMGV"
(Harmonic mean of the genotypic value),"RPGV"
(Relative performance of the genotypic values),"HMRPGV"
(Harmonic mean of the relative performance of the genotypic values), by callingblup_indexes()
internally"Pi_a", "Pi_f", "Pi_u"
(Superiority indexes for all, favorable and unfavorable environments, respectively, callingsuperiority()
internally)"Gai"
(Geometric adaptability index, callinggai()
internally)"S1"
(mean of the absolute rank differences of a genotype over the n environments),"S2"
(variance among the ranks over the k environments),"S3"
(sum of the absolute deviations),"S6"
(relative sum of squares of rank for each genotype), by callingHuehn()
internally"N1"
,"N2"
,"N3"
,"N4"
(Thennarasu"s statistics, callingThennarasu()
internally ).
References
Annicchiarico, P. 1992. Cultivar adaptation and recommendation from alfalfa trials in Northern Italy. Journal of Genetic & Breeding, 46:269-278
Ajay BC, Aravind J, Abdul Fiyaz R, Bera SK, Kumar N, Gangadhar K, Kona P (2018). “Modified AMMI Stability Index (MASI) for stability analysis.” ICAR-DGR Newsletter, 18, 4–5.
Ajay BC, Aravind J, Fiyaz RA, Kumar N, Lal C, Gangadhar K, Kona P, Dagla MC, Bera SK (2019). “Rectification of modified AMMI stability value (MASV).” Indian Journal of Genetics and Plant Breeding (The), 79, 726–731. https://www.isgpb.org/article/rectification-of-modified-ammi-stability-value-masv.
Annicchiarico P (1997). “Joint regression vs AMMI analysis of genotype-environment interactions for cereals in Italy.” Euphytica, 94(1), 53–62. doi:10.1023/A:1002954824178
Doring, T.F., and M. Reckling. 2018. Detecting global trends of cereal yield stability by adjusting the coefficient of variation. Eur. J. Agron. 99: 30-36. doi:10.1016/j.eja.2018.06.007
Doring, T.F., S. Knapp, and J.E. Cohen. 2015. Taylor's power law and the stability of crop yields. F. Crop. Res. 183: 294-302. doi:10.1016/j.fcr.2015.08.005
Eberhart, S.A., and W.A. Russell. 1966. Stability parameters for comparing Varieties. Crop Sci. 6:36-40. doi:10.2135/cropsci1966.0011183X000600010011x
Farshadfar E (2008) Incorporation of AMMI stability value and grain yield in a single non-parametric index (GSI) in bread wheat. Pakistan J Biol Sci 11:1791–1796. doi:10.3923/pjbs.2008.1791.1796
Fox, P.N., B. Skovmand, B.K. Thompson, H.J. Braun, and R. Cormier. 1990. Yield and adaptation of hexaploid spring triticale. Euphytica 47:57-64. doi:10.1007/BF00040364 .
Huehn, V.M. 1979. Beitrage zur erfassung der phanotypischen stabilitat. EDV Med. Biol. 10:112.
Jambhulkar NN, Rath NC, Bose LK, Subudhi HN, Biswajit M, Lipi D, Meher J (2017). “Stability analysis for grain yield in rice in demonstrations conducted during rabi season in India.” Oryza, 54(2), 236–240. doi:10.5958/2249-5266.2017.00030.3
Kang, M.S., and H.N. Pham. 1991. Simultaneous Selection for High Yielding and Stable Crop Genotypes. Agron. J. 83:161. doi:10.2134/agronj1991.00021962008300010037x
Lin, C.S., and M.R. Binns. 1988. A superiority measure of cultivar performance for cultivar x location data. Can. J. Plant Sci. 68:193-198. doi:10.4141/cjps88-018
Olivoto, T., A.D.C. L\'ucio, J.A.G. da silva, V.S. Marchioro, V.Q. de Souza, and E. Jost. 2019a. Mean performance and stability in multi-environment trials I: Combining features of AMMI and BLUP techniques. Agron. J. 111:2949-2960. doi:10.2134/agronj2019.03.0220
Mohammadi, R., & Amri, A. (2008). Comparison of parametric and non-parametric methods for selecting stable and adapted durum wheat genotypes in variable environments. Euphytica, 159(3), 419-432. doi:10.1007/s10681-007-9600-6
Shukla, G.K. 1972. Some statistical aspects of partitioning genotype-environmental components of variability. Heredity. 29:238-245. doi:10.1038/hdy.1972.87
Raju BMK (2002). “A study on AMMI model and its biplots.” Journal of the Indian Society of Agricultural Statistics, 55(3), 297–322.
Rao AR, Prabhakaran VT (2005). “Use of AMMI in simultaneous selection of genotypes for yield and stability.” Journal of the Indian Society of Agricultural Statistics, 59, 76–82.
Sneller CH, Kilgore-Norquest L, Dombek D (1997). “Repeatability of yield stability statistics in soybean.” Crop Science, 37(2), 383–390. doi:10.2135/cropsci1997.0011183X003700020013x
Thennarasu, K. 1995. On certain nonparametric procedures for studying genotype x environment interactions and yield stability. Ph.D. thesis. P.J. School, IARI, New Delhi, India.
Wricke, G. 1965. Zur berechnung der okovalenz bei sommerweizen und hafer. Z. Pflanzenzuchtg 52:127-138.
See also
acv()
, ammi_indexes()
, ecovalence()
, Fox()
, gai()
,
ge_reg()
, hmgv()
, hmrpgv()
, rpgv()
, Huehn()
, polar()
,
Shukla()
, superiority()
, Thennarasu()
, waas()
, waasb()
Author
Tiago Olivoto tiagoolivoto@gmail.com
Examples
# \donttest{
library(metan)
model <- ge_stats(data_ge, ENV, GEN, REP, GY)
#> Evaluating trait GY |============================================| 100% 00:00:21
get_model_data(model, "stats")
#> Class of the model: ge_stats
#> Variable extracted: stats
#> # A tibble: 10 × 44
#> var GEN Y CV ACV POLAR Var Shukla Wi_g Wi_f Wi_u Ecoval
#> <chr> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1 GY G1 2.60 35.2 34.1 0.0298 10.9 0.0280 84.4 89.2 81.1 1.22
#> 2 GY G10 2.47 42.3 38.6 0.136 14.2 0.244 59.2 64.6 54.4 7.96
#> 3 GY G2 2.74 34.0 35.2 0.0570 11.3 0.0861 82.8 95.3 75.6 3.03
#> 4 GY G3 2.96 29.9 33.8 0.0216 10.1 0.0121 104. 99.7 107. 0.725
#> 5 GY G4 2.64 31.4 31.0 -0.0537 8.93 0.0640 85.9 79.5 91.9 2.34
#> 6 GY G5 2.54 30.6 28.8 -0.119 7.82 0.0480 82.7 82.2 82.4 1.84
#> 7 GY G6 2.53 29.7 27.8 -0.147 7.34 0.0468 83.0 83.7 81.8 1.81
#> 8 GY G7 2.74 27.4 28.3 -0.133 7.33 0.122 83.9 77.6 93.4 4.16
#> 9 GY G8 3.00 30.4 35.1 0.0531 10.8 0.0712 98.8 90.5 107. 2.57
#> 10 GY G9 2.51 42.4 39.4 0.154 14.7 0.167 68.8 68.9 70.3 5.56
#> # … with 32 more variables: bij <dbl>, Sij <dbl>, R2 <dbl>, ASTAB <dbl>,
#> # ASI <dbl>, ASV <dbl>, AVAMGE <dbl>, DA <dbl>, DZ <dbl>, EV <dbl>, FA <dbl>,
#> # MASI <dbl>, MASV <dbl>, SIPC <dbl>, ZA <dbl>, WAAS <dbl>, WAASB <dbl>,
#> # HMGV <dbl>, RPGV <dbl>, HMRPGV <dbl>, Pi_a <dbl>, Pi_f <dbl>, Pi_u <dbl>,
#> # Gai <dbl>, S1 <dbl>, S2 <dbl>, S3 <dbl>, S6 <dbl>, N1 <dbl>, N2 <dbl>,
#> # N3 <dbl>, N4 <dbl>
# }