added clustering to main.R; add clustering heatmap in clustering.R