guessCellGrowthParams: Guess growth models parameters

Description Usage Arguments Details Value Author(s) Examples

View source: R/guessCellGrowthParams.R

Description

Guess initial parameters values for growth models

Usage

1
  guessCellGrowthParams(x, z, relative.height.at.lag = 0.1)

Arguments

x

numeric vector: time points

z

codenumeric vector: log(OD)

relative.height.at.lag

numeric scalar (see Details)

Details

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.

Value

A list with entries:

mu

numeric scalar: maximal growth rate parameter

l

numeric scalar: time lag parameter

z0

numeric scalar: minimal log(OD) parameter

zmax

numeric scalar: maximal log(OD) parameter

Author(s)

Julien Gagneur

Examples

1
2
3
4
5
6
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")

cellGrowth documentation built on Oct. 31, 2019, 8:38 a.m.