if (FALSE) {
library(agridat)
library(agriutilities)
data(besag.met)
dat <- besag.met
results <- check_design_met(
data = dat,
genotype = "gen",
trial = "county",
traits = c("yield"),
rep = "rep",
block = "block",
col = "col",
row = "row"
)
out <- single_trial_analysis(results, progress = FALSE)
met_results <- met_analysis(out, vcov = "fa2", progress = FALSE)
pp <- met_results$trial_effects
model <- met_results$met_models$yield
fa_summary(
model = model,
trial = "trial",
genotype = "genotype",
BLUEs_trial = pp,
mult_fa1 = -1,
mult_fa2 = -1,
filter_score = 1,
k_biplot = 10,
size_label_var = 3,
alpha_label_var = 0.5,
size_label_ind = 3,
alpha_label_ind = 0.8,
size_arrow = 0.2,
alpha_arrow = 0.1
)
}
Run the code above in your browser using DataLab