A R script to illustrate how to make box, bean, pirate, etc, plots

boxbeanpirateplots.R

# Some R code to demo boxplots, beanplots, pirateplots, etc

library(psych)
library(dplyr)
library(tidyr)
library(beanplot)
library(yarrr)

# Load the cushny data
data(cushny)

# Plot it using error.bars. Following example in psych package manual
error.bars(cushny[1:4],
           within=TRUE,
           ylab="Hours of sleep",
           xlab="Drug condition", 
           main="95% confidence of within subject effects")

# Reformat the data into a "long" and "tidy" format
Df <- select(cushny, Control, drug1, drug2L, drug2R) %>% 
  mutate(subject = factor(seq(length(Control)))) %>% 
  select(subject, everything()) %>% 
  gather(drug, sleep, Control:drug2R) %>% 
  arrange(subject)

# A box plot
boxplot(sleep ~ drug, 
        data=Df, 
        ylab='Hours of sleep', 
        xlab='Drug')

# Beanplot 
beanplot(sleep ~ drug, 
         data=Df, 
         ylab='Hours of sleep', 
         xlab='Drug',
         col='lightgray', 
         overallline = 'median')

# Pirateplot
pirateplot(formula = sleep ~ drug,
           data = Df,
           theme = 3,
           point.o = 1,
           ylab='Hours of sleep', 
           xlab='Drug')