7

我有一个带有布尔变量事件的数据框 d,指示某个事件是否在给定日期发生。我想创建一个新变量,指示距离最近的事件有多少观察(天)。

d=structure(list(date = structure(c(-365, -364, -363, -362, -361, 
-360, -359, -358, -357, -356, -355, -354, -353, -352, -351, -350, 
-349, -348, -347, -346), class = "Date"), event = c(TRUE, 
FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, 
FALSE, FALSE, FALSE, FALSE, TRUE, FALSE, FALSE, FALSE, FALSE, 
FALSE)), .Names = c("date", "event"), row.names = c(NA, 20L
), class = "data.frame")

有没有可以做到这一点的功能?

4

2 回答 2

5

就像是

apply(abs( sapply( which(d$event), "-", 1:nrow(d) )),1,min)

将@DWin 的答案概括为超过 2 个TRUE值。

于 2012-11-06T16:00:36.390 回答
4
> pmin( abs( sapply( which(d$event), "-", 1:nrow(d) )[,1] ) , 
        abs( sapply( which(d$event), "-", 1:nrow(d) )[,2] ) )
 [1] 0 1 2 3 4 5 6 7 6 5 4 3 2 1 0 1 2 3 4 5
于 2012-11-05T22:22:16.547 回答