Блог

Я пытаюсь создать конкурирующую модель выживания с риском, используя функцию crr (cmprsk) в R, и посредством предварительного анализа я хочу преобразовать две мои непрерывные переменные с помощью ограниченного преобразования кубического сплайна. Кто-нибудь знает способ сделать это, похожий на функцию rcs в пакете rms?

Вы можете использовать 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 (естественный кубический сплайн).