bhc                     Function to perform Bayesian Hierarchical
                        Clustering on a 2D array of discretised (i.e.
                        multinomial) data
