Binary data occurs in the special case where k=2. If there are m i observations in a subpopulation i, then the probability distribution of the number falling into the k categories y i = (y i1, y i2, .
The L matrix constructed to compute them is precisely the same as the one formed in PROC GLM. The LSMEANS statement is not available for multinomial distribution models for ordinal response data.