test(){ deriv(A1 = (-Cl*C -Cld*(C-C2) -(kon*C*R -koff*LR)*V)) deriv(A2 = Cld * (C-C2)) deriv(R = kin -kout*R -kon*C*R + koff*LR) deriv(LR = kon*C*R - koff*LR - keRL*LR) C = A1/V C2 = A2/Vt kin = R0*kout dosepoint(A1) sequence{R = R0} error(CEps = 0.1) #error(REps = 0.1) #error(LREps = 0.1) observe(CObs = C + C*CEps) #observe(RObs = R + R*REps) #observe(LRObs = LR + LR*LREps) fixef(V = c(, 0.05, ), Vt = c(, 0.1, ), Cl = c(, .001, ), Cld = c(, 0.003, ), R0 = c(, 12, )) fixef(kon = c(, 0.091, ), koff = c(, 0.0011, ), kout = c(, 0.0089, ), keRL = c(, .003, )) secondary(Kd = koff / kon, Km = (koff + keRL)/kon) }