Question

我正在对R中的公司网络进行分析，并试图将我的igraph结果导出到数据框中。

这是一个可重复的例子：

library(igraph)
sample <- data.frame(ID = 1:8, org_ID = c(5,4,1,2,2,2,5,7), mon = c("199801", "199802","199802","199802","199904","199912","200001", "200012"))

create.graphs <- function(df){
g <- graph.data.frame(d = df, directed = TRUE)
g <- simplify(g, remove.multiple = FALSE, remove.loops = TRUE)
E(g)$weight <- count_multiple(g)

#calculate global values
g$centrality <- centralization.degree(g)
#calculate local values
g$indegree <- degree(g, mode = "in",
                   loops = FALSE, normalized = FALSE)

return(g)
}

df.list <- split(sample, sample$mon)
g <- lapply(df.list, create.graphs)

如您所见，我有多个月的图表。我想将其导出为纵向数据，其中每行代表一个月（每个ID），每列代表相应的网络度量。

到目前为止，我已经设法创建了一个数据框，但不是如何通过图表列表运行它并将其放入合适的格式。另一个问题可能是图表具有不同数量的节点（一些节点大约有25个，其他节点大于40个节点），但理论上这应该被我的回归模型识别为缺失。

output <- data.frame(Centrality = g$`199801`$centrality,
        Indegree = g$`199801`$indegree)
output
summary(output)

我尝试为此编写类似于上述功能的功能，但遗憾的是无济于事。

提前感谢您阅读本文，非常感谢任何帮助

Answer 1

我想分享一下我是如何解决的（感谢Dave2e的建议）。

请注意，ci $ monat在原始数据中定义了我的时间段，因此每个时间点都有一行。

sumarTable <- data.frame(time = unique(ci$monat))

sumarTable$indegree <- lapply(g, function(x){x$indegree})
sumarTable$outdegree <- lapply(g, function(x){x$outdegree})
sumarTable$constraint <- lapply(g, function(x){x$constraint})

等

编辑：为了导出这些值，我不得不“压扁”列表：

sumarTable$indegree <- vapply(sumarTable$indegree, paste, collapse = ", ", character(1L))
sumarTable$outdegree <- vapply(sumarTable$outdegree, paste, collapse = ", ", character(1L))
sumarTable$constraint <- vapply(sumarTable$constraint, paste, collapse = ", ", character(1L))

从R

1 个答案: