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

posted by preyes323 – 2010-07-17 13:14 (5812 d 19:27 ago) – Posting: # 5648
Views: 16,341

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,654 posts in 4,992 threads, 1,571 registered users;
139 visitors (0 registered, 139 guests [including 14 identified bots]).
Forum time: 08:42 CEST (Europe/Vienna)

Always listen to experts.
They’ll tell you what can’t be done and why.
Then do it.    Robert A. Heinlein

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