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 <[email protected]> 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 > [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. -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ 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 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ _______________________________________________ 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.
