To Err is Human [General Statistics]
Hi d_labes,
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?
I see.
Thanks and Regards,
Paolo
❝ 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 [General Statistics]
- 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
- lme() answer and beyond ... Helmut 2010-07-14 13:47
- lme() answer and beyond ... d_labes 2010-07-14 10:41
- Mixed interest Helmut 2010-07-13 19:19
- Mixed interest d_labes 2010-07-13 15:11
- Phoenix/WinNonlin-experts around? Helmut 2010-07-13 13:42
- Phoenix/WinNonlin-experts around? d_labes 2010-07-13 12:15
- 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 Arr is Pirate d_labes 2010-07-13 15:47
- 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
- To Err is Humanpreyes323 2010-07-17 11:14
- To Err is Human d_labes 2010-07-16 10:42
- To Err is Human, but... Helmut 2010-07-13 13:45
- Data format ElMaestro 2011-02-06 11:54
- SAS System Viewer Helmut 2011-02-06 13:02
- SAS and/or Phoenix/WinNonlin-experts around? Helmut 2010-07-13 02:29
- Data format preyes323 2010-07-13 01:16
- Data format Helmut 2010-07-12 18:38
- Dataset preyes323 2010-07-12 18:20
- Dataset Helmut 2010-07-12 17:26