запись данных временного ряда в r

Question

Я пытаюсь перекодировать существующие данные со сверхурочной структурой.Мой набор данных выглядит следующим образом:

dput(z)

structure(list(democracy = c(0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 
1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L), year.x = 1967:2008, time = c(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, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 
41, 42)), .Names = c("democracy", "year.x", "time"), row.names = 176:217, class = "data.frame")

Так что я хочу создать новую переменную, скажем, time.democ, которая принимает значение ноль, если democracy==0, но снова начинает считать период времени, начиная сот 1, если democracy ==1, до democracy==0 снова.Я собираюсь сделать это для ряда стран, но я предполагаю, что обобщение достаточно просто с использованием ddply, если однажды я получу эту функцию правильно.Какие-либо предложения?

Я хотел бы получить это:

dput(z)

structure(list(democracy = c(0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 
1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L), year.x = 1967:2008, time = c(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, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 
41, 42), new.time = c(0, 0, 0, 0, 0, 0, 0, 1, 2, 0, 0, 0, 0, 
0, 0, 0, 0, 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)), .Names = c("democracy", 
"year.x", "time", "new.time"), row.names = 176:217, class = "data.frame")

Спасибо!

Answer 1

Вы можете использовать rle в сочетании с sequence, чтобы сделать это. rle выполняет кодирование длин серий, а sequence генерирует последовательности.

z$new.time <- sequence(rle(z$democracy)$lengths)
z$new.time[z$democracy==0] <- 0

head(z, 20)

    democracy year.x time new.time
176         0   1967    1        0
177         0   1968    2        0
178         0   1969    3        0
179         0   1970    4        0
180         0   1971    5        0
181         0   1972    6        0
182         1   1973    7        1
183         1   1974    8        2
184         1   1975    9        3
185         0   1976   10        0
186         0   1977   11        0
187         0   1978   12        0
188         0   1979   13        0
189         0   1980   14        0
190         0   1981   15        0
191         0   1982   16        0
192         1   1983   17        1
193         1   1984   18        2
194         1   1985   19        3
195         1   1986   20        4

Answer 2

Спасибо за ваши ответы. Я последовал вашим предложениям и в итоге написал функцию, чтобы я мог применить ее ко всем единицам в моем (продольном) наборе данных через ddply. Я публикую это, поскольку это может помочь некоторым другим, хотя я уверен, что есть более элегантные решения:

# is a long format data frame
new.time <- function(a){
    a <- a[order(a$year.x),]
    a$new.time <- sequence(rle(a$democracy)$lengths)-1
    a$new.time[a$democracy==0] <- 0
    return(a)
}

merged1 <- ddply(merged, .(country.x), new.time)