Details
The simulation is performed by simulating new response vectors from the fitted probabilities without
clustering, and comparing the maximized log likelihoods. The maximizations are performed by
profiling out the grouping factor. It is a very fast procedure, compared to glm, when the grouping
factor has many levels.