Discussion :

[khaoulagammoudi]

Bonjour,
J'ai mon tableau "ubiquity":

Code :

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

1
2
3
4
5
6
7
8
9
10
 
product	pays 1	pays 2	pays 3	pays 4
p1	1	1	1	1
p2	1	1	0	1
p3	0	0	1	0
p4	0	1	0	0
p5	0	1	0	1
p6	0	0	1	0
p7	1	1	1	1
p8	0	0	0	0

Les indices 1 et 0 indiquent si p est produit par un pays donné ou non.
Je voudrais calculer, pour tout deux produits (i et k), combien de pays (j) produisent ces deux produits.
Je voudrais aussi representer le resultat sous forme d'une matrice symetrique de produits comme suit:

Code :

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

1
2
3
4
5
6
7
 
  p1	p2	p3	p4	...
p1 	1	1	0	1
p2	6	0	1	2
p3	5	3	4	3
p4	0	1	0	1
...

Pour ce faire, j'ai essayé avec le code suivant:

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
 
> ubiquity <- read.csv('D:/ubiquity.csv', header = TRUE, sep = ";")
 >M<-matrix(nrow=nrow(ubiquity),*ncol=nrow(ubiquity))
> for (i in 1:nrow(ubiquity)){
+ for (k in i+1:nrow(ubiquity)){
+     t <- 0
+ for (j in 2:ncol(ubiquity)){
+ if(ubiquity[i,j]==ubiquity[k,j]){
+ t <- t+1 
+ }
+ }
+ M[i,k] <- t
+ }
+ }
Error in if (ubiquity[i, j] == ubiquity[k, j]) { : 
  valeur manquante là où TRUE / FALSE est requis

Je sais pas si la méthode que j'ai utilisé peut me donner le résultat souhaité et j'ai pas compris aussi l'erreur affiché sous R.
J'espère que vous pouvez m'aider.
Merci infiniment.

[DeepLearning]

Bonjour,

Je ne me suis pas basé sur le code que tu as fourni pour répondre à ton problème. Voici ce que j'ai fait, est-ce la réponse à ta question ?

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
#Nombre de lignes et de colonnes de la table générée pour l'exemple
nbrcol=4
nbrlin=8
 
#Valeurs aléatoires qui construiront la matrice ubiquity 
set.seed(777);x=sample(0:1,nbrcol*nbrlin,replace=TRUE)
 
#création de ubiquity pour que la table soit similaire à la tienne
ubiquity=matrix(x,ncol=nbrcol,nrow=nbrlin)
ubiquity<-as.data.frame(ubiquity)
rownames(ubiquity)<-paste("p",1:nrow(ubiquity),sep="")
colnames(ubiquity)<-paste("pays",1:ncol(ubiquity),sep=" ")
ubiquity
 
#Création de la matrice M que tu souhaites 
M<-matrix(nrow=nrow(ubiquity),ncol=nrow(ubiquity))
 
for(i in 1:nrow(ubiquity)){
  for(j in 1:nrow(ubiquity)){
    M[i,j]=sum((ubiquity[i,]!=0)&(ubiquity[j,]!=0))
  }
}

Évidement, toi tout ce qui t'intéresse c'est la création de la matrice M. J'ai utilisé une double boucle for pour la construire, ce genre de chose est à utiliser qu'en cas d’extrême urgence sur R à cause de la lenteur des boucles. Mais disons que je n'avais pas envie de me creuser la tête avec les fonction apply()

Donc voilà ce que ça fait à la fin

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
> ubiquity
   pays 1 pays 2 pays 3 pays 4
p1      1      1      0      0
p2      0      0      1      1
p3      0      1      0      0
p4      1      1      0      1
p5      1      1      0      1
p6      0      1      1      1
p7      0      1      1      1
p8      0      0      0      1
> M
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,]    2    0    1    2    2    1    1    0
[2,]    0    2    0    1    1    2    2    1
[3,]    1    0    1    1    1    1    1    0
[4,]    2    1    1    3    3    2    2    1
[5,]    2    1    1    3    3    2    2    1
[6,]    1    2    1    2    2    3    3    1
[7,]    1    2    1    2    2    3    3    1
[8,]    0    1    0    1    1    1    1    1

C'est bien ce que tu souhaites ?

[khaoulagammoudi]

Oui c'est ça ! merci bien