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. An instance of stanfit. integer. number of iterations used in each chain after subtracting warm-up samples. 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.