R/calculate_halflife.R

calculate_halflife <- function(times, y) {
  wmax <- which.max(y)
  wmin <- which(y < (max(y)/2))
  
  #if nothing meets the halflife condition...
  if(length(wmin[wmin > wmax]) == 0)
    return(NA)
  
  times[min(wmin[wmin > wmax])] - times[wmax]
}
wwiecek/httkgui documentation built on May 15, 2019, 6:31 p.m.