Model assessment for LDA based on the maximum abundance for each taxon.

Model assessment for LDA based on the maximum abundance for each taxon.

modelAssessment(
  ps,
  stan.fit,
  iterUse = 1000,
  ASVsIndexToPlot = c(1, 3, 10:14, 19:26, 36, 51:53, 148)
)

Arguments

ps

phyloseq object.
stan.fit

An instance of stanfit.
iterUse

integer. number of iterations used in each chain after subtracting warm-up samples.
ASVsIndexToPlot

An integer vector. Use to select ASVs to show the goodness of fit in histograms.

Value

A ggplot2 object. Histogram of data generated from the posterior estimates and the observed data.