Combinaisons + disjonctif complet

Version imprimable

16/03/2012, 16h30
A. D.

Combinaisons + disjonctif complet
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:

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:

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.
16/03/2012, 19h16
lilly74
Bonsoir,
pour l'instant, je n'ai trouvé que du bidouillage :roll:

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

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" :oops:

Bon week-end :)
Lilly

Bonsoir,

Voici une autre proposition:

Code:

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:

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:

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:

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é.

19/03/2012, 14h46
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.