plotar erro padrão usando "arrows"

Prezados, Como plotar o erro padrão em um gráfico de médias usando a função "arrows"? Agradeço desde logo a atenção. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei@gmail.com ; ocersm@lycos.com; rsmartins@usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

exemplo de dados seria bem-vindo (por exemplo, media vs. o q? grupos? um determinado valor de preditor?) set.seed(1) medias <- rnorm(10) grupos <- 1:10 dvpd <- seq(.5, 1, length.out=10) lsup <- medias+dvpd linf <- medias-dvpd plot(grupos, medias, ylim=c(min(linf), max(lsup))) arrows(grupos, medias, y1=lsup) arrows(grupos, medias, y1=linf) b

Benilton, Obrigado por usa resposta. Na realidade, eu quero plotar o erro padrão da média, e não o desvio. Não encontrei nenhum pacote que faça isso automáticamente, por isso estou recorrendo a lista. Uma amostra dos meus dados para facilitar: Trimestre peso_gonada 1 0,089422162 1 0,033689594 1 0,179681794 1 0,078145005 1 0,174349771 1 0,439640105 1 0,045500911 1 0,513336169 1 0,037032364 1 0,981350496 2 0,129142957 2 0,04385171 2 0,158729421 2 0,05851425 2 0,060010159 2 0,566533747 3 0,196643195 3 0,044722969 3 0,213180475 3 0,089670273 3 0,057524721 3 0,045667584 3 0,03180651 3 0,09469697 4 0,473498822 4 0,081268399 4 0,382127041 4 0,404346428 4 0,114257631 4 0,273045579 4 0,041819872 Agradeço desde logo pela ajuda! Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei@gmail.com ; ocersm@lycos.com; rsmartins@usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

erro-padrao, desvio-padrao... para os comandos que passei, nao faz diferenca o q os numeros gravados em 'dvpd' representam... o comando sera' o mesmo e vc so' precisa colocar "os valores certos" em dvpd.

Entendi, Benilton! Muito obrigado! R. Em 21 de outubro de 2011 14:00, Benilton Carvalho < beniltoncarvalho@gmail.com> escreveu:
erro-padrao, desvio-padrao... para os comandos que passei, nao faz diferenca o q os numeros gravados em 'dvpd' representam... o comando sera' o mesmo e vc so' precisa colocar "os valores certos" em dvpd. _______________________________________________ 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.
-- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei@gmail.com ; ocersm@lycos.com; rsmartins@usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

Dê cabeça eu lembro dessas gplots::plotmeans(), gplots::barplot2(), plotrix::plotCI(), Hmisc::xYplot(), latticeExtra::segplot(), paych::error.bars(), psych::error.bars.by() mas no R Graphical Manual (adicione aos favoritos do seu navedor) tem muuuuita coisa http://www.oga-lab.net/RGM2/search.php?query=error+bars À 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 ==========================================================================

Obrigado, Walmes. Já tinha escrito um pequeno script com as dicas do Benilton, mas depois testei o erro.bars.by do psych e ele faz exatamente a emsma coisa... Cheers, R. -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei@gmail.com ; ocersm@lycos.com; rsmartins@usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777
participantes (3)
-
Benilton Carvalho
-
Rodrigo Plei
-
Walmes Zeviani