Discussion :

[skadum]

Bonjour à tous,

J'ai une matrice numérique dont je souhaite pondérer (ou normaliser) toutes les valeurs par une ligne de celle-ci.
J'applique donc la formule:

Code :

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

1
2
3
4
##d'abord j'enregistre ma ligne 
ma_ligne = ma_matrice["ID",]
##et puis je divise ma matrice par les valeurs de la ligne
ma_matrice_pondérée = ma_matrice/ma_ligne

Jusque là pas de problème, je vérifie les ratios sur une ligne quelconque je trouve exactement les mêmes valeurs, tout va bien. Seulement si je vérifie la même ligne par laquelle j'ai pondérée, je ne trouve pas une série de 1 mais des valeurs autour de 1. Y aurait-il des histoires d'arrondis différents entre objets différents ? Bref je sèche mais je ne vois pas pourquoi cela fait ça surtout qu'en faisant: ma_matrice["ID",]/ma_matrice["ID",] , je trouve bien une série de 1.

Bref, est-ce que quelqu'un aurait déjà observé ce problème ?

Merci d'avance !

[A. D.]

Bonjour,

Tout d'abord, vous pourriez nous fournir un extrait de votre jeu de données afin qu'il nous soit plus facile de voir où se situe le problème.

Ensuite, est-ce que l'identifiant "ID" correspond bien au nom d'une ligne (et pas d'une colonne)?

Enfin, pour moi le calcul est faux car la division d'une matrice par un vecteur se fait (par défaut) par colonne !
En effet :

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
 
A<-matrix(c(1,2,5,2,4,8,5,10,15),byrow=TRUE,ncol=3)
A
     [,1] [,2] [,3]
[1,]    1    2    5
[2,]    2    4    8
[3,]    5   10   15
 
ma_ligne<-A[1,]
ma_ligne
[1] 1 2 5
 
mat_ponderee<-A/ma_ligne
mat_ponderee
     [,1] [,2] [,3]
[1,]    1    2    5
[2,]    1    2    4
[3,]    1    2    3

Bonne continuation


Cordialement,

A.D.

[skadum]

Bon j'ai finalement opté pour la boucle for, c'est beaucoup plus long mais au moins je suis sûr du résultat.
On m'a dit d'utiliser la fonction apply aussi mais pas testé encore, à suivre donc...

Merci pour la réponse A.D.

[A. D.]

Bonjour skadum,

Contente d'avoir pu aider
Tu pourrais peut-être poster le code que tu as écrit, ça pourrait servir à d'autres !

Sinon oui, je te conseille fortement de jeter un oeil du côté des fonctions de type "apply", c'est typiquement le genre de cas où elles peuvent être utiles. N'hésite pas à poser des questions à ce sujet si tu en as.

Bonne continuation


Cordialement,

A.D.

[lilly74]

Bon j'ai finalement opté pour la boucle for, c'est beaucoup plus long mais au moins je suis sûr du résultat.
On m'a dit d'utiliser la fonction apply aussi mais pas testé encore, à suivre donc...

Merci pour la réponse A.D.

Bonjour,
n'y a-t-il pas plus simple en passant par les transpositions?

Code :

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

1
2
3
4
5
6
7
8
A<-matrix(c(1,2,5,2,4,8,5,10,15),byrow=TRUE,ncol=3)
ma_ligne<-A[1,]
t(t(A)/ ma_ligne)
 
    [,1] [,2] [,3]
[1,]    1    1  1.0
[2,]    2    2  1.6
[3,]    5    5  3.0

N'hésitez pas à me dire si je suis à côté de la plaque

Bonne journée

[A. D.]

Bonjour,
n'y a-t-il pas plus simple en passant par les transpositions?

Bien vu ! Moi j'étais passée par un apply ET une transposition
Cela dit, c'est bien d'apprendre à se servir des fonctions de type "apply".

Bonne journée


Cordialement,

A.D.