0

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 番目のモデルは切片を正しく識別できません。ここで何が欠けていますか?

4

2 に答える 2