Euriann,
Considere o código abaixo.
library(doBy)
# Conjunto de ensaios em DBC em 6 locais (gerado por simulação).
da <- expand.grid(loc = gl(6, 1),
rpt = gl(3, 1),
blc = gl(4, 1),
trt = gl(5, 1))
da$prod <- rnorm(nrow(da))
# Análise
Hi everyone.
In my dataset I´ve got 36 genotypes in 5 locations in partially balanced
incomplete design 6x6, simple. We are analyzing productivity. I need to use
adjusted means to do the orders analysis.
TO obtain ANOVA individual to the 5 locations I run command:
mi <- lapply(split(da,
