内容简介:http://stackoverflow.com/questions/27520310/union-of-intersecting-vectors-in-a-list-in-r
我有一个向量列表如下.
data <- list(v1=c("a", "b", "c"), v2=c("g", "h", "k"), v3=c("c", "d"), v4=c("n", "a"), v5=c("h", "i"))
我正在努力实现以下目标
1)检查任何矢量是否相互相交.
2)如果发现相交向量,得到它们的联合.
所以想要的输出是
out <- list(v1=c("a", "b", "c", "d", "n"), v2=c("g", "h", "k", "i"))
我可以得到一组相交集的联合如下.
Reduce(union, list(data[[1]], data[[3]], data[[4]])) Reduce(union, list(data[[2]], data[[5]])
如何首先识别交叉向量?有没有办法将列表分成交叉向量组的列表?
更新
这是尝试使用data.table.获得所需的结果.但是对于这个 example 数据集中的大型列表来说仍然很慢.
datasets. data <- sapply(data, function(x) paste(x, collapse=", ")) data <- as.data.frame(data, stringsAsFactors = F) repeat { M <- nrow(data) data <- data.table( data , key = "data" ) data <- data[ , list(dataelement = unique(unlist(strsplit(data , ", " )))), by = list(data)] data <- data.table(data , key = "dataelement" ) data <- data[, list(data = paste0(sort(unique(unlist(strsplit(data, split=", ")))), collapse=", ")), by = "dataelement"] data$dataelement <- NULL data <- unique(data) N <- nrow(data) if (M == N) break } data <- strsplit(as.character(data$data) , "," )
这样就像一个图形问题,所以我喜欢用这个图书库,使用你的样本数据,你可以做
library(igraph) #build edgelist el <- do.call("rbind",lapply(data, embed, 2)) #make a graph gg <- graph.edgelist(el, directed=F) #partition the graph into disjoint sets split(V(gg)$name, clusters(gg)$membership) # $`1` # [1] "b" "a" "c" "d" "n" # # $`2` # [1] "h" "g" "k" "i"
我们可以查看结果
V(gg)$color=c("green","purple")[clusters(gg)$membership] plot(gg)
http://stackoverflow.com/questions/27520310/union-of-intersecting-vectors-in-a-list-in-r
以上就是本文的全部内容,希望本文的内容对大家的学习或者工作能带来一定的帮助,也希望大家多多支持 码农网
猜你喜欢:本站部分资源来源于网络,本站转载出于传递更多信息之目的,版权归原作者或者来源机构所有,如转载稿涉及版权问题,请联系我们。