Displaying Mean Value Labels on Boxplots: A Comprehensive Implementation Using R and ggplot2

Introduction and Problem Background

In data visualization, boxplots are a common chart type for displaying data distribution characteristics, intuitively showing median, quartiles, and outliers. However, in practical analysis, researchers often need to annotate the mean values of each group on boxplots for more precise comparison of central tendencies. Based on a typical Q&A from Stack Overflow, this article systematically explores how to effectively achieve this in the R environment, particularly using the ggplot2 package.

Data Preparation and Basic Visualization

We use the built-in R dataset PlantGrowth as an example, which contains weight measurements of plants under different treatment groups. First, load the necessary library and examine the data structure:

library(ggplot2)
data(PlantGrowth)
str(PlantGrowth)

Create a basic boxplot using ggplot2 and add mean points via the stat_summary function:

ggplot(data = PlantGrowth, aes(x = group, y = weight, fill = group)) + 
  geom_boxplot() + 
  stat_summary(fun = mean, colour = "darkred", geom = "point", 
               shape = 18, size = 3, show.legend = FALSE)

This code generates a visualization with three boxplots, where means are marked with red diamond points, but specific numerical values are not yet displayed.

Method 1: Calculating Means with Aggregate and Adding Labels

This is the method recommended by the accepted best answer. First, use the aggregate function to calculate means for each group:

means <- aggregate(weight ~ group, PlantGrowth, mean)
print(means)

After obtaining a data frame containing group and weight (mean), pass it as a new data layer to the geom_text function:

ggplot(data = PlantGrowth, aes(x = group, y = weight, fill = group)) + 
  geom_boxplot() + 
  stat_summary(fun = mean, colour = "darkred", geom = "point", 
               shape = 18, size = 3, show.legend = FALSE) + 
  geom_text(data = means, aes(label = round(weight, 2), y = weight + 0.08), 
            size = 4, colour = "blue")

Here, y = weight + 0.08 uses vertical offset to prevent label overlap with mean points, and round(weight, 2) formats the value to two decimal places. The core advantage of this method is high flexibility, allowing fine control over labels, such as adjusting color, size, and position.

Method 2: Directly Outputting Text with stat_summary

As a supplementary reference, another method involves using the stat_summary function to directly generate text labels without pre-calculating means:

ggplot(data = PlantGrowth, aes(x = group, y = weight, fill = group)) + 
  geom_boxplot() + 
  stat_summary(fun = mean, colour = "darkred", geom = "point", 
               shape = 18, size = 3, show.legend = FALSE) + 
  stat_summary(fun = mean, geom = "text", show.legend = FALSE, 
               aes(label = round(..y.., 1)), vjust = -0.7, colour = "red")

Here, ..y.. is an internal ggplot2 variable representing the computed mean; vjust = -0.7 controls vertical alignment. This method offers more concise code but has relatively limited customization options, making it suitable for rapid prototyping.

Technical Details and Optimization Suggestions

In practical applications, several points should be noted: First, label positions should be dynamically adjusted based on data range, e.g., using y = weight + 0.05 * max(PlantGrowth$weight) instead of fixed offsets. Second, if the data contains many groups, consider using geom_label instead of geom_text to enhance readability. Additionally, for non-normally distributed data, means may be affected by outliers, so combining median annotations can provide more comprehensive information.

Extended Applications and Conclusion

The methods introduced in this article are not limited to boxplots but can be extended to other statistical graphics, such as violin plots or dot plots. By combining data preprocessing with the dplyr package or formatting labels using scale functions, visualization effectiveness can be further improved. In summary, within the ggplot2 framework, through layer superposition and statistical transformations, complex data annotation needs can be efficiently met, providing strong support for scientific research and business analysis.