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@uol.com.br> escreveu: De: Mauro Sznelwar <sznelwar@uol.com.br> Assunto: Re: [R-br] comparação de substratos Para: r-br@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@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br