Well, the random statements give rise to the estimates for the gamma vector, so for any instance you would have to do something like:   /* Calculate the equivalent of a pred(blup) */ numblup = exp (X'beta + Z'gamma); prob_ind = numblup/(1+numblup); /* or calculate the equivalent of a red(noblup) */ numnoblup = exp(X'beta); prob_mean = numnoblup/(1+numnoblup); It appears that you get the beta's and gamma's from a dataset created with the outpost= option in the PROC BLIMM statement.  Without having one of those at hand, I am sort of stuck. I think that you could then calculate residuals by merging the calculated values back against the original dataset (or better SQL, where you can also calculate the residuals in a single call).   SteveDenham  ... View more

  Or if you want, comment out certain put statements using a macro variable which comments out lines when it holds an asterisk; %let debug=*;      data _null_;    input idno name $ startwght 3.;     put name @;    &debug      if startwght ne . then           put @15 startwght;       &debug                      else put;      datalines; 032 David 180 049 Amelia 145 126 Monica 219 Alan 210 ; run; ... View more