[R-br] Adicionado índices do teste de médias no barplot [Resolvido]
ASANTOS
alexandresantosbr em yahoo.com.br
Quarta Janeiro 4 10:08:17 BRST 2012
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 em gmail.com>
> *Para:* r-br em 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 em uol.com.br
> <mailto:sznelwar em 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 em listas.c3sl.ufpr.br <mailto: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.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br <mailto: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.
>
>
>
> _______________________________________________
> 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/20120104/292dfb8a/attachment.html>
Mais detalhes sobre a lista de discussão R-br