cph之后以R rms封装限制三次样条输出

Question

我正在R中开发COX回归模型。
我当前使用的模型如下

fh <- cph(S ~ rcs(MPV,4) + rcs(age,3) + BMI + smoking + hyperten + gender +
    rcs(FVCPP,3) + TLcoPP, x=TRUE, y=TRUE, surv=TRUE, time.inc=2*52) 

如果我想和这个人一起看

print(fh, latex = TRUE) 

对于coefs/SE/Wald，我得到3 MPV (MVP, MVP' and MVP'')，对于age (age, age')，我得到2。
有人可以向我解释这些输出是什么吗？即，我相信它们与我添加的受限三次样条曲线有关。

Answer 1

编写rcs(MPV,4)时，您定义了样条中要使用的结数。在这种情况下为4。类似地，rcs(age,3)定义了一个3节的样条。由于可识别性的限制，每个样条线会减去1个结。您可以认为这是为每个样条线定义一个截距。因此rcs(Age,3)是2个非线性基函数和一个截距的线性组合，而rcs(MPV,4)是3个非线性基函数和一个截距的线性组合，即，

和

在上述符号中，从打印语句中得到的是回归系数和，以及相应的标准误差，p值等。截距和通常设置为零，但是它们很重要，因为如果没有这些参数，模型拟合例程将不知道在y轴上的哪个位置约束样条曲线。

最后一点，您实际上可能对summary(fh)的输出更感兴趣。