partial AUCs [NCA / SHAM]

posted by PharmCat  – Russia, 2020-11-13 01:16 (252 d 15:03 ago) – Posting: # 22066
Views: 815

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:

Activity
 Admin contact
21,589 posts in 4,512 threads, 1,530 registered users;
online 21 (0 registered, 21 guests [including 5 identified bots]).
Forum time: Friday 17:19 CEST (Europe/Vienna)

Statistics is the grammar of science.    Karl Pearson

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