Nitpicker! [Two-Stage / GS Designs]

posted by ElMaestro  – Denmark, 2017-08-22 19:27 (2410 d 13:11 ago) – Posting: # 17733
Views: 27,002

Hi Hötzi,

❝ You have a point. IMHO, this question should be answered: yes!



And on basis of this story we will define Hötzi's universal law of BE:
ElMaestro is always right1.


By the way, since you are a speed devotee:
FindCritCV.Sec=function(Pwr.T, GMR, Nps, Is.St2=1, d=.00001, iniCV)
##the secant method has really nasty properties if the ini value is a little off.
{
  CV1=iniCV
  CV2=CV1*1.05
  y2=1
  while (abs(y2)>d)
  {
   y1=Power.calc(Nps, GMR, CV1, Is.St2)-Pwr.T
   y2=Power.calc(Nps, GMR, CV2, Is.St2)-Pwr.T
   CV3= CV2-y2*(CV2-CV1)/(y2-y1)
   CV1=CV2
   CV2=CV3
  }
  return(0.5*(CV1+CV2)) ##or play around
}

(...)
 for (Nps in 6:length(vecCritCVs))
  {
   if (Nps==6) iniCV=min(0.709*GMR-0.552, 0.709/GMR -0.552)
               else iniCV=vecCritCVs[Nps-1]
    vecCritCVs[Nps]<<-FindCritCV.Sec(0.8, GMR, Nps, 1, 1e-10, iniCV)
  }

(...)















Footnote 1:
Except when he isn't. Which also happens. With regrettable frequency...:-D

Pass or fail!
ElMaestro

Complete thread:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
93 visitors (0 registered, 93 guests [including 6 identified bots]).
Forum time: 07:39 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