pamk function

Partitioning around medoids with estimation of number of clusters