mbecVarianceStatsPlot: Plot Proportion of Variance for L(M)M

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

Description

Covariate-Variances as modeled by linear (mixed) models will be displayed as box-plots. It works with the output of 'mbecVarianceStats()' for methods 'lm' and 'lmm'. 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

mbecVarianceStatsPlot(variance.obj)

Arguments

variance.obj,

output of 'mbecVarianceStats' with method lm

Value

A ggplot2 box-plot object.

Examples

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

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