## 安装加载数据
if(!require(devtools)) install.packages("devtools")
devtools::install_github("kassambara/navdata")
install.packages(
  c("tidyverse", "igraph", "tidygraph", "ggraph")
)
library(pacman)
p_load("navdata","tidyverse", "igraph", "tidygraph", "ggraph")

library("navdata")
data("phone.call")
head(phone.call, 3)
## # A tibble: 3 x 3
##    source destination n.call
##              
## 1  France     Germany      9
## 2 Belgium      France      4
## 3  France       Spain      3

#  来源地去重
sources <- phone.call %>%
  distinct(source) %>%
  rename(label = source)
#目的地去重
destinations <- phone.call %>%
  distinct(destination) %>%
  rename(label = destination)
## 合并数据并添加一列索引
nodes <- full_join(sources, destinations, by = "label") 
nodes <- nodes %>%
  mutate(id = 1:nrow(nodes)) %>%
  select(id, everything())
head(nodes, 3)
## # A tibble: 3 x 2
##      id   label
##      
## 1     1  France
## 2     2 Belgium
## 3     3 Germany

# Rename the n.call column to weight
phone.call <- phone.call %>%
  rename(weight = n.call)
# (a) Join nodes id for source column
edges <- phone.call %>% 
  left_join(nodes, by = c("source" = "label")) %>% 
  rename(from = id)
# (b) Join nodes id for destination column
edges <- edges %>% 
  left_join(nodes, by = c("destination" = "label")) %>% 
  rename(to = id)
# (c) Select/keep only the columns from and to
edges <- select(edges, from, to, weight)
head(edges, 3)
## # A tibble: 3 x 3
##    from    to weight
##      
## 1     1     3      9
## 2     2     1      4
## 3     1     8      3

library(igraph)
net.igraph <- graph_from_data_frame(
  d = edges, vertices = nodes, 
  directed = TRUE
  )
  ## d 边列表，vertices 点列表，directer：有向或者无向
  ## network
  set.seed(123)
plot(net.igraph, edge.arrow.size = 0.2,
     layout = layout_with_graphopt)

library(tidygraph)
net.tidy <- tbl_graph(
  nodes = nodes, edges = edges, directed = TRUE
  )

library(ggraph)
ggraph(net.tidy, layout = "graphopt") + 
  geom_node_point(col = 'gold',size = 2) + # 点信息
  geom_edge_link(aes(width = weight), alpha = 0.8) +  # 边信息
  scale_edge_width(range = c(0.2, 2)) + # 控制粗细
  geom_node_text(aes(label = label), repel = TRUE) + # 增加节点的标签，reple避免节点重叠
  labs(edge_width = "phone.call") + # 图例标签
  theme_graph()

# Arc diagram
ggraph(net.tidy, layout = "linear") + 
  geom_edge_arc(aes(width = weight), alpha = 0.8) + 
  scale_edge_width(range = c(0.2, 2)) +
  geom_node_text(aes(label = label), repel = TRUE) +
  labs(edge_width = "Number of calls") +
  theme_graph()+
  theme(legend.position = "top")

# Coord diagram, circular
ggraph(net.tidy, layout = "linear", circular = TRUE) + 
  geom_edge_arc(aes(width = weight), alpha = 0.8) + 
  scale_edge_width(range = c(0.2, 2)) +
  geom_node_text(aes(label = label), repel = TRUE) +
  labs(edge_width = "Number of calls") +
  theme_graph()+
  theme(legend.position = "top")

library(tidyverse)
library("navdata")
data("phone.call2")
nodes <- phone.call2$nodes
edges <- phone.call2$edges
## 刚才数据集是phone.call，phone.call2可以直接提取点和边列表
##注意：该包中节点和边缘列表中的ID应该是从0开始的数字值。我们直接减1 即可
nodes_d3 <- mutate(nodes, id = id - 1)
edges_d3 <- mutate(edges, from = from - 1, to = to - 1)

library(networkD3)
forceNetwork(
  Links = edges_d3, Nodes = nodes_d3,bounded = T,
  Source = "from", Target = "to",      # so the network is directed.
  NodeID = "label", Group = "id", Value = "weight", 
  width = 700, height=500,legend = T,
  colourScale = JS("d3.scaleOrdinal(d3.schemeCategory10);"), ## 颜色模板
  opacity = 1, fontSize = 18, zoom = TRUE, opacityNoHover = 1  ## 显示标签信息
)
## 如果想自定义颜色新修改，直接括号中改为16进制颜色即可
 colourScale = JS("d3.scaleOrdinal([`#383867`, `#584c77`, `#e1eff7`,`#dd2c29`]);"),

sankeyNetwork(
  Links = edges_d3, Nodes = nodes_d3, 
  Source = "from", Target = "to", 
  NodeID = "label", Value = "weight", 
  fontSize = 16, unit = "Letter(s)")

library(visNetwork)
visNetwork(nodes, edges) %>%
  visLayout(randomSeed = 12)

visNetwork(nodes, edges) %>% 
  visIgraphLayout(layout = "layout_with_fr") %>% 
  visEdges(arrows = "middle") %>%
  visLayout(randomSeed = 1234)