I fear I am missing something.  So you want confidence bounds on 'result' at various times for future repetitions at the same random effects?  What about shifting to PROC GLIMMIX, where you can use the OUTPUT statement--   output out=out1 pred(blup ilink)=predicted lcl(blup ilink)=lower ucl(blup ilink)=upper;   This should give the blup's and confidence intervals on the original scale.   Steve Denham ... View more