Bioequivalence and Bioavailability Forum • bear v1.0.0 for R

bear v1.0.0 for R - first impressions [🇷 for BE/BA]

posted by Helmut – Vienna, Austria, 2008-07-14 06:12 (6183 d 04:50 ago) – Posting: # 2021
Views: 72,205

Dear Hsin-ya and Yung-jin!

Wonderful, thank you!

I had a quick view already, some first impressions below:

The plots for selection of data points in NCA should be in semilog scale rather than linear.
Maybe it’s possible to open only one graphics device and turn the history on – instead of opening a new device for every plot.
In sample size estimation you should consider rounding up to the next even number:

Example:
For 20% CV, theta 95%, power 80% bear comes up with:

--------------------------------------------------------------------------

                               <<Suggestion>>

--------------------------------------------------------------------------

             n>= 9.13  (sample size number for per sequence)

            Total sample size 18

--------------------------------------------------------------------------

I implemented Hauschke’s approximation in Excel (sorry, a quick shot) and obtained:
9.1763 / sequence = 18.3526 / total -> rounded up to next even number = 20
David Dubins’ FARTSSIE comes up with 19 / total
StudySize2.01: 9.1876 / sequence = 18.3752 / total
My R-code (modified from Patterson’s/Jones’ SAS) would give a total of 20 (power 83.47%).
With bear’s suggested n=18 I got power of only 79.1% (StudySize).

I failed in reading in the NCA test file. I tried it with tabulator and semicolon as field separators and got an error:

Error in readChar(con, 5L, useBytes = TRUE) : 

  cannot open connection

I would avoid the term "Tests of Hypothese using the Type I MS for SUBJECT(SEQUENCE) as an error term" and use "Tests of Hypothesis for SUBJECT(SEQUENCE) as an error term" instead. There were some discussions on the R-list and here also about Types I-III; don't fall into the trap of using SAS’ proprietary terminology…
IMHO a posteriori power is meaningless.

Two points for the wish-list:

CV_intra and CV_inter
intra-subject residuals – something everybody misses in WinNonlin…

Again, many thanks; I will continue testing – time allowing.

—
Dif-tor heh smusma 🖖🏼 Довге життя Україна!
Helmut Schütz

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Thread locked

Complete thread:

bear v1.0.0 - a data analytical tool for ABE in R yjlee168 2008-07-12 20:53
- bear v1.0.0 for R - first impressionsHelmut 2008-07-14 04:12
  - bear v1.0.0 for R - first impressions yjlee168 2008-07-14 10:26
    - bear v1.0.0 for R - first impressions Helmut 2008-07-14 14:02
      - bear v1.0.0 for R - first impressions yjlee168 2008-07-14 18:12
        
        bear v1.0.0 for R... Helmut 2008-07-14 18:56
        
        bear v1.0.0 for R... yjlee168 2008-07-14 19:30
        
        bear v1.0.0 for R... Helmut 2008-07-14 20:35
      - bear v1.0.0 for R - first impressions yjlee168 2008-07-15 09:44
- bear v1.0.0 - a data analytical tool for ABE in R martin 2008-07-14 12:18
  - bear v1.0.0 - a data analytical tool for ABE in R yjlee168 2008-07-14 19:02
    - bear v1.0.0 - a data analytical tool for ABE in R Aceto81 2008-07-15 10:07
      - bear v1.0.0 - a data analytical tool for ABE in R yjlee168 2008-07-16 07:22
        
        bear v1.0.0 - a data analytical tool for ABE in R Aceto81 2008-07-16 09:47
        
        bear v1.0.0 - a data analytical tool for ABE in R yjlee168 2008-07-16 21:32
      - lambda_z estimation yjlee168 2008-09-23 12:13
        
        lambda_z estimation d_labes 2008-09-23 14:41
        
        lambda_z estimation yjlee168 2008-09-24 22:41
        lambda_z estimation Aceto81 2008-09-25 10:56
        
        lambda_z estimation yjlee168 2008-09-25 21:00
        lambda_z estimation yjlee168 2008-09-26 08:03
        
        lambda_z estimation Aceto81 2008-09-26 10:25
        
        lambda_z estimation yjlee168 2008-09-28 00:08
        lambda_z estimation yjlee168 2008-09-29 13:32
        
        WinNonlin 5.2.1 vs. 6 beta Helmut 2008-09-29 15:54
        
        WinNonlin 5.2.1 vs. 6 beta yjlee168 2008-09-29 18:45
        WinNonlin 5.2.1 vs. 6 beta --> new finding yjlee168 2008-09-30 08:10
        
        Example 4 in SASophylistic d_labes 2008-09-30 09:43
        
        Example 4 in SASophylistic yjlee168 2008-09-30 12:23
        
        TTT method for lambda_z estimation yjlee168 2008-09-30 12:34
        
        TTT method for lambda_z estimation d_labes 2008-09-30 14:48
        
        TTT method for lambda_z estimation yjlee168 2008-09-30 20:28
        
        TTT method plus best fit combined d_labes 2008-10-01 08:43
        
        AIC or ARS as the best fit criterion? yjlee168 2008-10-02 12:40
        
        AIC or ARS as the best fit criterion? d_labes 2008-10-02 13:59
        
        TTT method for lambda_z estimation Aceto81 2008-09-30 15:11
        
        TTT method for lambda_z estimation yjlee168 2008-09-30 19:57
        
        TTT method for lambda_z estimation Aceto81 2008-10-01 14:18
        
        TTT method for lambda_z estimation yjlee168 2008-10-02 12:29
        
        TTT method for lambda_z estimation Aceto81 2008-10-03 15:33
        
        TTT method for lambda_z estimation yjlee168 2008-10-03 21:22
- bear v1.0.0 - a data analytical tool for ABE in R ElMaestro 2008-07-22 14:35
  - bear v1.0.0 - a data analytical tool for ABE in R yjlee168 2008-07-23 09:55
    - bear v1.0.0 - a data analytical tool for ABE in R ElMaestro 2008-07-24 10:04
      - bear v1.0.0 - a data analytical tool for ABE in R yjlee168 2008-07-25 20:38
        
        bear v1.0.0 - a data analytical tool for ABE in R ElMaestro 2008-07-28 08:42
    - bear v1.0.0 - a data analytical tool for ABE in R ElMaestro 2008-09-23 16:10
      - bear v1.0.0 - a data analytical tool for ABE in R yjlee168 2008-09-24 22:00
        
        bear v1.0.0 - a data analytical tool for ABE in R ElMaestro 2008-09-25 09:03
        
        bear v1.0.0 - a data analytical tool for ABE in R ElMaestro 2008-09-25 13:52
        
        bear v1.0.0 - a data analytical tool for ABE in R Helmut 2008-09-25 14:39
        
        bear v1.0.0 - a data analytical tool for ABE in R ElMaestro 2008-09-25 14:55
        
        bear v1.0.0 - a data analytical tool for ABE in R Helmut 2008-09-25 15:33
        
        GLM in R and the power to know d_labes 2008-09-25 16:45
        
        GLM in R and the power to know ElMaestro 2008-09-25 19:23
        
        nesting in R? yjlee168 2008-09-25 20:49
        
        nesting in R? ElMaestro 2008-09-25 22:01
        
        lm in R Helmut 2008-09-26 00:23
        
        lm in R ElMaestro 2008-09-26 09:15
        lm in R Aceto81 2008-09-26 12:26
  - Type III SS in balanced, crossover BE studies? yjlee168 2008-10-18 21:03
    - Type III SS in balanced, crossover BE studies? ElMaestro 2008-10-21 13:09
      - Type III SS in balanced, crossover BE studies? yjlee168 2008-10-21 13:19
        
        Only balanced, crossover BE studies? d_labes 2008-10-21 14:31
        
        Only balanced, crossover BE studies? yjlee168 2008-10-21 19:01
        
        Only balanced, crossover BE studies? ElMaestro 2008-10-24 11:26
        
        Thread locked Helmut 2008-10-24 11:58