R - кол-во упорядоченных последовательностей

Question

У меня есть проблема с последовательным подсчетом значения переменной.

Мои данные выглядят как этот

IDPERS value year divorce 
1 371102 married 1  0 
2 371102 divorced 2  1 
3 371102 divorced 3  1 
4 371102  NA 4  0 
5 371102 divorced 5  1 
6 371102 divorced 6  1 
7 371102  NA 7  0 
8 371102  NA 8  0 
9 371102 divorced 9  1 
10 371102 divorced 10  1 
11 371102 divorced 11  1 

То, что я хочу, чтобы подсчитать развод такой как

IDPERS value year divorce divorceEp 
1 371102 married 1  0   0 
2 371102 divorced 2  1   1 
3 371102 divorced 3  1   2 
4 371102  NA 4  0   2 
5 371102 divorced 5  1   3 
6 371102 divorced 6  1   4 
7 371102  NA 7  0   4 
8 371102  NA 8  0   4 
9 371102 divorced 9  1   5 
10 371102 divorced 10  1   6 
11 371102 divorced 11  1   7 

Любой ключ, как я мог рассчитывать правильно эту последовательность эпизодов?

library(dplyr) 
dta %>% mutate(1:sum(divorce)) 

Благодаря

dta = structure(list(IDPERS = c(371102, 371102, 371102, 371102, 371102, 
371102, 371102, 371102, 371102, 371102, 371102), value = c("married", 
"divorced", "divorced", "NA", "divorced", "divorced", "NA", "NA", 
"divorced", "divorced", "divorced"), year = 1:11, divorce = c(0, 
1, 1, 0, 1, 1, 0, 0, 1, 1, 1), divorceEp = c(0, 1, 2, 2, 3, 4, 
4, 4, 5, 6, 7)), class = "data.frame", .Names = c("IDPERS", "value", 
"year", "divorce", "divorceEp"), row.names = c(NA, -11L)) 

Answer 1

Вы можете использовать cumsum():

mutate(dta,divorceEp=cumsum(divorce)) 
    IDPERS value year divorce divorceEp 
1 371102 married 1  0   0 
2 371102 divorced 2  1   1 
3 371102 divorced 3  1   2 
4 371102  NA 4  0   2 
5 371102 divorced 5  1   3 
6 371102 divorced 6  1   4 
7 371102  NA 7  0   4 
8 371102  NA 8  0   4 
9 371102 divorced 9  1   5 
10 371102 divorced 10  1   6 
11 371102 divorced 11  1   7