例如以下一节,度数为2的一次样条曲线:
library(splines)
library(ISLR)
fit.spline <- lm(wage~bs(age, knots=c(42), degree=2), data=Wage)
summary(fit.spline)
我看到了一些出乎我的意料的估算。
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 57.349 3.950 14.518 < 2e-16 ***
bs(age, knots = c(42), degree = 2)1 59.511 5.786 10.285 < 2e-16 ***
bs(age, knots = c(42), degree = 2)2 65.722 4.076 16.122 < 2e-16 ***
bs(age, knots = c(42), degree = 2)3 37.170 9.722 3.823 0.000134 ***
有没有一种方法可以提取节点之前和之后的二次模型(以及其系数)? 也就是说,如何提取age = 42
切点之前和之后的两个二次模型?
使用summary(fit.spline)
可以得到系数,但是(据我所知),它们对于解释并不具有意义。