test(){
  deriv(Ce = Ke0*(C - Ce))
  cfMicro(A1, Ke)
  dosepoint(A1, = A1InfRate)
  C = A1 / V
  E = E0 - Imax * Ce^Gam / (IC50^Gam + Ce^Gam)
  error(EEps = 1)
  observe(EObs = E + EEps)
  stparm(V = tvV)
  stparm(Ke = tvKe)
  stparm(Ke0 = tvKe0 * exp(nKe0))
  stparm(IC50 = tvIC50 * exp(nIC50))
  stparm(Gam = tvGam * exp(nGam))
  stparm(E0 = tvE0 * exp(nE0))
  stparm(Imax = tvImax * exp(nImax))
  fixef(tvV(freeze) = c(, 1, ))
  fixef(tvKe(freeze) = c(, 1, ))
  fixef(tvKe0 = c(, 1, ))
  fixef(tvIC50 = c(, 1, ))
  fixef(tvGam = c(, 1, ))
  fixef(tvE0 = c(, 1, ))
  fixef(tvImax = c(, 1, ))
  secondary(Cl = tvKe*tvV)
  secondary(Ke_hl = log(2)/tvKe)
  ranef(diag(nKe0, nIC50, nGam, nE0, nImax) = c(1, 1, 1, 1, 1))
}
