r - 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（自然三次样条）函数的各种选项。