This section introduces the functionality of the base-R function boxplot. Note that for some data formats, the plot function with x being a factor variable will also work.
boxplot(len ~ supp*dose, data = ToothGrowth, col = c("white", "steelblue"), frame = FALSE)