当前位置: 首页>后端>正文

R语言可视化(二十九):聚类树图绘制

29. 聚类树图绘制

清除当前环境中的变量

rm(list=ls())

设置工作目录

setwd("C:/Users/Dell/Desktop/R_Plots/29dendrogram/")

使用dendrogram函数绘制聚类树图

# 查看内置示例数据
head(USArrests)
##            Murder Assault UrbanPop Rape
## Alabama      13.2     236       58 21.2
## Alaska       10.0     263       48 44.5
## Arizona       8.1     294       80 31.0
## Arkansas      8.8     190       50 19.5
## California    9.0     276       91 40.6
## Colorado      7.9     204       78 38.7

# 计算距离矩阵,默认method = "euclidean"计算欧氏距离
dists <- dist(USArrests,method = "euclidean") 
head(dists)
## [1]  37.17701  63.00833  46.92814  55.52477  41.93256 128.20694

# 进行层次聚类,method = "average"选择UPGMA聚类算法
hc <- hclust(dists, method = "ave")
hc
## 
## Call:
## hclust(d = dists, method = "ave")
## 
## Cluster method   : average 
## Distance         : euclidean 
## Number of objects: 50

# 将hclust对象转换为dendrogram对象
dend1 <- as.dendrogram(hc)
dend1
## 'dendrogram' with 2 branches and 50 members total, at height 152.314

# 绘制聚类树图,默认type = "rectangle"
plot(dend1, type = "rectangle", 
     ylab="Height",
     main="Cluster Dendrogram")
R语言可视化(二十九):聚类树图绘制,第1张
## "triangle" type and show inner nodes:
plot(dend1, 
     nodePar = list(pch = c(1,NA), cex = 1.2, lab.cex = 0.9),#设置节点的形状,大小和标签字体大小
     type = "triangle", center = TRUE)
R语言可视化(二十九):聚类树图绘制,第2张
plot(dend1, 
     edgePar = list(col = c("red","blue"), lty = 1:2),#设置节点边的颜色和线型
     dLeaf = 2, edge.root = TRUE)
R语言可视化(二十九):聚类树图绘制,第3张
plot(dend1, 
     nodePar = list(pch = 17:16, cex = 1.2:0.8, col = 2:3),
     horiz = TRUE)#水平放置聚类树
R语言可视化(二十九):聚类树图绘制,第4张
nP <- list(col = 3:2, cex = c(2.0, 0.8), pch =  21:22,
           bg =  c("light blue", "pink"),
           lab.cex = 0.8, lab.col = "tomato")
plot(dend1, 
     nodePar= nP, 
     edgePar = list(col = "gray", lwd = 2), 
     horiz = TRUE)
R语言可视化(二十九):聚类树图绘制,第5张
# plot dendrogram with some cuts
dend2 <- cut(dend1, h = 70)
dend2
## $upper
## 'dendrogram' with 2 branches and 4 members total, at height 152.314 
## 
## $lower
## $lower[[1]]
## 'dendrogram' with 2 branches and 2 members total, at height 38.52791 
## 
## $lower[[2]]
## 'dendrogram' with 2 branches and 14 members total, at height 44.28392 
## 
## $lower[[3]]
## 'dendrogram' with 2 branches and 14 members total, at height 44.83793 
## 
## $lower[[4]]
## 'dendrogram' with 2 branches and 20 members total, at height 54.74683

plot(dend2$upper, main = "Upper tree of cut at h=70")
R语言可视化(二十九):聚类树图绘制,第6张
##  dend2$lower is *NOT* a dendrogram, but a list of .. :
plot(dend2$lower[[1]], main = "First branch of lower tree with cut at h=70")
R语言可视化(二十九):聚类树图绘制,第7张
## "inner" and "leaf" edges in different type & color :
plot(dend2$lower[[2]], 
     nodePar = list(col = 1),   # non empty list
     edgePar = list(lty = 1:2, col = 2:1), 
     edge.root = TRUE)
R语言可视化(二十九):聚类树图绘制,第8张
plot(dend2$lower[[3]], 
     nodePar = list(col = 4), 
     horiz = TRUE, type = "tr")
R语言可视化(二十九):聚类树图绘制,第9张

使用ggdendro包绘制聚类树图

# 安装并加载所需的R包
#install.packages('ggdendro')
library(ggdendro)
library(ggplot2)

# 层次聚类
hc <- hclust(dist(USArrests), "ave")
hc
## 
## Call:
## hclust(d = dist(USArrests), method = "ave")
## 
## Cluster method   : average 
## Distance         : euclidean 
## Number of objects: 50

# Demonstrate plotting directly from object class hclust
# 使用ggdendrogram函数绘制聚类树
ggdendrogram(hc)
R语言可视化(二十九):聚类树图绘制,第10张
# 旋转90度
ggdendrogram(hc, rotate = TRUE)
R语言可视化(二十九):聚类树图绘制,第11张
# demonstrate converting hclust to dendro using dendro_data first
hcdata <- dendro_data(hc, type = "triangle")
hcdata
## $segments
##            x          y      xend      yend
## 1  17.994141 152.313999  4.765625 77.605024
## 2   4.765625  77.605024  1.500000 38.527912
## 3   1.500000  38.527912  1.000000  0.000000
## 4   1.500000  38.527912  2.000000  0.000000
## 5   4.765625  77.605024  8.031250 44.283922
## 6   8.031250  44.283922  3.875000 28.012211
## 7   3.875000  28.012211  3.000000  0.000000
## 8   3.875000  28.012211  4.750000 15.453120
## 9   4.750000  15.453120  4.000000  0.000000
## 10  4.750000  15.453120  5.500000 13.896043
## 11  5.500000  13.896043  5.000000  0.000000
## 12  5.500000  13.896043  6.000000  0.000000
## 13  8.031250  44.283922 12.187500 39.394633
## 14 12.187500  39.394633  9.625000 26.363428
## 15  9.625000  26.363428  7.750000 16.891499
## 16  7.750000  16.891499  7.000000  0.000000
## 17  7.750000  16.891499  8.500000 15.454449
## 18  8.500000  15.454449  8.000000  0.000000
## 19  8.500000  15.454449  9.000000  0.000000
## 20  9.625000  26.363428 11.500000 18.417331
## 21 11.500000  18.417331 10.500000  6.236986
## 22 10.500000   6.236986 10.000000  0.000000
## 23 10.500000   6.236986 11.000000  0.000000
## 24 11.500000  18.417331 12.500000 13.297368
## 25 12.500000  13.297368 12.000000  0.000000
## 26 12.500000  13.297368 13.000000  0.000000
## 27 12.187500  39.394633 14.750000 28.095803
## 28 14.750000  28.095803 14.000000  0.000000
## 29 14.750000  28.095803 15.500000 21.167192
## 30 15.500000  21.167192 15.000000  0.000000
## 31 15.500000  21.167192 16.000000  0.000000
## 32 17.994141 152.313999 31.222656 89.232093
## 33 31.222656  89.232093 23.796875 44.837933
## 34 23.796875  44.837933 20.343750 26.713777
## 35 20.343750  26.713777 17.937500 16.425489
## 36 17.937500  16.425489 17.000000  0.000000
## 37 17.937500  16.425489 18.875000 12.878100
## 38 18.875000  12.878100 18.000000  0.000000
## 39 18.875000  12.878100 19.750000 10.736739
## 40 19.750000  10.736739 19.000000  0.000000
## 41 19.750000  10.736739 20.500000  7.355270
## 42 20.500000   7.355270 20.000000  0.000000
## 43 20.500000   7.355270 21.000000  0.000000
## 44 20.343750  26.713777 22.750000 22.595978
## 45 22.750000  22.595978 22.000000  0.000000
## 46 22.750000  22.595978 23.500000 11.456439
## 47 23.500000  11.456439 23.000000  0.000000
## 48 23.500000  11.456439 24.000000  0.000000
## 49 23.796875  44.837933 27.250000 29.054195
## 50 27.250000  29.054195 25.750000 20.198479
## 51 25.750000  20.198479 25.000000  0.000000
## 52 25.750000  20.198479 26.500000 12.614278
## 53 26.500000  12.614278 26.000000  0.000000
## 54 26.500000  12.614278 27.000000  0.000000
## 55 27.250000  29.054195 28.750000 23.972143
## 56 28.750000  23.972143 28.000000  0.000000
## 57 28.750000  23.972143 29.500000 14.501034
## 58 29.500000  14.501034 29.000000  0.000000
## 59 29.500000  14.501034 30.000000  0.000000
## 60 31.222656  89.232093 38.648438 54.746831
## 61 38.648438  54.746831 34.437500 20.598507
## 62 34.437500  20.598507 31.875000 15.026107
## 63 31.875000  15.026107 31.000000  0.000000
## 64 31.875000  15.026107 32.750000 12.438692
## 65 32.750000  12.438692 32.000000  0.000000
## 66 32.750000  12.438692 33.500000  3.834058
## 67 33.500000   3.834058 33.000000  0.000000
## 68 33.500000   3.834058 34.000000  0.000000
## 69 34.437500  20.598507 37.000000 15.122897
## 70 37.000000  15.122897 35.500000  6.637771
## 71 35.500000   6.637771 35.000000  0.000000
## 72 35.500000   6.637771 36.000000  0.000000
## 73 37.000000  15.122897 38.500000 13.352260
## 74 38.500000  13.352260 37.500000  3.929377
## 75 37.500000   3.929377 37.000000  0.000000
## 76 37.500000   3.929377 38.000000  0.000000
## 77 38.500000  13.352260 39.500000  8.027453
## 78 39.500000   8.027453 39.000000  0.000000
## 79 39.500000   8.027453 40.000000  0.000000
## 80 38.648438  54.746831 42.859375 41.094765
## 81 42.859375  41.094765 41.000000  0.000000
## 82 42.859375  41.094765 44.718750 33.117815
## 83 44.718750  33.117815 42.750000 10.771175
## 84 42.750000  10.771175 42.000000  0.000000
## 85 42.750000  10.771175 43.500000  8.537564
## 86 43.500000   8.537564 43.000000  0.000000
## 87 43.500000   8.537564 44.000000  0.000000
## 88 44.718750  33.117815 46.687500 27.779904
## 89 46.687500  27.779904 45.500000 13.044922
## 90 45.500000  13.044922 45.000000  0.000000
## 91 45.500000  13.044922 46.000000  0.000000
## 92 46.687500  27.779904 47.875000 18.993398
## 93 47.875000  18.993398 47.000000  0.000000
## 94 47.875000  18.993398 48.750000 10.184218
## 95 48.750000  10.184218 48.000000  0.000000
## 96 48.750000  10.184218 49.500000  2.291288
## 97 49.500000   2.291288 49.000000  0.000000
## 98 49.500000   2.291288 50.000000  0.000000
## 
## $labels
##     x y          label
## 1   1 0        Florida
## 2   2 0 North Carolina
## 3   3 0     California
## 4   4 0       Maryland
## 5   5 0        Arizona
## 6   6 0     New Mexico
## 7   7 0       Delaware
## 8   8 0        Alabama
## 9   9 0      Louisiana
## 10 10 0       Illinois
## 11 11 0       New York
## 12 12 0       Michigan
## 13 13 0         Nevada
## 14 14 0         Alaska
## 15 15 0    Mississippi
## 16 16 0 South Carolina
## 17 17 0     Washington
## 18 18 0         Oregon
## 19 19 0        Wyoming
## 20 20 0       Oklahoma
## 21 21 0       Virginia
## 22 22 0   Rhode Island
## 23 23 0  Massachusetts
## 24 24 0     New Jersey
## 25 25 0       Missouri
## 26 26 0       Arkansas
## 27 27 0      Tennessee
## 28 28 0        Georgia
## 29 29 0       Colorado
## 30 30 0          Texas
## 31 31 0          Idaho
## 32 32 0       Nebraska
## 33 33 0       Kentucky
## 34 34 0        Montana
## 35 35 0           Ohio
## 36 36 0           Utah
## 37 37 0        Indiana
## 38 38 0         Kansas
## 39 39 0    Connecticut
## 40 40 0   Pennsylvania
## 41 41 0         Hawaii
## 42 42 0  West Virginia
## 43 43 0          Maine
## 44 44 0   South Dakota
## 45 45 0   North Dakota
## 46 46 0        Vermont
## 47 47 0      Minnesota
## 48 48 0      Wisconsin
## 49 49 0           Iowa
## 50 50 0  New Hampshire
## 
## $leaf_labels
## NULL
## 
## $class
## [1] "hclust"
## 
## attr(,"class")
## [1] "dendro"

ggdendrogram(hcdata, rotate = TRUE) + 
        labs(title = "Dendrogram in ggplot2")
R语言可视化(二十九):聚类树图绘制,第12张

使用ggraph包绘制聚类树图

# 安装并加载所需的R包
#install.packages("ggraph")
library(ggraph)
## Warning: package 'ggraph' was built under R version 3.6.3
library(igraph)
## 
## Attaching package: 'igraph'
## The following objects are masked from 'package:stats':
## 
##     decompose, spectrum
## The following object is masked from 'package:base':
## 
##     union
library(tidyverse)
## -- Attaching packages ------------------------------------- tidyverse 1.2.1 --
## √ tibble  2.1.3     √ purrr   0.3.2
## √ tidyr   1.1.2     √ dplyr   1.0.2
## √ readr   1.3.1     √ stringr 1.4.0
## √ tibble  2.1.3     √ forcats 0.4.0
## Warning: package 'tidyr' was built under R version 3.6.3
## Warning: package 'dplyr' was built under R version 3.6.3
## -- Conflicts ---------------------------------------- tidyverse_conflicts() --
## x dplyr::as_data_frame() masks tibble::as_data_frame(), igraph::as_data_frame()
## x purrr::compose()       masks igraph::compose()
## x tidyr::crossing()      masks igraph::crossing()
## x dplyr::filter()        masks stats::filter()
## x dplyr::groups()        masks igraph::groups()
## x dplyr::lag()           masks stats::lag()
## x purrr::simplify()      masks igraph::simplify()
library(RColorBrewer)
theme_set(theme_void())
## Warning: New theme missing the following elements: axis.title.x,
## axis.title.x.top, axis.title.y, axis.title.y.right, axis.text.x,
## axis.text.x.top, axis.text.y, axis.text.y.right, axis.ticks, axis.line,
## axis.line.x, axis.line.y, legend.background, legend.margin, legend.spacing,
## legend.spacing.x, legend.spacing.y, legend.key, legend.key.height,
## legend.key.width, legend.text.align, legend.title.align, legend.direction,
## legend.justification, legend.box.margin, legend.box.background,
## legend.box.spacing, panel.background, panel.border, panel.spacing.x,
## panel.spacing.y, panel.grid, panel.grid.minor, plot.background,
## strip.background, strip.placement, strip.text.x, strip.text.y

# 构建示例数据
# data: edge list
d1 <- data.frame(from="origin", to=paste("group", seq(1,7), sep=""))
head(d1)
##     from     to
## 1 origin group1
## 2 origin group2
## 3 origin group3
## 4 origin group4
## 5 origin group5
## 6 origin group6

d2 <- data.frame(from=rep(d1$to, each=7), to=paste("subgroup", seq(1,49), sep="_"))
head(d2)
##     from         to
## 1 group1 subgroup_1
## 2 group1 subgroup_2
## 3 group1 subgroup_3
## 4 group1 subgroup_4
## 5 group1 subgroup_5
## 6 group1 subgroup_6

edges <- rbind(d1, d2)
edges
##      from          to
## 1  origin      group1
## 2  origin      group2
## 3  origin      group3
## 4  origin      group4
## 5  origin      group5
## 6  origin      group6
## 7  origin      group7
## 8  group1  subgroup_1
## 9  group1  subgroup_2
## 10 group1  subgroup_3
## 11 group1  subgroup_4
## 12 group1  subgroup_5
## 13 group1  subgroup_6
## 14 group1  subgroup_7
## 15 group2  subgroup_8
## 16 group2  subgroup_9
## 17 group2 subgroup_10
## 18 group2 subgroup_11
## 19 group2 subgroup_12
## 20 group2 subgroup_13
## 21 group2 subgroup_14
## 22 group3 subgroup_15
## 23 group3 subgroup_16
## 24 group3 subgroup_17
## 25 group3 subgroup_18
## 26 group3 subgroup_19
## 27 group3 subgroup_20
## 28 group3 subgroup_21
## 29 group4 subgroup_22
## 30 group4 subgroup_23
## 31 group4 subgroup_24
## 32 group4 subgroup_25
## 33 group4 subgroup_26
## 34 group4 subgroup_27
## 35 group4 subgroup_28
## 36 group5 subgroup_29
## 37 group5 subgroup_30
## 38 group5 subgroup_31
## 39 group5 subgroup_32
## 40 group5 subgroup_33
## 41 group5 subgroup_34
## 42 group5 subgroup_35
## 43 group6 subgroup_36
## 44 group6 subgroup_37
## 45 group6 subgroup_38
## 46 group6 subgroup_39
## 47 group6 subgroup_40
## 48 group6 subgroup_41
## 49 group6 subgroup_42
## 50 group7 subgroup_43
## 51 group7 subgroup_44
## 52 group7 subgroup_45
## 53 group7 subgroup_46
## 54 group7 subgroup_47
## 55 group7 subgroup_48
## 56 group7 subgroup_49

# We can add a second data frame with information for each node!
name <- unique(c(as.character(edges$from), as.character(edges$to)))
vertices <- data.frame(
        name=name,
        group=c( rep(NA,8) ,  rep( paste("group", seq(1,7), sep=""), each=7)),
        cluster=sample(letters[1:4], length(name), replace=T),
        value=sample(seq(10,30), length(name), replace=T)
)
head(vertices)
##     name group cluster value
## 1 origin  <NA>       d    18
## 2 group1  <NA>       c    17
## 3 group2  <NA>       c    15
## 4 group3  <NA>       b    23
## 5 group4  <NA>       b    28
## 6 group5  <NA>       d    30

# Create a graph object
mygraph <- graph_from_data_frame( edges, vertices=vertices)
mygraph
## IGRAPH f42fbd3 DN-- 57 56 -- 
## + attr: name (v/c), group (v/c), cluster (v/c), value (v/n)
## + edges from f42fbd3 (vertex names):
##  [1] origin->group1      origin->group2      origin->group3     
##  [4] origin->group4      origin->group5      origin->group6     
##  [7] origin->group7      group1->subgroup_1  group1->subgroup_2 
## [10] group1->subgroup_3  group1->subgroup_4  group1->subgroup_5 
## [13] group1->subgroup_6  group1->subgroup_7  group2->subgroup_8 
## [16] group2->subgroup_9  group2->subgroup_10 group2->subgroup_11
## [19] group2->subgroup_12 group2->subgroup_13 group2->subgroup_14
## [22] group3->subgroup_15 group3->subgroup_16 group3->subgroup_17
## + ... omitted several edges

# 使用ggraph函数绘制聚类树图
ggraph(mygraph, layout = 'dendrogram') + 
        geom_edge_diagonal() 
R语言可视化(二十九):聚类树图绘制,第13张
# 绘制圆形的聚类树
ggraph(mygraph, layout = 'dendrogram', circular = TRUE) + 
        geom_edge_diagonal()
R语言可视化(二十九):聚类树图绘制,第14张
# 添加节点的标签,形状和信息
ggraph(mygraph, layout = 'dendrogram') + 
        geom_edge_diagonal() +
        geom_node_text(aes( label=name, filter=leaf, color=group) , angle=90 , hjust=1, nudge_y=-0.1) +
        geom_node_point(aes(filter=leaf, size=value, color=group) , alpha=0.6) +
        ylim(-.6, NA) +
        theme(legend.position="none")
R语言可视化(二十九):聚类树图绘制,第15张
# 构建测试数据集
# create a data frame giving the hierarchical structure of your individuals
d1=data.frame(from="origin", to=paste("group", seq(1,10), sep=""))
d2=data.frame(from=rep(d1$to, each=10), to=paste("subgroup", seq(1,100), sep="_"))
edges=rbind(d1, d2)

# create a vertices data.frame. One line per object of our hierarchy
vertices = data.frame(
        name = unique(c(as.character(edges$from), as.character(edges$to))) , 
        value = runif(111)
) 
# Let's add a column with the group of each name. It will be useful later to color points
vertices$group = edges$from[ match( vertices$name, edges$to ) ]

#Let's add information concerning the label we are going to add: angle, horizontal adjustement and potential flip
#calculate the ANGLE of the labels
vertices$id=NA
myleaves=which(is.na( match(vertices$name, edges$from) ))
nleaves=length(myleaves)
vertices$id[ myleaves ] = seq(1:nleaves)
vertices$angle= 90 - 360 * vertices$id / nleaves

# calculate the alignment of labels: right or left
# If I am on the left part of the plot, my labels have currently an angle < -90
vertices$hjust<-ifelse( vertices$angle < -90, 1, 0)

# flip angle BY to make them readable
vertices$angle<-ifelse(vertices$angle < -90, vertices$angle+180, vertices$angle)

# 查看测试数据
head(edges)
##     from     to
## 1 origin group1
## 2 origin group2
## 3 origin group3
## 4 origin group4
## 5 origin group5
## 6 origin group6

head(vertices)
##     name      value  group id angle hjust
## 1 origin 0.47071871   <NA> NA    NA    NA
## 2 group1 0.35606308 origin NA    NA    NA
## 3 group2 0.73638466 origin NA    NA    NA
## 4 group3 0.35833810 origin NA    NA    NA
## 5 group4 0.09394485 origin NA    NA    NA
## 6 group5 0.34061367 origin NA    NA    NA

# Create a graph object
mygraph <- graph_from_data_frame( edges, vertices=vertices )

# Make the plot
ggraph(mygraph, layout = 'dendrogram', circular = TRUE) + 
        geom_edge_diagonal(colour="grey") + #设置节点边的颜色
        # 设置节点的标签,字体大小,文本注释信息
        geom_node_text(aes(x = x*1.15, y=y*1.15, filter = leaf, label=name, angle = angle, hjust=hjust*0.4, colour=group), size=2.5, alpha=1) +
        # 设置节点的大小,颜色和透明度
        geom_node_point(aes(filter = leaf, x = x*1.07, y=y*1.07, colour=group, size=value, alpha=0.2)) +
        # 设置颜色的画板
        scale_colour_manual(values= rep( brewer.pal(9,"Paired") , 30)) +
        # 设置节点大小的范围
        scale_size_continuous( range = c(1,10) ) +
        theme_void() +
        theme(
                legend.position="none",
                plot.margin=unit(c(0,0,0,0),"cm"),
        ) +
        expand_limits(x = c(-1.3, 1.3), y = c(-1.3, 1.3))
R语言可视化(二十九):聚类树图绘制,第16张
sessionInfo()
## R version 3.6.0 (2019-04-26)
## Platform: x86_64-w64-mingw32/x64 (64-bit)
## Running under: Windows 10 x64 (build 18363)
## 
## Matrix products: default
## 
## locale:
## [1] LC_COLLATE=Chinese (Simplified)_China.936 
## [2] LC_CTYPE=Chinese (Simplified)_China.936   
## [3] LC_MONETARY=Chinese (Simplified)_China.936
## [4] LC_NUMERIC=C                              
## [5] LC_TIME=Chinese (Simplified)_China.936    
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
##  [1] RColorBrewer_1.1-2 forcats_0.4.0      stringr_1.4.0     
##  [4] dplyr_1.0.2        purrr_0.3.2        readr_1.3.1       
##  [7] tidyr_1.1.2        tibble_2.1.3       tidyverse_1.2.1   
## [10] igraph_1.2.4.1     ggraph_2.0.3       ggplot2_3.2.0     
## [13] ggdendro_0.1.22   
## 
## loaded via a namespace (and not attached):
##  [1] ggrepel_0.8.1      Rcpp_1.0.5         lubridate_1.7.4   
##  [4] lattice_0.20-38    assertthat_0.2.1   digest_0.6.20     
##  [7] ggforce_0.3.2      R6_2.4.0           cellranger_1.1.0  
## [10] backports_1.1.4    evaluate_0.14      httr_1.4.0        
## [13] pillar_1.4.2       rlang_0.4.7        lazyeval_0.2.2    
## [16] readxl_1.3.1       rstudioapi_0.10    rmarkdown_1.13    
## [19] labeling_0.3       polyclip_1.10-0    munsell_0.5.0     
## [22] broom_0.5.2        compiler_3.6.0     modelr_0.1.4      
## [25] xfun_0.8           pkgconfig_2.0.2    htmltools_0.3.6   
## [28] tidyselect_1.1.0   gridExtra_2.3      graphlayouts_0.7.0
## [31] viridisLite_0.3.0  crayon_1.3.4       withr_2.1.2       
## [34] MASS_7.3-51.4      grid_3.6.0         nlme_3.1-139      
## [37] jsonlite_1.6       gtable_0.3.0       lifecycle_0.2.0   
## [40] magrittr_1.5       scales_1.0.0       cli_1.1.0         
## [43] stringi_1.4.3      farver_1.1.0       viridis_0.5.1     
## [46] xml2_1.2.0         generics_0.0.2     vctrs_0.3.2       
## [49] tools_3.6.0        glue_1.4.2         tweenr_1.0.1      
## [52] hms_0.4.2          yaml_2.2.0         colorspace_1.4-1  
## [55] tidygraph_1.2.0    rvest_0.3.4        knitr_1.23        
## [58] haven_2.3.1

https://www.xamrdz.com/backend/3e21935111.html

相关文章: