partial AUCs [NCA / SHAM]

posted by PharmCat  – Russia, 2020-11-13 02:16 (1047 d 07:13 ago) – Posting: # 22066
Views: 1,924

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).

[image]

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

Complete thread:

UA Flag
Activity
 Admin contact
22,759 posts in 4,775 threads, 1,628 registered users;
18 visitors (0 registered, 18 guests [including 6 identified bots]).
Forum time: 10:30 CEST (Europe/Vienna)

Whenever a theory appears to you as the only possible one,
take this as a sign that you have neither understood the theory
nor the problem which it was intended to solve.    Karl R. Popper

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5