ste is NA when uprl does not cross HR=1

Oncostat · Jun 9, 2017 · c87265e · c87265e
1 parent eea1255
commit c87265e
Showing 1 changed file with 20 additions and 13 deletions.
diff --git a/R/predict.R b/R/predict.R
@@ -331,19 +331,26 @@ plot.predictSurrosurv <- function(x,
         ...
       )
       if (show.ste) {
-        points(STE[i],
-               0,
-               col = 2,
-               pch = '|',
-               font = 2)
-        text(
-          STE[i],
-          0,
-          col = 2,
-          font = 2,
-          adj = c(1.1, -.5),
-          labels = paste('STE =', round(exp(STE[i]), 2))
-        )
+        if (is.na(STE[i])) {
+          message(paste0(
+            '[', format.methodNames(x)[i], 
+            '] The STE does not exist: the upper limits ',
+            'of the prediction intervals never cross the line HR=1'))
+        } else {
+          points(STE[i],
+                 0,
+                 col = 2,
+                 pch = '|',
+                 font = 2)
+          text(
+            STE[i],
+            0,
+            col = 2,
+            font = 2,
+            adj = c(1.1, -.5),
+            labels = paste('STE =', round(exp(STE[i]), 2))
+          )
+        }
       }
       if (surro.stats) {
         mtext(