Calculating Treatment Means with Error Bars and p-Values in R Using ggplot2

Resolving Incorrect Results with ggplot2's scale_apply Function: A Known Issue and Possible Solutions