我需要找到连续几天我的数据达到阈值的位置。我正在寻找高于阈值的4个连续观察值。我想返回符合这些标准的系列的第一个观察点的位置。连续值高于阈值的位置和值
下面是一个例子的数据集:
eg = structure(list(t.date = structure(c(1L, 2L, 11L, 12L, 13L, 14L,
15L, 16L, 17L, 18L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L), .Label = c("4/30/11",
"5/1/11", "5/10/11", "5/11/11", "5/12/11", "5/13/11", "5/14/11",
"5/15/11", "5/16/11", "5/17/11", "5/2/11", "5/3/11", "5/4/11",
"5/5/11", "5/6/11", "5/7/11", "5/8/11", "5/9/11"), class = "factor"),
t.avg = c(4L, 4L, 5L, 6L, 10L, 18L, 18L, 18L, 18L, 12L, 10L,
10L, 8L, 8L, 9L, 10L, 6L, 5L)), .Names = c("date", "avg"
), row.names = c(NA, -18L), class = "data.frame")
我想去的地方平均满足标准(平均> 17 4天)的时间 一种方法:
eg$date %in% eg$date[which(eg$avg > 17)]
# [1] FALSE FALSE FALSE FALSE FALSE TRUE TRUE TRUE TRUE FALSE FALSE FALSE
# [13] FALSE FALSE FALSE FALSE FALSE FALSE
在这种情况下,我可以采取第一个案例TRUE
作为答案,但如果第二,第三或第四个不是TRUE
我需要在第一次约会中条件TRUE
:
eg$date[which(eg$avg > 17)]
# [1] 5/5/11 5/6/11 5/7/11 5/8/11
而且先观察该系列中的位置:
which(eg$avg > 17)
# [1] 6 7 8 9
我发现相关问题,但我一直无法弯曲我需要的方法。
非常感谢。使用基础R
难道不应该如果存在这样的功能,就是'rollmin'。或者也许rollapply(...,分,4) – 2014-11-22 00:13:23
好抓!我对这个字段名称是“平均”而感到困惑,并且一直跟着它。我编辑以反映您的评论。 – Jordan 2014-11-22 00:29:31
'rollmax'存在如此:'which(-rollmax(-eg $ avg,4)> 17)[1]' – 2014-11-22 13:53:00