r i = {{i}}
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") }
model$ANOVA model$statistics
model$means
if (ng >= 50) { model$groups[, -3] } else { model$groups }
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.