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))

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.

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