Hi,
I have IV and PO data from 3 animals. I was interested in fitting a simultaneous IV and PO model using Pheonix. In my source file I have added an ID for the type of dosing i.e. I sort the data using 1,2 ( all 3 animals that receive IV =1 and all three animals that receive PO = 2). Also, I have made a separate column for dose. My code is below.
test(){
deriv(Aa = - Ka * Aa - Km * Aa)
deriv(A1 = Ka * Aa - Ke * A1 - K12 * A1 + K21 * A2)
deriv(A2 = K12 * A1 - K21 * A2)
dosepoint(Aa, bioavail = (F))
dosepoint(A1)
C1 = A1 / V
error(CEps = 1)
observe(CObs1 = C1 * (1 + CEps))
stparm(Ka = tvKa * exp(nKa))
stparm(F = tvF)
stparm(V = tvV * exp(nV))
stparm(Ke = tvKe * exp(nKe))
stparm(K12 = tvK12 * exp(nK12))
stparm(K21 = tvK21 * exp(nK21))
fixef(tvKa= c(0, 0.021, ))
fixef(tvF = c(0, 0.3,1 ))
fixef(tvKm = c(0, 0.041, ))
fixef(tvV = c(, 135, ))
fixef(tvKe = c(0, 0.071, ))
fixef(tvK12 = c(0, 0.027,0.5 ))
fixef(tvK21 = c(0, 0.0028,0.5 ))
ranef(diag(nV, nKe, nKa, nK12, nK21) = c(1, 1, 1, 1, 1, 1, 1))
}
However, my model estimates a Ka and F for the IV group too.
Can anyone send me a test file for simultaneous IV and PO dosing and also let me know what I am missing in my code.
Thank you.