Caro Edson,
segue um exemplo

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


________________________________
 From: Edson Lira <[email protected]>
To: "[email protected]" <[email protected]>; Alisson Lucrecio 
<[email protected]> 
Sent: Monday, July 15, 2013 12:33 PM
Subject: Re: [R-br] Dados quantitativos em uma análise de variança
 


Alisson, contextualize o seu problema.

Variáveis quantitativas podem ser categorizadas e tratadas como fatores.

Depende do seu objetivo.
 
Edson Lira
Estatístico
Manaus-Amazonas


________________________________
 De: Alisson Lucrecio <[email protected]>
Para: "[email protected]" <[email protected]> 
Enviadas: Segunda-feira, 15 de Julho de 2013 10:48
Assunto: [R-br] Dados quantitativos em uma análise de variança
 


Caro Colegas r-br,

Bom dia.

Eu estou com uma duvida sobre dados quantitativos numa análise de variância se 
eu os considero como quantitativo ou como fator?

Obrigado.
 
Alisson Lucrécio da Costa
_______________________________________________
R-br mailing list
[email protected]
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.
_______________________________________________
R-br mailing list
[email protected]
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.

Responder a