r IGRAPHE-基于节点的公共连接和与第三个节点的相同类型的连接,在节点之间创建连接

r IGRAPHE-基于节点的公共连接和与第三个节点的相同类型的连接,在节点之间创建连接,r,igraph,R,Igraph,节点A和B当前未连接,如果它们满足两个条件,我想连接它们:1)它们都连接到同一个第三个节点;2) 与第三个节点的连接类型相同。假设A和B是同一个父亲的儿子,那么我想确定他们是兄弟姐妹。我如何指导r igraph创建此新领带?以下面的例子为例 edgelist <- read.table(text = " A C B C C D") graph <- graph.data.frame(edgelist, directed = F) E(graph)[1]$weight <- 2

节点A和B当前未连接,如果它们满足两个条件,我想连接它们:1)它们都连接到同一个第三个节点;2) 与第三个节点的连接类型相同。假设A和B是同一个父亲的儿子,那么我想确定他们是兄弟姐妹。我如何指导r igraph创建此新领带?以下面的例子为例

edgelist <- read.table(text = "
A C
B C
C D")
graph <- graph.data.frame(edgelist, directed = F)
E(graph)[1]$weight <- 2
E(graph)[2]$weight <- 2
E(graph)[3]$weight <- 1

IGRAPH 0dd6cf1 DNW- 4 3 -- 
+ attr: name (v/c), weight (e/n)
+ edges from 0dd6cf1 (vertex names):
[1] A->C B->C C->D

edgelist可能有一种更简单的方法,但我认为这可能会满足您的需求。我稍微扩展了您的示例,添加了多个案例

library(igraph)

options(stringsAsFactors = FALSE)

edgelist <- read.table(text = "
                       A C
                       B C
                       C D
                       D A
                       E C
                       D F
                       D G")

g <- graph.data.frame(edgelist, directed = F)

E(g)$weight <- c(2,2,1,1,2,2,2)

#plot graph, scaling the weights just to make them more obvious
plot(g, edge.width = E(g)$weight*3)
库(igraph)
选项(stringsAsFactors=FALSE)

edgelist您是否希望以传递方式建立这些连接?即,假设在您的示例中,有另一个节点E通过权重为2的链接连接到B。A和E一开始没有共享任何邻居,但在A连接到B后,A和E将共享B作为权重为2的邻居。A和E是否应该连接?因此在本例中,A和B之间的新连接将不具有重量2。(父母是2岁,我们可以说兄弟姐妹是1岁)。因此,A和E在A-B连接形成后不会连接。
#convert back to edgelist
el <- as.data.frame(get.edgelist(g))
el$weight <- E(g)$weight

ids <- unique(c(el$V1, el$V2))

#select sets of alters of all nodes that have edge weight equal to 2 
y <- lapply(ids, function(id) {

  x <- el[which(el$V1 == id | el$V2 == id),]
  x <- x[which(x$weight == 2),]
  alt_nodes <- setdiff(unique(c(x$V1, x$V2)), id)

})

#select sets that have 2 or more edges with weight 2
ly <- which(unlist(lapply(y, length))>= 2)

#add connections of weight 1 between all of these sets of nodes
res <- lapply(ly, function (i) {

  new_edge <- y[[i]]
  ne <- t(combn(new_edge,2))
  ne <- cbind(ne, rep(1, nrow(ne)))
  colnames(ne) <- names(el)
  el< <- rbind(el, ne)

})

#convert back to graph
g2  <-  graph.data.frame(el, directed  =  F)
E(g2)$weight <- el$weight

plot(g2, edge.width = as.numeric(E(g2)$weight)*3)