我可以在R studio的Cox回归中编写两个不同的事件时间吗?

时间:2018-09-12 12:58:26

标签: r regression analysis cox-regression

我下面有一个数据表。它包含有关176名特定肿瘤患者的信息。

一些患者复发(q $复发== 1),一些患者因此死于肿瘤(q $ death == 1)。 我知道从诊断之日起计算的复发时间(以月为单位)(q $ time.recur.months)和死亡的时间以月(q $ months.to.death)为单位。

q $ dfs代表经历复发或死亡的患者与未经历复发或死亡的患者:

q$dfs <- 1*I(q$recurrence==1 | q$death.by.meningioma=="yes") 

我想写一个Cox回归,但我想区分死亡和复发-因此要区分事件发生的时间。

总体思路是

cox <- coxph(Surv(TIME TO EVENTS, q$dfs) ~ 1 )

在q $ recurence == 1 的情况下,我希望活动时间为q $ time.recur.months ,但对于q $,我希望q $ months.to.death death.by.meningioma ==“是”

如果q $ recurrence == 1且q $ death.by.meningioma ==“ yes”,则活动时间应为q $ time.recur.months。

可以做到吗?

## My Data
q <- structure(list(recurrence = c(0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 
                               0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 
                               0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 
                               0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 
                               1L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, NA, 0L, 0L, 0L, 0L, 0L, 
                               0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 
                               0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 
                               0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 
                               0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 
                               0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 
                               0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 
                               0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L), time.recur.months = c(NA, NA, 
                                                                                      NA, NA, NA, NA, 92L, NA, NA, NA, 74L, NA, NA, NA, 2L, 8L, NA, 
                                                                                      NA, NA, NA, 58L, NA, NA, NA, NA, NA, 3L, NA, 4L, NA, NA, NA, 
                                                                                      NA, NA, NA, NA, NA, NA, NA, 39L, NA, NA, NA, NA, 15L, NA, NA, 
                                                                                      NA, NA, NA, NA, NA, NA, NA, NA, 12L, 56L, 57L, NA, NA, 49L, 17L, 
                                                                                      NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 5L, 
                                                                                      NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 9L, NA, 
                                                                                      89L, NA, NA, NA, 8L, 6L, 8L, 4L, NA, NA, NA, NA, NA, NA, NA, 
                                                                                      NA, NA, NA, NA, NA, 60L, NA, NA, 38L, NA, NA, NA, NA, NA, 90L, 
                                                                                      NA, 58L, 54L, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                      NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                      NA, NA, NA, NA, NA, NA, NA, NA, NA, 53L, NA, NA, 124L, NA, NA, 
                                                                                      NA, NA, NA, NA, 7L, NA), death = structure(c(1L, 1L, 2L, 1L, 
                                                                                                                                   1L, 1L, 1L, 3L, 2L, 3L, 1L, 1L, 1L, 1L, 3L, 3L, 1L, 1L, 2L, 3L, 
                                                                                                                                   2L, 1L, 1L, 1L, 2L, 1L, 3L, 1L, 3L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 
                                                                                                                                   1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 
                                                                                                                                   1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 
                                                                                                                                   1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 
                                                                                                                                   1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 3L, 3L, 
                                                                                                                                   3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
                                                                                                                                   1L, 1L, 1L, 3L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
                                                                                                                                   1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
                                                                                                                                   1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L, 1L, 
                                                                                                                                   1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 3L, 2L), .Label = c("", 
                                                                                                                                                                                               "no", "yes"), class = "factor"), days.to.death = c(NA, NA, 2143L, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, 21L, 1279L, 18L, NA, NA, NA, NA, 349L, 627L, 
                                                                                                                                                                                                                                                  NA, NA, 142L, 21L, 3168L, NA, NA, NA, 3598L, NA, 832L, NA, 1795L, 
                                                                                                                                                                                                                                                  NA, 491L, NA, NA, 3309L, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, NA, NA, NA, NA, 2076L, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, 2333L, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, 182L, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 654L, 1057L, 357L, 847L, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                  NA, 1160L, NA, NA, 3022L, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 3L, NA, 
                                                                                                                                                                                                                                                  NA, NA, NA, NA, NA, NA, NA, 2277L, NA, NA, 2679L, 1703L), months.to.death = c(NA, 
                                                                                                                                                                                                                                                                                                                                NA, 178.5833333, NA, NA, NA, NA, 1.75, 106.5833333, 1.5, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, 29.08333333, 52.25, NA, NA, 11.83333333, 1.75, 264, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, 299.8333333, NA, 69.33333333, NA, 149.5833333, NA, 
                                                                                                                                                                                                                                                                                                                                40.91666667, NA, NA, 275.75, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, NA, NA, NA, NA, NA, NA, 173, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, NA, NA, 194.4166667, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, NA, NA, NA, NA, 15.16666667, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 54.5, 
                                                                                                                                                                                                                                                                                                                                88.08333333, 29.75, 70.58333333, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 96.66666667, NA, NA, 
                                                                                                                                                                                                                                                                                                                                251.8333333, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.25, NA, NA, NA, 
                                                                                                                                                                                                                                                                                                                                NA, NA, NA, NA, NA, 189.75, NA, NA, 223.25, 141.9166667)), .Names = c("recurrence", 
                                                                                                                                                                                                                                                                                                                                                                                                      "time.recur.months", "death", "days.to.death", "months.to.death"
                                                                                                                                                                                                                                                                                                                                ), class = "data.frame", row.names = c(NA, -176L))

最佳 基督徒

0 个答案:

没有答案