[R-sig-ME] Predicted probabilites with CIs for multilevel logistic regression with prior weights

[Sam Crawley]

Hello all,

I am doing a multilevel binomial logistic regression using lme4, and the survey data I am using requires weights to be used. I would like to calculate various predicted probabilities with confidence intervals based on the estimated model. The predict function obviously doesn't give me standard errors, and the recommended method to get these is to use the bootMer function.

However, in my case, the weights provided are not integers, and the bootMer function exits with an error if the weights are not integers (I raised a GitHub issue about this, and was pointed to this list: https://github.com/lme4/lme4/issues/524 ).

So my question is, what is the best way to calculate the predicted probabilities (with confidence intervals) in my case?

I would appreciate any help you can give me, and I'm happy to provide more details if required.

Thanks,
Sam Crawley.