我正在尝试进行疾病模拟,我希望受感染的节点(is_infected>0)为红色,未感染的(is_infected=0)节点颜色为青色。
results <- results %>% mutate(
S = infected == 0,
E = 0 < infected & infected < 6,
I = infected >= 6 & infected <= 11,
R = infected > 11,
is_infected = infected > 0
)
net.layout.by.time <-
split(results, f = results$t%in% c(5, 7, 16, 40, 80)) %>%
lapply(FUN = right_join,
y = net.layout,
by = "id") %>%
bind_rows
net.layout.by.time <- split(results, f = results$t) %>%
lapply(FUN = right_join, y = net.layout, by = "id") %>%
bind_rows
net.layout.by.time %>%
filter(t %in% c(5, 6, 7, 20, 40)) %>%
ggplot(aes(xend = xend, yend = yend, x = x, y = y)) +
geom_edges(color = "lightgray") +
geom_nodes(aes(color = is_infected)) +
facet_wrap(~ t) +
theme_blank()