В R, принимая различия между значениями в двух строках на основе значений в других строках, с конкретными требованиями

Question

У меня есть несколько рядов (1000) данных с несколькими предметами, к которым у меня есть довольно простой вопрос, но я не совсем уверен, как мне ответить на него.Ниже приведен пример того, как выглядит мой фрейм данных.

     Subject  Time Current.State Transition.State Transition.Event L2 L1 C R2 R1 FoodCup
14       KM6   435            16               15                0  0  0 0  0  0       0
15       KM6   455            15                4                0  0  0 0  0  0       0
16       KM6   648             4                7                3  0  0 1  0  0       0
17       KM6   658             7               14                0  0  0 0  0  0       0
18       KM6   691            14               16                8  0  0 0  0  0       1
19       KM6   698            16                0                0  0  0 0  0  0       0
20       KM6   721            16               15                0  0  0 0  0  0       0
21       KM6   741            15                4                0  0  0 0  0  0       0
22       KM6   758             4                0                0  0  0 0  0  0       1
23       KM6   762             4                0                0  0  0 0  0  0       0
24       KM6   810             4                7                6  0  0 0  0  1       0
25       KM6   814             7                0                0  0  0 0  0  0       0
26       KM6   815             7                0                0  0  0 1  0  0       0
27       KM6   819             7                0                0  0  0 1  0  0       0
28       KM6   820             7               14                0  0  0 0  0  0       0
29       KM6   821            14                0                0  0  0 0  0  0       0
30       KM6   822            14                0                0  0  0 1  0  0       0
31       KM6   824            14                0                0  0  0 0  0  0       0
32       KM6   829            14                0                0  0  0 1  0  0       0
33       KM6   862            14               16                8  0  0 0  0  0       1
34       KM6   863            16                0                0  0  0 0  0  0       1

Задача, которую я ищу, состоит в следующем:

Сколько Time проходит между Transition.State == 7 иFoodCup == 1, но только для первого FoodCup == 1 с момента последнего Transition.State == 7.

Например, я хочу иметь возможность взять время в строке 16 (648), поскольку в этой строке есть Transition.State == 7, затем возьмите время (691) из строки 18, поскольку оно является первым FoodCup == 1 после строки 16, а затем получите значение этой разницы во времени (691-648).Но я хочу полностью игнорировать строку 22, даже если она имеет FoodCup == 1, поскольку над ней не было Transition.State == 7 "непосредственно".

Требуемым выходным значением будет новый data.frame с каждым случаем этого с разницей во времени на Subject (здесь не показано, но есть фактические значения df с)

Надеюсь, это имеет смысл.

Спасибо!

Answer 1

Может быть, что-то вроде следующего будет делать то, что вы хотите.

fun <- function(DF, want = 7){
  st <- which(DF[['Transition.State']] == want)
  fc <- which(DF[['FoodCup']] == 1)
  i <- findInterval(fc, st)
  i <- i[is.finite(i)][1]
  DF[['Time']][fc[i]] - DF[['Time']][st[i]]
}

fun(df1)
#[1] 43

fun(df1, 14)
#[1] 33

Данные в формате dput.

df1 <-
structure(list(Subject = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), 
.Label = "KM6", class = "factor"), Time = c(435L, 455L, 648L, 
658L, 691L, 698L, 721L, 741L, 
    758L, 762L, 810L, 814L, 815L, 819L, 820L, 821L, 822L, 824L, 
    829L, 862L, 863L), Current.State = c(16L, 15L, 4L, 7L, 14L, 
    16L, 16L, 15L, 4L, 4L, 4L, 7L, 7L, 7L, 7L, 14L, 14L, 14L, 
    14L, 14L, 16L), Transition.State = c(15L, 4L, 7L, 14L, 16L, 
    0L, 15L, 4L, 0L, 0L, 7L, 0L, 0L, 0L, 14L, 0L, 0L, 0L, 0L, 
    16L, 0L), Transition.Event = c(0L, 0L, 3L, 0L, 8L, 0L, 0L, 
    0L, 0L, 0L, 6L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 8L, 0L), 
    L2 = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), L1 = c(0L, 0L, 0L, 0L, 0L, 
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 
    0L), C = c(0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 
    1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L), R2 = c(0L, 0L, 0L, 0L, 
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 
    0L, 0L), R1 = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), FoodCup = c(0L, 
    0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 
    0L, 0L, 0L, 1L, 1L)), class = "data.frame", row.names = c("14", 
"15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", 
"26", "27", "28", "29", "30", "31", "32", "33", "34"))