A hierarchical modeling approach for clustering probability density functions