可视化应用实战案例：metacoder-相关进化树图的绘制

Posted 2022-02-04 文宇肃然

metacoder 包学习

安装和导入R包

#--选择安装cran或者github中的R包
# if(!require(metacoder))install.packages("metacoder")
# if(!require(metacoder))devtools::install_github("grunwaldlab/metacoder")

#--导入R包，开始学习
library(metacoder)
#---构造自己的数据
library(ggClusterNet)
library(tidyverse)
library(phyloseq)

学习数据结构

metacoder包使用的数据是tibble格式的数据框，一个OTU表格，包含原始OTU的count，并未抽平。

data(ps)
tax = ps %>%
  subset_taxa(Kingdom == "Bacteria") %>%
  filter_OTU_ps(100) %>%
  vegan_tax() %>%
  as.data.frame()

otu = ps %>%
  subset_taxa(Kingdom == "Bacteria") %>%
  filter_OTU_ps(100) %>%
  vegan_otu() %>% t() %>%
  as.data.frame()
head(otu)

##         KO1  KO2  KO3  KO4  KO5  KO6  OE1  OE2  OE3  OE4  OE5  OE6  WT1  WT2
## ASV_1  1113 1968  816 1372 1062 1087 1270 1637 1368  962 1247 1017 2345 2538
## ASV_28  253  303   54  43