[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