proc genmod data=resp; class id; model y=sm t t*t / link=logit dist=bin type3; repeated subject=id / type=un corrw; run; proc glimmix data=resp method=laplace; class id; model y=sm t / link=logit dist=bin s; random intercept t / type=un subject=id; run; proc glimmix data=resp method=laplace; class id; model y=sm t t*t / link=logit dist=bin s; random intercept / type=un subject=id; run; data resp; input s y1 y2 y3 y4 count; do i=1 to count; id+1; *increment id by 1 for every new record read into input buffer; y=y1; sm=s; t=7; output; y=y2; sm=s; t=8; output; y=y3; sm=s; t=9; output; y=y4; sm=s; t=10; output; end; drop s y1-y4 count; datalines; 0 0 0 0 0 237 0 0 0 0 1 10 0 0 0 1 0 15 0 0 0 1 1 4 0 0 1 0 0 16 0 0 1 0 1 2 0 0 1 1 0 7 0 0 1 1 1 3 0 1 0 0 0 24 0 1 0 0 1 3 0 1 0 1 0 3 0 1 0 1 1 2 0 1 1 0 0 6 0 1 1 0 1 2 0 1 1 1 0 5 0 1 1 1 1 11 1 0 0 0 0 118 1 0 0 0 1 6 1 0 0 1 0 8 1 0 0 1 1 2 1 0 1 0 0 11 1 0 1 0 1 1 1 0 1 1 0 6 1 0 1 1 1 4 1 1 0 0 0 7 1 1 0 0 1 3 1 1 0 1 0 3 1 1 0 1 1 1 1 1 1 0 0 4 1 1 1 0 1 2 1 1 1 1 0 4 1 1 1 1 1 7 ;