# NOT RUN {
library(metan)
# Considering p-value <= 0.05 to compute the WAAS
model <- waas(data_ge,
env = ENV,
gen = GEN,
rep = REP,
resp = GY)
# Declaring the number of axis to be used for computing WAAS
# and assigning a larger weight for the response variable when
# computing the WAASBY index.
model2 <- waas(data_ge,
env = ENV,
gen = GEN,
rep = REP,
resp = GY,
naxis = 3,
wresp = 60)
# Analyzing multiple variables (GY and HM) at the same time
# considering that smaller values of HM are better and higher
# values of GY are better, assigning a larger weight for the GY
# and a smaller weight for HM when computing WAASBY index.
model3 <- waas(data_ge,
env = ENV,
gen = GEN,
rep = REP,
resp = c(GY, HM),
mresp = c(100, 0),
wresp = c(60, 40))
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab