Bioequivalence and Bioavailability Forum

Hi Helmut,

I don't think that Certara (Phoenix) is not motivated to do better software. Look at LMM library - this is ancient as mammoth excrement and should be rewritten years ago (it is a big luck that we have much memory and CPU time, but this still doesn't work for a large number of observations).



So... In Julia this calculated fine:
Code:
df = CSV.file("pkpart.csv") |> DataFrame
pkds       = ClinicalTrialUtilities.pkimport(df; conc = :conc, time = :time)
pk         = ClinicalTrialUtilities.nca!(pkds, calcm = :luld)
println("AUCall: ", pk[1, :AUCall])
pall = pk[1, :AUCall]
ClinicalTrialUtilities.setdosetime!(pkds, ClinicalTrialUtilities.DoseTime(dose = 120, time = 0, tau = 4))
pk         = ClinicalTrialUtilities.nca!(pkds, calcm = :luld)
println("AUCtau 0 - 4: ", pk[1, :AUCtau])
p1 = pk[1, :AUCtau]
ClinicalTrialUtilities.setdosetime!(pkds, ClinicalTrialUtilities.DoseTime(dose = 120, time = 4, tau = 72-4))
pk         = ClinicalTrialUtilities.nca!(pkds, calcm = :luld)
println("AUCtau 4 - 72: ", pk[1, :AUCtau])
p2 = pk[1, :AUCtau]
println("Check pall ≈ p1 + p2: ", pall ≈ p1 + p2)


Out:
julia>
AUCall: 980.7754537633873
AUCtau 0 - 4: 54.04557219475543
AUCtau 4 - 72: 926.7298815686318
Check pall ≈ p1 + p2: true