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

posted by d_labes  – Berlin, Germany, 2009-10-07 10:46 (6089 d 06:32 ago) – Posting: # 4310
Views: 45,641

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,653 posts in 4,991 threads, 1,570 registered users;
209 visitors (0 registered, 209 guests [including 23 identified bots]).
Forum time: 17:18 CEST (Europe/Vienna)

The epistemological value of probability theory is based on the fact
that chance phenomena, considered collectively and on a grand scale,
create non-random regularity.    Andrey Kolmogorov

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