このリンクを使用して、カテゴリごとに色付きのラベルが付いた素敵な樹状図をプロットしています。2番目の答えは、以下のコードを使用するこのリンク( Rの水平デンドログラムのクラスターの周りのツリーカットと長方形)で私が見ているものです:
d <- dist(t(mat[,3:ncol(mat)]), method = "euclidean")
H.fit <- hclust(d, method="ward")
groups <- cutree(H.fit, k=16) # cut tree into clusters
hcdata<- dendro_data(H.fit, type="rectangle")
hcdata$labels <- merge(x = hcdata$labels, y = pm_gtex_comb, by.x = "label", by.y = "sample",all=TRUE)
ggplot() +
geom_segment(data=segment(hcdata), aes(x=x, y=y, xend=xend, yend=yend)) +
geom_text(data=label(hcdata), aes(x, y, label=label, hjust=0, color=cluster),
size=3) +
geom_rect(data=rect, aes(xmin=X1-.3, xmax=X2+.3, ymin=0, ymax=ymax),
color="red", fill=NA)+
geom_hline(yintercept=0.33, color="blue")+
coord_flip() + scale_y_reverse(expand=c(0.2, 0)) +
theme_dendro()
16個のクラスターがあり、145個のラベルがあるため、クラスターの一部を切り取って、それらのいくつかにのみ焦点を合わせ/切り取り/ズームインしたいので、いくつかのクラスターのみを表示できるようにします.何か方法はありますか?これは hclust オブジェクトにあります。これは、145 個のラベルで図が乱雑になるため、見栄えを良くするためだけのものです。ラベルごとに色を付けたいので、ggdendro はかなり適していると思います。
たとえば、このリンクでは、3) デンドログラムの拡大http://gastonsanchez.com/blog/how-to/2012/10/03/Dendrograms.htmlを見れば、