mbecRDAStatsPlot: Plot Proportion of Variance for pRDA

In buschlab/MBECS: Evaluation and correction of batch effects in microbiome data-sets

Description

Covariate-Variances as modeled by pRDA will be displayed as box-plots. It works with the output of 'mbecVarianceStats()' for the method 'rda'. Format of this output is a data.frame that contains a column for every model variable and as many rows as there are features (OTUs, Genes, ..). Multiple frames may be used as input by putting them into a list - IF the data.frames contain a column named 'type', this function will use 'facet_grid()' to display side-by-side panels to enable easy comparison.

Usage

mbecRDAStatsPlot(rda.obj)

Arguments

rda.obj,

list or single output of 'mbecVarianceStats' with method rda

Value

A ggplot2 box-plot object.

Examples

# This will return a paneled plot that shows results for three variance
# assessments.
df.var.rda <- mbecModelVariance(input.obj=dummy.mbec,
model.vars=c('group','batch'), method='rda', type='clr')
plot.rda <- mbecRDAStatsPlot(rda.obj=df.var.rda)

buschlab/MBECS documentation built on Jan. 21, 2022, 1:27 a.m.