[R-br] erro bar

[Raphael Saldanha]

Olá Geovane,

Seguindo o site abaixo, consegui fazer algo próximo com os seus dados. Você vai precisar dos pacotes plyr e ggplot2:

http://www.cookbook-r.com/Graphs/Plotting_means_and_error_bars_%28ggplot2%29/



#Função para criar um sumário estatístico dos dados (atente para o NC de 95%)

summarySE <- function(data=NULL, measurevar, groupvars=NULL, na.rm=FALSE, conf.interval=.95, .drop=TRUE) {
  require(plyr)
  
  # New version of length which can handle NA's: if na.rm==T, don't count them
  length2 <- function (x, na.rm=FALSE) {
    if (na.rm) sum(!is.na(x))
    else       length(x)
  }
  
  # This is does the summary; it's not easy to understand...
  datac <- ddply(data, groupvars, .drop=.drop,
                 .fun= function(xx, col, na.rm) {
                   c( N    = length2(xx[,col], na.rm=na.rm),
                      mean = mean   (xx[,col], na.rm=na.rm),
                      sd   = sd     (xx[,col], na.rm=na.rm)
                   )
                 },
                 measurevar,
                 na.rm
  )
  
  # Rename the "mean" column    
  datac <- rename(datac, c("mean"=measurevar))
  
  datac$se <- datac$sd / sqrt(datac$N)  # Calculate standard error of the mean
  
  # Confidence interval multiplier for standard error
  # Calculate t-statistic for confidence interval: 
  # e.g., if conf.interval is .95, use .975 (above/below), and use df=N-1
  ciMult <- qt(conf.interval/2 + .5, datac$N-1)
  datac$ci <- datac$se * ciMult
  
  return(datac)
}


# Ao aplicar a função seus dados, alguns NAs são retornados pois para algumas combinações de raça e sexo, há apenas um indivíduo
dadosSE <- summarySE(dados, measurevar="Peso", groupvars=c("Raça","Sexo"))

library(ggplot2)

# Error bars represent standard error of the mean
ggplot(dadosSE, aes(x=Raça, y=Peso, fill=Sexo)) + 
  geom_bar(position=position_dodge(), stat="identity") +
  geom_errorbar(aes(ymin=Peso-se, ymax=Peso+se),
                width=.2,                    # Width of the error bars
                position=position_dodge(.9))





---
Atenciosamente,Raphael Saldanha
rfsaldanha em outlook.com

Date: Wed, 19 Jun 2013 04:22:28 -0700
From: geovanecb em yahoo.com.br
To: r-br em listas.c3sl.ufpr.br
Subject: [R-br] erro bar



Olá pessoal tudo bem, como faço para construir o seguinte erro bar. Quero construir-lo da seguinte forma:

no eixo x  gostaria de ter a variável raça e no eixo y a variável raça e a legenda por sexo.




 
 
  Raça
  Peso
  Sexo
 
 
  Branco
  54
  M
 
 
  Branco
  89
  M
 
 
  Branco
  65
  M
 
 
  Negro
  89
  H
 
 
  Branco
  80
  H
 
 
  Negro
  83
  M
 
 
  Negro
  65
  H
 
 
  Negro
  69
  M
 
 
  Amarela 
  68
  H
 
 
  Amarela 
  98
  H
 
 
  Parda
  120
  H
 
 
  Parda
  150
  M
 
 
  Parda
  115
  M
 
 
  Indígena
  56
  M
 
 
  Indígena
  59
  M
 
 
  Branco
  88
  M
 



obrigado pessoal

_______________________________________________
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/20130619/91a20a02/attachment.html>