split cluster analysis into two functions; added code to produce files with cophenetic matrices