test(){
covariate(dose_grp) # taking amount of dose from input
D = dose_grp*exp(-Ke*t) # 1st order elimination of amount of dose
deriv(N = Kg * N - D * Kk * N) # Turnover model: linear stimulation of loss
sequence{N = 10000} # initialization of bacterial growth
error(NEps = 0.1) # initial estimate of residual error
observe(NObs = N + N*NEps) # multiplicative residual error model
fixef(Ke = c(, 1, )) # initial estimate of elimination rate constant
fixef(Kg = c(, 0.2, )) # ditto for growth rate constant
fixef(Kk = c(, 1, )) # ditto for kill rate constant
}