Discussion :

[mycaweb]

Bonjour, j'ai fais une petite routine qui permet de simuler 2 statistiques (enfin!). Et maintenant je cherche à calculer le coefficient de corrélation de Pearson. Pour cela j'ai repris ces lignes de codes :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
 
#include <math.h>
#define TINY 1.0e-20
 
void pearsn(float x[], float y[], unsigned long n, float *r, float *prob, float *z)
 
{
    float betai(float a, float b, float x);
    float erfcc(float x);
    unsigned long j;
    float yt, xt, t, df;
    float syy=0.0, sxy=0.0, sxx=0.0, ay=0.0, ax=0.0;
 
for (j=1;j<=n;j++)  {
    ax += x[j];
    ay += y[j];
    }
    ax /= n;
    ay /= n;
    for (j=1;j<=n;j++) {
        xt=x[j]-ax;
        yt=y[[j]-ay;
        sxx += xt*xt;
        syy += yt*yt;
        sxy += xt*yt;
    }
    *r=sxy/sqrt(sxx*syy);
    *z=0.5*log((1.0+(*r)+TINY)/(1.0-(*r)+TINY));
    df=n-2;
    t=(*r)*sqrt(df/((1.0-(*r)+TINY)*(1.0+(*r)+TINY)));
    *prob=betai(0.5*df, 0.5, df/(df+t*t));
}

En réalité mon problème est que je ne vois pas comment retourner le résultat de "*r" pour le reprendre plus loin dans mon programme, tout simplement . . .

[souviron34]

En réalité mon problème est que je ne vois pas comment retourner le résultat de "*r" pour le reprendre plus loin dans mon programme, tout simplement . . .

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
float x[100], y[100], r, prob, z ;
unsigned long n ;
..
pearsn( x, y, n, &r, &prob, &z);

[gwenBZH]

Même si je ne comprend pas ce que calcule ton code parce que les maths et moi, ça fait 3, je pense qu'un simple

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

return r;

qui retourne l'adresse de r (r étant un pointeur passé en paramètres de ta fonction) ou alors

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

return *r;

qui retourne la valeur qui se trouve à l'adresse r. (vu ce que tu demandes, c'est plutôt cette ligne d'ailleurs)

Seulement, si j'ai bien vu ton code, ni r, ni *r n'est modifié dans ta fonction, me trompe-je ?

(où alors j'ai rien compris à ta question)

edit: d'ailleurs j'ai pas compris la même chose que souviron34

[souviron34]

je pense qu'un simple

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

return r;

qui retourne l'adresse de r (r étant un pointeur passé en paramètres de ta fonction)

Puisque r est passé en paramètre, il ressort (on a passé son adresse)
Il est inutile, voire dangereux à très dangreux, d'en plus le retourner..

Seulement, si j'ai bien vu ton code, ni r, ni *r n'est modifié dans ta fonction, me trompe-je ?

t'as mal regardé

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

   *r=sxy/sqrt(sxx*syy);

[gwenBZH]

[QUOTE=souviron34;4271520]
Il est inutile, voire dangereux à très dangreux, d'en plus le retourner..
[QUOTE]

Oui, dans ce cas là, évidemment c'est inutile.

t'as mal regardé

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

   *r=sxy/sqrt(sxx*syy);

Arf, désolé . Et pourtant j'ai relu le code 3 fois pour voir où il était modifié.

*essuie ses lunettes, son écran et ouvre les yeux*

[diogene]

souviron34 :

Il est inutile, voire dangereux à très dangreux, d'en plus le retourner..

Inutile, soit, ça fait double emploi ! Mais dangereux, pourquoi ?

[jeroman]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
void pearsn(float x[], float y[], unsigned long n, float *r, float *prob, float *z)
{
    float betai(float a, float b, float x);
    float erfcc(float x);
 
/* ... */
}

Il ne faut pas déclarer des fonctions dans une définition de fonction.

[souviron34]

souviron34 :

Inutile, soit, ça fait double emploi ! Mais dangereux, pourquoi ?

parce que tu as 2 pointeurs pointant sur la même zone...

Alors la tentation est grande d'en libérer un et de se servir de l'autre.. ou réciproquement

Ou de faire des manips de pointeurs sur l'un et d'oublier l'autre...

Bref c'est extrêmement mélangeant dans la routine d'appel...

(le seul moyen "rationnel" de l'utiliser ainsi est d'affecter le retour au même pointeur que le paramètre, mais c'est pas clair, et ça continue à mélanger...)

Il ne faut pas déclarer des fonctions dans une définition de fonction.

D'une part ici ce n'est pas une définition mais un prototype.


Mais de toutes façons ne soyons pas aussi strict.. Cela dépend..

Un bon exemple est la fonction de tri utilisée dans qsort..

Dans beaucoup de bons exemples d'utilisation, tu trouveras le prototype (ce n'est pas la défintion) dans la routine d'appel à qsort..