1

我正在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) 

我得到 3 个coefs/SE/WaldetcMPV (MVP, MVP' and MVP'')和 2 个age (age, age').
有人可以向我解释这些输出是什么吗?即我相信它们与我添加的受限三次样条有关。

4

1 回答 1

0

当您编写 时rcs(MPV,4),您定义要在样条曲线中使用的结数;在本例中为 4。类似地,rcs(age,3)定义具有 3 个结的样条曲线。由于可识别性限制,从每个样条中减去 1 个结。您可以将其视为为每个样条定义一个截距。所以rcs(Age,3)是 2 个非线性基函数和一个截距rcs(MPV,4)的线性组合,而 3 个非线性基函数和一个截距的线性组合,即

在上面的符号中,您从 print 语句中得到的是回归系数,以及相应的标准误差、p 值等。截距通常设置为零,但它们很重要,因为没有它们,模型拟合例程如何不知道在 y 轴上的什么位置来约束样条曲线。

最后一点,您实际上可能对summary(fh).

于 2019-03-02T14:13:23.107 回答