1

我在 IntegerList 中有位置索引列表，我打算在给定阈值的情况下过滤它们，并且效果很好。但是，我想为每个 IntegerList 提取一个特定的过滤集以供进一步使用。我知道 myList 是嵌套列表，并且数据非常基于真实数据集进行模拟。有什么方法可以轻松优雅地检索想要的 IntegerList 吗？我怎样才能使它发生这种提取？

要运行迷你示例，需要以下库：

library(IRanges)
library(S4Vectors)

小例子：

myList <- list(f1=IntegerList(1,2,3,4,1,1,1,integer(0),1,2,4),
               f2=IntegerList(1,5,integer(0),integer(0),2,3,4,6,1,5,6),
               f3=IntegerList(1,4,6,7,2,3,3,7,2,5,7))

len <- Reduce('+', lapply(myList, lengths))
keepMe <- len >= length(myList)

我打算按如下方式过滤它们：

res.filt <- lapply(myList, function(elm) {
  ans <- list(keep=elm[keepMe], droped=elm[!keepMe])
  ans
})

我的粗略输出：

Keep.list <- list(f1.kp=res.filt$f1$keep, f2.kp=res.filt$f2$keep, f3.kp=res.filt$f3$keep)
Drop.list <- list(f1.dp=res.filt$f1$droped, f2.dp=res.filt$f2$droped, f3.dp=res.filt$f3$droped)

根据我粗略的输出，我怎样才能得到更优雅的输出？任何有效的方法来实现我的输出？谁能指出我该怎么做？或任何建议如何获得我的预期输出？提前致谢

4

1 回答 1

2

您用于过滤向量列表的思维过程/流程是合乎逻辑且非常理想的，但如果您使用以下内容，您可以将其收紧一点purrr：

library(purrr)

map(myList, lengths) %>% 
  reduce(`+`) %>% 
  map_lgl(`>=`, length(myList)) -> keep_me

keep_list <- map(myList, ~.[keep_me])
drop_list <- map(myList, ~.[!keep_me])
于 2016-10-01T15:37:49.790 回答