Discussion :

[Hutys]

Bonjour,

J'ai un vecteur x de longueur n dont les composantes sont des chaînes de caractère, et je souhaiterais créer une matrice mat de taille n x n telle que :

mat[i,j] =1 si x[i]==x[n]
et 0 sinon

Donc une matrice symétrique. J'ai décidé de ne remplir que la partie inférieure.
J'utilise le package Matrix, comme la matrice aura beaucoup de 0.

Voici le code que j'ai fait pour le moment :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
mat <- Matrix(0,nrow=length(x),ncol=length(x))
 
for (i in 1:length(x)){
  for (j in 1:length(x)){
    if (i>j)
    {
      mat[i,j] <- x[i]==x[j]
    }
  }
}

Je compte le faire pour plusieurs vecteurs x, et le souci est que n peut être très grand (entre 10 000 et 50 000). J'aimerais donc savoir s'il était possible d'optimiser ce code, ou s'il y avait un autre moyen de remplir la matrice sans boucle ?

Merci de votre aide !

[VonVelten]

Bonjour,

pour gagner du temps, il est déjà possible de se passer du test (i<j). Pour cela, on l'intègre avec une borne de boucle sur j dépendant de i :

Code R :

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

1
2
3
4
5
6
7
mat <- Matrix(0,nrow=length(x),ncol=length(x))
 
for (i in 2:length(x)){
  for (j in 1:(i-1)){
      mat[i,j] <- x[i]==x[j]
  }
}

cordialement

VV

[Hutys]

Merci de votre réponse, en effet ce code est plus rapide, mais malheureusement pas assez pour n grand.

J'ai donc essayé, pour chaque élément unique de x, de voir les différents indices où il apparaît dans mon vecteur, puis assigner à ma matrice, pour chaque combinaison de couple correspondante, la valeur 1.

Voici le code :

Code :

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

1
2
3
4
5
6
7
8
9
10
mat <- Matrix(0,nrow=length(x),ncol=length(x))
 
for (str in unique(x){
  if (length(which(x==str))>1){
    idx <- combn(which(x==str),2)
    for (k in 1:dim(idx)[2]){
      mat[idx[1,k],idx[2,k]] <- 1
    }
  }
}

Cela semble fonctionner correctement et est (beaucoup plus) rapide que ma proposition de base.

[faubry]

Créer la matrice par :

Code :

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

kronecker( x, t( x), "==")

A voir comment mettre le résultat dans une Matrix