Calculate the absorption profile of an in vivo concentration profile using the Loo-Riegelman approach. This is only valid for profiles that have been modelled with a two-compartments PK model.
The formula is Fabs(t)=(Cp(t)+Cperipheral(t)+K10*AUC0t(t))/(K10*AUC0inf) where K10=Cl/V and Cperipheral(t_n)=k12*Delta Cp*Delta t/2+k12/k21 * Cp(t_n-1)(1-exp(-k21*Delta t))+Cperipheral(t_n-1)*exp(-k21*Delta t)
In this implementation it is assumed that AUC0inf is the last AUC0t observed, meaning that Cp(t) has almost vanished in the last samples.
This protocol is much more accurate when the input Cp(t) is reinterpolated to a small time step like 0.5 minutes.
Reference: Leon Shargel, Susanna Wu-Pong, Andrew B.C. Yu. Applied Biopharmaceutics & Pharmacokinetics, 6e. McGraw Hill, 1999. Chap. 7
addSample(sampleName, t, y)¶
calculateCperipheral(t, Cp, k12, k21)¶