Use coord_cartesian() to hide the outlier instead of ylim(), so we don't subset.
Comment fix: exclude, not include, the outlier when calculating y limit.
Plot 4: R code + png