Les posteurs les plus actifs de la semaine
Aytan
les champs de la fonction "cov.rob" Vote_lcaples champs de la fonction "cov.rob" Voting_barles champs de la fonction "cov.rob" Vote_rcap 
c@ssoulet
les champs de la fonction "cov.rob" Vote_lcaples champs de la fonction "cov.rob" Voting_barles champs de la fonction "cov.rob" Vote_rcap 
Adddd
les champs de la fonction "cov.rob" Vote_lcaples champs de la fonction "cov.rob" Voting_barles champs de la fonction "cov.rob" Vote_rcap 
Eric Wajnberg
les champs de la fonction "cov.rob" Vote_lcaples champs de la fonction "cov.rob" Voting_barles champs de la fonction "cov.rob" Vote_rcap 
gg
les champs de la fonction "cov.rob" Vote_lcaples champs de la fonction "cov.rob" Voting_barles champs de la fonction "cov.rob" Vote_rcap 


les champs de la fonction "cov.rob"

Aller en bas

les champs de la fonction "cov.rob" Empty les champs de la fonction "cov.rob"

Message par sh le Ven 22 Avr 2011 - 8:30

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 .

sh

Nombre de messages : 24
Date d'inscription : 13/04/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

les champs de la fonction "cov.rob" Empty Re: les champs de la fonction "cov.rob"

Message par A.D. le Ven 22 Avr 2011 - 8:44

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 Smile


Cordialement,

A.D.

A.D.

Nombre de messages : 305
Age : 31
Localisation : Nantes
Date d'inscription : 02/12/2009

Voir le profil de l'utilisateur http://www.dacta.fr

Revenir en haut Aller en bas

les champs de la fonction "cov.rob" Empty Re: les champs de la fonction "cov.rob"

Message par sh le Sam 23 Avr 2011 - 21:38

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 Smile
SH


Dernière édition par sh le Ven 29 Avr 2011 - 14:06, édité 1 fois

sh

Nombre de messages : 24
Date d'inscription : 13/04/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

les champs de la fonction "cov.rob" Empty Re: les champs de la fonction "cov.rob"

Message par A.D. le Dim 24 Avr 2011 - 7:10

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 Smile


Cordialement,

A.D.

A.D.

Nombre de messages : 305
Age : 31
Localisation : Nantes
Date d'inscription : 02/12/2009

Voir le profil de l'utilisateur http://www.dacta.fr

Revenir en haut Aller en bas

les champs de la fonction "cov.rob" Empty Re: les champs de la fonction "cov.rob"

Message par sh le Dim 24 Avr 2011 - 8:15

merci pour votre précieuse aide.
s.H

sh

Nombre de messages : 24
Date d'inscription : 13/04/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

les champs de la fonction "cov.rob" Empty Re: les champs de la fonction "cov.rob"

Message par Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum