Dear R-helpers, Why the following calls result in a different output: library(splines) library(Design) x <- rnorm(10) ns(x, knots = 0, Boundary.knots = c(-1, 1)) rcs(x, c(-1, 0, 1)) Thanks in advance for any hints. Kind Regards, Moko Kosta