r - 为二分图检测 r 中的双团

Question

我正在尝试在 R 中重新创建 Biclique 社区方法 ( Lehmann, Schwartz, & Hansen, 2008 )，它依赖于 Ka,b biclique 的定义。下面的示例显示了两个相邻的 K2,2 bicliques - 第一个 clique 是 {A,B,1,2}，第二个 clique 是 {B,C,2,3}。我希望能够使用 R 识别这些集团，以便我可以将其应用于更广泛的数据集。

到目前为止，我已将我的尝试包含在 R 中，但我遇到了以下两个问题：

1. 如果我使用标准 walktrap.community 它识别社区但不允许集合 {B,2} 属于两个派系
2. 如果我使用更新的clique.community功能，这似乎无法识别派系或我不正确理解（或两者兼而有之）

示例代码：

library(igraph)

clique.community <- function(graph, k) {
  clq <- cliques(graph, min=k, max=k)
  edges <- c()
  for (i in seq_along(clq)) {
    for (j in seq_along(clq)) {
      if ( length(unique(c(clq[[i]], clq[[j]]))) == k+1 ) {
        edges <- c(edges, c(i,j))
      }
    }
  }
  clq.graph <- simplify(graph(edges))
  V(clq.graph)$name <- seq_len(vcount(clq.graph))
  comps <- decompose.graph(clq.graph)

  lapply(comps, function(x) {
    unique(unlist(clq[ V(x)$name ]))
  })
}

users <- c('A', 'A', 'B', 'B', 'B', 'C', 'C')
resources <- c(1, 2, 1, 2, 3, 2, 3)
cluster <- data.frame(users, resources)
matrix <- as.data.frame.matrix(table(cluster))

igraph <- graph.incidence(matrix)

clique.community(igraph, 2)
walktrap.community(igraph)

Answer 1

请注意，即使对于小（密集）图和 k,l 值，上述解决方案也会很快变得低效，因为它comb <- combn(vMode1, k)变得非常大。

在https://github.com/YupingLu/biclique正在开发的“biclique”包中可以找到更有效的解决方案。

Answer 2

我设法在Sisob 工作台中找到了一个脚本

computeBicliques <- function(graph, k, l) {

  vMode1 <- c()
  if (!is.null(V(graph)$type)) {

    vMode1 <- which(!V(graph)$type)
    vMode1 <- intersect(vMode1, which(degree(graph) >= l))
  }

  nb <- get.adjlist(graph)

  bicliques <- list()

  if (length(vMode1) >= k) {

    comb <- combn(vMode1, k)
    i <- 1
    sapply(1:ncol(comb), function(c) {

      commonNeighbours <- c()
      isFirst <- TRUE

      sapply(comb[,c], function(n) {

        if (isFirst) {

          isFirst <<- FALSE
          commonNeighbours <<- nb[[n]]
        } else {

          commonNeighbours <<- intersect(commonNeighbours, nb[[n]])
        }
      })

      if (length(commonNeighbours) >= l) {

        bicliques[[i]] <<- list(m1=comb[,c], m2=commonNeighbours)
      }

      i <<- i + 1
    })
  }
  bicliques
}