################################################### ### chunk number 1: ################################################### data(USArrests) h = hclust(dist(USArrests)) plot(h) ################################################### ### chunk number 2: ################################################### heatmap(as.matrix(USArrests)) ################################################### ### chunk number 3: ################################################### library(ctc) r2xcluster(USArrests,file='USArrests_xcluster.txt') ################################################### ### chunk number 4: ################################################### r2cluster(USArrests,file='USArrests_xcluster.txt') ################################################### ### chunk number 5: ################################################### hr = hclust(dist(USArrests)) hc = hclust(dist(t(USArrests))) ################################################### ### chunk number 6: ################################################### write(hc2Newick(hr),file='hclust.newick') ################################################### ### chunk number 7: ################################################### r2atr(hc,file="cluster.atr") r2gtr(hr,file="cluster.gtr") r2cdt(hr,hc,USArrests ,file="cluster.cdt") ################################################### ### chunk number 8: ################################################### hclust2treeview(USArrests,file="cluster.cdt")