Reputation: 121
I'm trying to add errorsbars to a graph I created:
As you can see the error bars are not in the right location and I can't figure out why not. The following is the code I used:
df_midtermclean$gender <- factor(df_midtermclean$ï..gender, labels = c('Male', 'Female'))
#reshaping from wide to long format
data_long <- melt(df_midtermclean,
# ID variables - all the variables to keep but not split apart on
id.vars=c("gender"),
# The source columns
measure.vars=c("sSkills", "sPerform", "sComplex", "sIQ" ),
# Name of the destination column that will identify the original
# column that the measurement came from
variable.name="measures",
value.name="score")
p2 <- ggplot(data_long, aes(measures, score)) + stat_summary(fun.y = mean, geom = "bar", aes(fill = gender), position = "dodge")
p2 <- p2 + stat_summary(geom = "errorbar", fun.data = mean_cl_normal, position = position_dodge(.9), width = .1)
p2 <- p2 + labs(x = "Continuous Measures", y = "Participant Scores", fill = "Gender")
p2
Data.
data_long <-
structure(list(gender = structure(c(1L, 1L, 1L,
1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L,
2L, 2L, 1L, 1L, 1L), .Label = c("Male", "Female"),
class = "factor"), measures = structure(c(1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L), .Label = c("sSkills",
"sPerform", "sComplex", "sIQ"), class = "factor"),
score = c(4.75, 3, 3.5, 3.75, 3, 3.25, 3.25, 3,
2.75, 2.75, 2.5, 2.75, 3.75, 3.5, 3.75, 3.75,
3.25, 3.25, 4, 3.25)), row.names = c(NA, 20L),
class = "data.frame")
Upvotes: 0
Views: 760
Reputation: 76673
I believe that something like the following will do what you want.
agg <- aggregate(score ~ gender + measures, data_long, function(x) c(mu = mean(x), se = sd(x)))
agg <- cbind(agg[-ncol(agg)], agg[[ncol(agg)]])
p2 <- ggplot(agg, aes(measures, mu, fill = gender, group = gender)) +
geom_bar(stat = "identity", position = "dodge") +
geom_errorbar(aes(ymin = mu - se, ymax = mu + se),
position = position_dodge(0.9),
width = 0.1)
p2
Upvotes: 0
Reputation: 419
I think the second line of your code should be of the following format:
geom_errorbar(aes(ymin = lower, ymax = upper), position = dodge, width = 0.25)
I am not sure what stat_summary is in your code, can you explain that ? However, you need to explicitly mention the lower and upper bounds of the error bars, I don't think that is included in your code.
Edit: try tweaking the position=position=position_dodge(0.9), along with the size and widths.
Upvotes: 0
Reputation: 1004
Unfortunatly your code is unreproducible ! However, I may be found a solution on http://www.sthda.com/english/wiki/ggplot2-error-bars-quick-start-guide-r-software-and-data-visualization
library(ggplot2)
df <- ToothGrowth
df$dose <- as.factor(df$dose)
# Function to calculate the mean and the standard deviation
data_summary <- function(data, varname, groupnames){
require(plyr)
summary_func <- function(x, col){
c(mean = mean(x[[col]], na.rm=TRUE),
sd = sd(x[[col]], na.rm=TRUE))
}
data_sum <- ddply(data, groupnames, .fun=summary_func, varname)
data_sum <- rename(data_sum, c("mean" = varname))
return(data_sum)
}
# New dataframe
df2 <- data_summary(ToothGrowth, varname="len",
groupnames=c("supp", "dose"))
df2$dose <- as.factor(df2$dose)
ggplot(df2, aes(x=dose, y=len, fill=supp)) +
geom_bar(stat="identity", color="black", position=position_dodge()) +
geom_errorbar(aes(ymin=len-sd, ymax=len+sd), width=.2, position=position_dodge(.9))
Upvotes: 1