如何将与三次样条拟合的 lm() 的输出转换为回归方程

机器算法验证 r 样条
2022-01-23 09:31:23

我有一些代码和输出,我想构建一个模型。我不知道如何使用此输出构建模型:

 require("splines")
 x   <- c(0.2,   0.23,   0.26,   0.29,   0.33,   0.46,    0.53 )
 y   <- c(0.211, 0.2026, 0.2034, 0.2167, 0.2177, 0.19225, 0.182)
 fit <- lm(y ~ ns(x,3))
 summary(fit)

请注意,ns()生成自然三次样条的 B 样条基矩阵。因此,该模型y针对x使用三个自由度的 B 样条进行回归。这种模型的方程会是什么样子?

1个回答
require(rms)
f <- ols(y ~ rcs(x, 3))  # 2 d.f. for x
Function(f)  # represent fitted function in simplest R form
latex(f)     # typeset algebraic representation of fit

rcs“受限三次样条”是自然样条的另一种表示。