Kudos! [Power / Sample Size]

posted by ElMaestro  – Denmark, 2017-02-05 20:15 (1450 d 09:46 ago) – Posting: # 17016
Views: 21,326

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:

Activity
 Admin contact
21,312 posts in 4,445 threads, 1,489 registered users;
online 6 (0 registered, 6 guests [including 2 identified bots]).
Forum time: Tuesday 06:02 UTC (Europe/Vienna)

Any one who considers arithmetical methods
of producing random digits is, of course,
in a state of sin.    John von Neumann

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