test <- function(time, parms) {
with(as.list(parms), {
y <- (K * y0) / (y0 + (K - y0) * exp(-mumax * time)) + y_shift
return(as.matrix(data.frame(time=time, y=y)))
})
}
mygrowthmodel <- growthmodel(test, c("y0", "mumax", "K", "y_shift"))
Run the code above in your browser using DataLab