Após um post anterior (Rótulo e dendrograma de folhas coloridas em r) Eu tenho uma questão de acompanhamento.

Minhas perguntas são semelhantes ao post mencionado, mas eu imagino que isso pode ser feito usando o macaco (por exemplo,plot(as.phylo(fit), type="fan", labelCol) como tem mais tipo de filogenia.

As perguntas postadas foram:

Como posso mostrar os códigos de grupo no rótulo da folha (em vez do número da amostra)?

Desejo atribuir uma cor a cada grupo de código e colorir o rótulo da folha de acordo com ele (pode acontecer de eles não estarem no mesmo clade e por isso eu possa encontrar mais informações)?

E o exemplo de código é:

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))

## make unique rownames (equal rownames are not allowed)
rownames(sample) <- make.unique(groupCodes)

colorCodes <- c(A="red", B="green", C="blue", D="yellow")


## perform clustering
distSamples <- dist(sample)
hc <- hclust(distSamples)

## function to set label color
labelCol <- function(x) {
  if (is.leaf(x)) {
    ## fetch label
    label <- attr(x, "label")
    code <- substr(label, 1, 1)
    ## use the following line to reset the label to one letter code
    # attr(x, "label") <- code
    attr(x, "nodePar") <- list(lab.col=colorCodes[code])
  }
  return(x)
}

## apply labelCol on all nodes of the dendrogram
d <- dendrapply(as.dendrogram(hc), labelCol)

plot(d)