[R-br] comparação de substratos

Emmanuel Arnhold emmanuelarnhold em yahoo.com.br
Sexta Abril 29 08:31:40 BRT 2011


Olá Mauro. Copiei a tabela para um editor de texto e de lá para o R. Depois colei os comandos no R e não ocorreu erro. Não sei o que esta ocorrendo no seu caso. Talves você não criou os objetos com o grau de liberdade e quadrado médio do erro. Se quiser pode usar:
 
df<-df.residual(resultado)
MSerror<-deviance(resultado)/df
 
Ok.


--- Em sex, 29/4/11, Mauro Sznelwar <sznelwar em uol.com.br> escreveu:


De: Mauro Sznelwar <sznelwar em uol.com.br>
Assunto: Re: [R-br] comparação de substratos
Para: r-br em listas.c3sl.ufpr.br
Data: Sexta-feira, 29 de Abril de 2011, 1:38





Estava tentando rodar os comandos, e deu problema a partir deste:
with(subset(y, Linhagens == "L1"), HSD.test(Peso,Proteína, df,  MSError))
Erro em qtukey(p, nranges, nmeans, df, lower.tail, log.p) : 
  Argumento não-numérico para função matemática







Se entendi bem, o exemplo abaixo pode ajudar. 

  






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) 

  

resultado=aov(Peso~Proteína+Linhagens+Proteína*Linhagens ) 

  

anova(resultado) 

  

require(agricolae) 

  

cv.model(resultado) 

  

 df=6 

MSError=0.004167 

  

  

#Melhor teor de proteína dentro da linhagem 1 

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

  

#Melhor teor de proteína dentro da linhagem 2 

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

  

#Melhor linhagem dentro do teor de proteína de 14% 

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

  

#Melhor linhagem dentro do teor de proteína de 16% 

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

  

#Melhor linhagem dentro do teor de proteína de 18% 

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

  






Proteína 


Linhagem 


D.M.S.(5%) Linha 



L1 


L2 



14% 


1,65 cB 


2,25 aA 


0,1580 



16% 


1,95 bB 


2,25 aA 



18% 


2,35 aA 


2,30 aA 



D.M.S.(5%) Coluna 


0,1981 


  

Médias seguidas de uma mesma letra maiúscula (em uma mesma linha) e minúscula (em uma mesma coluna), são estatisticamente iguais pelo teste de Tukey ao nível de 5% de probabilidade. 

  






 
 


-----Anexo incorporado-----


_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110429/179301a8/attachment-0001.html>


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