[R-br] Teste de Tukey em experimento fatorial

[Ailton Andrade de Oliveira]

Uma outra sugestão é neste caso não usar o cbind. Ele está 
"desfatorizando" A, B e C.

exp4=data.frame(A,B,C,y)
exp4aov=aov(y~A*B*C,data=exp4)
TukeyHSD(exp4aov)

 >Não conferi tudo, mas creio que o que está incorreto é esta parte:
 >
 >exp3=data.frame(cbind(A,B,C,y))
 >exp3
 >exp3$A=as.factor(A)
 >exp3$B=as.factor(B)
 >exp3$C=as.factor(C)
 >exp3aov=aov(y~A*B*C,data=exp3)
 >summary(exp3aov)
 >TukeyHSD(exp3aov)
 >
 >On Mon, Dec 16, 2013 at 9:32 AM, Ari Clecius <ari072000 em 
yahoo.com.br> wrote:
 >> Prezados, como faço para realizar o teste de Tukey neste experimento?
 >>
 >> Att,
 >>
 >> Ari Clecius Alves de Lima
 >> Engenheiro Químico
 >> Me. Engenharia Civil
 >> (085)87687786
 >> (085)33669042
 >>
 >> # Experimento factorial completo
 >> library(qualityTools)
 >> exp1 = facDesign( k = 3 , replicates=2)
 >> names ( exp1 ) = c ( " pH " , " Volume da amostra" , "Volume de 
eluição" )
 >> units(exp1)=c("", "mL", "mL")
 >> lows ( exp1 ) = c (3, 250,2)
 >> highs ( exp1 ) = c (7.5,500,4)
 >> summary( exp1 )
 >> y=c(301194.3, 163197.05, 236720.8, 67201.2, 142970.3, 145672.55, 
346889.9,
 >> 1470375.4, 156877.15, 74782.75, 129978, 144724.8, 67715.85, 74768.9,
 >> 149722.25, 168859.55)
 >> response(exp1)=data.frame(y)
 >> summary( exp1 )
 >> lm.1=lm(y~A*B*C,data=exp1)
 >> summary(lm.1)
 >>
 >> exp2=data.frame(exp1)
 >> A=c(exp2$A)
 >> B=c(exp2$B)
 >> C=c(exp2$C)
 >> y=(exp2$y)
 >> exp3=data.frame(cbind(A,B,C,y))
 >> exp3
 >> A=as.factor(A)
 >> B=as.factor(B)
 >> C=as.factor(C)
 >> exp3aov=aov(y~A*B*C,data=exp3)
 >> TukeyHSD(exp3aov)
 >>
 >> Error in TukeyHSD.aov(exp3aov) : nenhum fator no modelo ajustado
 >> In addition: Warning messages:
 >> 1: In replications(paste("~", xx), data = mf) : non-factors ignored: A
 >> 2: In replications(paste("~", xx), data = mf) : non-factors ignored: B
 >> 3: In replications(paste("~", xx), data = mf) : non-factors ignored: C
 >> 4: In replications(paste("~", xx), data = mf) : non-factors ignored: 
A, B
 >> 5: In replications(paste("~", xx), data = mf) : non-factors ignored: 
A, C
 >> 6: In replications(paste("~", xx), data = mf) : non-factors ignored: 
B, C
 >> 7: In replications(paste("~", xx), data = mf) :
 >>   non-factors ignored: A, B, C
 >>
 >> _______________________________________________
 >> 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.
 >
 >
 >
 >--
 >///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
 >Jose Claudio Faria
 >Estatistica
 >UESC/DCET/Brasil
 >joseclaudio.faria at gmail.com
 >Telefones:
 >55(73)3680.5545 - UESC
 >55(73)9100.7351 - TIM
 >55(73)8817.6159 - OI
 >///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

-- 
Um abraço,
Ailton Andrade de Oliveira