Optimiser remplissage d'une matrice sous condition

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:

---

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 !

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:

---

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

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:

---

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.

Créer la matrice par :

Code:

---

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

---

A voir comment mettre le résultat dans une Matrix