[R-br] comparação de substratos

[Emmanuel Arnhold]

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 em gmail.com> escreveu:


De: Gilson Sanchez <gilsonsch em gmail.com>
Assunto: [R-br] comparação de substratos
Para: "Grupo R-Br" <r-br em 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 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/20110428/fa6d1367/attachment.html>