Kudos! [Power / Sample Size]

posted by ElMaestro  – Denmark, 2017-02-05 22:15 (2608 d 17:59 ago) – Posting: # 17016
Views: 25,399

Hi Hötzi,

❝ Quick & extremely dirty:

library(PowerTOST)

❝ 100*power.TOST(alpha=0.5, CV=0.3, theta0=0.95, n=40, theta1=0.8, theta2=1)-1.25

❝ # [1] 76.57964


I can't follow this, alpha 0.5? Substraction of 1.25% somehow from ... what?

Instead, I propose something along the lines of:

pwrDK222BE=function(alpha, CV, theta0, theta1, theta2, n)
{
  ##power the standard way is P(GMR is below theta2)-P(GMR is below theta1)
  ##and thus we just plug in upper or lower= 1 to get us what we need.
  P= power.TOST(alpha=alpha, CV=CV,
       theta0=theta0, n=n, theta1=theta1, theta2=theta2)
  Plow= power.TOST(alpha=alpha, CV=CV,
        theta0=theta0, n=n, theta1=theta1, theta2=1.0) 
  Phi  =power.TOST(alpha=alpha, CV=CV,
       theta0=theta0, n=n, theta1=1.0, theta2=theta2)
  return(P-Plow-Phi)
}

pwrDK222BE(0.05, 0.3, 0.95, 0.8, 1.25, 40)


[1] 0.7652256

Pass or fail!
ElMaestro

Complete thread:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
119 visitors (0 registered, 119 guests [including 6 identified bots]).
Forum time: 16:15 CET (Europe/Vienna)

With four parameters I can fit an elephant,
and with five I can make him wiggle his trunk.    John von Neumann

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