Как создать зависимую от времени зависимость от использования в `coxph` модели в пакете выживания в r? - PullRequest
Новая
       17

Как создать зависимую от времени зависимость от использования в `coxph` модели в пакете выживания в r?

0 голосов
/ 04 марта 2019

У меня есть такой набор данных: 

  id       status  time event   
  <chr>     <dbl> <dbl> <time>  
1 22016890      1   459 425 days
2 22042190      1    83  NA days
3 18108887      1   373  NA days
4 17706799      1    33  NA days
5 17039253      1   531  NA days
6 23223163      1   230  NA days
7 15848065      1     1  NA days

, и мне нужно создать изменяющуюся по времени версию event.Как мне это сделать?

Я пытался последовать примеру 3.2 из этой виньетки , используя survival::tmerge, но мне не удалось.

Набор данных: 

nd <- structure(list(id = c("22016890", "22042190", "18108887", "17706799", 
"17039253", "23223163", "15848065"), status = c(1, 1, 1, 1, 1, 
1, 1), time = c(459, 83, 373, 33, 531, 230, 1), event = structure(c(425, 
NA, NA, NA, NA, NA, NA), class = "difftime", units = "days")), row.names = c(NA, 
-7L), class = c("tbl_df", "tbl", "data.frame"))

Как мне нужны мои данные: 

# A tibble: 8 x 5
  id       start   end status event
  <chr>    <dbl> <dbl>  <dbl> <dbl>
1 22016890     0   425      1     0
2 22016890   425   459      1     1
3 22042190     0    83      1     0
4 18108887     0   373      1     0
5 17706799     0    33      1     0
6 17039253     0   531      1     0
7 23223163     0   230      1     0
8 15848065     0     1      1     0

1 Ответ

0 голосов
/ 06 марта 2019

Мне удалось сделать это, используя: 

library(survival)
simpledata <- tmerge(data1 = nd[, c(1:3)], data2 = nd, id = id, tstop = time)
simpledata <- tmerge(simpledata, nd, id=id, infect = tdc(event))
simpledata
        id status time tstart tstop infect
1 22016890      1  459      0   425      0
2 22016890      1  459    425   459      1
3 22042190      1   83      0    83      0
4 18108887      1  373      0   373      0
5 17706799      1   33      0    33      0
6 17039253      1  531      0   531      0
7 23223163      1  230      0   230      0
8 15848065      1    1      0     1      0
...