Hello,

I’m fitting a multinomial logit mode to an outcome with 4 categories, and I’m estimating K-1 intercepts. I’m also specifying p[K]=0 to restrict the 4th category to a fixed value. The outcome categories are labeled 1-4 in stan data and it looks like Stan uses the last K=4 category as a reference. I get three coefficients for alpha - a1, a2, a3. How are these ordered - for example is a1 an intercept for K=1, or is a1 the next closest level to the reference, which would be K=3?

Related, is there any reference material on how Stan handles categorical outcomes/reference values and how/if one can modify that?