Dear all,
I am trying to write a PK model for inverse gaussian absorption following subcutaneous administration of a drug (see Csajka C, Drover D, Verotta D. The use of a sum of inverse Gaussian functions to describe the absorption profile of drugs exhibiting complex absorption. Pharm Res. 2005;22(8):1227-35). The model text is below. When running the model, I consistently get an error "Model execution failed. Initial fcntheta1 value NaN". I tried different initial estimates and different algorithms without success. Any suggestion?
Thank you,
Bruno
test(){
input = (MAT/(6.283*CV^2*t^3))^0.5*exp(-((t-MAT)^2)/(2*CV^2*MAT*t))
deriv(Aa = - input * Aa)
deriv(A1 = input * Aa - Cl * C - Cl2 * (C - C2))
deriv(A2 = Cl2 * (C - C2))
dosepoint(Aa)
C = A1 / V
C2 = A2 / V2
error(CEps = 0.1)
observe(CObs = C * (1 + CEps))
stparm(MAT = tvMAT * exp(nMAT))
stparm(CV = tvCV * exp(nCV))
stparm(V = tvV * exp(nV))
stparm(V2 = tvV2 * exp(nV2))
stparm(Cl = tvCl * exp(nCl))
stparm(Cl2 = tvCl2 * exp(nCl2))
fixef(tvMAT = c(, 1, ))
fixef(tvCV = c(, 1, ))
fixef(tvV = c(, 1, ))
fixef(tvV2 = c(, 1, ))
fixef(tvCl = c(, 1, ))
fixef(tvCl2 = c(, 1, ))
ranef(diag(nV, nCl, nMAT, nCV, nV2, nCl2) = c(1, 1, 1, 1, 1, 1))
}
override test(){
fixef(tvMAT = c(,5,))
fixef(tvCV = c(,0.1,))
fixef(tvV = c(,2876.17166666667,))
fixef(tvV2 = c(,11299.7133333333,))
fixef(tvCl = c(,19.7714833333333,))
fixef(tvCl2 = c(,80.6368833333333,))
}