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

[Gilenio Borges Fernandes]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130614/9a000d7a/attachment.html>