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. --- Em qui, 28/4/11, Gilson Sanchez <gilsonsch@gmail.com> escreveu: De: Gilson Sanchez <gilsonsch@gmail.com> Assunto: [R-br] comparação de substratos Para: "Grupo R-Br" <r-br@listas.c3sl.ufpr.br> Data: Quinta-feira, 28 de Abril de 2011, 14:13 Pessoal mais uma vez precissando ajuda de vc's Estou com um experimento de 3 especies de mudas com 6 tratamentos (substratos) e 3 repetições (exemplo abaixo) eu estou fazendo uma comparação dos substratos por espécie individualmete, mas eu não consigo fazer a compração entre os tratamentos com as espécies, ou seja eu quero comparar trata1 entre as espe1, espe2, e espe3; trat2 entre as esp1, esp2 e esp3 e assim por diante. algem pode ajudar. assi eu teria uma tabela com Tukey na linha compração dos substratos por éspecie indivualmente e na coluna comparação dos substratos entre as espécies. obrigado ############################### ### ANALIZANDO VARIAVEL nf1 ### ############################### #ANÁLISE DE VARIÂNCIA anova.nf<-aov(nf~trat, data = dados) anova(anova.nf) #COMPRAÇÃO DAS MÉDIAS TUKEY 5% require(laercio) LTukey(anova.nf,"trat") # TESTE DE TUKEY ### TESTES DE HOMOGENEIDADE DE VARIÂNCIAS bartlett.test(nf~trat) ### TESTE DE NORMALIDADE shapiro.test(residuals(anova.nf)) Dados: especie trat nf dc alt pts... esp1 1 2 3 4 5 esp1 1 2 3 4 5 esp1 1 2 3 4 5 esp1 2 2 3 4 5 esp1 2 2 3 4 5 esp1 2 2 3 4 5 ... esp1 6 2 3 4 5 esp2 1 2 3 4 5 esp2 1 2 3 4 5 esp2 1 2 3 4 5 esp2 2 2 3 4 5 esp2 2 2 3 4 5 esp2 2 2 3 4 5 ... esp2 6 2 3 4 5 esp3 1 2 3 4 5 esp3 1 2 3 4 5 esp3 1 2 3 4 5 esp3 2 2 3 4 5 esp3 2 2 3 4 5 esp3 2 2 3 4 5 ... esp1 6 2 3 4 5 -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows. -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br