colors <- c("#F8766D", "#7CAE00", "#00BFC4", "#C77CFF") library(ggplot2) ggplot(df, aes(x = categorical_FOD_FODs, y = Frequency, fill = Condition_FODs)) + geom_bar(stat = "identity", position = position_dodge()) + scale_fill_manual(values = colors, labels = c("Verbal", "Vm", "Manual", "Mv")) + labs( x = "FODs (direction)", y = "Frequency of FODs", title = "FODs values (ascending/descending) in four condition", fill = "Conditions (renamed)" ) + theme_minimal() + theme( plot.title = element_text(size = 14, face = "bold"), axis.title = element_text(size = 12, face = "bold"), axis.text = element_text(size = 10), legend.title = element_text(size = 10), legend.text = element_text(size = 10), panel.border = element_rect(color = "black", fill = NA, size = 1) ) I want you to improve this plot to show the result more clear than now. So please add error bars or confidence intervals to show the variability in the data.

# I will add the confidence intervals in my plot to show the variability in the data.

ggplot(df, aes(x = categorical_FOD_FODs, y = Frequency_FODs, fill = Condition_FODs)) +   
  geom_bar(stat = "identity", position = position_dodge()) + 
  geom_errorbar(aes(ymin = Frequency_FODs - sd, ymax = Frequency_FODs + sd), position = position_dodge(.9), width = .25) +  
  scale_fill_manual(values = colors, 
                    labels = c("Verbal", "Vm", "Manual", "Mv")) +   
  labs(     
    x = "FODs (direction)",     
    y = "Frequency of FODs",  
    title = "FODs values (ascending/descending) in four condition",
    fill = "Conditions (renamed)"   
  ) +   
  theme_minimal() +
  theme(