Discussion :

[benoitbarnouin]

Bonjour,

Je débute sur R et je ne réussi pas à trouver de solution à mon problème :

J'ai une matrice avec 16 variables (b1, b2, b3, ... b16) et 54 lignes (je n'en met ici que 2 pour illustrer). Je souhaite créer de nouvelles variables sous la forme b1/b2 puis b1+b2 et enfin (b1/b2)+b3.
Sauf que je souhaite générer toutes les possibilités pour mes 16 variables (donc plusieurs centaines de possibilités.

Voilà à quoi ressemble mes données :

b1	b2	b3	b4	b5	b6	b7	b8	b9	b10	b11	b12	b13	b14	b15	b16
0.486	0.379	0.943	0.021	0.736	0.159	0.320	0.803	0.002	0.123	0.665	0.653	0.312	0.468	0.786	0.051
0.331	0.495	0.782	0.174	0.695	0.478	0.951	0.367	0.265	0.688	0.079	0.718	0.325	0.069	0.478	0.782

J'arriver manuellement à générer les nouvelles variables pour tous les ratios possibles par exemple, mais cela me force à rédiger 256 calculs puisqu'il existe 256 ratios possible avec 16 variables.

Je souhaite donc automatiser ce traitement et il me semble qu'il faut utiliser une boucle for pour y arriver. j'ai également regardé la fonction mutate() mais je n'ai pas réussi à comprendre comment l'appliquer à ma situation.

Merci pour votre aide,

Cordialement,

Benoit

[tototode]

Bonjour,

Tout dépend des combinaisons que tu cherches à établir. Pour les sommes, b1+b2 c'est la même chose que b2+b1 donc on n'a pas besoin de ces doublons ?
Si oui alors tu peux faire :

Code :

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

1
2
3
# si tes données sont dans un objet tab
somme <- combn(16, 2, function(x) rowSums(tab[,x]))
colnames(somme) <- combn(16, 2, function(x) paste(x, collapse = "_"))

Par contre pour les ratios est-ce qu'il te faut tous les ratios possibles ? b1/b2 et b2/b1 ?

Code :

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

1
2
3
ratio <- combn(16, 2, function(x) tab[,x[1]]/tab[,x[2]])
ratio <- cbind(ratio, 1/ratio) # pour avoir les deux cas de figures
colnames(ratio) <- c(colnames(somme), sub("(^\\d+)_(\\d+$)", "\\2_\\1", colnames(somme), perl = T))

Après pour ce qui est de ratios divisés par des variables il faut préciser un peu plus quelles sont les combinaisons que tu souhaites.
cdlt