Debutant bloqué pour calcul manuel de correlation (TP débutant)

Bonjour à tous, je débute en R et lors d'un TP visant à nous faire comprendre quelles formules entrent en compte lors des utilisations des commandes (exemple : cov(), mean(), cor()) afin que nous voyons les differences induites par R par rapport aux formules de base.

Nous devons donc comparer les version "formulée" et les version "commandes R" (je parles des résultats) afin d'étudier les variations tout se passe bien jusqu'a la dernière question, celle de cor()...

Mon problème est le suivant :

j'utilise cars, voici mes étapes :

Code:

---

1 2 3 4

data(cars) attach(car) x <- speed /* je met la colonne de la vitesse dans x, et celle de la distance de freinage correspondante dans y */ y <- dist

---

lorsque j’essaye de rentrer manuellement la formule de la corrélation suivante :

Pièce jointe 422140

(que l'on me demande de mettre une variable r) je tape la ligne suivante

Code:

---

r = (sum(length(x)-mean(x))*sum(length(y)*mean(y)))/sqrt(sum(length(x)-(mean(x))^2))*sqrt(sum(length(y)-(mean(y))^2))

---

je ne sais pas si mes parenthèses sont male mises (j'ai essayé plus d'une heure, plusieurs combinaisons) mais il en retourne toujours l'output suivant :
Pièce jointe 422144

si vous savez ou j'me plante, je suis preneur, je commence a me démotivé (d'autant plus je sais tres bien qu'on a inventé cor() pour cette raison justement) !

Bien cordialement,

Pierre

Bonjour,

Première astuce : découper votre code pour tester pas à pas plutôt que d'exécuter directement une grosse ligne de code comme ça ;)

Quelques pistes :

• Y a-t-il des NA (valeurs manquantes) dans vos données x et y ?
• Le résultat de la soustraction passée dans la racine carrée ("sqrt") est-il / peut-il être négatif ? Si oui, cela va coincer...

HTH ! :D


Cordialement,


A.D.

Es-tu vraiment sûr que ton code correspond à la formule ? Début de réponse : que vient faire la fonction length dans le code ?

Bonjour,

Il y a des erreurs dans la formule, le length je ne comprends pas , est ce que sum(x*y) = sum(x)*sum(y) ?

je te laisse réfléchir un peu si t'arrives pas, tu peux demander la formule.

Bon courage.

:bravo:

1) Il y a tellement d'erreurs basiques dans le code de Aztharos que je lui conseillerais d'abord de relire son code et, éventuellement, d’approfondir en général les règles de la programmation avant de lui donner la solution.

2) Pour réagir à l'objectif de l'exercice "Bonjour à tous, je débute en R et lors d'un TP visant à nous faire comprendre quelles formules entrent en compte lors des utilisations des commandes (exemple : cov(), mean(), cor()) afin que nous voyons les differences induites par R par rapport aux formules de base", je suis en désaccord avec l'intitulé. R n'induit pas des différences par rapport aux formules de base, R est un outil statistique donc propose les formulations de base pour l'estimation statistique de ces valeurs à partir d'un échantillon de taille finie, pour une population supposée infinie (ou de très nombreux ordres de grandeurs en taille par rapport à celle de l'échantillon) et non les formulations de base pour une population finie en incluant alors tous les individus de la population. J'espère donc que Aztharos a simplifié l'intitulé car autrement, je le considère quant à moi, comme erroné.