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

posted by preyes323 – 2010-07-17 13:14 (5417 d 14:00 ago) – Posting: # 5648
Views: 13,851

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 :cool:.)


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. :ponder:


I see. :-)

Thanks and Regards,
Paolo

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,669 registered users;
169 visitors (0 registered, 169 guests [including 13 identified bots]).
Forum time: 03:14 CEST (Europe/Vienna)

By all means let’s be open-minded, but not so open-minded
that our brains drop out.    Richard Dawkins

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