[R-br] Digest R-br, volume 75, assunto 8

Domingo Março 12 20:13:37 BRT 2017

Obrigado Walmes,

Obrigado pela indicação. Mas acabei optando por realizar o cálculo dos quadrados médios combinados e graus de liberdade combinados manualmente mesmo, era isso que eu estava na dúvida... e de fato deu resultado diferente do método que estava de exemplo no pacote agricolae.

Acredito eu que está certo como fiz.

Se a alguém interessar, ou se disponibilizar de conferir, mando os dados e comandos:

# Análise de um experimento em faixas com interação significativa:

dados<- structure(list(local = c(1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4,
1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1,
2, 3, 4), trat = c(0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60,
0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60, 0, 0, 0, 0, 30, 30,
30, 30, 60, 60, 60, 60), bloco = c(1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3,
3, 3, 3, 3, 3, 3), yield = c(150, 163, 164, 171, 188, 174, 182,
182, 220, 170, 170, 182, 147, 141, 133, 140, 157, 166, 144, 143,
173, 165, 152, 170, 146, 140, 149, 150, 172, 158, 142, 160, 204,
163, 150, 175)), .Names = c("local", "trat", "bloco", "yield"
), row.names = c(NA, 36L), class = "data.frame")

model<-with(dados,strip.plot(bloco, trat, local, yield))

# Se interação não significativa:

comparison1<-with(dados,HSD.test(yield,trat,model$gl.a,model$Ea));comparison1

comparison2<-with(dados,HSD.test(yield,local,model$gl.b,model$Eb)); comparison2

# Se tiver interação significativa:

QMea<- model$Ea ;Dfea<-model$gl.a
QMeb<- model$Eb; Dfeb<-model$gl.b
QMec<- model$Ec ; Dfec<-model$gl.c

a <-3; b <-4

#Erro medio a e c
QMeac <-(QMea+(b-1)*QMec)/b

##Graus de liberdade a e c
Glac<-((QMea+(b-1)*QMec)^2) / ( ((QMea^2)/Dfea)  + ((((b-1)*QMec)^2)/Dfec) )

#Erro medio b e c
QMebc<-(QMeb+(a-1)*QMec)/a

##Graus de liberdade a e c
Glbc<-((QMeb+(a-1)*QMec)^2) / ( ((QMeb^2)/Dfeb)  + ((((a-1)*QMec)^2)/Dfec) )

#---------------------------------------------------------------------------------
# Comparação de médias por Tukey:

## Tratamentos dentro de locais:

L1<-subset(dados,dados$local=="1")
t1<-HSD.test(L1$yield,L1$trat,Glac,QMeac)
t1

L2<-subset(dados,dados$local=="2")
t2<-HSD.test(L2$yield,L2$trat,Glac,QMeac)
t2

L3<-subset(dados,dados$local=="3")
t3<-HSD.test(L3$yield,L3$trat,Glac,QMeac)
t3

L4<-subset(dados,dados$local=="4")
t4<-HSD.test(L4$yield,L4$trat,Glac,QMeac)
t4

## Locais dentro de tratamentos:

A1<-subset(dados,dados$trat=="0")
ta1<-HSD.test(A1$yield,A1$local,Glbc,QMebc)
ta1

A2<-subset(dados,dados$trat=="30")
ta2<-HSD.test(A2$yield,A2$local,Glbc,QMebc)
ta2

A3<-subset(dados,dados$trat=="60")
ta3<-HSD.test(A3$yield,A3$local,Glbc,QMebc)
ta3

# FIM ?

Abraços!

..................
Me. Eng. Agr. Maurício Sangiogo
TEL:(55) 96822030

________________________________
De: R-br <r-br-bounces em listas.c3sl.ufpr.br> em nome de r-br-request em listas.c3sl.ufpr.br <r-br-request em listas.c3sl.ufpr.br>
Enviado: domingo, 12 de março de 2017 12:00
Para: r-br em listas.c3sl.ufpr.br
Assunto: Digest R-br, volume 75, assunto 8

Enviar submissões para a lista de discussão R-br para
        r-br em listas.c3sl.ufpr.br

Para se cadastrar ou descadastrar via WWW, visite o endereço
        https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br>
listas.inf.ufpr.br
A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português).

ou, via email, envie uma mensagem com a palavra 'help' no assunto ou
corpo da mensagem para
        r-br-request em listas.c3sl.ufpr.br

Você poderá entrar em contato com a pessoa que gerencia a lista pelo
endereço
        r-br-owner em listas.c3sl.ufpr.br

Quando responder, por favor edite sua linha Assunto assim ela será
mais específica que "Re: Contents of R-br digest..."

Tópicos de Hoje:

   1. Re: Desdobramento de um experimento em faixas (Walmes Zeviani)

----------------------------------------------------------------------

Message: 1
Date: Sat, 11 Mar 2017 17:44:03 -0300
From: Walmes Zeviani <walmeszeviani em gmail.com>
To: Maurício Sangiogo <ms_sangiogo em hotmail.com>, a lista Brasileira
        oficial de discussão do programa R.  <r-br em listas.c3sl.ufpr.br>
Subject: Re: [R-br] Desdobramento de um experimento em faixas
Message-ID:
        <CAFU=EkYuWZnFRExk8eNhcLSq3u1=OG0Yg=ZQUNeGoxQQc8obzQ em mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Verifique se essa postagem do histórico da lista serve:
http://r-br.2285057.n4.nabble.com/R-br-Anova-em-faixas-td4655899.html. Caso
R-br - [R-br] Anova em faixas<http://r-br.2285057.n4.nabble.com/R-br-Anova-em-faixas-td4655899.html>
r-br.2285057.n4.nabble.com
[R-br] Anova em faixas. Srs, Preciso de uma ajuda. Preciso fazer uma ANOVA de um experimento com três fatores, porém o problema está no delineamento que foi em faixas. Tenho o quadro da análise...

não funcione, envie no corpo da mensagem os resultados que você tem como
sendo os certos e o resultado que você está conseguindo. Não esqueça de
disponibilizar acesso aos dados usados no corpo da mensagem (use dput()).

À disposição.
Walmes.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170311/02ec2ff1/attachment-0001.html>

------------------------------

Subject: Legenda do Digest

_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br>
listas.inf.ufpr.br
A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português).

------------------------------

Fim da Digest R-br, volume 75, assunto 8
****************************************
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170312/f17a8af0/attachment.html>