教学手册:学习手册
MaAsLin2
# BiocManager::install("Maaslin2",force = TRUE)
library(Maaslin2)
# 用的是相对丰度,行名为-ID行样本,列为细菌
input_data = system.file("extdata", "HMP2_taxonomy.tsv", package="Maaslin2") # The abundance table file
input_data
# 行名为-ID行样本
input_metadata = system.file("extdata", "HMP2_metadata.tsv", package="Maaslin2") # The metadata table file
input_metadata
df_input_data = read.table(file = input_data,
header = TRUE,