My z-value was imprecise. Note below for those who are interested.

Thanks to Matthew, Jason and Peter for off-list replies.

data one; input exposure disease count; cards; 1 1 25 1 0 50 0 1 25 0 0 75 ; run;

data two; input #1 @5 a #2 @5 b #3 @5 c #4 @5 d; or=(a*d)/(b*c); lb=exp(log(or)-probit(0.975)*(1/a+1/b+1/c+1/d)**0.5); ub=exp(log(or)+probit(0.975)*(1/a+1/b+1/c+1/d)**0.5); cards; 1 1 25 1 0 50 0 1 25 0 0 75 ; run;

proc freq data=one order=data noprint; tables exposure*disease / cmh; output out=one1 (keep=_mhor_ l_mhor u_mhor) cmh; weight count; run;

data _null_; merge two one1; put a= b= c= d= /; put or= _mhor_= /; put @5 lb= @25 ub= /; put l_mhor= u_mhor= /; run;

Kevin

------------------------ CONFIDENTIALITY NOTICE: This e-mail transmission, and any documents, files or previous e-mail messages attached to it may contain proprietary, privileged or confidential information. If you are not an intended recipient, or a person responsible for delivering it to the intended recipient, you are hereby notified that any disclosure, copying, distribution or use of any of the information contained in or attached to this transmission is STRICTLY PROHIBITED. If you have received this transmission in error, please immediately notify me by reply e-mail and destroy the original transmission and its attachments without saving them in any manner. -----------------------