In reply to wudneh ketema.

I do not have SAS IML code for fitting mixed models. I use PROC MIXED or GLIMMIX.