ANOVA in R using the function lapply ans terms
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, f=da$local), aov, formula=producao ~ repeticao/bloco + cultivar) lapply(mi, summary)#not adjusted for blocks But to transform cultivar adjusted for blocks (to obtain adjusted means) I need use command terms(formula,keep.order=TRUE) I DO NOT KNOW HOW TO USE THIS IN THE CORRECT ORDER. Something like this? m0 <- lapply(split(da, f=da$local), aov(terms(producao ~ repeticao/bloco + cultivar,keep.order = TRUE),da)) lapply(m0, summary)#adjusted to blocks Thanks Euriann
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 conjunta. m0 <- aov(terms(prod ~ loc/rpt/blc + trt, keep.order = TRUE), data = da) anova(m0) # Médias amostrais. aggregate(prod ~ trt, data = da, FUN = mean) # Médias ajustadas. LSmeans(m0, effect = "trt") # Análises por local. fits <- lapply(split(da, f = da$loc), FUN = aov, formula = terms(prod ~ rpt/blc + trt, keep.order = TRUE)) lapply(fits, anova) # Médias amostrais por local. aggregate(prod ~ trt + loc, data = da, FUN = mean) # Médias ajustadas. lapply(fits, FUN = LSmeans, effect = "trt") # Para o delineamento usado nesse exemplo, as médias ajustadas coincidem # com as médias amostrais. À disposição. Walmes.
participantes (2)
-
Euriann Yamamoto -
Walmes Zeviani