Tidygraph и igraph - построить график из несоответствия данных
Я могу построить объект графа в igraph из двух датафреймов без проблем. Когда я пытаюсь сделать то же самое в tidygraph, я получаю ошибки. Позвольте мне продемонстрировать. Сначала я загружаю свои исходные данные (данные с доски объявлений):
library(dplyr)
library(tidyr)
library(tidygraph)
library(lubridate)
library(iterpc)
library(igraph)
df <- data.frame(author_id = c(2,4,8,16,4,8,2,256,512,8),
topic_id = c(101,101,101,101,301,301,501,501,501,501),
time = as.POSIXct(c("2011-08-16 20:20:11", "2011-08-16 21:10:00", "2011-08-17 06:30:10",
"2011-08-17 10:08:32", "2011-08-20 22:23:01","2011-08-20 23:03:03",
"2011-08-25 17:05:01", "2011-08-25 19:15:10", "2011-08-25 20:07:11",
"2011-08-25 23:59:59")),
vendor = as.logical(c("FALSE", "FALSE", "TRUE", "FALSE", "FALSE",
"TRUE", "FALSE", "FALSE", "FALSE", "TRUE")))
Затем я создаю уникальный список узлов (люди, которые публикуют вещи на доске объявлений):
node <- df %>% distinct(author_id, vendor) %>% rename(id = author_id) %>% mutate(vendor = as.numeric(vendor))
Затем мой крайний список (люди, связанные через дискуссионную ветку (тему)):
edge <- df %>%
group_by(topic_id) %>%
do(data.frame(getall(iterpc(table(.$author_id), 2, replace =TRUE)))) %>%
filter(X1 != X2) %>% rename(from = X1, to = X2) %>% select(to, from, topic_id)
С помощью igraph я могу создать этот граф объект:
test_net <- graph_from_data_frame(d = edge, directed = F, vertices = node)
plot(test_net)
Это выглядит хорошо. Теперь я пытаюсь сделать то же самое с Tidygraph:
tidy_net <- tbl_graph(nodes = node, edges = edge, directed = F)
Error in add_vertices(gr, nrow(nodes) - gorder(gr)) : At type_indexededgelist.c:369 : cannot add negative number of vertices, Invalid value
Хлоп! Однако, когда я импортирую объект igraph в тидограф:
tidy_net <- as_tbl_graph(test_net)
plot(tidy_net)
Все работает! Что здесь происходит? Пожалуйста помоги.
1 ответ
Я думаю, потому что ваши узлы id
и края to
а также from
являются числовыми, это предполагает, что должны быть узлы для каждого целого числа между min(node$id)
(2) и max(node$id)
(512). Вы можете обойти это, принуждая их к персонажам. Кроме того, ваш iterpc
Команда не работает для меня должным образом, поэтому я преобразовал его в tidyr
версия расширения ваших данных.
node <-
df %>%
distinct(author_id, vendor) %>%
rename(id = author_id) %>%
mutate(vendor = as.numeric(vendor)) %>%
mutate(id = as.character(id))
edge <-
df %>%
group_by(topic_id) %>%
expand(topic_id, from = author_id, to = author_id) %>%
filter(from < to) %>%
select(to, from, topic_id) %>%
mutate_at(vars(to, from), as.character)
tidy_net <- tbl_graph(nodes = node, edges = edge, directed = F)
plot(tidy_net)