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?