mtsi: Multi-trait stability index

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/mtsi.R

Description

[Stable]

Computes the multi-trait stability index proposed by Olivoto et al. (2019)

Usage

1
2
3
4
5
6
7
8
mtsi(
  .data,
  index = "waasby",
  ideotype = NULL,
  SI = 15,
  mineval = 1,
  verbose = TRUE
)

Arguments

.data

An object of class waasb or waas.

index

If index = 'waasby' (default) both stability and mean performance are considered. If index = 'waasb' the multi-trait index will be computed considering the stability of genotypes only. More details can be seen in waasb() and waas() functions.

ideotype

A vector of length nvar where nvar is the number of variables used to plan the ideotype. Use 'h' to indicate the traits in which higher values are desired or 'l' to indicate the variables in which lower values are desired. For example, ideotype = c("h, h, h, h, l") will consider that the ideotype has higher values for the first four traits and lower values for the last trait.

SI

An integer (0-100). The selection intensity in percentage of the total number of genotypes.

mineval

The minimum value so that an eigenvector is retained in the factor analysis.

verbose

If verbose = TRUE (Default) then some results are shown in the console.

Value

An object of class mtsi with the following items:

Author(s)

Tiago Olivoto tiagoolivoto@gmail.com

References

Olivoto, T., A.D.C. L\'ucio, J.A.G. da silva, B.G. Sari, and M.I. Diel. 2019. Mean performance and stability in multi-environment trials II: Selection based on multiple traits. Agron. J. 111:2961-2969. doi: 10.2134/agronj2019.03.0220

See Also

mgidi(), waasb(), get_model_data()

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
library(metan)
# Based on stability only, for both GY and HM, higher is better
mtsi_model <-
    waasb(data_ge,
    env = ENV,
    gen = GEN,
    rep = REP,
    resp = c(GY, HM))
mtsi_index <-
    mtsi(mtsi_model, index = 'waasb')


# Based on mean performance and stability (using pipe operator %>%)
# GY: higher is better
# HM: lower is better

mtsi_index2 <-
 data_ge %>%
 waasb(ENV, GEN, REP,
       resp = c(GY, HM),
       mresp = c("h, l")) %>%
 mtsi()

metan documentation built on Nov. 10, 2021, 9:11 a.m.