Colegas, estou iniciando o aprendizado com a plataforma R. Alguém tem dicas, ou sumário dos principais comandos e funções para iniciantes? Obrigada! Cláudia christina, ouvinte mestrado Epidemiologia/saude coletiva/IMS/UERJ Em 1 de abril de 2011 11:29, Emmanuel Arnhold <emmanuelarnhold@yahoo.com.br>escreveu:
# 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*
* *
* *
* *
* *
* *
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br