[R-br] Desdobramento da interação via glht:multcomp

Emmanuel Arnhold emmanuelarnhold em yahoo.com.br
Sexta Abril 1 11:29:34 BRT 2011


# Colegas ... no exemplo abaixo tento realizar uma análise do desdobramento (via glht:multcomp) do fator Proteína dentro do nível L1 e L2 do fator Linhagens.
 
#Quem tiver um tempo e interesse gostaria que desse uma olhada e me desse uma opinião ou uma forma mais fácil de fazer o desdobramento
 





Proteína

Linhagens

Peso


p1

L1

1.6


p1

L1

1.7


p1

L2

2.3


p1

L2

2.2


p2

L1

1.9


p2

L1

2


p2

L2

2.2


p2

L2

2.3


p3

L1

2.3


p3

L1

2.4


p3

L2

2.3


p3

L2

2.3

y=read.table("clipboard", h=T)
 
attach(y)
 
m=aov(Peso~Proteína+Linhagens+Proteína*Linhagens )
anova(m)
 
 # análise via agricolae
 
require(agricolae)
df<-df.residual(m)MSError<-deviance(m)/ df.residual(m)
 
 
#Melhor teor de proteína dentro da linhagem 1

with(subset(y, Linhagens == "L1"), HSD.test(Peso,Proteína, df,  MSError, group=F))

#Melhor teor de proteína dentro da linhagem 2

with(subset(y, Linhagens == "L2"), HSD.test(Peso,Proteína, df,  MSError, group=F))
 
 
 
 
 
 
 
 
# análise via multcomp
 
require(multcomp)
 
a=subset(y, Linhagens == "L1")
b=subset(y, Linhagens == "L2")
 
##Melhor teor de proteína dentro da linhagem 1, OK!!!!
 
y1=with(a, glht (m,  linfct = mcp(Proteína = "Tukey" )))
summary(y1)
 
#Melhor teor de proteína dentro da linhagem 2, RODA DENTRO DA LINHAGEM 1?????

y2=with(b, glht (m,  linfct = mcp(Proteína ="Tukey")))
summary(y2)
 
 
 
#Criei um sistema que segue abaixo. Consiste em trocar os códigos dos níveis do fator Linhagens
 
Linhagens=c('L2', 'L2', 'L1', 'L1', 'L2', 'L2', 'L1', 'L1', 'L2', 'L2', 'L1', 'L1')
Linhagens=as.factor(Linhagens)
y=data.frame(Proteína, Linhagens, Peso)
 
a=subset(y, Linhagens == "L1")
m2=aov(Peso~Proteína+Linhagens+Proteína*Linhagens)
 
 
#Melhor teor de proteína dentro da linhagem 2, OK!!! Agora temos o resultado correto!!!
y1=with(a, glht (m2,  linfct = mcp(Proteína = "Tukey" )))
summary(y1)
 
 
Perguntas.
 
Alguém tem uma solução mais fácil? Prática?
 
Como trocar o código dos níveis dos fatores de maneira mais fácil?
 
 
 
Obrigado desde já por qualquer colaboração.
 
Emmanuel Arnhold
 
 
 
 
 
 
 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110401/4be9cbc6/attachment.html>


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