17

私のサンプルには5つのグループコードがあります(サンプル名/種/などとして機能しますが、繰り返します)。

したがって、ヘルプが素晴らしいという2つの問題があります。

  • (サンプル番号の代わりに) リーフ ラベルにグループ コードを表示するにはどうすればよいですか?

  • 各コード グループに色を割り当て、それに応じてリーフ ラベルに色を付けたいと思います (それらが同じクレードに属さず、それによって詳細情報を見つけることができる場合があります)。

私のスクリプトでそうすることは可能ですか (ape または ggdendro):

sample<-read.table("C:/.../DOutput.txt", header=F, sep="")
groupCodes <- sample[,1]
sample2<-sample[,2:100] 
d <- dist(sample2, method = "euclidean")  
fit <- hclust(d, method="ward")
plot(as.phylo(fit), type="fan") 
ggdendrogram(fit, theme_dendro=FALSE)  

read.table を置き換えるランダムなデータフレーム:

sample = data.frame(matrix(floor(abs(rnorm(20000)*100)),ncol=200))
groupCodes <- c(rep("A",25), rep("B",25), rep("C",25), rep("D",25)) # fixed error
sample2 <- data.frame(cbind(groupCodes), sample) 
4

2 に答える 2