[R-br] rect.hclust()

[Mauricio Figueira]

Oi, estou tentando delimitar os grupos no Cluster analysis, porém a função rect.hclust() extrapola os limites dos retângulos, deixando-os abertos.

Os dados estão disponíveis na íntegra no endereço https://www.dropbox.com/s/8vxymeq8wkvhg1l/dados.csv

O script que estou usando segue abaixo.

require(vegan)

dados<-read.csv("dados.csv",h=T, row.names=1, sep=";")
dados
dis<-vegdist (dados, method = "euclidean")

dis
range(dis)
clust<- hclust(dis, "average")
plot(clust, main="Análise de Agrupamento", sub="Ligação pela média", 
xlab="Locais e orgãos", ylab="Dissimilaridade (Distância Euclidiana)")
cor(dis, cophenetic(clust))
rect.hclust(clust, k=2, border="red")#desenhar dois retângulos
x <- rect.hclust(clust, h=40)#cortar na altura de 40

Se alguém puder me ajudar, desde já agradeço.

Obrigado Augusto Ribas pelas sugestões de postagem.

Maurício Figueira

Mestrando em Engenharia Florestal (UFSM)
Herbário do Departamento de Ciências Florestais
 (HDCF)
Currículo: http://lattes.cnpq.br/1022600405818850
Tel.:(55) 9663-4920 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121122/d6a6fdf5/attachment.html>