对连续重复的奇数行求和并保留第一行,删除其余行。

3

我遇到了一个问题 - 如何对连续的重复奇数行求和,并删除除第一行之外的所有行。我已经搞清楚如何对连续的重复行求和并删除除第一行之外的所有行(链接:https://stackoverflow.com/a/32588960/11323232)。但是在这个项目中,我想要对连续的重复奇数行进行求和,而不是所有连续的重复行。

 ia<-c(1,1,2,NA,2,1,1,1,1,2,1,2)
 time<-c(4.5,2.4,3.6,1.5,1.2,4.9,6.4,4.4, 4.7, 7.3,2.3, 4.3)
 a<-as.data.frame(cbind(ia, time))

  a
   ia time
1   1  4.5
2   1  2.4
3   2  3.6
5   2  1.2
6   1  4.9
7   1  6.4
8   1  4.4
9   1  4.7
10  2  7.3
11  1  2.3
12  2  4.3

to 

 a
   ia time
1   1  6.9
3   2  3.6
5   2  1.2
6   1  20.4
10  2  7.3
11  1  2.3
12  2  4.3

如何编辑以下代码以达到我的目标,即将连续重复的奇数行相加,并删除除第一行以外的所有行?
 result <- a %>%
 filter(na.locf(ia) == na.locf(ia, fromLast = TRUE)) %>%
 mutate(ia = na.locf(ia)) %>%
 mutate(change = ia != lag(ia, default = FALSE)) %>%
 group_by(group = cumsum(change), ia) %>%
 # this part
 summarise(time = sum(time))

1
感谢tmfmnk和Patrik_P。我还有另一个问题。如果“time”与列表长度相同,我该怎么办?例如:time <- list(c(4.5,2), 2.4,3.6,1.5,1.2,4.9,6.4,4.4, 4.7, 7.3,2.3, 4.3)。 - hees
2个回答

1
可能的一种方法是:
a %>%
 group_by(grp = with(rle(ia), rep(seq_along(lengths), lengths))) %>%
 mutate(grp2 = ia %/% 2 == 0,
        time = sum(time)) %>%
 filter(!grp2 | (grp2 & row_number() == 1)) %>%
 ungroup() %>%
 select(-grp, -grp2)

      ia  time
  <dbl> <dbl>
1     1   6.9
2     2   3.6
3     2   1.2
4     1  20.4
5     2   7.3
6     1   2.3
7     2   4.3

0
你可以尝试使用 data.table 实现以下代码:
library(data.table)
ia <- c(1,1,2,NA,2,1,1,1,1,2,1,2)
time <- c(4.5,2.4,3.6,1.5,1.2,4.9,6.4,4.4, 4.7, 7.3,2.3, 4.3)
a <- data.table(ia, time)
a[, sum(time), by=.(ia, rleid(!ia %% 2 == 0))]

提供

##   ia rleid   V1
##1:  1     1  6.9
##2:  2     2  3.6
##3: NA     3  1.5
##4:  2     4  1.2
##5:  1     5 20.4
##6:  2     6  7.3
##7:  1     7  2.3
##8:  2     8  4.3

太好了!谢谢。嗨,patrik_P。如果“time”的长度与列表相同,我该怎么办?例如,time <- list(c(4.5,2), 2.4,3.6,1.5,1.2,4.9,6.4,4.4, 4.7, 7.3,2.3, 4.3)。 - hees
如果您将列表分配为data.table中的列,则可以执行以下操作:time2 <- list(c(4.5,2), 2.4,3.6,1.5,1.2,4.9,6.4,4.4, 4.7, 7.3,2.3, 4.3); a2 <- data.table(ia, time2); a2[, .(SUMtime = sum(unlist(time2))), by=.(ia, rleid(!ia %% 2 == 0))] - Patrik_P
再次感谢。可能有些误导。您能否看一下以下问题?我已经发布了一个新的问题,描述了更新请求。https://dev59.com/0bPma4cB1Zd3GeqPlBhC - hees

网页内容由stack overflow 提供, 点击上面的
可以查看英文原文,
原文链接