[R-br] Teste de Tukey em experimento fatorial

Tiago Souza Marçal tiagosouzamarcal em hotmail.com
Segunda Dezembro 16 12:50:28 BRST 2013


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 em yahoo.com.br
To: r-br em 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 LimaEngenheiro QuímicoMe. Engenharia Civil(085)87687786(085)33669042
# Experimento factorial completolibrary(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))exp3A=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 ajustadoIn addition: Warning messages:1: In replications(paste("~", xx), data = mf) : non-factors ignored: A2: In replications(paste("~", xx), data = mf) : non-factors ignored: B3: In replications(paste("~", xx), data = mf) : non-factors ignored: C4: In replications(paste("~", xx), data = mf) : non-factors ignored: A,
 B5: In replications(paste("~", xx), data = mf) : non-factors ignored: A, C6: In replications(paste("~", xx), data = mf) : non-factors ignored: B, C7: 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. 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131216/5b768269/attachment.html>


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