partial AUCs [NCA / SHAM]

posted by PharmCat  – Russia, 2020-11-13 01:16 (14 d 13:30 ago) – Posting: # 22066
Views: 233

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,214 posts in 4,427 threads, 1,481 registered users;
online 19 (0 registered, 19 guests [including 14 identified bots]).
Forum time: Friday 14:47 CET (Europe/Vienna)

Biostatistician. One who has neither the intellect for mathematics
nor the commitment for medicine but likes to dabble in both.    Stephen Senn

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