# NOT RUN {
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(100, 0)) %>%
mtsi()
# }
Run the code above in your browser using DataLab