[R-br] Teste de Tukey em experimento fatorial

[Ari Clecius]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131216/8c9253ef/attachment.html>