HGLMs for analysis of correlated non-normal data