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

posted by d_labes  – Berlin, Germany, 2009-10-07 10:46 (6096 d 01:40 ago) – Posting: # 4310
Views: 45,704

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,654 posts in 4,992 threads, 1,571 registered users;
138 visitors (0 registered, 138 guests [including 17 identified bots]).
Forum time: 12:27 CEST (Europe/Vienna)

Always listen to experts.
They’ll tell you what can’t be done and why.
Then do it.    Robert A. Heinlein

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