Экспортируйте подогнанные сплайны регрессии (построенные 'bs' или 'ns') как кусочно-полиномы

Question

Возьмем, например, следующий одноузловой, второй степени, сплайн:

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) дает коэффициенты, но (насколько я понимаю) они не имеют смысла для интерпретации.

Answer 1

Меня постоянно просили заключить идею в моем первоначальном ответе в удобную для пользователя функцию, способную репараметризовать подобранную линейную или обобщенную линейную модель с bs или же ns срок. В конце концов я выкатил небольшой пакет R SplinesUtils по адресу https://github.com/ZheyuanLi/SplinesUtils (с руководством по пакету PDF-версии). Вы можете установить его через

## make you have `devtools` package avaiable
devtools::install_github("ZheyuanLi/SplinesUtils")

Функция для использования здесь RegBsplineAsPiecePoly,

library(SplinesUtils)

library(splines)
library(ISLR)
fit.spline <- lm(wage ~ bs(age, knots=c(42), degree=2), data = Wage)

ans1 <- RegBsplineAsPiecePoly(fit.spline, "bs(age, knots = c(42), degree = 2)")
ans1
#2 piecewise polynomials of degree 2 are constructed!
#Use 'summary' to export all of them.
#The first 2 are printed below.
#8.2e-15 + 4.96 * (x - 18) + 0.0991 * (x - 18) ^ 2
#61.9 + 0.2 * (x - 42) + 0.0224 * (x - 42) ^ 2

## coefficients as a matrix
ans1$PiecePoly$coef
#              [,1]        [,2]
#[1,]  8.204641e-15 61.91542748
#[2,]  4.959286e+00  0.20033307
#[3,] -9.914485e-02 -0.02240887

## knots
ans1$knots
#[1] 18 42 80

По умолчанию функция параметризована кусочно-полиномами в сдвинутом виде (см. ?PiecePoly). Вы можете установить shift = FALSE для несмещенной версии.

ans2 <- RegBsplineAsPiecePoly(fit.spline, "bs(age, knots = c(42), degree = 2)",
                              shift = FALSE)
ans2
#2 piecewise polynomials of degree 2 are constructed!
#Use 'summary' to export all of them.
#The first 2 are printed below.
#-121 + 8.53 * x + 0.0991 * x ^ 2
#14 + 2.08 * x + 0.0224 * x ^ 2

## coefficients as a matrix
ans2$PiecePoly$coef
#              [,1]        [,2]
#[1,] -121.39007747 13.97219046
#[2,]    8.52850050  2.08267822
#[3,]   -0.09914485 -0.02240887

Вы можете предсказать сплайны с predict,

xg <- 18:80
yg1 <- predict(ans1, xg)  ## use shifted form
yg2 <- predict(ans2, xg)  ## use non-shifted form
all.equal(yg1, yg2)
#[1] TRUE

Но так как в модели имеется перехват, прогнозируемые значения будут отличаться от прогнозирования модели по перехвату.

yh <- predict(fit.spline, data.frame(age = xg))
intercept <- coef(fit.spline)[[1]]
all.equal(yh, yg1 + intercept, check.attributes = FALSE)
#[1] TRUE

Пакет имеет summary, print, plot, predict а также solve методы для класса "PiecePoly". Исследуйте пакет для больше.