聚类2元整数向量

Question

假设我是一组2元素的整数向量，例如，这个矩阵中的每一行：

my.mat = matrix(data = c(1,4,2,5,7,2,3,6,3,8), ncol = 2, byrow=T)
> my.mat
     [,1] [,2]
[1,]    1    4
[2,]    2    5
[3,]    7    2
[4,]    3    6
[5,]    3    8

或者另一种选择是：

my.list = list(c(1,4), c(2,5), c(7,2), c(3,6), c(3,8))

我正在寻找一个函数，它将这些向量相互交叉，并返回一个list，其中每个元素都是所有相交向量的值。不与任何其他向量相交的向量的值也应作为列表元素返回。

简单地说，这应该返回以下列表：

list(c(1,4), c(2,5,7,2), c(3,6,3,8))

如您所见，向量c(1,4)不与任何其他向量相交，因此它作为单个元素列表返回。另外两个列表元素是交叉的。

Answer 1

对于这种类型的问题，我喜欢使用igraph包，因为很容易将行看作图中的边，并且您正在查找所有连接的块。

my.mat = matrix(data = c(1,4,2,5,7,2,3,6,3,8), ncol = 2, byrow=T)

library(igraph)
class(my.mat)<-"character"
gg <- graph.edgelist(my.mat, directed=T)
vset <- clusters(gg)$membership

xx<-sapply(unique(vset), function(s) {
    es <- E(gg)[from(V(gg)[vset==s])]
    as.numeric(V(gg)[t(get.edges(gg,es))]$name)
})

这会把你的矩阵变成一个边缘列表(只有当你有不连续的数字时，才需要转换到字符，所以它更一般地起作用)。这使得这张图

​

​

然后，clusters()函数将它们划分为子组。然后，我遍历每个子群，提取每个子群边缘的顶点名称。

最后，xx的值将是

[[1]]
[1] 1 4

[[2]]
[1] 2 5 7 2

[[3]]
[1] 3 6 3 8

所以这是一个连接值的列表。