Below you can see the code with which the map was created
library(ggplot2)
library(plotly)
data("ToothGrowth")
ToothGrowth$dose <- as.factor(ToothGrowth$dose)
tl_vs_st <- ggplot(aes(x = dose, y = len), data = ToothGrowth)+ facet_grid(~ supp)
tl_vs_st <- tl_vs_st + geom_boxplot(aes(fill = dose)) + theme_bw()
tl_vs_st <- tl_vs_st + ggtitle("Tooth length versus method of supplement")
tl_vs_st <- tl_vs_st + xlab("Supplement type") + ylab("Tooth Length")
tl_vs_st <- tl_vs_st + theme(plot.title = element_text(face = "bold",vjust = 3, hjust = 0.5))
plot <- ggplotly(tl_vs_st)
plot