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

Muito obrigado a todos, vou tentar usar o ExpDes. Att,   Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)87687786 (085)33669042 Em Segunda-feira, 16 de Dezembro de 2013 11:50, Tiago Souza Marçal <tiagosouzamarcal@hotmail.com> escreveu: Uma opção mais simples seria utilizar o package "ExpDes". #################################################################   Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES   Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas   #################################################################    ________________________________ Date: Mon, 16 Dec 2013 03:32:41 -0800 From: ari072000@yahoo.com.br To: r-br@listas.c3sl.ufpr.br Subject: [R-br] Teste de Tukey em experimento fatorial 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@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. _______________________________________________ R-br mailing list R-br@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.