2

I have a multilevel logistic regression (Rasch), fit in JAGS. Specifically,

$\text{logit} (P(\text{Win})) = \alpha_p + \gamma_w + \delta_{p,w}$

The prior on $\alpha_p$ is partially pooled, $\alpha_p \sim N(\mu_\alpha,\tau_\alpha)$, and the rest of the priors are normal with mean 0. The model "works", i.e., converges over multiple chains, with sensible Gelman Rhat values.

Now, for interpretation, I'd like to factor out the intercept.

$\text{logit} (P(\text{Win})) = \mu_0 + \alpha_p + \gamma_w + \delta_{p,w}$

If I leave the prior on $\alpha_p$ as is, the model does not converge, presumably because $\mu_0$ is not identifiable with respect to $\mu_\alpha$.

The model does converge if I change the prior to $\alpha_p \sim N(0,\tau_\alpha)$. But if I do that, do I still have partial pooling in my model? That is, do the $\alpha_p$ for those $p$ who have few observations still borrow strength from the other $\alpha_p$?

Community
  • 1
Jack Tanner
  • 4,552
  • 3
  • 27
  • 39
  • *Sidenote:* As pointed out earlier, there's no interaction term between item and person parameters in a Rasch model (in other words, items are assumed to be equally discriminative). There was a good overview on fitting common IRT models with BUGS (most code should run smoothly under JAGS) in the [Journal of Statistical Software](http://www.jstatsoft.org/v36/c01/). There was also a nice workshop during one of the annual ISOQOL conference: [The Use of R and WinBUGS in Fitting Item Response Theory Models](http://idecide.mskcc.org/stats/multi/irt_winbugs.php). – chl Aug 06 '12 at 09:22

0 Answers0