Performance of a mixed effects logistic regression model for binary outcomes with unequal cluster size