Above is what the G matrix dimension and G matrix estimate is. But there is something wrong as I mentioned in the first post. ... View more

@binfang wrote: This means there should only be 3 rows of estimate in the covariance output, but there are 4 rows instead, which is not what I expected since one level should be used as reference.    I don't think this is correct, I think there are covariance estimates for every level of your 4-level categorical variable. But honestly, it has been a very long time since I did PROC GLIMMIX with DIST=MULTINOMIAL. Maybe someone knows. ... View more