data one; input study $ trt $ y n; cards; solomon celeb 18 1356 solomon placebo 3 679 pfizer celeb 2 285 pfizer placebo 0 140 presap celeb 9 933 presap placebo 3 628 proc print; * proc logistic;class study trt/param=ref; * model y/n=study trt; run; %include 'K:\expand.binomial.data.txt'; %expand(dataset=one,successes=y,trials=n); proc print; * proc logistic;class study trt/param=ref; * model y=study trt; run;