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. O script que estou usando segue abaixo. #dissimilaridade EUCLIDIANA 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. 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
Não da para reproduzir seus exemplo, vc não forneceu nada chamado dados. Não esqueça de avisar também os pacotes que esta usando. vegdist é do pacote vegan certo? http://www.leg.ufpr.br/doku.php/software:rbr?&#guia_de_postagem_da_r-br_em_c... Da uma olhada no guia de postagem, tente usar o comando dput(dados) e poste aqui, pelo menos parte dos seus dados ou dados criados de exemplo. Assim será possivel reproduzir sua figura e assim tentar ajudar. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/<http://augustoribas.heliohost.org> Lattes: http://lattes.cnpq.br/7355685961127056
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
participantes (2)
-
Augusto Ribas -
Mauricio Figueira