Chciałbym zobaczyć boxplots kombinacji czynników iKazano mi użyć kraty za to. Próbowałem tego i wygląda tak:

 Ale teraz chciałbym również dodać statystyki ANOVA do każdej z grup. Prawdopodobnie statystyki powinny wyświetlać wartość p w każdym panelu (w kolorze białym poniżej np. „Australia”). Jak to zrobić w sieci? Zauważ, że wcale nie nalegam na kraty ...

Przykładowy kod:

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(...)
})

Kod do wykonywania ANOVA dla jednej z grup i dowyodrębnić wartość p czy to jest:

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