¿Cómo hacer un boxplot de datos agrupados parecido al de la figura?

Quisiera poder hacer un boxplot como el de la figura en R estadístico preferiblemente. Tengo 5 tratamientos (Variedades de arroz) a las cuales se le midió la altura de 20 plantas de cada variedad cada cierto tiempo. En el ejemplo se especifica el siguiente código pero no lo entiendo.

p <- plot_ly(ggplot2::diamonds, x = ~cut, y = ~price, color = ~clarity, type = "box") %>%
  layout(boxmode = "group")
# Create a shareable link to your chart
# Set up API credentials: https://plot.ly/r/getting-started
chart_link = api_create(p, filename="box/group")
chart_link