Dear colleague
The 30 minutes (assuming you use minute foir time) shoukld be put as the duration of the extravascular infusion.
I did it with the interface ands the final code is below:
test(){
deriv(Aa = - (Aa * K_dosing_to_lungs))
deriv(Alungs = (Aa * K_dosing_to_lungs)- (Alungs * klung_plasma))
deriv(A1 = (Alungs * klung_plasma)- (Q_C2_C3 * C)- (Q_C2_C4 * C)- (Q_C2_C5 * (C - C5)))
deriv(Aliver = (Q_C2_C3 * C)- (CL_Cliver_A0liver * Cliver))
deriv(Akidney = (Q_C2_C4 * C)- (CL_Ckidney_A0kidney * Ckidney))
deriv(A5 = (Q_C2_C5 * (C - C5)))
urinecpt(A0liver = (CL_Cliver_A0liver * Cliver))
urinecpt(A0kidney = (CL_Ckidney_A0kidney * Ckidney))
dosepoint(Aa, duration = (30), idosevar = AaDose, infdosevar = AaInfDose, infratevar = AaInfRate)
Clung = Alungs / Vlungs
C = A1 / V
Cliver = Aliver / Vliver
Ckidney = Akidney / Vkidney
C5 = A5 / V5
error(CEpslung = 0.1)
observe(CObslung = Clung * (1 + CEpslung))
error(CEpsliver = 0.1)
observe(CObsliver = Cliver * (1 + CEpsliver))
error(CEpskidney = 0.1)
observe(CObskidney = Ckidney * (1 + CEpskidney))
error(CEpsplasma = 0.1)
observe(CObsplasma = C * (1 + CEpsplasma))
deriv(E = Kin * (1 + Emax * C / (C + EC50)) - Kout * E)
sequence{E = (Kin) / (Kout)}
error(EEps = 1)
observe(EObs = E + EEps)
stparm(Vlungs = tvVlungs * exp(nVlungs))
stparm(V = tvV * exp(nV))
stparm(Vliver = tvVliver * exp(nVliver))
stparm(Vkidney = tvVkidney * exp(nVkidney))
stparm(V5 = tvV5 * exp(nV5))
stparm(K_dosing_to_lungs = tvK_dosing_to_lungs * exp(nK_dosing_to_lungs))
stparm(klung_plasma = tvklung_plasma * exp(nklung_plasma))
stparm(Q_C2_C3 = tvQ_C2_C3 * exp(nQ_C2_C3))
stparm(Q_C2_C4 = tvQ_C2_C4 * exp(nQ_C2_C4))
stparm(Q_C2_C5 = tvQ_C2_C5 * exp(nQ_C2_C5))
stparm(CL_Cliver_A0liver = tvCL_Cliver_A0liver * exp(nCL_Cliver_A0liver))
stparm(CL_Ckidney_A0kidney = tvCL_Ckidney_A0kidney * exp(nCL_Ckidney_A0kidney))
stparm(Emax = tvEmax * exp(nEmax))
stparm(EC50 = tvEC50 * exp(nEC50))
stparm(Kin = tvKin * exp(nKin))
stparm(Kout = tvKout * exp(nKout))
fixef(tvVlungs = c(, 1, ))
fixef(tvV = c(, 1, ))
fixef(tvVliver = c(, 1, ))
fixef(tvVkidney = c(, 1, ))
fixef(tvV5 = c(, 1, ))
fixef(tvK_dosing_to_lungs = c(, 1, ))
fixef(tvklung_plasma = c(, 1, ))
fixef(tvQ_C2_C3 = c(, 1, ))
fixef(tvQ_C2_C4 = c(, 1, ))
fixef(tvQ_C2_C5 = c(, 1, ))
fixef(tvCL_Cliver_A0liver = c(, 1, ))
fixef(tvCL_Ckidney_A0kidney = c(, 1, ))
fixef(tvEmax = c(, 1, ))
fixef(tvEC50 = c(, 1, ))
fixef(tvKin = c(, 1, ))
fixef(tvKout = c(, 1, ))
ranef(diag(nV5, nQ_C2_C3, nQ_C2_C5, nK_dosing_to_lungs, nVlungs, nklung_plasma, nVliver, nV, nVkidney, nCL_Cliver_A0liver, nCL_Ckidney_A0kidney, nEmax, nEC50, nKin, nKout, nQ_C2_C4) = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1))
}
About how to find initial estimates, I need the data set to help you. There are different techniques like pseudo NCA in conjunction with the initial estimate tab and information from the biology but without data there is no one way to explain all this.
Please do something to your data in such a way I can get those and helping you fitting the model to the data.
You can change the headers names, the units. This usually is enough.
The reason I need the data is because then I can work on a strategy to find initial estimates and to upgrade the model if necessary.
Best Regards
​Serge
​poppharm@gmail.com just in case if you need it one day