Bioequivalence and Bioavailability Forum

To Err is Human [General Sta­tis­tics]

posted by preyes323 – 2010-07-17 13:14 (5813 d 15:34 ago) – Posting: # 5648
Views: 16,346

Hi d_labes,

❝ I have calculated using (is theta in your code)

❝ theta2 = ((log(1.25)/0.2936)2=(0.76)2

❝ (BTW: Yes nitpickers, of course I know it's really

❝ (log(1.25)/sqrt(log(0.3*0.3+1)))^2=(0.7601283...)^2 .)

Yes, I included theta in my code. From what I see our main difference might be in the formula for the crit bound. The formula I used is the one from the FDA progesterone guideline. Did you use the same?

data SABE.upper_cmax;
   merge SABE.dglm1_cmax SABE.iglm2_cmax;
   theta=((log(1.25))/0.25)**2;
   y=-theta*s2wr;
   boundy=y*dfd/cinv(0.95,dfd);
   sWR=sqrt(s2wr);
   critbound=(x+y)+sqrt(((boundx-x)**2)+((boundy-y)**2));
run;

❝ Dear Paolo: Sorry but this question is plaguing us as well and is source of our continuing discussion here.

❝ Especially we don't know if one can trust in the results of the FDA Proc Mixed code in case of the extra-reference design your data came from.

I see.

Thanks and Regards,
Paolo

Complete thread:

• SABE Reference Variability and CI Computation preyes323 2010-07-12 17:12 
  • Dataset Helmut 2010-07-12 17:26
    • Dataset preyes323 2010-07-12 18:20
      • Data format Helmut 2010-07-12 18:38
        • Data format preyes323 2010-07-13 01:16
          • SAS and/or Phoenix/WinNonlin-experts around? Helmut 2010-07-13 02:29
            • Phoenix/WinNonlin-experts around? d_labes 2010-07-13 12:15
              • Phoenix/WinNonlin-experts around? Helmut 2010-07-13 13:42
                • Mixed interest d_labes 2010-07-13 15:11
                  • Mixed interest Helmut 2010-07-13 19:19
                    • lme() answer and beyond ... d_labes 2010-07-14 10:41
                      • lme() answer and beyond ... Helmut 2010-07-14 13:47
                        • AIC, BIC and that all ... d_labes 2010-07-15 12:01
          • To Err is Human d_labes 2010-07-13 11:51
            • To Err is Human, but... Helmut 2010-07-13 13:45
              • ... to Arr is Pirate d_labes 2010-07-13 15:47
                • Brilliant page!!! ElMaestro 2010-07-14 21:03
            • To Err is Human preyes323 2010-07-14 14:45
              • Heads up! Helmut 2010-07-14 15:17
              • To Err is Teacher d_labes 2010-07-15 11:05
            • To Err is Human preyes323 2010-07-16 02:28
              • To Err is Human d_labes 2010-07-16 10:42
                • To Err is Humanpreyes323 2010-07-17 11:14
                  • Regulatory constants d_labes 2010-07-19 10:00
          • Data format ElMaestro 2011-02-06 11:54
            • SAS System Viewer Helmut 2011-02-06 13:02