Hi all,    I keep getting the following error messages when I run my code: "R-side random effects are not supported for the multinomial distribution." and "Least square means are not supported for the multinomial distribution"   My data set has milk yield from from 32 cows over 21 days (DIM = days in milk i.e. the number of days they have been lactating). Treatment is if cows were given hay or no hay. Cows were blocked by parity and expected calving date. My code is:   PROC GLIMMIX DATA = haytrial; CLASS treatment DIM treatment block; MODEL milkyield=treatment*DIM/DDFM=kr2; RANDOM block; RANDOM intercept/ SUBJECT=cow(block); RANDOM DIM/RESIDUAL SUBJECT=cow(block) TYPE=ar(1); LSMEANS treatment*DIM; RUN;   I am completely new to SAS so any help is much appreciated! ... View more