Initial sample size guess for the Potvin methods [Two-Stage / GS Designs]

posted by ElMaestro  – Belgium?, 2017-08-19 16:14 (1193 d 03:09 ago) – Posting: # 17711
Views: 25,058

(edited by ElMaestro on 2017-08-19 17:14)

Hi Hötzi,

» Error in Power.calc(nps1, GMR, CV, Is.St2 = 1) :
»   could not find function "Power.calc"
» Did I specify the arguments wrongly?

That's why I included the remark:
##equation in Potvin et al.
I thought you'd be plugging in your own, but be that as it may.

An implementation of the function is here:

for (i in 1:5000)
    vecQT1[i]=qt(p=1-alpha1, df=i)
    vecQT2[i]=qt(p=1-alpha2, df=i)

Power.calc=function(Nps, GMR, CV, Is.St2=0)
  if (Is.St2==T) df=df-1
  ##cat("1-a=", 1-alpha, "df=",df,"\n")
  if (Is.St2==0) q=vecQT1[df]
    else         q=vecQT2[df]
  pw=pt( (nom1/den) -q, df) - pt( (nom2/den) +q, df)
  if (pw<0) pw=0.01


Edit: Forgot to write: For classical B the alphas are (0.0294, 0.0294) and for C they are (0.05, 0.0294) though that isn't related to the topic. Second, Hötzi, you gave a little table which I am not very sure how to read. Note this isn't about the sample size calculation itself, this is about the initial guess used for the sample size calculation by up/down iteration.

I could be wrong, but...

Best regards,

No, of course you do not need to audit your CRO if it was inspected in 1968 by the agency of Crabongostan.

Complete thread:

 Admin contact
21,206 posts in 4,425 threads, 1,481 registered users;
online 2 (0 registered, 2 guests [including 2 identified bots]).
Forum time: Tuesday 18:24 CET (Europe/Vienna)

All we know about the world teaches us that the effects of A and B
are always different—in some decimal place—for any A and B.
Thus asking “are the effects different?” is foolish.    John W. Tukey

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