test(){ covariate(Cp) deriv(Ce = Ke0 * (Cp - Ce)) E = Emax * (Ce^Gam / (EC50^Gam + Ce^Gam)) error(EEps = 1) observe(EObs = E + EEps) fixef(Ke0 = c(, 1, )) fixef(EC50 = c(, 1000, )) fixef(Emax = c(, 100, )) fixef(Gam = c(,2, )) }