ggplot2 :: geom_box（）：如何汇总（x，fill）对？

Question

我可以轻松地从钻石数据集创建分组框图：

dd = diamonds[diamonds$color %in% c("E","J"),]
ggplot(dd, aes(x=cut, y=carat)) + 
    geom_boxplot(aes(fill=color), outlier.shape=NA)

# note: below does *not* work for `geom_text()`
ggplot(dd, aes(x=cut, y=carat, fill=color)) + 
    geom_boxplot(outlier.shape=NA)

如何标记x变量定义的每个fill组？

我正在寻找在stat_summary()对之间使用x的解决方案来计算，例如

下每个红色和绿色对之间的p值

这大致对应this question，但区别在于数据按fill分组，我想标记这些组。