Hi, there. I am looking for a package to fit the following binomial-normal model Y_ij ~ Binomial (N_ij, P_ij) Logit(P_ij) = \beta_0i+\beta_1*X+e_ij \beta_0i ~ N(\beta_0,\sigma_b^2) e_ij ~ N(0,\sigma^2) This model has two variance components, one random effect \beta_0i and one error e_ij Can lme4 be used for this? I would appreciate a lot if some example code can be provided. Thanks. Yulei [[alternative HTML version deleted]]