Hi Pharsight team,
I've come across a new one for me. I'm attempting to run a target mediated drug disposition model and originally ran the FO method and then ran the ELS. The ELS took a long time to run but after about 20 hours convergence resulted in a plausible estimates for theta. however since standard errors were not computed for the theta's (and we know that publications like to see these in the data we submit) I elected to attempt to run a different algorithmn to see if it made a difference. Unfortunately I received this error:
FATAL ERROR IN QRPEM/IMPEM!
Model not suitable for QRPEM analysis
Possibly non linear covariate model or
some other unimplemented feature
I've also received other errors that simply state the run aborted due to a "Fortran" error...please explain!
The model I am using is not that complicated as far as textual models go. ive modified a TMDD model proposed by Mager and Jusko in 2002 for a mixed absorption and paramaterized it for clearance as opposed to microconstants.
Is there some workarond, or is my model (below) somehow flawed for the QRPEM approach?
thanks
Elliot
test(){
deriv(A1 = - (Cl * C) + (Aa * Ka)- (Cl2 * (C - C2)) + (Alymph * dur)-(kon*C*Rf)+(koff*DR) )
deriv(Aa = - (Aa * Ka))
deriv(A2 = (Cl2 * (C - C2)))
deriv(Alymph = - (Alymph * dur))
deriv(DR=(C*kon*Rf)-(koff+kint)*DR)
Rf=Rmax-DR
C = A1 / V
dosepoint(Aa, bioavail = (F), idosevar = AaDose, infdosevar = AaInfDose, infratevar = AaInfRate)
C2 = A2 / V2
error(CEps = 0.50)
observe(CObs = C * (1 + CEps))
dosepoint(Alymph, tlag = (Tlag2), bioavail = (1-F), duration = (dur), idosevar = AlymphDose, infdosevar = AlymphInfDose, infratevar = AlymphInfRate)
stparm(V = tvV * exp(nV))
stparm(Cl = tvCl * exp(nCl))
stparm(Ka = tvKa * exp(nKa))
stparm(V2 = tvV2 * exp(nV2))
stparm(Cl2 = tvCl2 * exp(nCl2))
stparm(F = ilogit(tvF + nF))
stparm(dur = tvdur * exp(ndur))
stparm(Tlag2 = tvTlag2 * exp(nTlag2))
stparm(kint = tvkint*exp(nkint))
stparm(kon = tvkon*exp(nkon))
stparm(koff = tvkoff*exp(nkoff))
stparm(Rmax = tvRmax*exp(nRmax))