Discussion :

[doc funky]

Bonjour,

je ne saurais même pas donner un nom approprié à ma manipulation donc dsl si il y a un doublon...

j'ai une colonne 'booléen' et j'aimerai créer une 2e colonne qui identifie un groupe jusqu'au changement d'état:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
c1     c2
0     1
0     1
0     1
1     2
1     2
0     3
0     3
0     3
0     3
1     4
1     4
1     4

Avec mon reflexe de codeur de type 'vba' je vais tourner en boucle sur la colonne et i=1 'while' c1 == 0 puis i=i+1 puis c2=i 'while' c1=1 etc...

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
 
i <- 0
l <- 1
while (l<nrow(gpe_data)) {
 
  i <- i+1
 
  while (gpe_data[l,1] == 0) {
    gpe_data[l,2] <- i
    l <- l+1
  }
 
  i <- i+1
 
  while (gpe_data[l,1] == 1) {
    gpe_data[l,2] <- i
    l <- l+1
  }
 
 
}

bon ben c'est bien ça marche mais j'ai a peu près 1 500 000 lignes, c'est donc très long a exécuter et surtout pas la meilleure approche en R. Auriez-vous une idée plus approprié?


Merci pour votre aide. (j’espère que c'est clair et lisible...)

[Sengar]

Bonjour,

Code :

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

c2<-cumsum(c(1,abs(diff(c1))))

[doc funky]

Bonjour Sengar,

merci pour ta réponse (brève mais efficace). Peux-tu expliquer en détail ce que fait cette fonction stp? ça fonctionne mais je ne sais pas comment...


Allé, encore merci, j'ai ce qu'il faut en tout cas.

[Sengar]

diff() permet de faire le differentiel d'un vecteur.

En clair : 1 1 0 1 0 0 devient 0 -1 1 -1 0

abs() prend la valeur absolue : du coup, abs(diff()) sur ton vecteur prend forcément ou la valeur 0 (si il n'y a pas eu de changement d'état) ou 1 (s'il vient d'y a avoir un changement d'état).
0 -1 1 -1 0 devient 0 1 1 1 0


Je rajoute un élement 1 grâce à c(1,) en première position de ce vecteur pour dire qu'on partira de l'état numérotté 1.
0 1 1 1 0 devient 1 0 1 1 1 0


Une fois que tu as repéré les instants de changements d'état, tu peux utiliser cumsum() qui fait la somme cumulée d'un vecteur :
1 0 1 1 1 0 devient 1 1 2 3 4 4, ce qui est le résultat attendu.

[doc funky]

génial, je vais méditer la dessus... MERCI!