qqplot

Bonjour,
Quelqu'un saurait-il faire un qqplot avec un logiciel autre que SAS et R ?
Mon problème : je ne trouve pas le même qqplot avec R et SAS !
Les données :
108.9163237
109.739369
110.9739369
105.4421769
107.4829932
109.5238095
109.1397849
110.7526882
107.5268817
110.3603604
110.3603604
109.009009
108.7121212
110.9848485
108.7121212
109.6666667
109.6666667
107.6666667
110.9739369
111.3854595
110.5624143
117.6870748
111.5646259
109.5238095
113.9784946
115.5913978
115.5913978
113.0630631
110.3603604
106.3063063
108.7121212
112.1212121
113.2575758
112.6666667
117.6666667
131.6666667
110.6995885
108.6419753
113.1687243
108.1632653
108.1632653
106.122449
111.2903226
109.6774194
109.6774194
112.1621622
108.1081081
112.1621622
110.2272727
109.0909091
110.2272727
111
112
112

Qqplot R :
22501
22502

Réponses

  • Salut Steven,

    Voici le résultat de Matlab :
    22503
  • Pourquoi ne cherches-tu pas dans la doc ce que font R ou SAS quand tu leur demandes un qqplot ?
  • Merci Kuja.
    TT>> y'a pas 36 manières de faire un qqplot, si ?
  • @SN : visiblement ils ne tracent pas la même droite.
  • Je ne comprend pas ta démarche. Tu constates une différence et tu dis ensuite qu'il n'y a qu'une manière de le faire. J'ai dû raté un truc.
  • Je te demandais s'il y a d'autres manières. La doc SAS ne dit rien de spécial, donc c'est censé être le qqplot classique.
  • Mais pourquoi demander s'il y a plusieurs manières quand tu pars de la constatation qu'il y a plusieurs manières !? A moins que tu ne soupçonnes un bug quelque part ? Tu as essayé avec d'autres données ? Tu as calculé à la main l'équation de ta droite préférée ? Tu as regardé le code (je ne sais plus si c'est lisible) ? Quelle est cette droite qui serait classique dans les qqplots ?
  • Salut Steven,

    j'arrive à retrouver tes deux graphiques avec SAS. Apparemment le deuxième est tracé suivant les quantiles théoriques d'ordre $\frac{i}{n+1}$ et le premier d'ordre $m(i) = \begin{cases} 1 - m(n) & i = 1\\ \\
    \dfrac{i - 0.3175}{n + 0.365} & i = 2, 3, \ldots, n-1\\ \\
    0.5^{1/n} & i = n.\end{cases}$.

    Pour ce qui est de la droite tracée, je peux juste te dire que SAS estime l'ordonnée à l'origine et la pente par les estimations de la moyenne et de la variance de l'échantillon.

    Cordialement,
    bd

    Edit : Il y a effectivement plusieurs méthodes de construction d'un qq-plot
  • Apparemment c'est SAS qui donne le bon, car voici la fonction de répartition théorique de la gaussienne de moyenne mu-chapeau et d’écart-type sigma-chapeau et la fonction de répartition empirique de la série de données

    22504
    xx.jpg 23.6K
  • Voilà comment R trace la droite, d'où vient cette idée ??
    ‘qqline’ adds a line to a normal quantile-quantile plot which passes through the first and third quartiles.
  • Après un rapide calcul et étant donnée le peu de poids que semblent avoir les quelques observations aberrantes de la fin du graphique, je dirais que SAS estime les coefficient de la droite qu'il trace grâce au estimateur des moments de la loi normale et R estime les coefficients de la droite qu'il trace par régression.

    Edit : Raté
  • Content de voir que cela devient un peu plus scientifique par ici ;)
  • J'ai toujours grandement préféré la comparaison de la fonction de répartition théorique avec l'empirique, comme dans mon dernier graphique, pluôt que le qqplot ; et suite à cette aventure je préfère encore plus !
  • Voilà comment obtenir le même graphique que SAS avec R :
    mu <- mean(x)
    sigma <- sd(x)
    qqnorm(x)
    abline(mu,sigma)
    
  • Steven Neutral écrivait : http://www.les-mathematiques.net/phorum/read.php?13,728798,728798#msg-728798

    Peut-on avoir le code R utilisé ?
    Merci ^^

    [Inutile de répéter un message illisible ! Un lien suffit. AD]
Connectez-vous ou Inscrivez-vous pour répondre.