*Prime minister approval study; data Data1; do ID=1 to 794; approve=1; time=0; output; approve=1; time=1; output; end; do ID=795 to 944; approve=1; time=0; output; approve=0; time=1; output; end; do ID=945 to 1030; approve=0; time=0; output; approve=1; time=1; output; end; do ID=1031 to 1600; approve=0; time=0; output; approve=0; time=1; output; end; proc logistic data=Data1; strata ID; model approve(event='1')=time; run; *Diabetes case control study; data Data1; do ID=1 to 9; case=1; diab=1; output; case=0; diab=1; output; end; do ID=10 to 25; case=1; diab=0; output; case=0; diab=1; output; end; do ID=26 to 62; case=1; diab=1; output; case=0; diab=0; output; end; do ID=63 to 144; case=1; diab=0; output; case=0; diab=0; output; end; proc logistic data=Data1; strata ID; model case(event='1')=diab; run; *Endometrial Cancer case control study; data d1; do ID=1 to 63; do Outcome = 1 to 0 by -1; input Gall Hyper @@; output; end; end; datalines; 0 0 0 0 0 0 0 0 0 1 0 1 0 0 1 0 1 0 0 1 0 1 0 0 1 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 1 0 0 1 1 0 1 0 1 0 0 1 0 1 0 0 0 0 1 1 0 0 1 1 0 0 0 1 0 1 0 0 0 0 1 1 0 1 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 1 0 0 0 1 1 0 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 0 0 0 1 1 1 1 0 0 0 1 0 1 0 0 0 1 0 1 0 1 0 1 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 1 1 0 0 0 1 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 1 0 1 0 1 0 0 0 1 0 0 1 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 1 1 0 0 1 0 1 0 0 1 0 0 1 0 0 0 ; proc logistic data=d1; strata ID; model outcome(event='1')= Gall Hyper; run;