6

我有发明人在专利方面的合作数据。每个发明人都是一个节点，每条边代表两个发明人合作的专利。一些专利的发明人超过 2 个，因此一些专利用多个边表示。

我想对至少一位发明人位于博伊西但并非所有发明人都位于博伊西的专利进行子图化。其他专利和发明人需要排除在选择之外。

例如：

gg <- graph.atlas(711)
V(gg)$name <- 1:7
V(gg)$city <- c("BOISE","NEW YORK","NEW YORK","BOISE","BOISE","LA","LA")
V(gg)$color <- ifelse(V(gg)$city=="BOISE", "orange","yellow")
gg<-delete.edges(gg, E(gg, P=c(1,2,2,3,2,7,7,6,7,3,3,4,3,5,4,5,5,6,6,1))) 
gg <- add.edges(gg,c(1,4,4,5,5,1),attr=list(patent=1))
gg <- add.edges(gg,c(7,5,5,4,4,7),attr=list(patent=2))
gg <- add.edges(gg,c(7,3,3,5,5,7),attr=list(patent=3))
gg <- add.edges(gg,c(2,7,7,6,6,2),attr=list(patent=4))
gg <- add.edges(gg,c(6,4),attr=list(patent=5))
plot(gg, edge.label=E(gg)$patent)

产生：

网络示例

从这个网络中，我只想对所有出现在专利 2、3、5 边缘的节点进行子图绘制。

在此示例中，节点 1 不应出现在子图中。此外，还应排除与专利 #1 有关的从节点 5 到节点 4 的边。

我一直在为这个问题苦苦挣扎一段时间。这可能吗？

4

1 回答 1

7

这个怎么样

#final all patent names
patents <- unique(edge.attributes(gg)$patent)

#check if criteria are met for patent
okpatents <- sapply(patents, function(p){
    cities <- V(gg)[inc(E(gg)[patent==p])]$city
    nc <- sum(cities=="BOISE")
    return(nc>0 & nc < length(cities))
})

#extract subgraph
gs <- subgraph.edges(gg, E(gg)[patent %in% patents[okpatents]])

#verify
plot(gs, edge.label=E(gs)$patent)

在此处输入图像描述

PS。非常好的可重现示例；）

于 2014-08-26T14:18:59.327 回答