dada2 (version 1.0.3)

plotComplementarySubstitutions: Plot Substitution Pairs from DADA Result

Description

This is similar to original DADA article, Figure 6.

Usage

plotComplementarySubstitutions(dadaOut, facetByGrp = TRUE)

Arguments

dadaOut

(Required). A dada-class object.

facetByGrp

(Optional). Default TRUE. Whether to plot all substitution groups together in one panel or separately on a grid of panels with a linear model fit.

Value

A ggplot2 object. Will be rendered to default device if printed, or can be stored and further modified. See ggsave for additional options.

Examples

derep1 = derepFastq(system.file("extdata", "sam1F.fastq.gz", package="dada2"), verbose = TRUE)
dada1 <- dada(derep1, err = inflateErr(tperr1, 2), selfConsist = TRUE) 
plotComplementarySubstitutions(dada1)