nci = scan("nci.data") nci=matrix(nci,nrow=6830,byrow=T) dist(nci) # this makes a 6830x6830 matrix of distances between genes. I wasn't able to cluster this on my mac. Let me know if you have success. nci = t(nci) # this is 64 x 6830 hclust(dist(nci,method="average")) # as one possibility; this reproduces Figure 14.12 type=scan("ncinames.Rdata", what="character") rownames(nci) = type plot(hclust(dist(nci),method="average"),cex=0.6)