gráfico de médias com desvio determinado pelo usuário
Prezados, Gostraia de fazer um gráfico mostrando a média de porcentagem larvas transportadas por localidade, tendo como a linha de desvio o intervalo de confianca 95% (± 95% IC) da média dos mínimos quadrados (least square means) É possível fazer isso usando os recursos gráficos do R? A base de dados que tenho é assim: Local propLarvasTransp 1 0.75 1 0.37 1 0.42 1 0.25 2 0 2 0.22 2 0.15 2 0.05 3 0.11 3 0.14 3 0.07 3 0.16 . . . Encontrei em livros e na internet como plotar todos os tipos de desvios, menos esse bendito intervalo de confiança. Agradeço qualquer ajuda desde logo. 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
Rodrigo, Você pode usar coisas prontas no R, como http://addictedtor.free.fr/graphiques/graphcode.php?graph=54 http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=72 Ou pode construir com os recursos gráficos do R, de onde é possível ter mais controle, veja da <- data.frame(trat=gl(5,8)) da$y <- as.numeric(da$trat)+rnorm(nrow(da)) m0 <- lm(y~trat, da) new <- data.frame(trat=levels(da$trat)) new$pred <- predict(m0, newdata=new, interval="confidence") str(new) ylim <- c(0, max(new$pred)*1.05) bp <- barplot(new$pred[,1], ylim=ylim) arrows(bp, new$pred[,2], bp, new$pred[,3], code=3, angle=90) box() À 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 ==========================================================================
Obrigadão, Walmes! Vou tentar os gráficos prontos e também tentar adaptar o seu pequeno script para meus dados. 1 abç! 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
participantes (2)
-
Rodrigo Plei -
Walmes Zeviani