Probleme pour calcul de distance mahalanobis

[Varin Thibault]

Bonjour, j'ai besoin de calculer une matrice de distance avec la fonction mahalanobis du package stats.
J'ai donc effectue les etapes suivantes (pour tester la commande) :

mat<-matrix(sample(c(1:1000),size=1000,replace=F),100,100)
mahalanobis(mat, colMeans(mat), cov(mat))

R m affiche le message d erreur suivant :

Error in solve.default(cov, ...) :
Lapack routine dgesv: system is exactly singular

Est-ce que par hasard quelqu un comprend la raison de ce message?

Merci d'avance!
Cordialement
Thibault

[François Bonnot]

Bonjour,,

Lapack routine dgesv: system is exactly singular

Est-ce que par hasard quelqu un comprend la raison de ce message?

Oui : le système est exactement singulier :)

cov(mat) ne peut être inversible puisque mat n'est pas de plein rang car les mêmes colonnes sont répétées.

Donc (par exemple) :

Code : Tout sélectionner

mat<-matrix(sample(c(1:1000),size=1000,replace=F),100,10)

FB

[Varin Thibault]

Bonjour, merci pour la reponse,
autant pour moi...
Je vais essayer avec des donnees reelles pour voir si tout se passe bien!
Merci!!!
Thibault