Nitpicker! [Two-Stage / GS Designs]

posted by ElMaestro  – Denmark, 2017-08-22 19:27 (2734 d 16:11 ago) – Posting: # 17733
Views: 28,801

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
23,381 posts in 4,914 threads, 1,662 registered users;
38 visitors (1 registered, 37 guests [including 9 identified bots]).
Forum time: 10:39 CET (Europe/Vienna)

Science is built up with facts, as a house is with stones.
But a collection of facts is no more a science
than a heap of stones is a house.    Henri Poincaré

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