[R-br] identificação de grupos em um cluster
Humberto Hazin
hghazin em hotmail.com
Sexta Março 2 20:37:06 BRT 2012
Olá Antonio,
A melhor solução que encontrei para resolver esse problema foi (deve
haver outra forma melhor, mais funcionou):
library(vegan)
data(dune)
dis <- vegdist(dune)
cluc <- hclust(dis, "complete")
plot(cluc)
rect.hclust(cluc, 3)
grp <- cutree(cluc, 3)
grp
namesI<-c("I")
namesII<-c("II")
namesIII<-c("III")
legend(0.4,0.98,namesIII,cex=1.5,bty="n",bg ="white",col="white")
legend(5.2,0.98,namesI,cex=1.5,bty="n",bg ='white')
legend(14.2,0.98,namesII,cex=1.5,bty="n",bg ='white')
Um abraço
Humberto Hazin
Em 3/2/2012 8:02 AM, Antonio Silva escreveu:
> Olá
>
> Gostaria de saber se há como indicar no dendrograma, na altura do
> corte, o número referente ao grupo, como indicado em cutree.
>
> Observei que em alguns dendrogramas que os grupos selecionados de
> acordo com um dado o nível de corte não seguem necessariamente uma
> ordem ascendente (grupo 3 pode não ser o terceiro da esquerda para
> direita).
>
> Seria algo semelhante ao da figura
> https://www.crops.org/images/publications/cs/42/5/1584f4.jpeg
>
> Por exemplo exemplo:
>
> library(vegan)
> data(dune)
> dis <- vegdist(dune)
> cluc <- hclust(dis, "complete")
> plot(cluc)
> rect.hclust(cluc, 3)
> grp <- cutree(cluc, 3)
> grp
>
> Notem que o grupo com os objetos 17, 19, 11 e 18, posicionado a
> esquerda no dedrograma é o grupo de número 3, e o grupo com os objetos
> 1,2, ..., 7, representado no meio do dendrograma é o de número 1.
>
> Eu gostaria que, no nível de corte, aparecesse no dendrograma, ao lado
> das barras verticais de cada cluster, a indicação 3, 1, 2.
>
> Agradeço qualquer dica. Abraços
>
> Antônio
>
> --
> Antônio Olinto Ávila da Silva
> Biólogo / Oceanógrafo
> Instituto de Pesca
> São Paulo, Brasil
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120302/c4518352/attachment.html>
Mais detalhes sobre a lista de discussão R-br