bcluster function

Wrapper function for b-cluster analysis