Boxplot with quantile whiskers

myBoxplot.R

myboxplot <- function(datalist, quantiles = c(.1, .9), names, ...){
  
  
  bp <- boxplot(datalist, plot = FALSE)
  
  quants <- sapply(datalist, function(df){
    
    quantile(df, quantiles, na.rm = T)
    
  })
  
  bp$stats[1,] <- quants[1,]
  bp$stats[5,] <- quants[2,]
  bp$names <- names
  
  bxp(bp, ...)
    
}