Precision of PowerTOST [Power / Sample Size]

posted by Helmut Homepage – Vienna, Austria, 2017-12-09 18:00 (2322 d 16:58 ago) – Posting: # 18049
Views: 13,864

Hi mittyri,

❝ I tried to figure out the precision of PowerTOST estimations with the following code:

  scABELdata$power[j] <- power.scABEL(CV=0.8, n=54, theta0=0.95,

                                      design="2x3x3", nsims=10000,

                                      setseed=F)

❝ please correct if it is wrong.


I think so.
You got the right answer to a wrong question. For theta0=0.95 and targetpower=0.8 sample sizes are 57 for design="2x3x3" and 56 for design="2x2x3". Therefore, with 54 power will be <0.8 for both designs:

library(PowerTOST)
power.scABEL(CV=0.8, theta0=0.95, n=57, design="2x3x3")
# [1] 0.82231
power.scABEL(CV=0.8, theta0=0.95, n=54, design="2x3x3")
# [1] 0.79917
power.scABEL(CV=0.8, theta0=0.95, n=56, design="2x2x3")
# [1] 0.81545
power.scABEL(CV=0.8, theta0=0.95, n=54, design="2x2x3")
# [1] 0.79883


❝ […] I would say that sample mean 80% in PowerTOST is a real unfortune :cool:


On the contrary; fortunately works as designed with the right sample sizes. :-D

[image] [image]

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

UA Flag
Activity
 Admin contact
22,987 posts in 4,824 threads, 1,663 registered users;
90 visitors (1 registered, 89 guests [including 7 identified bots]).
Forum time: 11:58 CEST (Europe/Vienna)

The only way to comprehend what mathematicians mean by Infinity
is to contemplate the extent of human stupidity.    Voltaire

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