My multinomial logistic model is written something like this.

1) How can I check if this model makes use of the efficient samplers in the glm module?

2) If it does not, is there a way to rewrite this model in JAGS to takes advantage of the glm module?

# K is number of categories # o is index of observations # exp.eta is the exponentiated linear component of a GLM # pr is the probability of the observation generating a response in category k # y is the observed response for (k in 1:K) { exp.eta[o,k] <- exp(inprod(X[o,],beta)) pr[o,k] <- exp.eta[o,k] / sum(exp.eta[o,1:K]) } y[o] ~ dcat(pr[o,1:K])