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

1
2
3
4
5
6
7
8
9
10
#--选择安装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)

1
2
3
4
5
6
7
8
9
10
11
12
13
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)

1
2
3
##         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