crr模型的受限三次样条（cmprsk）

Question

我试图使用R中的crr函数（cmprsk）建立竞争性风险生存模型，并通过初步分析，我想使用受限三次样条变换来变换两个连续变量。有人知道类似于rms软件包中的rcs函数的方法吗？

Answer 1

生成splines时可以使用model.matrix包

例如。其中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样条）或ns（自然三次样条）功能可用的各种选项。