É 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@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2013/6/14 Gilenio Borges Fernandes <gilenio@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@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.