Mahalanobis para o dendograma!

Pessoal como fazer para construir um dendrograma com a matriz de distancia de Mahalanobis? Já que Mahalanobis não esta na função dist! Já dei uma pesquisada na lista mesmo e no google, mas não consegui! x <- matrix(rnorm(100*3), ncol = 3)D2 <- mahalanobis(x, colMeans(x), cov(x))hc <- hclust(D2,method="ward.D") # erro in if (is.na(n) || n > 65536L) stop("size cannot be NA nor exceed 65536") : valor ausente onde TRUE/FALSE necessário André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

André, essa função "mahalanobis" não retorna uma matriz, retorna somente um vetor. será que não é isso, link abaixo, que vc quer? http://stats.stackexchange.com/questions/65705/pairwise-mahalanobis-distance Em 24/06/2015 15:41, Andre Oliveira escreveu:
Pessoal como fazer para construir um dendrograma com a matriz de distancia de Mahalanobis? Já que Mahalanobis não esta na função dist! Já dei uma pesquisada na lista mesmo e no google, mas não consegui!
x <- matrix(rnorm(100*3), ncol = 3) D2 <- mahalanobis(x, colMeans(x), cov(x)) hc <- hclust(D2,method="ward.D") # erro in if (is.na(n) || n > 65536L) stop("size cannot be NA nor exceed 65536") : valor ausente onde TRUE/FALSE necessário
André Oliveira Souza.
--- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

Andre, tenta assim install.packages("StatMatch") library(StatMatch) mahalanobis.dist(iris[1:6,1:4]) Em 24 de junho de 2015 18:54, David Feitosa <contato@davidfeitosa.com> escreveu:
Desconheço essa "mahalanobis", mas eu precisei criar um dendograma uma vez para plotar um gráfico e criei baseado num list e na documentação.
Em 24 de junho de 2015 16:28, Cleber Borges <klebyn@yahoo.com.br> escreveu:
mahalanobis
Atenciosamente,
David F.
_______________________________________________ R-br mailing list R-br@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.
-- Bruno C. Vidigal

Pois é Bruno eu já havia pesquisado aqui e achei este pacote também, mas comparando com outros não batem matriz que procuro! x1 <- c(131.37, 132.37, 134.47, 135.50, 136.17) x2 <- c(133.60, 132.70, 133.80, 132.30, 130.33)x3 <- c(99.17, 99.07, 96.03, 94.53, 93.50)x4 <- c(50.53, 50.23, 50.57, 51.97, 51.37)x <- cbind(x1, x2, x3, x4)require(StatMatch)D2=mahalanobis.dist(x)D2d2=as.matrix(D2) d2d=as.dist(d2) d 1 2 3 4 2 2.828427 3 2.828427 2.828427 4 2.828427 2.828427 2.828427 5 2.828427 2.828427 2.828427 2.828427 require(biotools) x1 <- c(131.37, 132.37, 134.47, 135.50, 136.17)x2 <- c(133.60, 132.70, 133.80, 132.30, 130.33)x3 <- c(99.17, 99.07, 96.03, 94.53, 93.50)x4 <- c(50.53, 50.23, 50.57, 51.97, 51.37)x <- cbind(x1, x2, x3, x4)Cov <- matrix(c(21.112,0.038,0.078,2.01, 0.038,23.486,5.2,2.844,0.078,5.2,24.18,1.134, 2.01,2.844,1.134,10.154), 4, 4) # Não entendi como foi calculada tal matriz Cov tentei var(x) e não bateu com Cov! D2=D2.dist(x, Cov)D2 1 2 3 4 2 0.09103867 3 0.90451337 0.73086040 4 1.88407643 1.59779871 0.44281866 5 2.69996452 2.17905736 0.91121818 0.21984612 André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES Em Quarta-feira, 24 de Junho de 2015 21:07, bruno vidigal <caetano.vidigal@gmail.com> escreveu: Andre, tenta assim install.packages("StatMatch")library(StatMatch) mahalanobis.dist(iris[1:6,1:4]) Em 24 de junho de 2015 18:54, David Feitosa <contato@davidfeitosa.com> escreveu: Desconheço essa "mahalanobis",mas eu precisei criar um dendograma uma vez para plotar um gráfico ecriei baseado num list e na documentação. Em 24 de junho de 2015 16:28, Cleber Borges <klebyn@yahoo.com.br> escreveu: mahalanobis Atenciosamente, David F. _______________________________________________ R-br mailing list R-br@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. -- Bruno C. Vidigal _______________________________________________ R-br mailing list R-br@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.
participantes (4)
-
Andre Oliveira
-
bruno vidigal
-
Cleber Borges
-
David Feitosa