Bioequivalence and Bioavailability Forum • Power of TOST (R package MBESS)

Power of TOST (R package MBESS) [🇷 for BE/BA]

posted by Helmut – Vienna, Austria, 2009-10-07 04:18 (5733 d 05:16 ago) – Posting: # 4309
Views: 34,006

Dear D Labes!

❝ Have a look at the R package MBESS. If I had seen it right there is a function power.equivalence.md that will give you some spare time.

Quick tests

require(MBESS)

n  <- c(8,12,18,24,30,40,60) # total sample size

df <- n-2

power.equivalence.md.plot(.05, logscale=TRUE, .8, 1.25, .20,

                          n, df, 'Diletti, Figure 1c')

In the text output power is the theoretical value 0.05 (=α) at the borders of the acceptance range.

Now for ElMaestro's candidates:

n <- 6

power.equivalence.md(.05, logscale=TRUE, 0.8, 1.25, 0.95, 0.65, n, n-2)

[1] 0.001048352

n <- 38

power.equivalence.md(.05, logscale=TRUE, 0.8, 1.25, 0.95, 0.65, n, n-2)

[1] 0.01392505

Yet some other values for our collection. ;-)

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

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

Complete thread:

non-central t, power, R: ElMaestro 2009-10-03 19:05
- non-central t, power, R: Helmut 2009-10-03 21:45
  - non-central t, power, R: ElMaestro 2009-10-04 21:51
    - non-central t, power, R: Helmut 2009-10-04 23:58
- non-central t is not the question d_labes 2009-10-05 12:58
  - SAS -> R -> SAS Helmut 2009-10-05 14:53
    - Lost in translation d_labes 2009-10-06 14:31
      - Lost in translation Helmut 2009-10-06 14:42
        
        Macacus cynocephalus d_labes 2009-10-06 15:14
        
        Lost in translation Helmut 2009-10-06 21:46
  - non-central t is not the question ElMaestro 2009-10-05 17:32
    - Power of TOST in R d_labes 2009-10-06 16:45
      - Power of TOST (R package MBESS)Helmut 2009-10-07 02:18
        
        Power of TOST (R package MBESS) d_labes 2009-10-07 08:46
        
        Power of TOST (R package MBESS) Helmut 2009-10-07 12:38
      - The Power at limits d_labes 2009-10-07 11:32
        
        Power at limits 2 d_labes 2009-10-07 15:15
        
        package MBESS Helmut 2009-10-07 22:36
        
        Package MBESS and Power curiosity d_labes 2009-10-15 14:57
        
        Package MBESS and Power curiosity ElMaestro 2009-10-17 12:26
        
        Apfelstrudel and EFG 2.01 Helmut 2009-10-17 13:44
        EFG and another power curiosity d_labes 2009-10-20 11:35
        
        EFG and another power curiosity ElMaestro 2009-11-23 19:03
        
        EFG and another power curiosity Helmut 2009-11-23 19:17
- non-central t, power, R: yjlee168 2009-10-07 00:08
  - PASS Helmut 2009-10-07 01:16
    - PASS yjlee168 2009-10-12 23:17
- R code for power ElMaestro 2009-10-09 11:32
  - Belgian beer ? Ohlbe 2009-10-09 11:50
  - Another one - R code for power of TROST d_labes 2009-10-09 13:50