[R-br] Teste de Tukey em experimento fatorial

Fernando Souza nandodesouza em gmail.com
Segunda Dezembro 16 17:27:17 BRST 2013


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


Mais detalhes sobre a lista de discussão R-br