r i = {{i}}

{{i+1}}. Analysis for trait r traits[i]

y <- dfr[, traits[i]]

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

{{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
}


AGROFIMS/hagrofims documentation built on May 6, 2020, 7:43 p.m.