可以用类似这样的语句:
m.bsp1 <- lm(y ~ bs(x, df = 6, knots = c(-2.5, -1, 2)))
这里调用的bs()是splines包里的函数,多项式样条B样条基。
函数形式是——
bs(x, df = NULL, knots = NULL, degree = 3, intercept = FALSE, Boundary.knots = range(x))
有一个文档写的特别详细,可以参考一下:
http://web.as.uky.edu/statistics/users/pbreheny/764-F11/notes/9-27.pdf