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

posted by d_labes  – Berlin, Germany, 2009-10-07 10:46 (5736 d 08:27 ago) – Posting: # 4310
Views: 34,054

Dear Helmut!

Quick tests


Oh actually quick! But dirty :-P .
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 :-D .
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.

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,672 registered users;
57 visitors (0 registered, 57 guests [including 13 identified bots]).
Forum time: 19:14 CEST (Europe/Vienna)

Medical researches can be divided into two sorts:
those who think that meta is better and those
who believe that pooling is fooling.    Stephen Senn

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