Data Analysis Using Hierarchical Generalized Linear Models with R