Hierarchical modeling for clustered data