如何在ggplot中按Tukey显著性字母匹配箱线图与标签颜色？

我想在ggplot2或ggboxplot中，让箱线图和对应的Tukey显著性字母标签使用相同的颜色，但不知道怎么自动或更优雅地实现（比如用terrain.colors这类调色函数）。目前只能手动设置颜色来达到预期效果：箱线图和对应Tukey显著性字母标签颜色一致。以下是基于Stack Overflow已有答案编写的代码，希望能得到帮助：

library(plyr)
library(ggplot2)
library(multcompView)

set.seed(0)
lev <- gl(3, 10)
y <- c(rnorm(10), rnorm(10) + 0.1, rnorm(10) + 3)
d <- data.frame(lev=lev, y=y)

a <- aov(y~lev, data=d)
tHSD <- TukeyHSD(a, ordered = FALSE, conf.level = 0.95)

generate_label_df <- function(HSD, flev){
  # Extract labels and factor levels from Tukey post-hoc 
  Tukey.levels <- HSD[[flev]][,4]
  Tukey.labels <- multcompLetters(Tukey.levels)['Letters']
  plot.labels <- names(Tukey.labels[['Letters']])
  
  # Get highest quantile for Tukey's 5 number summary and add a bit of space to buffer between    
  # upper quantile and label placement
  boxplot.df <- ddply(d, flev, function (x) max(fivenum(x$y)) + 0.2)
  
  # Create a data frame out of the factor levels and Tukey's homogenous group letters
  plot.levels <- data.frame(plot.labels, labels = Tukey.labels[['Letters']],
                            stringsAsFactors = FALSE)
  
  # Merge it with the labels
  labels.df <- merge(plot.levels, boxplot.df, by.x = 'plot.labels', by.y = flev, sort = FALSE)
  
  return(labels.df)
}

#Generate ggplot

ggplot(d, aes(x=lev, y=y)) + geom_boxplot(fill = c("green", "green", "orange")) +
   geom_text(data = generate_label_df(tHSD, 'lev'), colour = c("green","orange", "green"), aes(x = plot.labels, y = V1, label = labels )) +
   scale_colour_manual(values=c("green", "green", "orange"))

要实现箱线图与标签颜色的自动匹配，核心是给每个因子水平绑定对应颜色，让填充色和文字色共享同一映射规则，无需手动逐个指定。以下是优化后的实现：

library(plyr)
library(ggplot2)
library(multcompView)

set.seed(0)
lev <- gl(3, 10)
y <- c(rnorm(10), rnorm(10) + 0.1, rnorm(10) + 3)
d <- data.frame(lev=lev, y=y)

# 1. 自动生成颜色映射：按因子水平数量生成对应颜色
factor_levels <- levels(d$lev)
color_map <- terrain.colors(length(factor_levels))
names(color_map) <- factor_levels # 绑定颜色到因子水平

a <- aov(y~lev, data=d)
tHSD <- TukeyHSD(a, ordered = FALSE, conf.level = 0.95)

# 2. 修改标签生成函数，加入颜色信息
generate_label_df <- function(HSD, flev, color_map){
  Tukey.levels <- HSD[[flev]][,4]
  Tukey.labels <- multcompLetters(Tukey.levels)['Letters']
  plot.labels <- names(Tukey.labels[['Letters']])
  
  boxplot.df <- ddply(d, flev, function (x) max(fivenum(x$y)) + 0.2)
  
  # 新增color列，关联对应因子水平的颜色
  plot.levels <- data.frame(plot.labels, 
                            labels = Tukey.labels[['Letters']],
                            color = color_map[plot.labels],
                            stringsAsFactors = FALSE)
  
  labels.df <- merge(plot.levels, boxplot.df, by.x = 'plot.labels', by.y = flev, sort = FALSE)
  return(labels.df)
}

# 3. 生成带颜色信息的标签数据框
labels_df <- generate_label_df(tHSD, 'lev', color_map)

# 4. 绘制图形，自动关联颜色
ggplot(d, aes(x=lev, y=y)) + 
  geom_boxplot(aes(fill = lev)) + # 按因子水平映射填充色
  geom_text(data = labels_df, aes(x = plot.labels, y = V1, label = labels, color = plot.labels)) + # 按因子水平映射文字色
  scale_fill_manual(values = color_map) + # 应用自定义颜色集
  scale_color_manual(values = color_map) + # 文字色与填充色保持一致
  theme_minimal()

关键优化点 #

• 自动调色：用terrain.colors（或viridis、RColorBrewer等调色函数）根据因子水平数量自动生成颜色，无需手动输入颜色值；
• 颜色绑定：将颜色与因子水平一一绑定，确保箱线图和标签的颜色始终对应；
• 扩展性强：后续新增因子水平时，无需修改颜色相关代码，调色函数会自动适配新的水平数量。

内容的提问来源于stack exchange，提问作者Darwin PC