在R中的向量中查找连续值

时间:2014-07-19 06:21:48

标签: r

如果我有这样的矢量:

dat <- c(1,2,3,4,5,19,20,21,56,80,81,92)

如何将其分解为列表:

[[1]]
1 2 3 4 5

[[2]]
19 20 21

[[3]]
56

[[4]]
80 81

[[5]]
92

2 个答案:

答案 0 :(得分:11)

只需将splitdiff结合使用:

> split(dat, cumsum(c(1, diff(dat) != 1)))
$`1`
[1] 1 2 3 4 5

$`2`
[1] 19 20 21

$`3`
[1] 56

$`4`
[1] 80 81

$`5`
[1] 92

不完全是你要求的,但&#34; R.utils&#34;包有几个相关的有趣功能:

library(R.utils)
seqToIntervals(dat)
#      from to
# [1,]    1  5
# [2,]   19 21
# [3,]   56 56
# [4,]   80 81
# [5,]   92 92
seqToHumanReadable(dat)
# [1] "1-5, 19-21, 56, 80-81, 92"

答案 1 :(得分:1)

我认为Robert Krzyzanowski是正确的。所以这里有一个整合,涉及将矢量放入一个tibble(数据框)。

library(tidyverse)
# library(dplyr)
# library(tidyr)

df <- c(1,2,3,4,5,19,20,21,56,80,81,92) %>%
  tibble(dat = .)

# using lag()
df %>%
  group_by(seq_id = cumsum(dat != lag(dat) + 1 | is.na(dat != lag(dat) + 1)) %>%
  nest()

# using diff()
df %>%
  group_by(seq_id = cumsum(c(1, diff(dat)) != 1)) %>%
  nest()

当然,您不需要将结果组嵌套到列表列中,而是可以执行某种摘要操作。