Bioequivalence and Bioavailability Forum

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

posted by d_labes  – Berlin, Germany, 2009-10-07 10:46 (5733 d 06:38 ago) – Posting: # 4310
Views: 34,030

Dear Helmut!

❝ Quick tests

Oh actually quick! But dirty .
I assume that your 5th parameter is meant as CV.
But the function requires sigma, I think.

With that modification

library(MBESS)
CV <- 0.65
sigma<-sqrt(log(1+(CV)**2))

n <- 6
power.equivalence.md(0.05,TRUE, 0.8, 1.25, 0.95, sigma, n, n-2)
[1] 0.001620383

n <- 38
power.equivalence.md(0.05,TRUE, 0.8, 1.25, 0.95, sigma, n, n-2)
[1] 0.04509401

The answer for n=38 seems to me very sasophylistic home-brewed magic and Q-functionalistic .
Not so astonishing if you look at the author of that pice of R-software.

PS: I have just seen that the documentation suggest that CV can be used in case of log. But this seems erroneous.

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