r i = {{i}}

{{i+1}}. Analysis for variable r vars[i]

y <- dfr[, vars[i]]

if (ng >= 50) {
  model <- pepa::PBIB.test.mod(block, geno, rep, y, k, method)
} else {
  model <- agricolae::PBIB.test(block, geno, rep, y, k, method)
}

{{i+1}}.1. ANOVA

model$ANOVA
model$statistics

{{i+1}}.2. Adjusted means

model$means

{{i+1}}.3. LSD test

if (ng >= 50) {
  model$groups[, -3]
} else {
  model$groups
}


reyzaguirre/pepa documentation built on March 29, 2025, 9:56 p.m.