我正在使用" CompetingRiskFrailty"运行 CompetingRiskFrailtySurvfitCreat e功能。在R中打包。这个例子效果很好。但是当我尝试自己的数据时,它就陷入了困境。我不知道它是否仍在运行,需要更多时间还是无限循环?有没有办法找出发生了什么?
提前致谢!!
我的数据中有前20个。我总共获得了3034个障碍
X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 X11 X12 X13 X14 X15
1 101 140 0 0 0 1 1 0 0 1 1 0 0 1 0
2 401 421 0 0 1 1 1 0 0 1 1 0 0 0 0
3 601 519 0 0 0 0 1 0 0 0 1 1 0 0 0
4 602 349 0 0 0 0 0 0 0 1 1 1 0 0 0
5 701 111 0 1 0 0 0 1 NaN 1 0 0 0 0 0
6 701 488 0 1 0 0 0 1 NaN 1 0 0 0 0 0
7 701 488 0 1 0 0 0 1 NaN 1 0 0 0 0 0
8 702 614 0 0 0 0 0 0 0 1 1 0 0 0 0
9 702 614 0 0 0 1 1 1 0 1 0 0 0 0 0
10 1001 362 0 0 1 1 1 0 0 1 1 0 0 0 0
11 1001 362 0 0 1 0 0 0 1 1 1 1 0 0 0
12 1001 362 0 0 1 1 1 0 0 1 1 1 0 0 0
13 1002 204 0 0 1 0 1 1 0 1 0 1 0 0 0
14 1002 204 0 0 1 0 1 1 0 1 0 1 0 0 0
15 1002 204 0 0 1 1 1 0 1 1 1 1 0 0 0
16 1002 204 0 0 1 1 1 1 0 1 0 1 0 0 0
17 1101 51 0 0 0 1 1 0 0 1 1 1 0 0 0
18 1201 789 0 0 1 1 1 0 0 1 1 0 0 1 0
19 1202 328 0 0 1 0 0 1 1 1 0 0 0 1 0
20 1202 328 0 0 1 0 0 0 0 0 1 0 0 1 0
以下是我的代码:
ID <- D_u$X1
surv.time <- D_u$X2
status <- cbind(D_u$X3, D_u$X4, D_u$X5)
my.Surv <- CompetingRiskFrailtySurv(ID=ID,surv.time=surv.time,status=status)
my.Survfit <- CompetingRiskFrailtySurvfitCreate(my.Surv ~ X6 + X7 + X8 + X9 + X10 + X11 + X12 + X13 + X14 + X15, data = D_u , risk.names=c("N","N without S","D without N"))