I use the graphical tool to build a model for the simutaneously fitting of both IV and SC data and calculation of PK parameters including bioavailability. Here is the test for the model:
test(){
deriv(A1 = - (Cl * C)- (Cl2 * (C - C2)) + (Aa * Ka))
urinecpt(A0 = (Cl * C))
deriv(A2 = (Cl2 * (C - C2)))
deriv(Aa = - (Aa * Ka))
C = A1 / V
dosepoint(A1, idosevar = A1Dose, infdosevar = A1InfDose, infratevar = A1InfRate)
C2 = A2 / V2
error(CEps = 1)
observe(CObs = C * (1 + CEps))
dosepoint(Aa, bioavail = (F), idosevar = AaDose, infdosevar = AaInfDose, infratevar = AaInfRate)
stparm(V = tvV * exp(nV))
stparm(Cl = tvCl * exp(nCl))
stparm(V2 = tvV2 * exp(nV2))
stparm(Cl2 = tvCl2 * exp(nCl2))
stparm(Ka = tvKa * exp(nKa))
stparm(F = tvF * exp(nF))
fixef(tvV = c(, 2936, ))
fixef(tvCl = c(, 238, ))
fixef(tvV2 = c(, 2852, ))
fixef(tvCl2 = c(, 859, ))
fixef(tvKa = c(, 1.28, ))
fixef(tvF = c(, 0.77, ))
ranef(diag(nV, nV2, nCl, nCl2, nKa, nF) = c(1, 1, 1, 1, 1, 1))
}
However, the bioavailability result does not make sense. It is over 100%. I know the BA is about 80% using the NCA analysis. A file with text and graph model is attached. Do you see any errors in the model text and graph?
Thanks for the help!
Wendy
/extranet/media/kunena/attachments/legacy/images/model_image.jpg