5

私の樹形図はひどく醜く、読めない寸前で、通常は次のようになります。

ここに画像の説明を入力

library(TraMineR)
library(cluster)
data(biofam)
lab <- c("P","L","M","LM","C","LC","LMC","D")
biofam.seq <- seqdef(biofam[1:500,10:25], states=lab)

ccost <- seqsubm(biofam.seq, method = "CONSTANT", cval = 2, with.missing=TRUE)
sequences.OM <- seqdist(biofam.seq, method = "OM", norm= TRUE, sm = ccost,     
with.missing=TRUE)

clusterward <- agnes(sequences.OM, diss = TRUE, method = "ward")
plot(clusterward, which.plots = 2)

私が作成したいのは、次のようなものです。つまり、ラベルのサイズを慎重に制御して実際に表示できる丸い樹形図を意味します。

ここに画像の説明を入力

Rでこれを達成するにはどうすればよいですか?

4

1 に答える 1

9

次の解決策は最適ではないかもしれませんが、試してみる価値があります。

library(ape)
CL1 <- as.hclust(clusterward)
CL2 <- as.phylo(CL1)
plot(CL2, type="fan", cex=0.5)

ここに画像の説明を入力

主な問題は明らかにオブジェクトがまだ多すぎるため、ラベルが多すぎるという事実です。ラベルをオフにするには、引数を使用しますshow.tip.label=FALSE。次のようにして、デバイス全体を占有するマージンを取り除くこともできますno.margin=TRUE

plot(CL2, type="fan", show.tip.label=FALSE, no.margin=TRUE)

ここに画像の説明を入力

于 2014-01-28T16:14:38.607 に答える