[R-br] ANOVA in R using the function lapply ans terms

[Euriann Yamamoto]

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




-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20171024/1271f1d5/attachment.html>