MaAsLin2用户手册

MaAsLin2是MaAsLin的下一代。

MaAsLin2是用于发现临床元数据和微生物元组学特征之间的多变量关联的综合R包。MaAsLin2依赖于一般线性模型来适应大多数现代流行病学研究设计,包括多种分析方法(支持多协变量和重复测量)、过滤、规范化和转换选项,以便为您的特定研究定制分析。

如果您使用MaAsLin2软件,请引用我们的手稿:Mallick et al.(2020+)。“人口规模元组学研究中的多变量关联”(准备中)。

如果您有问题,请发邮件给谷歌组MaAsLin用户


描述

MaAsLin2在人群规模的流行病学研究中发现了微生物组多组学特征和复杂元数据之间的关联。该软件包括多种分析方法(支持多个协变量和重复度量)、过滤、规范化和转换选项,以定制您的特定研究的分析。

需求

MaAsLin2是一个R包,可以在命令行上运行,也可以作为R函数运行。

安装

MaAsLin2可以从命令行运行,也可以作为R函数运行。如果只从命令行运行,则不需要安装MaAsLin2包,但需要安装MaAsLin2依赖项。

从命令行

  1. 下载源代码:MaAsLin2.tar.gz
  2. 解压下载:
    • $ tar xzf maaslin2.tar.gz
  3. 安装Bioconductor依赖项edgeR和metagenomeSeq。
  4. 安装CRAN依赖项:
    • $ R - q - e”install.packages (c(‘lmerTest’,‘pbapply’,‘车’,‘dplyr’,“素食主义者”,“化学计量学”,‘ggplot2’,‘pheatmap’,‘哈’,‘日志’,‘data.table’,‘MuMIn’,‘glmmTMB’,‘质量’,‘cplm’,‘pscl’),回购= ' http://cran.r-project.org ')”
  5. 安装MaAsLin2包(只有r,如果作为r函数运行必须安装):
    • $ R CMD安装maaslin2

从R

先安装Bioconductor,再安装Maaslin2

如果(!install.packages("BiocManager") BiocManager::install("Maaslin2")

如何跑步

MaAsLin2可以从命令行运行,也可以作为R函数运行。两个方法都需要相同的参数,具有相同的选项,并使用相同的默认设置。

输入文件

MaAsLin2需要两个输入文件。

  1. 数据(或特性)文件
    • 该文件是用制表符分隔的。
    • 将特征格式化为列,将示例格式化为行。
    • 这种形式的转置也是可以的。
    • 该文件中的可能特征包括分类法或基因。
  2. 元数据文件
    • 该文件是用制表符分隔的。
    • 将特征格式化为列,将示例格式化为行。
    • 这种形式的转置也是可以的。
    • 此文件中可能的元数据包括性别或年龄。

数据文件可以包含元数据文件中不包含的样例(以及相反的情况)。对于这两种情况,没有包含在这两个文件中的样本将从分析中删除。另外,示例在两个文件中的顺序也不需要相同。

注意:如果将MaAsLin2作为函数运行,则数据和元数据输入可以为类型data.frame而不是文件的路径。

输出文件

MaAsLin2生成两种类型的输出文件:数据和可视化。

  1. 数据输出文件
    • all_results.tsv
      • 这包括与返回的data.frame相同的数据。
      • 该文件包含按q-value递增顺序排列的所有结果。
      • 第一列是元数据和特性名称。
      • 接下来的两列是来自模型的值和系数。
      • 下一列是模型的标准差。
      • N列是数据点的总数。
      • N.not.zero列为非零数据点的总和。
      • 计算的pvalue是倒数第二列。
      • qvalue是用p.adjust用修正法。
    • significant_results.tsv
      • 该文件是第一个文件中结果的子集。
      • 它只包括与q值<=到阈值的关联。
    • residuals.rds
      • 该文件包含每个特征的残差数据帧。
    • fitted.rds
      • 该文件包含一个数据帧,其中包含每个特征的拟合值。
    • ranef.rds
      • 该文件包含为每个特征提取随机效果的数据帧(当指定随机效果时)。
    • maaslin2.log
      • 该文件包含运行的所有日志信息。
      • 它包括所有设置、警告、错误和运行的步骤。
  2. 可视化输出文件
    • heatmap.pdf
      • 该文件包含重要关联的热图。
    • [a - z和0 - 9]+ . pdf
      • 为每个重要的关联生成一个图。
      • 散点图用于连续元数据。
      • 箱形图用于分类数据。
      • 绘制的数据点经过归一化、滤波和变换。

运行演示

示例输入文件可以在本月/ extdataMaAsLin2源文件的文件夹。提供的文件是由HMP2数据生成的,可从该数据下载https://ibdmdb.org/

HMP2_taxonomy.tsv:是一个以选项卡分隔的文件,列为物种,行为样本。它是分类文件的一个子集,所以它只包括所有样本的物种丰度。

HMP2_metadata.tsv:是一个以制表符分隔的文件,其中示例为行,元数据为列。它是元数据文件的一个子集,因此它只包含一些字段。

命令行

Maaslin2美元。R——transform=AST——fixed_effects="diagnosis,dysbiosisnonIBD,dysbiosisUC,dysbiosisCD,抗生素,年龄"——random_effects="site,subject"——normalization=NONE——standardization =FALSE inst/extdata/HMP2_taxonomy. R——transform=AST——fixed_effects="diagnosis,dysbiosisnonIBD,dysbiosisUC,dysbiosisCD,抗生素,年龄"tsv本月/ extdata / HMP2_metadata。tsv demo_output

  • 确保提供MaAsLin2可执行文件的完整路径(即/R/ MaAsLin2 .R)。
  • 在demo命令中:
    • HMP2_taxonomy.tsv是数据(或特性)文件的路径吗
    • HMP2_metadata.tsv是元数据文件的路径吗
    • demo_output是写入输出的文件夹的路径吗

在R

library(Maaslin2) input_data <- system. input_data。文件(“extdata”、“HMP2_taxonomy。tsv', package="Maaslin2") input_metadata <-system. tsv'文件(“extdata”、“HMP2_metadata。tsv', package=" NONE ") fit_data <- Maaslin2(input_data, input_metadata, 'demo_output', transform =" AST", fixed_effects = c('diagnosis', 'dysbiosisnonIBD','dysbiosisUC','dysbiosisCD', '抗生素','年龄'),random_effects = c('site', 'subject'), reference =" diagnosis,nonIBD", normalization = 'NONE', standardization = FALSE)
# 22-04-26 17:00:41 INFO::写入函数参数到日志文件## 22-04-26 17:00:41 INFO::验证所选选项是否有效## 22-04-26 17:00:41 INFO::确定输入文件的格式## 22-04-26 17:00:41 INFO::输入格式是数据样本作为行,元数据样本作为行## 22-04-26 17:00:41 INFO::随机效果的公式:expr ~ (1 | site) + (1 | subject) ## 22-04-26 17:00:41 INFO::固定效应公式:expr ~诊断+ dysbiosisnonIBD + dysbiosisUC + dysbiosisCD +抗生素+年龄## 22-04-26 17:00:41 INFO::根据最小丰度和最小流行率筛选数据## 22-04-26 17:00:41 INFO::数据中的总样本:1595 ## 22-04-26 17:00:41 INFO::需要的最小样本和最小丰度对于一个不被过滤的特征:159.500000 ## 22-04-26 17:00:41 INFO::总过滤特征:0 ## 22-04-26 17:00:41 INFO::来自丰度和流行过滤的过滤特征名:## 22-04-26 17:00:41 INFO::来自方差过滤的过滤特征名:0 ## 22-04-26 17:00:41 INFO::运行选定的归一化方法:NONE ## 22-04-26 17:00:41 INFO::绕过z-score应用程序到元数据## 22-04-26 17:00:41 INFO::运行选定的转换方法:运行选定的分析方法:拟合模型到特征号1,bifidobacterium . teenens
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号2,双歧杆菌。信息:拟合模型到特征号3,双歧杆菌
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号4,双歧杆菌。假atenulatum
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的功能号5,Collinsella。拟合模型的特征号6,拟杆菌科
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号7,拟杆菌。信息:拟合模型的特征号8,拟杆菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号9,拟杆菌
## boundary (singular) fit:参见help('isSingular')
##警告:模型未能收敛于1个负特征值:-5.6e+00
信息:拟合模型的特征号10,拟杆菌。拟合模型的特征号11,拟杆菌。finegoldii
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号12,拟杆菌。脆弱
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号13,拟杆菌类
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号14,拟杆菌
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号15,拟杆菌。信息:拟合模型的特征号16,拟杆菌科
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号17,拟杆菌。信息:拟合模型的特征号18,拟杆菌。thetaiotaomicron
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型特征号19,Bacteroides.uniformis
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号20,拟杆菌类
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的功能编号21,拟杆菌。木聚糖isolvens
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号22,拟杆菌。INFO::拟合模型到功能号23,barnesiella . testinihominis ## 22-04-26 17:00:46 INFO::拟合模型到功能号24,Coprobacter。信息:拟合模型的特征号25,臭杆菌。内脏
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号26,副杆菌类
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号27,副杆菌类。信息::拟合模型的特征号28,副杆菌类。信息::拟合模型的特征号29,parabacteroides .未分类
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的功能号30,Paraprevotella。信息::适合模型的功能号31,Paraprevotella。信息:拟合模型的特征号32,普雷沃氏菌.copri
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的功能号33,Alistipes.finegoldii
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的功能号34,Alistipes.onderdonkii
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的功能号35,Alistipes.putredinis
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的功能编号36,Alistipes.shahii
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的功能号37,Alistipes。信息:拟合模型的特征号38,链球菌。信息::拟合模型的特征号39,梭状芽胞杆菌
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的特征编号40,梭状芽孢杆菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征编号41,梭状芽胞杆菌
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的功能号码42,梭状芽孢杆菌
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的功能编号43,梭状芽孢杆菌。信息:拟合模型的特征号44,梭状芽孢杆菌。信息:拟合模型的特征号45,梭状芽孢杆菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号46,黄酮ifractor.plautii
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征编号47,真菌性细菌。信息:拟合模型的特征编号48,真菌性细菌。信息:拟合模型的特征号49,真菌性细菌。信息::拟合模型的特征号50,真菌性细菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征编号为51,Eubacterium.sp.3.1.31
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征编号52,真细菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号码53,瘤胃球菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号54,瘤胃球菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号55,瘤胃球菌。信息::适合模型的特征号56,Coprococcus.comes ## 22-04-26 17:00:51信息::适合模型的特征号57,Dorea。信息::拟合模型的特征号58,lachnospiraceae . bacteria .1.1.57 faa
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号59,lachnospiraceae . bacteria .3.1.46 faa
## boundary (singular) fit:参见help('isSingular')
信息::适合的模型特征号60,roseburia .人族
## boundary (singular) fit:参见help('isSingular')
## 22-04-26 17:00:51 INFO::拟合模型到功能号61,Roseburia.inulinivorans
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的功能号63,罗斯伯里亚。未分类
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的功能号64,Oscillibacter。信息:拟合模型的特征号65,peptostreptococcaceae . nonname .unclassified
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号66,Faecalibacterium。信息::拟合模型的特征号67,瘤胃球菌.bromii
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号68,瘤胃球菌。信息::拟合模型的特征号69,瘤胃球菌。信息:拟合模型的特征号70,亚doligranulum。信息:拟合模型的特征号71,Coprobacillus。信息:拟合模型的特征号码72,酸氨基球菌。未分类
##在checkConv(attr(opt, " derived "), opt$par, ctrl = control$checkConv,: ## Model failed to convergence with max|grad| = 0.00291214 (tol = 0.002, component 1)
信息::适合模型的功能号码73,Dialister.invisus
## boundary (singular) fit:参见help('isSingular')
##警告:模型未能收敛于1个负特征值:-2.1e+02
信息::拟合模型的特征号74,细孔菌
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的特征号75,Veillonella.dispar
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号76,绒毛
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的功能号码77,Veillonella。信息::拟合模型到特征号78,burkholderials . bacteria .1.1.47
## boundary (singular) fit:参见help('isSingular')
信息::适合模型的功能号80,Sutterella。信息:拟合模型的特征号81,Bilophila。信息:拟合模型的特征号82,大肠杆菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的特征号83,埃希氏菌。信息:拟合模型的特征号84,克雷伯菌肺炎
## boundary (singular) fit:参见help('isSingular')
信息:拟合模型的特征号85,副流感嗜血杆菌
## boundary (singular) fit:参见help('isSingular')
信息::拟合模型的功能号86,Akkermansia.muciniphila
## boundary (singular) fit:参见help('isSingular')
##警告:模型未能收敛于1负特征值:-2.2e+02
信息::拟合模型的特征号87,c2likevirv。## 22-04-26 17:00:55 INFO::写入残差到文件demo_output/residuals。rds ## 22-04-26 17:00:55 INFO::写入适合的值到文件demo_output/ fits。将提取的随机效果写入demo_output/ranef文件。rds ## 22-04-26 17:00:55 INFO::写入所有结果到文件(按增加q-values的顺序):demo_output/all_results。将重要结果(小于或等于阈值0.250000的结果)写入文件(通过增加q-values的顺序):demo_output/significant。tsv ## 2015-04-26 17:00:55 INFO::写入重要结果的热图到文件:demo_output/heatmap.pdf
##警告在xtfrm.data.frame(x):不能xtfrm数据帧
## 22-04-26 17:00:56 INFO::绘制从最重要到最不重要的关联,按元数据分组## 22-04-26 17:00:56 INFO::绘制元数据编号1的数据,dysbiosisCD ## 22-04-26 17:00:56 INFO::为分类数据,dysbiosisCD vs Faecalibacterium.prausnitzii的箱图
创建分类数据箱线图,dysbiosisCD vs Bacteroides.uniformis
创建分类数据箱线图,dysbiosisCD vs Eubacterium.rectale
创建分类数据箱线图,dysbiosisCD vs alistips .putredinis
创建分类数据箱线图,dysbiosisCD vs subdoligranum .unclassified
创建分类数据箱线图,dysbiosisCD vs .大肠杆菌
创建分类数据箱线图,dysbiosisCD vs Bacteroides.vulgatus
创建分类数据箱线图,dysbiosisCD vs . clostridium . clostridioformme
创建分类数据箱线图,dysbiosisCD vs克雷伯菌肺炎
创建分类数据箱线图,dysbiosisCD vs . Clostridium.hathewayi
创建分类数据箱线图,dysbiosisCD vs alistips .shahii
创建分类数据箱线图,dysbiosisCD vs Ruminococcus.obeum
创建分类数据箱线图,dysbiosisCD vs Roseburia.inulinivorans
创建分类数据箱线图,dysbiosisCD vs Bacteroides.thetaiotaomicron
创建分类数据箱线图,dysbiosisCD vs copprococcus .comes
创建分类数据箱线图,dysbiosisCD vs Veillonella.unclassified
创建分类数据箱线图,dysbiosisCD vs lachnospiraceae . bacteria .3.1.46 faa
创建分类数据箱线图,dysbiosisCD vs Sutterella.wadsworthensis
创建分类数据箱线图,dysbiosisCD vs臭气菌
创建分类数据箱线图,dysbiosisCD vs parabacteroides . disasonis
创建分类数据箱线图,dysbiosisCD vs Roseburia.hominis
创建分类数据箱线图,dysbiosisCD vs burkholderials . bacteria .1.1.47
创建分类数据箱线图,dysbiosisCD vs Escherichia.unclassified
创建分类数据箱线图,dysbiosisCD vs Bilophila.unclassified
创建分类数据箱线图,dysbiosisCD vs Dorea.longicatena
创建分类数据箱线图,dysbiosisCD vs aliistips .unclassified
创建分类数据箱线图,dysbiosisCD vs . Clostridium.symbiosum
创建分类数据箱线图,dysbiosisCD vs . Eubacterium.hallii
创建分类数据箱线图,dysbiosisCD vs parasutella . secrementhominis
创建分类数据箱线图,dysbiosisCD vs拟杆菌类
创建分类数据箱线图,dysbiosisCD vs alistips .finegoldii
创建分类数据箱线图,dysbiosisCD vs Eubacterium.ventriosum
创建分类数据箱线图,dysbiosisCD vs . Clostridium.nexile
创建分类数据箱线图,dysbiosisCD vs . Clostridium.leptum
创建分类数据箱线图,dysbiosisCD vs Veillonella.parvula
创建分类数据箱线图,dysbiosisCD vs拟杆菌类。卵形
创建分类数据箱线图,dysbiosisCD vs oscillibacteria .unclassified
创建分类数据箱线图,dysbiosisCD vs Bacteroides.finegoldii
创建分类数据箱线图,dysbiosisCD vs Ruminococcus.bromii
创建分类数据箱线图,dysbiosisCD vs Coprobacillus.unclassified
创建分类数据箱线图,dysbiosisCD vs Parabacteroides.merdae
创建分类数据箱线图,dysbiosisCD vs Bacteroides.eggerthii
创建分类数据箱线图,dysbiosisCD vs Veillonella.dispar
创建分类数据箱线图,dysbiosisCD vs Collinsella.aerofaciens
创建分类数据箱线图,dysbiosisCD vs拟杆菌类
创建分类数据箱线图,dysbiosisCD vs Ruminococcus.gnavus
创建分类数据箱线图,dysbiosisCD vs Paraprevotella.clara
创建分类数据箱线图,dysbiosisCD vs拟杆菌类
创建分类数据箱线图,dysbiosisCD vs . Ruminococcus.lactaris
创建分类数据箱线图,dysbiosisCD vs Paraprevotella.unclassified
创建分类数据箱线图,dysbiosisCD vs Ruminococcus.torques
创建分类数据箱线图,dysbiosisCD vs c2likevvirus .unclassified
创建分类数据箱线图,dysbiosisCD vs .faecis
创建分类数据箱线图,dysbiosisCD vs双歧杆菌
创建分类数据箱线图,dysbiosisCD vs Eubacterium.siraeum
创建分类数据箱图,dysbiosisCD vs Dialister.invisus
创建分类数据箱线图,dysbiosisCD vs parabacteroid .unclassified
为分类数据创建箱线图,dysbiosisUC vs subdoligranum .unclassified
创建分类数据箱线图,dysbiosisUC vs .普氏杆菌
创建分类数据箱线图,dysbiosisUC vs . Bacteroides.fragilis
创建分类数据箱线图,dysbiosisUC vs拟杆菌类
创建分类数据箱线图,dysbiosisUC vs Faecalibacterium.prausnitzii
创建分类数据箱线图,dysbiosisUC vs oscillibacteria .unclassified
创建分类数据箱线图,dysbiosisUC vs Bacteroides.uniformis
创建分类数据箱线图,dysbiosisUC vs . Eubacterium.hallii
创建分类数据箱线图,dysbiosisUC vs拟杆菌
创建分类数据箱线图,dysbiosisUC vs克雷伯菌肺炎
创建分类数据箱线图,dysbiosisUC vs . Coprobacillus.unclassified
创建分类数据箱线图,dysbiosisUC vs拟杆菌类
创建分类数据箱线图,dysbiosisUC vs Eubacterium.rectale
创建分类数据箱线图,dysbiosisUC vs . Ruminococcus.gnavus
创建分类数据箱线图,dysbiosisUC vs拟杆菌类。卵形
创建分类数据箱线图,dysbiosisUC vs lachnospiraceae . bacteria .3.1.46 faa
创建分类数据箱线图,dysbiosisUC vs Eubacterium.siraeum
创建分类数据箱线图,dysbiosisUC vs alistips .putredinis
创建分类数据箱线图,dysbiosisUC vs Ruminococcus.obeum
创建分类数据箱线图,dysbiosisUC vs .埃希氏大肠杆菌
创建分类数据箱线图,dysbiosisUC vs . Bacteroides.eggerthii
创建分类数据箱线图,dysbiosisUC vs . Clostridium.leptum
创建分类数据箱线图,dysbiosisUC vs barnesiella . testinihominis
创建分类数据箱线图,dysbiosisUC vs Alistipes.shahii
创建分类数据箱线图,dysbiosisUC vs Collinsella.aerofaciens
创建分类数据箱线图,dysbiosisUC vs . Clostridium.nexile
创建分类数据箱线图,dysbiosisUC vs拟杆菌类
创建分类数据箱线图,dysbiosisUC vs Parabacteroides.goldsteinii
创建分类数据箱线图,dysbiosisUC vs Eubacterium.ventriosum
创建分类数据箱线图,dysbiosisUC vs copprococcus .comes
创建分类数据箱线图,dysbiosisUC vs双歧杆菌
创建分类数据箱线图,dysbiosisUC vs . Veillonella.atypica
为分类数据,dysbiosisnonIBD vs Faecalibacterium.prausnitzii创建箱线图
创建分类数据箱线图,dysbiosisnonIBD vs Eubacterium.rectale
创建分类数据箱线图,dysbiosisnonIBD vs .普氏杆菌
创建分类数据箱线图,dysbiosisnonIBD vs Eubacterium.sp.3.1.31
创建分类数据箱线图,dysbiosisnonIBD vs Roseburia.hominis
创建分类数据箱线图,dysbiosisnonIBD vs臭气杆菌.内脏
创建分类数据箱线图,dysbiosisnonIBD vs Bilophila.unclassified
创建分类数据箱线图,ibd与大肠埃希菌
创建分类数据箱线图,dysbiosisnonIBD vs subdoligranum .unclassified
创建分类数据箱线图,dysbiosisnonIBD vs . Clostridium.leptum
创建分类数据箱线图,dysbiosisnonIBD vs . Bacteroides.fragilis
创建分类数据箱线图,dysbiosisnonIBD vs克雷伯氏菌肺炎
创建分类数据箱线图,dysbiosisnonIBD vs Ruminococcus.torques
创建分类数据箱线图,dysbiosisnonIBD vs . clostridium . clostridioformme
创建分类数据箱线图,dysbiosisnonIBD vs alistips .putredinis
创建分类数据箱线图,dysbiosisnonIBD vs Coprobacillus.unclassified
创建分类数据箱线图,dysbiosisnonIBD vs Eubacterium.eligens
创建分类数据箱线图,dysbiosisnonIBD vs Veillonella.dispar
创建分类数据箱线图,dysbiosisnonIBD vs Bacteroides.eggerthii
创建分类数据箱线图,dysbiosisnonIBD vs lachnospiraceae .3.1.46 faa
创建分类数据箱线图,dysbiosisnonIBD vs Ruminococcus.obeum
创建分类数据箱线图,dysbiosisnonIBD vs Ruminococcus.bromii
创建分类数据箱线图,dysbiosisnonIBD vs parasutella .排泄物门菌
创建分类数据箱线图,dysbiosisnonIBD vs alistips .shahii
创建分类数据箱线图,dysbiosisnonIBD vs barnesiella . testinihominis
创建分类数据箱线图,dysbiosisnonIBD vs burkholderians . bacteria .1.1.47
创建分类数据箱线图,dysbiosisnonIBD vs Paraprevotella.unclassified
为分类数据创建箱线图,抗生素vs .大肠杆菌
创建分类数据箱线图,抗生素vs .inulinivorans
创建分类数据箱线图,抗生素vs .月桂花
创建分类数据箱线图,抗生素vs . bacteroidales . bacteria .ph8
创建分类数据箱线图,抗生素vs人玫瑰buria
创建分类数据箱图,抗生素vs Dialister.invisus
创建分类数据箱线图,抗生素vs .克雷伯氏菌肺炎
创建分类数据箱线图,抗生素vs真菌性细菌
创建分类数据箱线图,抗生素vs . Sutterella.wadsworthensis
创建分类数据箱线图,抗生素vs .thetaiotaomicron
创建分类数据箱线图,抗生素vs .溃疡性瘤胃球菌
创建分类数据箱线图,抗生素vs .finegoldii
创建分类数据箱线图,抗生素vs .onderdonkii
创建分类数据箱线图,抗生素vs . Dorea.longicatena
创建分类数据箱线图,抗生素vs .假atenulatum双歧杆菌
创建分类数据箱线图,抗生素vs .双歧杆菌
创建分类数据箱线图,抗生素vs .bolteae
创建分类数据箱线图,抗生素vs .hathewayi
创建分类数据箱线图,抗生素vs . Faecalibacterium.prausnitzii
创建分类数据箱线图,抗生素vs .溴氏瘤胃球菌
创建分类数据箱线图,抗生素vs真菌性细菌
创建分类数据箱线图,抗生素vs .瘤胃球菌
创建分类数据箱线图,抗生素vs脆弱拟杆菌
创建分类数据箱线图,抗生素vs .未分类大肠杆菌
创建分类数据箱线图,抗生素vs .非分类菌
创建分类数据箱线图,抗生素vs . Eubacterium.sp.3.1.31
创建分类数据箱线图,抗生素vs腐霉杆菌
创建分类数据箱线图,抗生素vs细小细孔菌
创建分类数据箱线图,抗生素vs拟杆菌类
创建分类数据箱线图,抗生素vs副杆菌类
创建分类数据箱线图,抗生素vs副杆菌类
创建分类数据箱线图,抗生素vs .粪门菌
创建分类数据箱线图,抗生素vs .双歧杆菌
创建分类数据箱线图,抗生素vs . lachnospiraceae . bacteria .1.1.57 faa
创建分类数据箱线图,抗生素vs Akkermansia.muciniphila
创建分类数据箱线图,抗生素vs .共生梭菌
创建分类数据箱线图,抗生素vs . Paraprevotella.clara
创建分类数据箱线图,抗生素vs .粪便拟杆菌
为连续数据创建散点图,年龄vs . Clostridium.clostridioforme
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
创建散点图连续数据,年龄vs副流感嗜血杆菌
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs双歧杆菌
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs subdoligranum .unclassified
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs拟杆菌类
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs Faecalibacterium.prausnitzii
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
创建散点图的连续数据,年龄vs .共生梭菌
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
创建散点图的连续数据,年龄vs .bromii
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs Akkermansia.muciniphila
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs Veillonella.dispar
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs Collinsella.aerofaciens
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs . roseburgia . testinalis
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs拟杆菌类。thetaiotaomicron
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs Dialister.invisus
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs酸胺球菌。未分类
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs Veillonella.unclassified
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
创建散点图的连续数据,年龄vs双歧杆菌
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs aliistipes .unclassified
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为连续数据创建散点图,年龄vs bacteroidales . bacteria .ph8
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
使用公式'y ~ x'创建geom_smooth()
##警告:删除19行包含丢失的值(geom_point)。
为分类数据创建箱线图,诊断与clostridium . clostridioformme
创建分类数据箱线图,诊断与梭状芽胞杆菌
创建分类数据箱线图,诊断与Alistipes.putredinis
创建分类数据箱线图,诊断vs oscillibacteria .unclassified
创建分类数据箱线图,诊断vs aliistipes .shahii
创建分类数据箱线图,诊断vs Collinsella.aerofaciens
创建分类数据箱线图,诊断vs Akkermansia.muciniphila
创建分类数据箱线图,诊断vs .未分类coprobacillus
创建分类数据箱线图,诊断vs副杆菌类
创建分类数据箱线图,诊断vs .bromii瘤胃球菌
创建分类数据箱线图,诊断与瘤胃球菌。骨痂
创建分类数据箱线图,诊断vs . Sutterella.wadsworthensis
创建分类数据箱线图,诊断vs .溴氏瘤胃球菌
创建分类数据箱线图,诊断vs Akkermansia.muciniphila
创建分类数据箱线图,诊断vs瘤胃球菌
创建分类数据箱线图,诊断vs bacteroidales . bacteria .ph8
创建分类数据箱线图,诊断vs双歧杆菌。青少年
创建分类数据箱线图,诊断vs脆弱拟杆菌类
创建分类数据箱线图,诊断vs Alistipes.finegoldii
创建分类数据箱线图,诊断vs梭状芽孢杆菌
创建分类数据箱线图,诊断vs . Sutterella.wadsworthensis
创建分类数据箱线图,诊断vs .未分类埃希氏病
创建分类数据箱线图,诊断vs乳状瘤胃球菌
创建分类数据箱线图,诊断vs未分类Bilophila.unclassified
创建分类数据箱线图,诊断vs胃球菌
创建分类数据箱线图,诊断vs subdoligranum .unclassified
创建分类数据箱线图,诊断vs .未分类coprobacillus
创建分类数据箱线图,诊断与副杆菌类。未分类

会话信息

使用下面的命令可以显示在R中运行demo的会话信息。

sessionInfo ()
## R版本4.2.0 RC (22-04-19 r82224) ##平台:x86_64-pc-linux-gnu(64位)##运行在:Ubuntu 20.04.4 LTS ## ##矩阵产品:default ## BLAS: /home/biocbuild/bbs-3.15-bio /R/lib/libRblas. ##因此## LAPACK: /home/biocbuild/bbs-3.15-bio /R/lib/libRlapack。因此## ## locale: ## [1] LC_CTYPE=en_US。UTF-8 LC_NUMERIC= c# [3] LC_TIME=en_GB LC_COLLATE= c# [5] LC_MONETARY=en_US。utf - 8 LC_MESSAGES = en_US。UTF-8 ## [7] LC_PAPER=en_US。UTF-8 LC_NAME= c# [9] LC_ADDRESS=C LC_TELEPHONE= c# [11] LC_MEASUREMENT=en_US。UTF-8 LC_IDENTIFICATION=C ## ##附加的基本包:## [1]stats graphics grDevices utils datasets methods base ## ##其他附加的包:## [1]Maaslin2_1.10.0 ## ##通过命名空间加载(并没有附加):# [1] Rcpp_1.0.8.3 mvtnorm_1.1-3 lattice_0.20-45 ## [4] assertthat_0.2.1 digest_0.6.29 utf8_1.2.2 ## [7] R6_2.5.1 pcaPP_2.0-1 evaluate_0.15 ## [10] highrr_0.9 ggplot_3.3.5 pillar_1.7.0 ## [13] rlang_1.0.2 minqa_1.2.4 data.table_1.14.2 ## [13] nloptr_2.0.0 jquerylib_0.1.4 Matrix_1.4-1 ## [19] hash_2.2.6.2 rmarkdown_2.14 labeling_0.4.2 ## [25] splines_4.2.0 lme4_1.1-29 string_1 .4.0 ## [25] pheatmap_1.0.12 munsell_0.5.0 compiler_4.2.0 ## [28] numd_2016.8 -1.1 xfun_0.30 pkgconfig_2.0.3 ## [10]# [40] dplyr_1.0.8 withr_2.5.0 MASS_7.3-57 ## [46] gtable_0.3.0 lifecycle_1.0.1 DBI_1.1.2 ## [49] magrittr_2.0.3 pbapply_1.5-0 cli_3. 1.0 ## [55] getopt_1.20.3 lpsymphony_1.24.0 robustbase_0.95-0 ## [58] optparse_1.7.1 bslib_0.3.1 ellipsis_0.3.2 ## [61]generics_0.1.2 vctrs_0.4.1 boot_1.3-28 ## [64] RColorBrewer_1.1-3 tools_4.2.0 glue_1.6.2 ## [67] DEoptimR_1.0-11 purrr_0.3.4 parallel_4.2.0 ## [70] fastmap_1.1.0 yaml_2.3.5 colorspace_2.0-3 ## [73] knitr_1.38 sass_0.4.1

选项

运行MaAsLin2 help打印选项列表和默认设置。

Maaslin2美元。R -help用法:./R/Maaslin2.R . help选项<数据。tsv > < metadata.tsv >

选项:-h, -help显示帮助信息并退出

-b MIN_VARIANCE,——MIN_VARIANCE =MIN_VARIANCE保留方差大于[默认值:0.0]的特征,——max_significant = max_significant显著性的q值阈值[默认值:0]-a MIN_ABUNDANCE,——MIN_ABUNDANCE =MIN_ABUNDANCE每个特征的最小丰度[默认值:0]-b MIN_VARIANCE,——MIN_VARIANCE =MIN_VARIANCE0.25] -n NORMALIZATION,——NORMALIZATION =NORMALIZATION应用的规范化方法[Default: TSS][选项:TSS, CLR, CSS, NONE, TMM] -t TRANSFORM,——TRANSFORM =TRANSFORM应用的转换[Default: LOG][选项:LOG, LOGIT, AST, NONE] -m ANALYSIS_METHOD,——ANALYSIS_METHOD =ANALYSIS_METHOD应用的分析方法[默认:LM][选项:LM, CPLM, NEGBIN, ZINB] -r RANDOM_EFFECTS,——RANDOM_EFFECTS =RANDOM_EFFECTS模型的随机效果,多个效果用逗号分隔[默认值:none] -f FIXED_EFFECTS,——FIXED_EFFECTS =FIXED_EFFECTS模型的固定效果,多个效果用逗号分隔[默认值:all] -c CORRECTION,——CORRECTION =CORRECTION计算q值的校正方法[默认值:none]-o PLOT_SCATTER,——PLOT_SCATTER =PLOT_SCATTER使连续的元数据处于相同的比例上[默认值:TRUE] -l PLOT_HEATMAP,——PLOT_HEATMAP =PLOT_HEATMAP为显著关联生成一个热图[默认值:TRUE] -i HEATMAP_FIRST_N,——HEATMAP_FIRST_N =HEATMAP_FIRST_N在热图中,绘制具有显著关联的前N个特征[默认值:TRUE] -o PLOT_SCATTER,——PLOT_SCATTER =PLOT_SCATTER为显著关联生成一个散点图[默认值:TRUE]-e CORES,——CORES =CORES并行运行的R个进程的数量[默认值:1]-d REFERENCE,——REFERENCE =REFERENCE作为引用的因子,对于有两个以上级别的变量,以'variable, REFERENCE '分号分隔的字符串形式提供[默认值:NA]

故障排除

  1. 问:当我从命令行运行时,我看到了错误Maaslin2。R:命令未找到.我该如何解决这个问题?
    • 答:在运行Maaslin2.R时提供可执行文件的完整路径。
  2. 问:当我作为一个函数运行时,我看到了错误库错误(Maaslin2):没有名为Maaslin2的包.我该如何解决这个问题?
    • 答:安装R包,然后尝试再次加载库。
  3. 问:当我试图安装R包时,我看到有关依赖项没有安装的错误。为什么会这样?
    • 答:安装R包不会自动安装MaAsLin2所需的包。请先安装依赖项,然后再安装MaAsLin2 R包。