2013-01-13 81 views
0

我正在尝试通过我的发布网络将我的协作网络化。我用igraph,它很棒。然而,由于我的顶点(代表我在网络中的那个顶点)到所有共同作者的边缘都有一个漂亮的图表。我想从我的顶点删除一些只通过另一个作者相关的作者。基本上是作者,我只是共同创作。无论如何,我已经确定了这些顶点,并且我知道我的顶点。现在我无法找到一种方法去除仅将这一组边连接到我的边。更一般地说,如何从两组矢量中去除边,例如V(g)[a]和V(g)[b]?r igraph:从两组矢量中删除边缘

感谢,

这里有一个例子:

au1 <- c('deb', 'art', 'deb', 'seb', 'deb', 'deb', 'mar', 'mar', 'joy', 'deb') 
au2 <- c('art', 'deb', 'soy', 'deb', 'joy', 'ani', 'deb', 'deb', 'nem', 'mar') 
au3 <- c('mar', 'lio', 'mil', 'mar', 'ani', 'lul', 'nem', 'art', 'deb', 'tat') 


tata <- data.frame(au1, au2, au3) 
xaulist2 <- levels(factor(unlist(tata[,]))) 
xaulist <- levels(as.factor(xaulist2)) 
xaulist_att <- c(rep('prime', 2), 'main', 'second', 'second', rep('prime', 3), 'second', rep('prime', 3)) 
au_att <- data.frame(au_name=xaulist, level=xaulist_att) 

# matrix list preparation 
tutu <- matrix(NA, nrow=length(xaulist), ncol=dim(tata)[1]) # row are authors and col are papers 
for (i in 1:length(xaulist)) 
{ 
    for (j in 1:dim(tata)[1]) 
    { 
    ifelse('TRUE' %in% as.character(tata[j,]==xaulist[i]), tutu[i,j] <- 1, tutu[i,j] <- 0) 
    } 
} 
tutu[is.na(tutu)] <- 0 

tutu[tutu>=1] <- 1 # change it to a Boolean matrix 
termMatrix <- tutu %*% t(tutu) 

# build a graph from the above matrix 
g <- graph.adjacency(termMatrix, weighted=T, mode = 'undirected') 
g <- simplify(g) # remove loops 
V(g)$label <- xaulist # set labels of vertices 
V(g)$degree <- degree(g) # set degrees of vertices 
V(g)[xaulist_att=='second']$color <- 'red' 
V(g)[xaulist_att=='main']$color <- 'blue' 
set.seed(112) # set seed to make the layout reproducible 
plot(g) 

所以,问题是你如何从有属性“第二”到具有属性的“主”的一个作者删除边缘,即是从红色到蓝色的?

再次感谢,

回答

5

下面是做到这一点的一种方式,它认为这是非常具有可读性:

g2 <- delete.edges(g, E(g)[ V(g)[xaulist_att == 'second'] %--% 
          V(g)[xaulist_att == 'main' ] ]) 

## plot the results 
coords <- layout.auto(g) 
layout(rbind(1:2)) 
plot(g, layout=coords, main="g") 
plot(g2, layout=coords, main="g2") 

详情请参见?iterators

+0

yeahhhhhhh!有用。谢谢, – jejuba

1

没有重复的例子,我不知道我有明白你想要什么。但是,这应该可以帮助您:

##deletes from g edge 1->2. 
g <- delete.edges(g, E(g, P=c(1,2)))