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)