ggplot(data, aes(y = y)) + # Create ggplot without outliers geom_boxplot(outlier.shape = NA) + coord_cartesian(ylim = quantile(data$y, c(0.1, 0.9)))