splines::ns()
Rの関数に問題があるようです。
簡単なダミー問題を作成しました
dat <- data.frame(t <- seq(0, 6, .01),
x <- rnorm(length(t), sd = 1),
y <- 5 + t - x^2 + rnorm(length(t), sd = .33))
lm(y ~ t + I(x^2), data = dat)
library(splines)
lm(y ~ t + ns(x, knots = c(0), Boundary.knots = c(-3, 3)), data = dat)
最初のモデルは正常に機能しますが、2 番目のモデルは切片を正しく識別できません。ここで何が欠けていますか?