Discussion :

[A. D.]

Bonjour à tous,

Je vais essayer d'exposer mon problème le plus clairement possible...
Je suis à la recherche d'une fonction R qui me permettrait de passer de la matrice P à la matrice Z suivantes :

Code :

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

1
2
3
4
5
6
7
8
9
P<-t(combn(4,2))
P
     [,1] [,2]
[1,]    1    2
[2,]    1    3
[3,]    1    4
[4,]    2    3
[5,]    2    4
[6,]    3    4

Code :

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

1
2
3
4
5
6
7
8
9
 
Z
     [,1] [,2] [3,] [4,]
[1,]   1    1    0    0
[2,]   1    0    1    0
[3,]   1    0    0    1
[4,]   0    1    1    0
[5,]   0    1    0    1
[6,]   0    0    1    1

Ou bien qui me permettrait d'obtenir directement la matrice Z...


Pour le moment j'ai réussi à obtenir ce que je souhaite en passant la matrice P au format disjonctif complet (fontion "acm.disjonctif" du package "ade4" par exemple) et avec une bidouille pour sommer les colonnes ayant le même nom.

Voilà, si vous souhaitez que je donne des précisions, n'hésitez pas.
Merci d'avance


Cordialement,

A.D.

[lilly74]

Bonsoir,
pour l'instant, je n'ai trouvé que du bidouillage

Voici quand même le code, si cela peut vous "inspirer":

Code :

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

1
2
3
4
5
6
7
for (i in 1:4){
 ind <- which(P == i)
 test <- which(ind > 6)
 ind[test] <- ifelse(is.na(test), test, ind[test] %% 6)
 ind[ind == 0] <- 6 # car 12 modulo 6 vaut 0...
 mat[ind , i] <- 1
}

Oui, oui, je n'ai pas levé la boucle "for"

Bon week-end
Lilly

[hoccha]

Bonsoir,

Voici une autre proposition:

Code :

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

1
2
3
4
5
6
7
8
9
# Cette fonction produit un vecteur de taille "lengthVect" où les éléments en position "pos" sont à "1".
fct <- function(pos,lengthVect){
   res <- rep(0,lengthVect)
   for(ii in pos) res[ii] <- 1
   return(res)
}
 
P<-t(combn(4,2))
Z<- t(apply(P,1,fct,lengthVect=4))

Et voici des exemples de résultats:

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
 
P<-t(combn(4,2))
Z<- t(apply(P,1,fct,lengthVect=4))
> print(P)
     [,1] [,2]
[1,]    1    2
[2,]    1    3
[3,]    1    4
[4,]    2    3
[5,]    2    4
[6,]    3    4
> print(Z)
     [,1] [,2] [,3] [,4]
[1,]    1    1    0    0
[2,]    1    0    1    0
[3,]    1    0    0    1
[4,]    0    1    1    0
[5,]    0    1    0    1
[6,]    0    0    1    1

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
 
P<-t(combn(5,3))
Z<- t(apply(P,1,fct,lengthVect=5))
> print(P)
      [,1] [,2] [,3]
 [1,]    1    2    3
 [2,]    1    2    4
 [3,]    1    2    5
 [4,]    1    3    4
 [5,]    1    3    5
 [6,]    1    4    5
 [7,]    2    3    4
 [8,]    2    3    5
 [9,]    2    4    5
[10,]    3    4    5
> print(Z)
      [,1] [,2] [,3] [,4] [,5]
 [1,]    1    1    1    0    0
 [2,]    1    1    0    1    0
 [3,]    1    1    0    0    1
 [4,]    1    0    1    1    0
 [5,]    1    0    1    0    1
 [6,]    1    0    0    1    1
 [7,]    0    1    1    1    0
 [8,]    0    1    1    0    1
 [9,]    0    1    0    1    1
[10,]    0    0    1    1    1

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
> P<-t(combn(3,2))
> Z<- t(apply(P,1,fct,lengthVect=3))
> print(P)
     [,1] [,2]
[1,]    1    2
[2,]    1    3
[3,]    2    3
> print(Z)
     [,1] [,2] [,3]
[1,]    1    1    0
[2,]    1    0    1
[3,]    0    1    1

J'espère que je vous ai aidé.

[A. D.]

Merci à tous les deux pour vos réponses !

Ce sont effectivement de bonnes alternatives à ce que j'avais fait (où je me complique quand même un peu la vie en passant par la fonction "acm.disjonctif" du package "ade4").
Je me demandais juste si une telle fonction existait déjà dans un des nombreux packages de R, mais peut-être pas...

Encore merci !


Cordialement,

A.D.