Boxplot em relação a dois fatores usando ggplot2 em R

Eu sou muito novo no R e em todos os pacotes no R. Olhei para a documentação do ggplot2, mas não consegui encontrar isso. Eu quero um gráfico de caixa de variável boxthisem relação a dois fatores f1e f2. Suponha que ambos f1e f2são variáveis ​​fatoriais e cada um deles use dois valores e boxthisé uma variável contínua. Eu quero ficar 4 boxplots em um gráfico, cada um correspondendo a uma combinação das possíveis combinações que f1e f2pode tomar. Eu acho que usando a funcionalidade básica em R, isso pode ser feito por

> boxplot(boxthis ~ f1 * f2 , data = datasetname) 

Agradecemos antecipadamente por qualquer ajuda.

Posso pensar em duas maneiras de conseguir isso:

1. Crie todas as combinações de f1e f2fora da ggplotfunção

library(ggplot2)

df <- data.frame(f1=factor(rbinom(100, 1, 0.45), label=c("m","w")), 
                 f2=factor(rbinom(100, 1, 0.45), label=c("young","old")),
                 boxthis=rnorm(100))

df$f1f2 <- interaction(df$f1, df$f2)

ggplot(aes(y = boxthis, x = f1f2), data = df) + geom_boxplot()

2. use cor / preenchimento / etc.

ggplot(aes(y = boxthis, x = f2, fill = f1), data = df) + geom_boxplot()

Aqui está outra maneira:

ggplot(datasetname) + boxplot(aes(x=as.factor(f1),y=boxthis)) + facet_wrap(~as.factor(f2), nrow=5)

Fornece algo como o seguinte: