Bom dia pessoal, Conforme rotina do Walmes eu consegui solucionar o problema, obrigado pela ajuda de todos, segue exemplo com solução: require(gplots) #tratmentos no tempo1 trat1<-rep(1:3,100) resp1<-rpois(300,25) med1<-tapply(resp1, trat1,mean)##Media dos tratamentos ep1<-tapply(resp1,trat1,sd)/sqrt(tapply(resp1,trat1,length))###Erros padrao da media #tratmentos no tempo2 trat2<-rep(1:3,100) resp2<-rpois(300,10) med2<-tapply(resp2, trat2,mean, na.rm = TRUE)##Media dos tratamentos ep2<-tapply(resp2,trat2,sd, na.rm = TRUE)/sqrt(tapply(resp2,trat2,length))###Erros padrao da media #tratmentos no tempo3 trat3<-rep(1:3,100) resp3<-rpois(300,5) med3<-tapply(resp3, trat3,mean)##Media dos tratamentos ep3<-tapply(resp3,trat3,sd)/sqrt(tapply(resp3,trat3,length))###Erros padrao da media ### # Montando as tabelas para cada tempo: a<- med1 b<- med2 c<-med3 a2<-ep1 b2<-ep2 c2<-ep3 # Unindo as tabelas TAB <- cbind(a,b,c) TAB2<-cbind(a2,b2,c2) # Criando o barplot múltiplo mp <- barplot2(TAB, beside = TRUE, axisnames = FALSE, plot.ci=T, ci.u=TAB+TAB2,ci.l=TAB-TAB2, legend.text=T, ylim=c(0,30),col =c("black","grey","white")) # Adicionando os tempos mtext(1, at = colMeans(mp), text = rep(c("Tempo1", "Tempo2", "Tempo3"), 3),line = 1, cex = 0.75) ##Adicionado os indices do teste de médias text(mp, TAB, label=paste(format(TAB, dig=4), c("a","a","a","a","a","a","a","a")), pos=3,cex=0.85) -- Alexandre DOS SANTOS Engenheiro Florestal, Msc. Laboratório de Entomologia Florestal Departamento de Entomologia Universidade Federal de Lavras Caixa Postal 3037 37200-000 - Lavras - Minas Gerais - Brasil Tel: +55 35 92230304 Em 03-01-2012 21:39, Alexandre Santos escreveu:
Obrigado Walmes, Vou tentar fazer pelo seu link, após resolver coloco um código reproduzível com a solução. E desculpa Mauro o pacote é o gplots. *Alexandre dos Santos *Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122
------------------------------------------------------------------------ *De:* Zé Augusto <jabnaz@gmail.com> *Para:* r-br@listas.c3sl.ufpr.br *Enviadas:* Terça-feira, 3 de Janeiro de 2012 21:15 *Assunto:* Re: [R-br] Adicionado índices do teste de médias no barplot
Mauro
O que existe é o pacote "gplots". Veja em http://crantastic.org/packages/gplots
José Augusto
Em 3 de janeiro de 2012 21:01, Mauro Sznelwar <sznelwar@uol.com.br <mailto:sznelwar@uol.com.br>> escreveu:
Estava tentando reproduzir este código, mas não consegui achar esta bibliotéca gplot, diz que não existe no R14.0!
### require(gplot) # Criando o barplot múltiplo mp <- barplot2(TAB, beside = TRUE, axisnames = FALSE, plot.ci <http://plot.ci>=T, ci.u=TAB+TAB2,ci.l=TAB-TAB2, ylab="Number of Bradysia difformis", legend.text=T) # Adicionando os meses mtext(1, at = colMeans(mp), text = rep(c("April", "May", "June", "July","August", "September","October"), 7),line = 1, cex = 0.75)
Porém gostaria de saber se existe algum comando que consigo especificar um vetor contendo os índices e faze-los aparecer sobre as barras do erro padrão, ao invés de especificar índice por índice como text().
Obrigado,
-- Alexandre DOS SANTOS Engenheiro Florestal, Msc. Laboratório de Entomologia Florestal Departamento de Entomologia Universidade Federal de Lavras Caixa Postal 3037 37200-000 - Lavras - Minas Gerais - Brasil Tel: +55 35 92230304
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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.
_______________________________________________ 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.