Kudos! [Power / Sample Size]

posted by ElMaestro  – Denmark, 2017-02-05 22:15 (3076 d 03:08 ago) – Posting: # 17016
Views: 29,577

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
23,428 posts in 4,929 threads, 1,691 registered users;
34 visitors (0 registered, 34 guests [including 16 identified bots]).
Forum time: 02:24 CEST (Europe/Vienna)

To know that we know what we know,
and to know that we do not know what we do not know,
that is true knowledge.    Nicolaus Copernicus

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