library(SiZer) x <- 1:10 y <- c(100, 200, 300, 350, 400, 500, 600, 800, 900, 600) model.pwl = piecewise.linear(x, y, middle = 1, CI = FALSE, bootstrap.samples = 1000, sig.level = 0.05) plot(model.pwl)