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

[Edson Lira]

Dá uma olhada nos posts da lista que há uns dias atrás, essa discussão estava em pauta. Mas acho que é possível. 

[. ]'s.
Edson Lira
Estatístico
Ma-Am

Em 23/07/2013, às 10:34, Alisson Lucrecio <alissonluc em yahoo.com.br> escreveu:

> 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)
> 
> Obrigado.
>  
> Alisson Lucrécio da Costa
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130723/96c646a7/attachment.html>