les champs de la fonction "cov.rob"

bonjour,
Dans le cadre de mon travail , je dois calculer des distance robuste et pour cela j'utilise:

library(MASS)
cov.rob(x, cor = FALSE, quantile.used= ((n + p + 1) / 2)), method="mcd")

et j'obtiens un message d'erreur .
alors comment faire pour remplir correctement les champs de cette fonction??
autre question: pour une analyse multivariées correcte doit-je faire d'abord une autre analyse discriminante??
merci à tous d'avoir pris le temps de lire ma question .

Bonjour,

Quel est le message d'erreur renvoyé par R? Pouvez-vous mettre un petit extrait de vos données pour que nous puissions éventuellement tester la fonction "cov.rob"?

Personellement je ne connais pas la méthode du calcul des distances robustes, mais dans l'aide de R concernant la fonction "cov.rob", pour le paramètre "quantile.used", il est mis :
quantile.used = floor ( (n + p + 1) / 2 ) .
Peut-être que votre problème vient de l'absence de "floor" dans votre commande?...

Bonne continuation


Cordialement,

A.D.

salut ,
j'ai enfin réussi à trouver les distances robustes (le seuil c'est la limite qu'on ne veut pas dépasser par exemple on tolère que au plus 5% de la population soit hors norme )

Code:


seuil=sqrt(qchisq(0.95,p))
abline(h=seuil)

mcd.est=cov.mcd(X)

D.robuste=sqrt(mahalanobis(X,mcd.est$center,mcd.est$cov))
plot(D.robuste,xlab="Indice",ylab="DM") 
abline(h=seuil)


j'espère que cela aidera d'autres personnes
SH

Bonjour,

Comme je vous l'ai dit dans mon précédent message, je ne connais pas la méthode des distances robustes, je ne peux donc pas vous aider concernant la procédure à suivre.

Cependant, concernant votre deuxième question :

sh a écrit: quel code pourra me donner la liste des observations dont D.robuste > h (D.robuste:distance et h une valeur numérique) .

Si "D.robuste$distance" est le vecteur qui contient vos observations (numériques), alors pour trouver tous les éléments de ce vecteur qui sont supérieurs à h, vous pouvez faire comme ceci :

Code:

 D.robuste$distance [ which ( D.robuste$distance > h ) ]

Bonne continuation


Cordialement,

A.D.

merci pour votre précieuse aide.
s.H