将序列中的值分配给一组连续的行，使某些行为空

Question

我试图将几个连续的行分组(并赋予它们相同的值)，同时将一些行留空(当不满足某个条件时)。

我的数据是位置（xy坐标）、测量日期/时间以及测量之间的时间跨度。有些简化，它们看起来是这样的：

ID   X     Y      Time    Span
1    3445  7671   0:00    -
2    3312  7677   4:00    4
3    3309  7680   12:00   8
4    3299  7681   16:00   4
5    3243  7655   20:00   4
6    3222  7612   4:00    8
7    3260  7633   0:00    4
8    3254  7641   8:00    8
9    3230  7612   0:00    16
10   3203  7656   4:00    4
11   3202  7678   8:00    4
12   3159  7609   20:00   12
...

我想为在4小时内测量的每个位置序列分配一个值，并使我的数据看起来像这样：

ID   X     Y      Time    Span  Sequence
1    3445  7671   0:00    -     -
2    3312  7677   4:00    4     1
3    3309  7680   12:00   8     NA
4    3299  7681   16:00   4     2
5    3243  7655   20:00   4     2
6    3222  7612   4:00    8     NA
7    3260  7633   0:00    4     3
8    3254  7641   8:00    8     NA
9    3230  7612   0:00    16    NA
10   3203  7656   4:00    4     4
11   3202  7678   8:00    4     4
12   3159  7609   20:00   12    NA

我尝试了几种带有循环“for”加上“ifelse”条件的算法，例如：

Sequence <- for (i in 1:max(ID)) {
ifelse (Span <= 4, i+1, "NA")
}

没有任何运气。我知道我的尝试是不正确的，但是我的编程技能非常基础，我在网络上没有发现任何类似的问题。

任何想法都将不胜感激！

Answer 1

count = 0
for (i in 1:max(ID)) {
      Sequence[i] = ifelse(Span[i] <= 4, count <- count+1, NA)
}

Answer 2

这是使用rle和rep的另一种选择。我们将假设您的＜code＞数据。帧被命名为“test”。

首先，初始化“序列”列，用NA填充。

test$Sequence <- NA

其次，指定要匹配的条件，在本例中为 test$Span == 4。

x <- test$Span == 4

第三，使用 rle 输出（长度和值）的组合来获取序列中每个新运行的发生次数。

spanSeq <- rle(x)$lengths[rle(x)$values == TRUE]

最后，使用 rep，并将 times 参数设置为在步骤 3 中获得的结果。根据 test$Span == 4 匹配的索引，对 test$Sequence 的所需值进行子集化，并将它们替换为新序列。

test$Sequence[x] <- rep(seq_along(spanSeq), times = spanSeq)
test
#    ID    X    Y  Time Span Sequence
# 1   1 3445 7671  0:00    -       NA
# 2   2 3312 7677  4:00    4        1
# 3   3 3309 7680 12:00    8       NA
# 4   4 3299 7681 16:00    4        2
# 5   5 3243 7655 20:00    4        2
# 6   6 3222 7612  4:00    8       NA
# 7   7 3260 7633  0:00    4        3
# 8   8 3254 7641  8:00    8       NA
# 9   9 3230 7612  0:00   16       NA
# 10 10 3203 7656  4:00    4        4
# 11 11 3202 7678  8:00    4        4
# 12 12 3159 7609 20:00   12       NA

一旦您了解了所涉及的步骤，您还可以直接使用within（）来完成此操作。下面将给出相同的结果：

within(test, {
  Sequence <- NA
  spanSeq <- rle(Span == 4)$lengths[rle(Span == 4)$values == TRUE]
  Sequence[Span == 4] <- rep(seq_along(spanSeq), times = spanSeq)
  rm(spanSeq)
})

Answer 3

这是一个长长的一行：

ifelse(x <- DF$Span == 4, cumsum(c(head(x, 1), tail(x, -1) - head(x, -1) == 1)), NA)
# [1] NA  1 NA  2  2 NA  3 NA NA  4  4 NA

解释：

< li> x是一个TRUE/FALSE向量，表示< code>Span是< code>4。 < li> tail(x，-1)是编写< code>x[2:length(x)]的一种安全方式 < li> head(x，-1)是编写< code > x[1:(length(x)-1)] 的一种安全方式 < li> tail(x，-1) - head(x，-1) == 1是一个真/假向量，表示我们从< code>Span！= 4到< code>Span == 4。 < li >由于上面的向量比< code>x短一个元素，所以我在它前面添加了< code>head(x，1)。< code>head(x，1)是编写< code>x[1]的一种安全方式。 < li >然后我取< code>cumsum，这样它就把向量TRUE/FALSE转换成一个递增整数的向量:其中< code>Span从< code >跳转！=4到< code>==4它增加1，否则保持不变。 < li >所有内容都包装在< code>ifelse中，因此您只能看到< code>x为真的数字，即< code>Span == 4。