[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