ProbBreed employs Bayesian statistics to analyse multi-environment trials’ data, and uses its outputs to estimate the marginal and pairwise probabilities of superior performance and superior stability of the genotypes, as well as their conditional probability of superior performance. The method is thoroughly described at https://doi.org/10.1007/s00122-022-04041-y.
You can install the CRAN version of ProbBreed
using the following command:
Alternatively, you can install the development version of ProbBreed
from GitHub with:
A basic workflow using the available data is:
library(ProbBreed)
mod = bayes_met(data = maize,
gen = "Hybrid",
loc = "Location",
repl = c("Rep", "Block"),
year = NULL,
reg = "Region",
res.het = F,
trait = "GY",
iter = 6000, cores = 4, chains = 4)
outs = extr_outs(data = maize, trait = "GY", model = mod,
probs = c(0.05, 0.95),
check.stan.diag = FALSE,
verbose = TRUE)
results = prob_sup(data = maize, trait = "GY", gen = "Hybrid", loc = "Location",
mod.output = outs, reg = 'Region', year = NULL, int = .2,
increase = TRUE, save.df = FALSE, interactive = FALSE,
verbose = TRUE)