test(){ covariate(C) deriv(Ce = Ke0*(C - Ce)) #cfMicro(A1, Cl / V) #dosepoint(A1) #C = A1 / V E = E0 * (1 - Ce / (IC50 + Ce)) error(EEps = 1) observe(EObs = E + EEps) #stparm(V = tvV) #stparm(Cl = tvCl) stparm(Ke0 = tvKe0 * exp(nKe0)) stparm(IC50 = tvIC50 * exp(nIC50)) stparm(E0 = tvE0 * exp(nE0)) #fixef(tvV(freeze) = c(, 1, )) #fixef(tvCl(freeze) = c(, 1, )) fixef(tvKe0 = c(, 0.04, )) fixef(tvIC50 = c(, 40, )) fixef(tvE0(freeze) = c(, 100, )) ranef(diag(nIC50, nE0, nKe0) = c(1, 1, 1)) }