guessCellGrowthParams {cellGrowth} | R Documentation |
Guess initial parameters values for growth models
guessCellGrowthParams(x, z, relative.height.at.lag = 0.1)
x |
|
z |
codenumeric vector: log(OD) |
relative.height.at.lag |
|
The relative.height.at.lag parameter should be close to the relative height of the point, where the curve reaches its maximal slope. If the fitting fails, try to set this parameter to a different value.
A list with entries:
mu |
|
l |
|
z0 |
|
zmax |
|
Julien Gagneur
x <- 1:1000 z <- gompertz(x, mu=0.01, l=200, z0=1, zmax=5)+rnorm(length(x),mean=0,sd=0.25) guess <- guessCellGrowthParams(x,z,relative.height.at.lag=0.5) fit <- nls(z~gompertz(x,mu,l,z0,zmax),start=guess) plot(x,z) lines(x,predict(fit,x),lwd=2,col="red")