r - 各グループのボックスプロットに追加の統計をプロットする方法は?

Question

因子の組み合わせの箱ひげ図を見たいのですが、そのために格子を使用するように言われました。私はそれを試してみましたが、次のようになります。

しかし、ここで、各グループに ANOVA 統計も追加したいと思います。おそらく、統計は各パネルに p 値を表示する必要があります (たとえば、「オーストラリア」の下の白い部分)。格子でこれを行う方法は？私は格子をまったく主張していないことに注意してください...

コード例:

set.seed(123)
n <- 300
country <- sample(c("Europe", "Africa", "Asia", "Australia"), n, replace = TRUE)
type <- sample(c("city", "river", "village"), n, replace = TRUE)
month <- sample(c("may", "june", "july"), n, replace = TRUE)
x <- rnorm(n)
df <- data.frame(x, country, type, month)

bwplot(x ~ type|country+month, data = df, panel=function(...) {
    panel.abline(h=0, col="green")
    panel.bwplot(...)
})

グループの 1 つに対して ANOVA を実行し、p 値を抽出するコードは次のとおりです。

model <- aov(x ~ type, data = df[df$country == 'Africa' & df$month == 'may',])
p_value <- summary(model)[[1]][["Pr(>F)"]][2]