Hierarchical Generalized Linear Models In R