Em 16/12/2013 08:32, Ari Clecius escreveu:
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.
Você pode utilizar a função TukeyHSD(lm.1,c("A","B","C"), essa função
pertenção ao pacote stats (ou seja ela vem instalada com o R)
---
Este email está limpo de vírus e malwares porque a proteção do avast! Antivírus
está ativa.
http://www.avast.com
_______________________________________________
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.
