Hey guys,

Really hoping you can help. I have a data set with dengue incidence over time for all 33 departments of Colombia. I have created a simple line graph showing the mean annual incidence for all departments (2007-2018) using a data frame that has only 12 observations and 2 variables (incidence and year).

However I want to add error bars to this graph based on the annual incidence of each department, so drawing from a different data frame. Is this possible?

The code i'm using does not come up with any errors but doesn't do anything either:

ggplot(data=mean_all_dept,aes(x=Year, y=log10(CI100k)))+

geom_line(color="purple")+

geom_point()+

labs(title = "Mean annual dengue incidence for 33 departments")+

labs(y="Cumulative Incidence (per 100,000)", x="Year")

##line plot of CI for all departments against time with error bars

library(ggplot2)

df <- mean_all_dept

df$Year <- as.factor(df$Year)

head(df)

mean_summary <- function(meanz, CI100k, Year, Dept) {

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(meanz, Year, Dept, .fun=summary_func,

CI100k)

data_sum <- rename(data_sum, c("Mean" = CI100k))

return(data_sum)

}

##Summarise the data

df3 <- mean_summary(meanz, varname="CI100k",

groupnames=c("Dept", "Year"))

head(df3)

p <- ggplot(df3, aes(x=Year, y=log10(CI100k)))+

geom_line(color="purple")+

geom_point+

geom_errorbar(aes(ymin=CI100k-sd, ymax=CI100k+sd), width=.2,

position=position_dodge(0.05))

p+labs(title = "Mean annual dengue incidence for 33 departments", y="Cumulative Incidence (per 100,000)", x="Year")

Thank you so much for any help!!