我是靠谱客的博主 忧伤龙猫,最近开发中收集的这篇文章主要介绍可视化应用实战案例:绘制交互式+pdf+png等多格式桑基图,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

写在前面

桑基图挺好看的,用于观察不同门类之间的从属关系,并且绘制很漂亮的结构图,当然可以用于很多个地方。这里我们用微生物组数据的phyloseq对象,很快很方便的为大家构建一个桑基图。所以如果你有phyloseq对象直接来试试吧。

实战

导入需要R包

library(ggClusterNet)
library(tidyverse)
library(phyloseq)
library(tidyverse)
library(viridis)
library(patchwork)
library(networkD3)
data(ps)

我们将微生物数据按照属水平合并,然后取前五十个丰度最高的属,然后去除Unassigned”) %>%  vegan_tax() %>%的微生物。注意这几个函数存在于ggClusterNet中,注意下载安装。

tax = ps %>% 
  ggClusterNet::tax_glom_wt(ranks = 6) %>%
  filter_OTU_ps(50) %>%
  subset_taxa(Genus  != "Unassigned") %>%
  vegan_tax() %>%
  as.data.frame()
head(tax)

构建连接,一个源—-一个目标。标记一下不同分类等级的标签。

id2 = c("k","p","c","o","f","g")
dat = NULL
for (i in 1:5) {
  dat <- tax[,c(i,i+1)] %>% distinct(.keep_all = TRUE) 
  colnames(dat) = c("so

最后

以上就是忧伤龙猫为你收集整理的可视化应用实战案例:绘制交互式+pdf+png等多格式桑基图的全部内容,希望文章能够帮你解决可视化应用实战案例:绘制交互式+pdf+png等多格式桑基图所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(60)

评论列表共有 0 条评论

立即
投稿
返回
顶部