[R-br] Analise de variância fator com níveis quantitativos

[Alisson Lucrecio]

Caro Colegas r-br,

Boa tarde.

Nesse exemplo abaixo é certo transformar um dado quantitativo em fator para depois fazer a analise de variância?  Obrigado

Alisson Lucrécio da Costa


adi <- expand.grid(cult=gl(1,5,la=LETTERS[1]), fert=101)
fat <- expand.grid(cult=gl(5,5,la=LETTERS[2:6]), fert=seq(50,150,25))
da <- rbind(adi, fat)

theta <- c(c(-194.29, -197.26, -197.85, -203.03, -190.20, -190.45),
           c(9.1797, 8.2686, 8.6437, 9.3438, 8.8773, 8.1872),
           c(-0.03382, -0.03479, -0.03632, -0.03341, -0.03597, -0.03675))

X <- model.matrix(~-1+cult/(fert+I(fert^2)), data=da)
da$eta <- X%*%matrix(theta)
da$y <- da$eta+rnorm(nrow(da),0,30)

require(lattice)
xyplot(y~fert|cult, data=da, type=c("p","a"))

da$Fert <- factor(da$fert) #??????
levels(da$Fert)
levels(da$cult)

m0 <- aov(y~cult*Fert, data=da)
anova(m0)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130715/bf98fef0/attachment.html>