[R] extract rownames by cluster order
라벨:
Informatics
이메일로 전송BlogThis!Facebook에서 공유
data = data.matrix(data)
distance = dist(data)
cluster = hclust(distance, method="ward")
dendrogram = as.dendrogram(cluster)
Rowv = rowMeans(data, na.rm = T)
dendrogram = reorder(dendrogram, Rowv)
## Produce the heatmap from the calculated dendrogram.
## Don't allow it to re-order rows because we have already re-ordered them above.
reorderfun = function(d,w) { d }
png("heatmap.png", res=150, height=22,width=17,units="in")
heatmap(data,col=colors,breaks=breaks,scale="none",Colv=NA,Rowv=dendrogram,labRow=NA, reorderfun=reorderfun)
dev.off()
## Re-order the original data using the computed dendrogram
rowInd = rev(order.dendrogram(dendrogram))
di = dim(data)
nc = di[2L]
nr = di[1L]
colInd = 1L:nc
data_ordered <- data-blogger-escaped-col.names="T)" data-blogger-escaped-colind="" data-blogger-escaped-data="" data-blogger-escaped-data_ordered="" data-blogger-escaped-pre="" data-blogger-escaped-quote="F," data-blogger-escaped-row.names="T," data-blogger-escaped-rowind="" data-blogger-escaped-rows.txt="" data-blogger-escaped-sep="\t" data-blogger-escaped-write.table="">
이메일로 전송BlogThis!Facebook에서 공유
라벨:
Informatics
-