如果连续因子与R一起编辑/合并行

时间:2016-07-15 08:10:14

标签: r date merge factors

我有这个样本

    data <- structure(list(id = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 
2L, 2L), .Label = c("a", "b"), class = "factor"), minTime = structure(c(1L, 
2L, 3L, 4L, 5L, 7L, 6L, 8L, 9L, 10L), .Label = c("2014-06-06 07:39:50", 
"2014-06-07 02:24:32", "2014-06-07 15:14:29", "2014-06-07 2:29", 
"2014-06-08 5:40", "2014-06-18 17:54:42", "2014-06-18 2:45", 
"2014-06-19 02:37:53", "2014-06-19 19:15", "2014-06-19 22:15"
), class = "factor"), maxTime = structure(c(1L, 3L, 4L, 2L, 5L, 
6L, 7L, 8L, 9L, 10L), .Label = c("2014-06-07 01:41:31", "2014-06-07 10:01", 
"2014-06-07 14:44:08", "2014-06-07 22:31:02", "2014-06-08 5:50", 
"2014-06-18 2:50", "2014-06-19 01:49:05", "2014-06-19 18:51:36", 
"2014-06-19 20:15", "2014-06-19 23:15"), class = "factor"), duration.minutes = c(NA, 
740L, 437L, 452L, NA, NA, 474L, 974L, 4062L, 353L), event = structure(c(1L, 
4L, 4L, 2L, 2L, 1L, 4L, 4L, 3L, 4L), .Label = c("enter", "exit", 
"stop", "trip"), class = "factor")), .Names = c("id", "minTime", 
"maxTime", "duration.minutes", "event"), class = "data.frame", row.names = c(NA, 
-10L))

我想合并活动&#34;旅行&#34;如果它们对于每个id是连续的。

如果已经合并了id a和b的连续出行:

  • 编辑了maxTime以获得旅行的全部内容
  • duration.minutes也被编辑。

例如,此输出:

dput(output.wanted)
structure(list(id = structure(c(1L, 1L, 1L, 2L, 2L, 1L, 2L), .Label = c("a", 
"b"), class = "factor"), minTime = structure(c(5L, 6L, 7L, 2L, 
1L, 3L, 4L), .Label = c("18.6.2014 17:54", "18.6.2014 2:45", 
"19.6.2014 19:15", "19.6.2014 22:15", "6.6.2014 7:39", "7.6.2014 2:24", 
"7.6.2014 2:29"), class = "factor"), maxTime = structure(c(5L, 
7L, 6L, 1L, 2L, 3L, 4L), .Label = c("18.6.2014 2:50", "19.6.2014 18:51", 
"19.6.2014 20:15", "19.6.2014 23:15", "7.6.2014 1:41", "7.6.2014 10:01", 
"7.6.2014 22:31"), class = "factor"), duration.minutes = c(NA, 
1177L, 452L, NA, 1448L, 4062L, 353L), event = structure(c(1L, 
4L, 2L, 1L, 4L, 3L, 4L), .Label = c("enter", "exit", "stop", 
"trip"), class = "factor")), .Names = c("id", "minTime", "maxTime", 
"duration.minutes", "event"), class = "data.frame", row.names = c(NA, 
-7L))

关键是我想要合并所有事件(旅行或停止),如果它们是相同的id连续。但如果停止或旅行是唯一的,它会保持不变

我一直试图用group_by和mutate来做,但我有点失落......

- 我找到了解决方案,但事件stops

仍存在一个小问题

以下是一个示例:

> dput(total)
structure(list(Ship = c(205482000, 205482000, 205482000, 205482000, 
205482000, 205482000, 205482000, 205482000, 205482000, 205482000, 
205482000, 205482000), minTime = structure(c(1401570241, 1401969219, 
1401981860, 1402052108, 1402768362, 1402772602, 1402841443, 1402855773, 
1403056361, 1403278916, 1403290856, 1403367735), class = c("POSIXct", 
"POSIXt")), maxTime = structure(c(1401966639, 1401980399, 1402051849, 
1402056430, 1402772313, 1402839873, 1402852433, 1403052550, 1403276355, 
1403289496, 1403367596, 1403371285), class = c("POSIXct", "POSIXt"
)), duration.minutes = structure(c(6607, 186, 1166, NA, NA, 1121, 
183, 3280, 3667, 176, 1279, NA), class = "difftime", units = "mins"), 
    event = c("stop", "trip", "trip", "exit", "enter", "trip", 
    "trip", "stop", "stop", "trip", "trip", "exit"), dist.sailed = c(NA, 
    50254.2034817555, 349194.108518887, NA, NA, 347816.081064252, 
    50035.8859874946, NA, NA, 49982.687612038, 351978.737678528, 
    NA)), .Names = c("Ship", "minTime", "maxTime", "duration.minutes", 
"event", "dist.sailed"), class = "data.frame", row.names = c(4L, 
5L, 6L, 2L, 1L, 7L, 8L, 9L, 10L, 11L, 12L, 3L))

以下代码未生成停止的完整持续时间,但仅生成NA:

 total <- total %>% 
    group_by(Ship) %>% 
    mutate(new_id = data.table::rleid(event)) %>% 
    group_by(event, new_id, Ship) %>% 
    mutate(duration.minutes = ifelse(event == 'trip', sum(duration.minutes), duration.minutes), maxTime = tail(maxTime, 1))%>% 
    mutate(duration.minutes = ifelse(event == 'stop', sum(duration.minutes), duration.minutes), maxTime = tail(maxTime, 1))%>% 
    mutate(dist.sailed = ifelse(event == 'trip', sum(dist.sailed), dist.sailed), dist.sailed = tail(dist.sailed, 1)) %>% 
    filter(!duplicated(duration.minutes)) %>% 
    select(-new_id)

我将mutate(duration.minutes = ifelse(event == 'stop', sum(duration.minutes), duration.minutes), maxTime = tail(maxTime, 1))%>%添加到@Sotos的代码中

2 个答案:

答案 0 :(得分:2)

这有点乱,

library(dplyr)
data %>% 
   group_by(id) %>% 
   mutate(new_id = data.table::rleid(event)) %>% 
   group_by(event, new_id, id) %>% 
   mutate(duration.minutes = ifelse(event == 'trip', sum(duration.minutes), duration.minutes), maxTime = tail(maxTime, 1)) %>% 
   filter(!duplicated(duration.minutes)) %>% 
   select(-new_id) 

#new_id     id             minTime             maxTime duration.minutes  event
#   <int> <fctr>              <fctr>              <fctr>            <int> <fctr>
#1      1      a 2014-06-06 07:39:50 2014-06-07 01:41:31               NA  enter
#2      2      a 2014-06-07 02:24:32 2014-06-07 22:31:02             1177   trip
#3      3      a     2014-06-07 2:29     2014-06-08 5:50              452   exit
#4      3      a     2014-06-08 5:40     2014-06-08 5:50               NA   exit
#5      1      b     2014-06-18 2:45     2014-06-18 2:50               NA  enter
#6      2      b 2014-06-18 17:54:42 2014-06-19 18:51:36             1448   trip
#7      3      b    2014-06-19 19:15    2014-06-19 20:15             4062   stop
#8      4      b    2014-06-19 22:15    2014-06-19 23:15              353   trip

如果我们为stop而不是trip运行相同的代码,我们会得到以下结果

#new_id      Ship             minTime             maxTime duration.minutes event dist.sailed
#    <int>     <dbl>              <time>              <time>            <dbl> <chr>       <dbl>
#1       1 205482000 2014-06-01 00:04:01 2014-06-05 14:10:39             6607  stop          NA
#2       2 205482000 2014-06-05 14:53:39 2014-06-06 13:50:49              186  trip    50254.20
#3       2 205482000 2014-06-05 18:24:20 2014-06-06 13:50:49             1166  trip   349194.11
#4       3 205482000 2014-06-06 13:55:08 2014-06-06 15:07:10               NA  exit          NA
#5       4 205482000 2014-06-14 20:52:42 2014-06-14 21:58:33               NA enter          NA
#6       5 205482000 2014-06-14 22:03:22 2014-06-15 20:13:53             1121  trip   347816.08
#7       5 205482000 2014-06-15 17:10:43 2014-06-15 20:13:53              183  trip    50035.89
#8       6 205482000 2014-06-15 21:09:33 2014-06-20 17:59:15             6947  stop          NA
#9       7 205482000 2014-06-20 18:41:56 2014-06-21 19:19:56              176  trip    49982.69
#10      7 205482000 2014-06-20 22:00:56 2014-06-21 19:19:56             1279  trip   351978.74
#11      8 205482000 2014-06-21 19:22:15 2014-06-21 20:21:25               NA  exit          NA

答案 1 :(得分:0)

如果您使用mutate包,则必须使用some.df %>% group_by(event) %>% summarise(sum.maxTime = sum(maxTime), sum.duration = sum(duration.minutes)) ,而不是Mutate以某种方式合并数据。像

summarise

{{1}}添加或更改data.frame中的列,但结果将始终与原始行具有相同的行数。 {{1}}之后的行号将是数据集中组的数量。

希望这有帮助