Violin plots

library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.1     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(ggthemes)

control <- data.frame(time = rnorm(10, mean = 7,  sd = 0.6),
                      group = factor("ctrl"))
drug_a  <- data.frame(time = rnorm( 8, mean = 9,  sd = 0.6),
                      group = factor("drgA"))
drug_b  <- data.frame(time = rnorm( 9, mean = 7,  sd = 0.6),
                      group = factor("drgB"))
drug_c  <- data.frame(time = rnorm( 7, mean = 7,  sd = 0.6),
                      group = factor("drgC"))
drug_d  <- data.frame(time = rnorm( 8, mean = 11, sd = 0.6),
                      group = factor("drgD"))
  
d <- rbind(control, drug_a, drug_b, drug_c, drug_d)

plot(time ~ group, data = d)

p <- ggplot(d, aes(x = group, y = time))
# basic violin plot
p + geom_violin()

# add points
p + geom_violin() + geom_jitter(height = 0, width = 0.1)

# no trim
p + geom_violin(trim=FALSE) + geom_jitter(height = 0, width = 0.1)

# add color
p + geom_violin(aes(fill = group), trim=FALSE) + geom_jitter(height = 0, width = 0.1) 

# add economist-style color
p + geom_violin(aes(fill = group), trim=FALSE) + geom_jitter(height = 0, width = 0.1) + scale_fill_economist()

# add boxplot
 p + geom_violin(aes(fill = group), trim=FALSE) + geom_boxplot(width = 0.1) + geom_jitter(height = 0, width = 0.1) + scale_fill_economist()