improve root search intervals

philchalmers · Dec 4, 2024 · 6459a49 · 6459a49
1 parent eb6dc4d
commit 6459a49
Showing 1 changed file with 4 additions and 6 deletions.
diff --git a/R/PBA.R b/R/PBA.R
@@ -338,13 +338,11 @@ PBA <- function(f, interval, ..., p = .6,
                                                  probs = c(.05, .95)),
                                max.interval=interval,median=med, CI=predCI)
         predCIs_root[1L] <- SimSolveUniroot(SimMod=SimMod, b=predCIs[2L],
-                                            interval=quantile(medhistory[medhistory != 0],
-                                                              probs = c(.05, .95)),
-                                            max.interval=interval,median=med, CI=predCI)[1L]
+                                            interval=c(min(medhistory[medhistory != 0]), predCIs[1]),
+                                            max.interval=c(interval[1], predCIs[1]), median=med, CI=predCI)[1L]
         predCIs_root[2L] <- SimSolveUniroot(SimMod=SimMod, b=predCIs[3L],
-                                            interval=quantile(medhistory[medhistory != 0],
-                                                              probs = c(.05, .95)),
-                                            max.interval=interval,median=med, CI=predCI)[1L]
+                                            interval=c(predCIs[1], max(medhistory[medhistory != 0])),
+                                            max.interval=c(predCIs[1], interval[2]), median=med, CI=predCI)[1L]
     }
     if(verbose)
         cat("\n")