Извлечение списка ребер с условиями из объекта igraph в R
Я работаю с неориентированным объектом igraph, состоящим из узлов разных типов (например, мужчины M желтого цвета и женщины F оранжевого цвета):
g <- graph.atlas(711)
V(g)$name <- 1:7
V(g)$gender <- c("M","F","M","M","M","F","F")
V(g)$color <- ifelse(V(g)$gender=="F", "orange","yellow")
g<-delete.edges(g, E(g, P=c(1,2,2,3,2,7,7,6,7,3,3,4,3,5,4,5,5,6,6,1)))
g<-add.edges(g,c(1,4,4,5,5,1,4,7,7,3,3,5,5,7,2,7,7,6,6,2,6,4))
plot(g)
Я хотел бы извлечь список ребер, состоящий из ребер, соединяющих узлы разных типов (мужчин и женщин):
edgelist <- rbind(c(3,7),
c(4,6),
c(4,7),
c(5,7))
assortativity
использует долю ребер, соединяющих вершины типа M и F, но я не знаю, чтобы извлечь из этих ребер простоту.
get.edgelist
возвращает только весь список ребер без возможности установки условий.
2 ответа
Решение
Вы можете использовать %--%
селектор, чтобы найти ребра, которые соединяют мужские узлы с женскими узлами. Например
E(g)[V(g)[gender=="M"] %--% V(g)[gender=="F"]]
V(g)[gender=="M"]
находит все "мужские" узлы и V(g)[gender=="F"]
находит все женские узлы и %--%
находит все ребра между двумя наборами.
edges = get.edgelist(g)
edges[rowSums(apply(edges, 2, function(x) get.vertex.attribute(g, "gender", x)) == "M") == 1,]
# [,1] [,2]
#[1,] 4 7
#[2,] 3 7
#[3,] 5 7
#[4,] 4 6