Improving the Hosmer-Lemeshow goodness-of-fit test in large models with replicated Bernoulli trials