Function for ISNI computation when the longitudinal/clustered binary outcome follows a GLMM.
Useful links