Bayesian modeling and model–based decision analysis 203
data theta1; set postout;
keep iteration theta1_1-theta1_37; run;
data theta2; set postout;
keep iteration theta2_1-theta2_37; run;
data u1; set postout;
keep iteration u1_1-u1_37; run;
proc transpose data=theta1 out=theta1t(rename=(col1=theta1));
by iteration; run;
proc transpose data=theta2 out=theta2t(rename=(col1=theta2));
by iteration; run;
proc transpose data=u1 out=u1t(rename=(col1=u1));
by iteration; run;
data pred;
merge postout (keep=iteration emax ec50 e0 var_y) theta1t theta2t u1t;
by iteration;
cl=exp(theta1);
v=exp(theta2);
do time=10 to 480 by 10;
tinf=240;
t1=1;
if time>tinf then t1=0;
t2=tinf*(1-t1)+t1*time;
do dose=0.1 to 5 by 0.2;
pk=(dose/cl)*(1-exp(-cl*t2/v))*exp(-cl*(1-t1)*(time-tinf)/v); ...