Como definir um modelo para esse experimento?
Boa tarde! Pessoal, tenho um experimento em DIC, com 6 repetições, em esquema de parcelas subdivididas, onde tenho nas parcelas um fatorial (2 variedades X 2 vinhaças) e nas subparcelas aplicou-se aluminio (4 níveis: 60, 45, 35 e 20%). Minha dúvida: como defino no R um modelo inicial para analisar isso? Seria algo talvez parecido com isso: aov(y ~ variedade*vinhaça*aluminio + Error(parcelas/subparcelas)) ? Desde já, muito obrigadaSimone D. Sartorio
Se eu entendi correto, minhas sugestões são: da <- expand.grid(var=gl(2,1), vin=gl(2,1), alu=gl(4,1), r=1:6) da$ue <- with(da, interaction(var, vin, r)) da$y <- rnorm(nrow(da)) str(da) require(nlme) m0 <- aov(y~var*vin*alu+Error(ue), data=da) summary(m0) m1 <- lme(y~var*vin*alu, random=~1|ue, data=da) anova(m1) No entanto, eu reforço minha preferência pelo modelo misto e alguns pontos são: * no caso balanceado são equivalentes; * no caso desbalanceado, o modelo misto não requer balanceamento portanto continua viável; * existem mais métodos disponíveis para objetos de classe "lme" do que para "aovlist", como por exemplo, comparações múltiplas de hipóteses. À disposição. Walmes.
participantes (2)
-
Simone D. Sartorio -
walmes .