[R-br] Teste de Tukey em experimento fatorial

[Jose Claudio Faria]

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
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\