假设我有以下数据。框架,其中pos是位置坐标。我已经包括了一个变量thresh,其中val大于给定的阈值t。
set.seed(123)
n <- 20
t <- 0
DF <- data.frame(pos = seq(from = 0, by = 0.3, length.out = n),
val = sample(-2:5, size = n, replace = TRUE))
DF$thresh <- DF$val > t
DF
## pos val thresh
## 1 0.0 0 FALSE
## 2 0.3 4 TRUE
## 3 0.6 1 TRUE
## 4 0.9 5 TRUE
## 5 1.2 5 TRUE
## 6 1.5 -2 FALSE
## 7 1.8 2 TRUE
## 8 2.1 5 TRUE
## 9 2.4 2 TRUE
## 10 2.7 1 TRUE
## 11 3.0 5 TRUE
## 12 3.3 1 TRUE
## 13 3.6 3 TRUE
## 14 3.9 2 TRUE
## 15 4.2 -2 FALSE
## 16 4.5 5 TRUE
## 17 4.8 -1 FALSE
## 18 5.1 -2 FALSE
## 19 5.4 0 FALSE
## 20 5.7 5 TRUE
我怎么能得到区域坐标,其中val是积极的,即在上面的例子:
0.3 - 1.2,
1.8 - 3.9,
4.5 - 4.5,
5.7 - 5.7
我曾想过拆分数据。按阈值帧,然后从每个数据的第一行和最后一行访问pos。框架列表元素,但这只会将所有TRUE和FALSE子集组合在一起。有没有办法根据真值将thresh变量转换为字符,并丢弃假值?
split(DF, DF$thresh) # not what I want
## $`FALSE`
## pos val thresh
## 1 0.0 0 FALSE
## 6 1.5 -2 FALSE
## 15 4.2 -2 FALSE
## 17 4.8 -1 FALSE
## 18 5.1 -2 FALSE
## 19 5.4 0 FALSE
##
## $`TRUE`
## pos val thresh
## 2 0.3 4 TRUE
## 3 0.6 1 TRUE
## 4 0.9 5 TRUE
## 5 1.2 5 TRUE
## 7 1.8 2 TRUE
## 8 2.1 5 TRUE
## 9 2.4 2 TRUE
## 10 2.7 1 TRUE
## 11 3.0 5 TRUE
## 12 3.3 1 TRUE
## 13 3.6 3 TRUE
## 14 3.9 2 TRUE
## 16 4.5 5 TRUE
## 20 5.7 5 TRUE
我尝试过的另一件笨重的事情是cumsum
,但它同样包含错误的行:
split(DF, cumsum(DF$thresh == 0)) # not what I want but close to it...
## $`1`
## pos val thresh
## 1 0.0 0 FALSE
## 2 0.3 4 TRUE
## 3 0.6 1 TRUE
## 4 0.9 5 TRUE
## 5 1.2 5 TRUE
##
## $`2`
## pos val thresh
## 6 1.5 -2 FALSE
## 7 1.8 2 TRUE
## 8 2.1 5 TRUE
## 9 2.4 2 TRUE
## 10 2.7 1 TRUE
## 11 3.0 5 TRUE
## 12 3.3 1 TRUE
## 13 3.6 3 TRUE
## 14 3.9 2 TRUE
##
## $`3`
## pos val thresh
## 15 4.2 -2 FALSE
## 16 4.5 5 TRUE
##
## $`4`
## pos val thresh
## 17 4.8 -1 FALSE
##
## $`5`
## pos val thresh
## 18 5.1 -2 FALSE
##
## $`6`
## pos val thresh
## 19 5.4 0 FALSE
## 20 5.7 5 TRUE
这里有一个带有数据的选项。表
。我们使用rleid
创建一个分组变量,根据“thresh”和split
对“pos”进行子集划分。
DT <- setDT(DF)[,pos[thresh] ,.(gr=rleid(thresh))]
split(DT$V1, DT$gr)
#$`2`
#[1] 0.3 0.6 0.9 1.2
#$`4`
#[1] 1.8 2.1 2.4 2.7 3.0 3.3 3.6 3.9
#$`6`
#[1] 4.5
#$`8`
#[1] 5.7
或者我们可以使用rle
frombase R
创建分组变量,然后基于该变量进行split
gr <- inverse.rle(within.list(rle(DF$thresh), values <- seq_along(values)))
with(DF, split(pos[thresh], gr[thresh]))
或者正如@thelatemail所提到的,cumsum
也可以在使用“thresh”进行分组后用于分组。
with(DF, split(pos[thresh],cumsum(!thresh)[thresh]))