Bioequivalence and Bioavailability Forum • possible solutions?

possible solutions? [🇷 for BE/BA]

posted by yjlee168 – Kaohsiung, Taiwan, 2016-02-02 22:30 (3293 d 19:30 ago) – Posting: # 15921
Views: 25,356
(edited on 2016-02-03 08:46)

Dear Astea and mittyri,

I think that I have found the solution for your datasets. This time I keep all user's input data unchanged. It used to assign a 'NA' to any conc.=0 except the first data point (time =0). Calculations for AUC and AUMC are fine for all-linear trapezoidal method. If linear-up/log-down trapezoidal method is applied, the following situations will switch to all-linear: (1) C_i=C_i-1; (2) C_i>C_i-1; and (3) C_i=0 or C_i-1=0 (new) for the trapezoidal calculation between the data point of (t_i-1, C_i-1) and (t_i, C_i).

So for Astea's dataset (all linear),I got

    

 subj  time conc  AUC(0-t) AUMC(0-t)

    1  0.00    0     0.000     0.000

    1  0.25    0     0.000     0.000

    1  0.50  125    15.625     7.812

    1  0.75  567   102.125    68.781

    1  1.00  932   289.500   238.438

    1  1.50 1343   858.250   975.062

    1  2.00 1739  1628.750  2348.188

    1  3.00 1604  3300.250  6493.188

    1  4.00 1460  4832.250 11819.188

    1  8.00  797  9346.250 36251.188

    1 12.00  383 11706.250 58195.188

    1 24.00   72 14436.250 96139.188

And with same Astea's dataset (linear-up/log-down),

    

 subj  time conc  AUC(0-t) AUMC(0-t)

    1  0.00    0     0.000     0.000

    1  0.25    0     0.000     0.000

    1  0.50  125    15.625     7.812

    1  0.75  567   102.125    68.781

    1  1.00  932   289.500   238.438

    1  1.50 1343   858.250   975.062

    1  2.00 1739  1628.750  2348.188

    1  3.00 1604  3299.341  6513.416

    1  4.00 1460  4830.212 11859.468

    1  8.00  797  9211.243 37267.003

    1 12.00  383 11471.007 59317.527

    1 24.00   72 13703.908 95940.683

for mittyri's dataset (all linear),

    

 subj  time conc  AUC(0-t) AUMC(0-t)

    1  0.00    0     0.000     0.000

    1  0.25    0     0.000     0.000

    1  0.50  125    15.625     7.812

    1  0.75    0    31.250    15.625

    1  1.00  932   147.750   132.125

    1  1.50 1343   716.500   868.750

    1  2.00 1739  1487.000  2241.875

    1  3.00 1604  3158.500  6386.875

    1  4.00 1460  4690.500 11712.875

    1  8.00  797  9204.500 36144.875

    1 12.00  383 11564.500 58088.875

    1 24.00   72 14294.500 96032.875

The only thing remained unsolved is how to do a semi-log plot for these zero conc.? It will be skipped if included for the plot. Then we will see a broken line since log(0) cannot be plotted. Any suggestion for this? Thanks in advanced.

—
All the best,
-- Yung-jin Lee
bear v2.9.2:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here

Complete thread:

bear v2.7.0 released yjlee168 2015-12-07 09:18 [🇷 for BE/BA]
- Suggestion ElMaestro 2015-12-07 10:12
  - Suggestion yjlee168 2015-12-07 10:55
    - Suggestion Helmut 2015-12-07 14:01
      - Suggestion yjlee168 2015-12-07 18:51
- bear v2.7.0 released Astea 2015-12-10 16:36
  - if Clast is zero in bear yjlee168 2015-12-10 17:23
    - if Clast is zero in bear Astea 2015-12-10 20:50
      - if Clast is zero in bear yjlee168 2015-12-11 03:28
        
        if Clast is zero in bear Astea 2015-12-11 09:39
        
        if Clast is zero in bear yjlee168 2015-12-11 17:16
    - bear v2.7.1 yjlee168 2015-12-15 10:21
      - bear v2.7.1 Astea 2015-12-15 19:32
        
        BLLOQ ≠ 0 with measured drug plasma conc. yjlee168 2015-12-15 20:17
        
        BLLOQ ≠ 0 with measured drug plasma conc. Helmut 2015-12-16 01:45
        
        BLLOQ ≠ 0 with measured drug plasma conc. yjlee168 2015-12-16 11:38
        
        BLLOQ ≠ 0 ∀ t > 0 Helmut 2015-12-16 12:40
        
        BLLOQ ∋ { 0, ≈0, BLLOQ{...} ∀ t > 0} yjlee168 2015-12-16 18:22
        
        BLLOQ at first two time ponts, AUC calculation Astea 2016-01-30 17:11
        
        zero-zero with AUC calculation in bear yjlee168 2016-01-30 19:00
        
        BQL rules in bear mittyri 2016-02-01 22:16
        
        BQL rules in bear yjlee168 2016-02-02 10:59
        possible solutions?yjlee168 2016-02-02 21:30
        zero (0), missing (NA), or BLLOQ (0 or NA) in bear yjlee168 2016-02-03 09:04
        
        zero (0), missing (NA), or BLLOQ (0 or NA) in bear mahmoud-teaima 2016-02-08 21:28
        
        Please avoid Full Quotes (“TOFU”)! Ohlbe 2016-02-09 00:51
        
        error message in bear mahmoud-teaima 2016-02-09 04:55
        
        please check your dataset file yjlee168 2016-02-09 06:40
        
        please check your dataset file mahmoud-teaima 2016-02-09 15:00
        
        please check your dataset file yjlee168 2016-02-09 15:38
        
        A possible R bug? yjlee168 2016-02-09 19:14