Накопленная заболеваемость в статистике R - PullRequest
0 голосов
/ 19 октября 2019
library(cmprsk)   
library(MASS)
data(Melanoma)
Melanoma$event <- NA
Melanoma$event[Melanoma$status == 2] <- 0
Melanoma$event[Melanoma$status == 1] <- 1
Melanoma$event[Melanoma$status == 3] <- 2

cif=cuminc(ftime=Melanoma$time, fstatus=Melanoma$event,cencode=0)


cifd <- data.frame(
  cause = factor(rep(c(2, 1),c(length(cif[[1]]$time), length(cif[[2]]$time)))),
  time  = c(cif[[1]]$time, cif[[2]]$time),
  cif   = c(cif[[1]]$est,  cif[[2]]$est)) 

этот приведенный выше пример кода оценивает совокупную заболеваемость для конкурирующего риска (cif), а затем создает фрейм данных (cifd) для хранения совокупных значений заболеваемости.

мои вопросы -

1. Является ли риск более полного изучения события 1 = 0,3387175 при наличии конкурирующего события 2? Как указано в cif?

2, почему cif дублируется каждый раз в списке cif [[1]] $ time для события 1 и события 2 в cif [[2]] $ time?

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...