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)