Alex Reinhart
February 27, 2014
\[ g(\mathbb{E}[Y_{ij}]) = \beta_0 + X_{ij} \beta + U_{ij} + \epsilon_{ij} \]
\( g \): link function:
CID.Gibbs(rep.mat[,c("user1","user2")],
components=list(LSM(2)),
burnin=800, thin=10,
draws=2000)