Hello,
I am trying to fit a lag absorption model, but I encountered an error. Since I have successfully implemented lag absorption in the past, I revisited a previous lag model that had run successfully and found that the same error occurs. Below is a snapshot of my dose control block and the stack trace of the error.
Could you please look into this or suggest a solution if this issue has been resolved before?
` @dosecontrol begin
lags = (; Oral = tvlag_or * exp(η[6])
end
@dynamics begin
Oral' = -Ka_or*Oral
Central' = Ka_or*Oral - (CL/Vc)*Central
end
Can you try removing the variability on the lag parameter and run the model. It is well documented that having BSV on lag parameters result in unstable models that provide different results each time.
I think I found the issue, or at least what I think may be causing this error.
When I explicitly define the random effects (e.g., ωcl ∈ RealDomain(lower = 0.0001),…) rather than the PDiagDomain(x) set up, the lag absorption model ran fine.
I know I have tried both ways in the past and it worked. Perhaps this is a bug because random effect was not even included on the tvlag, or the new Pumas update does not allow this anymore.