比如下面的一个结,二度,样条:
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)产量系数,但(据我的理解)它们对解释没有意义.