Diagrama de caja con respecto a dos factores usando ggplot2 en R

Soy muy nuevo en R y en cualquier paquete en R. Miré la documentación de ggplot2 pero no pude encontrar esto. Quiero un diagrama de caja de variable boxthiscon respecto a dos factores f1y f2. Que se supone que tanto f1y f2son variables de factor y cada uno de ellos tiene dos valores y boxthises una variable continua. Quiero conseguir 4 diagramas de caja en un gráfico, cada uno correspondiente a una combinación de las posibles combinaciones que f1y f2puede tomar. Creo que usando la funcionalidad básica en R, esto puede hacerse por

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

Gracias de antemano por cualquier ayuda.

Se me ocurren dos formas de lograr esto:

1. Cree todas las combinaciones de f1y f2fuera de la ggplotfunción

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 color / relleno / etc.

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

Aquí hay otra forma:

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

Da algo como lo siguiente: