r i = {{i}}
r traits[i]
y <- data[, traits[i]] n_geno <- length(unique(data[, geno])) data[, block] <- factor(data[, block]) data[, geno] <- factor(data[, geno]) data[, rep] <- factor(data[, rep]) if(n_geno>=50){ model <- pepa::PBIB.test.mod(data[, block], data[, geno], data[, rep], y, k, method = "VC") } else { model <- agricolae::PBIB.test(data[, block], data[, geno], data[, rep], y, k, method = "VC") }
model$ANOVA model$statistics
model$means
if(n_geno>= 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.