ASA Connect

 View Only

R survivalROC function with TP's > 1

  • 1.  R survivalROC function with TP's > 1

    Posted 05-22-2017 11:31

    Dear All,

    While using the R survivalROC function on my data I got TP's that were greater than 1. Please advise where the problem may lie.

    I attached my printout for your convenience:

     

    $cut.values

    [1]       -Inf -2.0392992 -1.6639958 -1.4763441 -1.2886924 -1.1010407 -0.9133890 -0.7257373 -0.5380856 -0.3504339

    [11] -0.1627822  0.0248695  0.2125212  0.4001729  0.5878246  0.7754763  0.9631280  1.1507797  1.3384314  1.5260831

    [21]  1.7137348  1.9013865  2.0890382  2.8396451  3.2149485

     

    $TP

    [1] 1.00000000 1.00373900 0.99018323 0.99829624 0.98992139 1.00005297 1.01080076 1.00710510 1.01966745 1.01071444

    [11] 0.98560138 0.94548627 0.91418073 0.75007299 0.70314998 0.62260339 0.50940277 0.45280246 0.39620216 0.33960185

    [21] 0.16980092 0.11320062 0.05660031 0.00000000 0.00000000

     

    $FP

    [1] 1.00000000 0.98368249 0.91081649 0.87800974 0.81904904 0.75508366 0.70625802 0.66133830 0.61202195 0.52260485

    [11] 0.33041327 0.28003606 0.15074461 0.08797951 0.07005610 0.06122550 0.06122550 0.06122550 0.04591913 0.03061275

    [21] 0.03061275 0.03061275 0.01530638 0.01530638 0.00000000

     

    $predict.time

    [1] 36

     

    $Survival

    [1] 0.7871356

     

    $AUC

    [1] 0.9242652

    Thank You,

    Benjy



    ------------------------------
    Benjamin A. Levinson, PhD
    NYU School of Medicine
    ------------------------------