我有非常广泛的网络数据
df <- data.frame(name= c("Sue", "Rob", "Jacob", "Terry"),
f1 = c("Rob", "Sue", "Rob", "Jacob"),
f2 = c(NA, "Terry", "Terry", NA),
f3 = c(NA, "Jacob", NA, NA),
bf1 = c(1,1,1,0),
bf2 = c(NA, 1, 0, NA),
bf3 = c(NA, 0, NA, NA) )
变量f1
到f3
指的是朋友Sue,Rob,Jacob和Terry。变量bf1
到bf3
是指他们是否将其视为最好的朋友。
我现在想要根据此数据框创建一个R igraph
的图表项,并包含所有信息。到目前为止,我有:
adj_list <- subset(df, select = c(name, f1:f3))
m <- as.matrix(adj_list)
el1 <- cbind(m[, 1], c(m[, -1]))
el1 <- el1[complete.cases(el1),]
el1[2,] <- as.numeric(factor(c(el1[1,])))
g1 <- graph_from_edgelist(el1)
但我不知道如何将次要信息导入图形对象。有任何想法吗?
答案 0 :(得分:1)
一种方法是将数据重新整形为长格式。
# Get indices of columns to group together
friendsIndices <- with(dat, grep("^f", names(dat))) # strings beginning with 'f`
bffIndices <- with(dat, grep("^bf", names(dat))) # strings beginning with 'bf`
整形
ldat <- reshape(dat,
varying=list(friendsIndices, bffIndices),
direction="long", idvar="name")
# Tidy up long data -
# remove time column
# remove rows with `NA`
ldat <- ldat[!is.na(ldat$f1), -match("time", names(ldat))]
ldat
# name f1 bf1
# Sue.1 Sue Rob 1
# Rob.1 Rob Sue 1
# Jacob.1 Jacob Rob 1
# Terry.1 Terry Jacob 0
# Rob.2 Rob Terry 1
# Jacob.2 Jacob Terry 0
# Rob.3 Rob Jacob 0
现在可以将数据读入igraph
library(igraph)
g <- graph_from_data_frame(ldat)
bf
属性存储在边缘属性E(g)$bf1
中。
然后,您可以在基于此属性进行绘图时添加边缘颜色或权重。
plot(g,
edge.width=2*(E(g)$bf1+1),
edge.color=c("red", "blue")[E(g)$bf1+1],
edge.curved=TRUE)