- Fixed missing detail about def for count data in documentation for bmdMA

- Fixed bmd.edfct for BC model for cases where maximum is reached for very small dose values.

Author: JensBaalkilde

R/bmd.edfct.R

@@ -213,7 +213,17 @@ bmd.edfct <- function(object){
             expVal <- exp(parmVec[1]*(log(dose)-log(parmVec[4])))
             parmVec[5]*(1+expVal*(1-parmVec[1]))-(parmVec[3]-parmVec[2])*expVal*parmVec[1]/dose
           }
-          maxAt <- uniroot(helpEqn, interval)$root
+          maxAt <- try(uniroot(helpEqn, interval)$root, TRUE)
+          
+          # solution if maxAt fails
+          if(inherits(maxAt, "try-error")){
+            find_maxAt_tries <- 1
+            while(inherits(maxAt, "try-error") & (find_maxAt_tries<10)){
+              interval[1] <- interval[1]/1e2
+              maxAt <- try(uniroot(helpEqn, interval)$root, TRUE)
+              find_maxAt_tries <- find_maxAt_tries+1
+            }
+          }
 
           eqn <- function(dose) {tempVal*(1+exp(parmVec[1]*(log(dose)-log(parmVec[4]))))-(1+parmVec[5]*dose/(parmVec[3]-parmVec[2]))}
           EDp <- uniroot(eqn, lower = maxAt, upper = upper)$root

R/bmdMA.R

@@ -57,10 +57,10 @@
 #' dose-response models), 2 SD for "hybridSD" background and 0.9 for
 #' "hybridpercentile"
 #' @param def character string specifying the definition of the benchmark dose
-#' to use in the calculations. "excess" , "additional" and "point" are for
-#' binomial response whereas "relative", "extra", "added", "hybridExc" (excess
+#' to use in the calculations. "excess", "additional" and "point" are for
+#' binomial response. "relative", "extra", "added", "hybridExc" (excess
 #' hybrid), "hybridAdd" (additional hybrid), and "point" are for continuous
-#' response
+#' response. "relative", "extra", and "point" are for count response data.
 #' 
 #' "excess" - BMR is defined as: BMR = (f(BMD) - p0)/(1 - p0).  Works for
 #' binomial response. BMR should be between 0 and 1.