simdata <- read.table("hclust_in.txt") #read in similarity data dissim <- function(x)(c(36-x)) #convert to dissimilarity data dissimdata <- as.dist(apply(simdata[1:20], c(1,2), dissim)) names = simdata[,21] #read in object names clus <- hclust(dissimdata, "ward") #hierarchical clustering, ward method plot(clus, labels=names) #plot tree