R- CRR模型(CMPRSK)的限制立方样条

Question

我正在尝试使用R中的CRR函数（CMPRSK）在R中进行竞争风险存活模型，并且通过初步分析，我想通过限制的立方样条来改变我的两个连续变量转型。有人知道一种类似于RMS软件包中RCS函数的方法吗？

Answer 1

生成model.matrix

时，您可以使用splines软件包

例如。其中cov3和cov4是您想要转换的连续变量：

covariates <- model.matrix(~cov1 + cov2 + bs(cov3, ...) + ns(cov4, ...), data = df)[,-1]
model <- with(df, crr(ftime, fstatus, covariates, failcode = 1, cencode = 0))

...可以包括bs（B-Spline）或ns（天然立方样条）功能的各种选项。