Posthoc labels on anova boxplot in R

Если у меня есть некоторые данные и я провожу ANOVA и пост-хок тесты, как мне сделать боксплот, который добавляет пост-хок классификацию автоматически, а не редактирует рисунок вне R?

Например, вот некоторые данные для начала:

install.packages("reshape", dependencies=T)
library(reshape)

x <- rnorm(30)
y <- rnorm(30)+1
z <- rnorm(30)+0.5

data.1 <- data.frame(x, y, z)
data.2 <- melt(data.1)

Вот код для запуска простого одностороннего ANOVA и всех незапланированных пост-хок тестов сравнения:

linear.model <- lm(value~variable, data=data.2)
anova(linear.model)

# Analysis of Variance Table
# Response: value
#           Df Sum Sq Mean Sq F value   Pr(>F)   
# variable   2 10.942  5.4710  5.8628 0.004087 **
# Residuals 87 81.185  0.9332     

TukeyHSD(aov(linear.model))

# Tukey multiple comparisons of means
# 95% family-wise confidence level
# Fit: aov(formula = linear.model)
# $variable
          # diff        lwr        upr     p adj
# y-x  0.8344105  0.2396705 1.42915051 0.0034468
# z-x  0.2593612 -0.3353788 0.85410126 0.5539050
# z-y -0.5750493 -1.1697893 0.01969078 0.0602975

На данный момент я хочу отнести x к группе "a", y к группе "b" и z к группе "a,b". Я могу сделать боксплот, но как аннотировать его буквами?

boxplot(value~variable, data=data.2)