Precision of PowerTOST [Power / Sample Size]

posted by mittyri  – Russia, 2017-12-08 23:56 (2688 d 12:51 ago) – Posting: # 18048
Views: 15,075

Dear Detlew,

thank you for clarification
I tried to figure out the precision of PowerTOST estimations with the following code:
library("ggplot2")
library("PowerTOST")
reps <- 1E6
scABELdata <- data.frame(rep = 1:reps)
for(i in 1:reps){
  set.seed(i)
  scABELdata$power[i] <- power.scABEL(CV=0.8, n=54, theta0=0.95, design="2x3x3", nsims=10000, setseed=F)
}

ggplot(scABELdata, aes(power))+
  geom_density(fill = 2, alpha = 0.3)+
  theme_bw() +
  ggtitle(sprintf("%d reps: Mean is %.4f, SD is %.4f; %.2f %% obs are less than 0.8", reps, mean(scABELdata$power), sd(scABELdata$power), sum(scABELdata$power<0.8)/length(scABELdata$power)*100))


please correct if it is wrong.
The resulted plot:
[image]

BAsed on the mean and se calculated from sd and n=20 I would say that sample mean 80% in PowerTOST is a real unfortune :cool:

❝ BTW: Both Laszlos recommend to use PowerTOST if they were asked :thumb up:.

Fully agree and thanks again for this wonderful package.

Kind regards,
Mittyri

Complete thread:

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

Do not put your faith in what statistics say until you have carefully
considered what they do not say.    William W. Watt

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