[R-br] Teste de hipoteses em análise de grupo de experimentos

walmes . walmeszeviani em gmail.com
Sábado Junho 15 12:04:26 BRT 2013


É isso?

resp <- c(109.30,124.00,100.00,87.45,115.12,93.10,95.60,101.00,98.50,
          68.00,47.00,78.00,118.00,141.60,117.70,112.30,91.30,125.65,115.00,

109.00,97.50,104.00,78.00,87.00,117.70,127.20,116.50,108.00,156.70,
          93.10,83.00,111.50,114.00,128.00,65.00,84.00,122.50,114.60,129.00,
          108.00,99.75,135.40,111.00,140.00,83.50,89.00,88.50,60.00,126.60,

112.00,123.60,103.20,113.00,158.50,130.00,119.00,87.50,74.00,82.00,67.00)
trat <- gl(5,12,label=paste("N",1:5,sep=""))
exp <- rep(gl(4,3,label=paste("E",1:4,sep="")),5)
bloco <- rep(gl(3,1,label=paste("bloco",1:3,sep="")),20)
dados <- data.frame(trat=trat,exp=exp,bloco=bloco,resp=resp)

xtabs(~trat+bloco+exp, data=dados)

mod <- aov(resp ~ trat * exp + exp/bloco, data = dados)

# análise individuais
mi <- lapply(split(dados, f=dados$exp), aov, formula=resp~bloco+trat)
lapply(mi, summary)

# análise conjunta
mod <- aov(resp~exp/bloco+trat+Error(exp:trat), data=dados)
summary(mod)

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes em ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


2013/6/14 Gilenio Borges Fernandes <gilenio em ufba.br>

> Gostaria de contar com ajuda para os seguinte problema - analise conjunta
> de 4 experimentos em blocos casualizados completos com 3 repetições
>
> resp=c(109.30,124.00,100.00,87.45,115.12,93.10,95.60,101.00,98.50,
>        68.00,47.00,78.00,118.00,141.60,117.70,112.30,91.30,125.65,115.00,
>        109.00,97.50,104.00,78.00,87.00,117.70,127.20,116.50,108.00,156.70,
>        93.10,83.00,111.50,114.00,128.00,65.00,84.00,122.50,114.60,129.00,
>        108.00,99.75,135.40,111.00,140.00,83.50,89.00,88.50,60.00,126.60,
>
>  112.00,123.60,103.20,113.00,158.50,130.00,119.00,87.50,74.00,82.00,67.00)
> trat=gl(5,12,label=paste("N",1:5,sep=""))
> exp<-rep(gl(4,3,label=paste("E",1:4,sep="")),5)
> bloco<-rep(gl(3,1,label=paste("bloco",1:3,sep="")),20)
> dados=data.frame(trat=trat,exp=exp,bloco=bloco,resp=resp)
> mod = aov(resp ~ trat * exp + exp/bloco, data = dados)
> summary(mod)
>
>
> A declaração do modelo (mod) com a função aov testa tudo com o resíduo
> (32 gl). Quero testar tratamento (trat) e experimento (exp) com a a
> interação trat*exp e
> esta com o resíduo (32 gl). Como declarar os efeitos no modelo?
>
>
> --
> Gilenio Borges Fernandes
> Professor Associado III
> Universidade Federal da Bahia
> Instituto de Matemática
> Departamento de Estatística
> Av. Adhemar de Barros, s/n – Ondina.
> 40.170-110 - Salvador - BA, Brasil
> Tel.: (071)3283-6280/6336  Fax:  (071)3283-6276
> URL: http://lattes.cnpq.br/6764860618464860
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130615/c653c802/attachment.html>


Mais detalhes sobre a lista de discussão R-br