Bioequivalence and Bioavailability Forum

Gentleman,

Could you please help me out with my sas code for the 95% upper confidence bound? I am in the process of writing the sas code for a 3-way partial rep study (1T, Reference Rep), I followed the codes from FDA's Progesterone guidance but I am stuck at the 95% upper confidence bound computation. How do you get SAS to compute using the (boundx, x) from IGLM2 and (boundy, y) from DGLM1?I got the correct data for IGLM2. DGLM1 is where I have issue, I got 2 sets of (boundy and y), 1 set from MODEL and 1 set from ERROR. The critbound computation should take (boundx, x) from IGLM2 and (boundy, y) from ERROR from DGLM1, correct? For some reason SAS took (boundy, y) from MODEL from DGLM1 along with (boundx, x) from IGLM2. It should've taken (boundy, y) from ERROR from DGLM1.

Here is my code (partial)...
.
.
proc glm data=scavbe;
class sequence;
model ilat=sequence/clparm alpha=0.1;
estimate 'average' intercept 1 sequence 0.3333333333 0.3333333333 0.3333333333;
ods output overallanova=iglm1;
ods output Estimates=iglm2;
ods output NObs=iglm3;
title1 'scaled average BE';
run;

data iglm2a; set iglm2;
pointest=exp(estimate);
x=estimate**2 - StdErr**2;
boundx=(max((abs(LowerCL)),(abs(UpperCL))))**2;
run;

proc glm data=scavbe;
class sequence;
model dlat=sequence;
ods output overallanova=dglm1;
ods output NObs=dglm3;
title1 'scaled average BE';
run;

data dglm1a; set dglm1;
dfd=df;
s2wr=ms/2;

data UCL;
merge iglm2a dglm1a;
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;

Thank you.

Jag

---

Edit: Category changed. BTW, if you are always addressing “Gentleman” (why singular?) you will decrease your chances getting answers from the ladies here. I know four female biostatisticians here personally. [Helmut]