Oi Gilson, Não tenho certeza se entendi bem. Tente isso aqui: with(seus_dados, tapply(sua_variável_de_interesse_a_ser_observada, list(substrato, espécie), mean)) Abraço!! Em 28 de abril de 2011 10:13, Gilson Sanchez <gilsonsch@gmail.com> escreveu:

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.

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

-- Nory Daniel de Carvalho Erazo Técnico em Sistemática Vegetal/Herbário Coordenação de Pesquisas em Botânica-CPBO Instituto Nacional de Pesquisas da Amazônia-INPA Av. André Araújo, 2936, Aleixo, CEP 69060-001, Manaus-AM Fone: (92) 3643-3125 Celular: (92) 8196-4472 E-mail: nory@inpa.gov.br norydaniel@yahoo.com.br norydaniel@hotmail.com norydaniel@gmail.com

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