# Example dataset
data("data", package = "seliNDRIx", envir = environment())
traits <- c("tmy", "py", "fatyield")
fixed <- c("farm", "soc", "poc")
random <- c("sire")
economic_values <- c(1, 0.85, 0.65)
results <- mixed_si(data = data, traits = traits,
fixed = fixed, random = random, economic_values = economic_values)
Run the code above in your browser using DataLab