clip.R

data <- -10:10
lims <- range(data)
x <- sample(data, 50, replace = TRUE)

par(mfrow=c(1, 5))

plot(x, ylim=lims, main="initial data")

# only positive
plot(pmax(x, 0), ylim=lims, main="clip bottom")

# only negative
plot(pmin(x, 0), ylim=lims, main="clip top")

# clip 1-3
plot(pmax(pmin(x,3), 1), ylim=lims, main="nested pmax & pmin")

clip <- function(x, lower, upper) {
  pmax(pmin(x, upper), lower)
}

plot(clip(x, 1, 3), ylim=lims, main="clip function")