R/bootfit.R

In RemkoDuursma/fitplc: Fit Hydraulic Vulnerability Curves

# This works for this object only
#' @importFrom stats residuals
#' @importFrom stats update
#' @importFrom stats nls
#' @importFrom stats predict
#' @importFrom stats quantile
bootfit <- function(fit, n=999, maxnbad=50, Data, startList, weights=NULL){
  

  
  nrow <- length(residuals(fit))
  
  refits <- list()
  ndone <- 0
  nbad <- 0
  
  while(ndone < n && nbad < maxnbad){
    bootsample <- sample(nrow, replace=TRUE)
    
    tryfit <- try(update(fit, subset=bootsample, data=Data, start=startList, weights=weights), silent=TRUE)

    if(!inherits(tryfit, "try-error")){
      ndone <- ndone + 1
      nbad <- 0
      refits[[ndone]] <- tryfit  
    } else {
      nbad <- nbad + 1
    }
  
  if(nbad == maxnbad)warning("Too many consecutive fits have failed.")
  }
  
  
  coefs <- t(sapply(refits,coef))
  return(coefs)
}