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