Adicionado índices do teste de médias no barplot
Boa tarde pessoal, Estou fazendo um barplot com a função barplot2 do pacote gplot, porém eu gostaria de colocar os índices do teste de média do Tukey sobre as barras do erro padrão, para isso fiz: ### require(gplot) # Criando o barplot múltiplo mp <- barplot2(TAB, beside = TRUE, axisnames = FALSE, 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
Alexandre, Não entendi a pergunta e o seu código não é reproduzível. Veja se este post é útil http://ridiculas.wordpress.com/2011/05/21/fatorial-com-desdobramento-da-inte... e verifique o conteúdo desse script que tem no final a confecção de gráficos com médias usando a barplot2() http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/cap14ancova-iso.R À 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 twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
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=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
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>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=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 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.
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> 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=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
>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.
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.
participantes (5)
-
Alexandre Santos -
ASANTOS
-
Mauro Sznelwar -
Walmes Zeviani -
Zé Augusto