R/fit_logistic.R

In jeffkimbrel/platereadR: A package for dealing with Biolog and plate read data

#' Fit a single sample with logistic equation
#'
#' @param x Time vector
#' @param y Growth vector
#'
#' @export

fit_logistic = function(x, y) {
  require("growthrates")
  lower = c(y0 = -1, mumax = 0, K = 0)
  upper = c(y0 = 1,  mumax = 10, K = 10)
  p <- c(y0 = min(y), mumax = 0.02, K = max(y))

  d = fit_growthmodel(FUN = grow_logistic, p = p, time = x, y = y, lower = lower, upper = upper)
  return(d)
}