View source: R/fitDistribution.R
Fit a Zero-Inflated Gaussian (ZIG) distribution for each taxon of the count
data. The model estimation procedure is performed by metagenomeSeq
fitZig function without assuming the presence of any group in
the samples (design matrix equal to a column of ones.)
fitZIG(object, assay_name = "counts", verbose = TRUE)
a phyloseq object, a TreeSummarizedExperiment object, or a matrix of counts.
the name of the assay to extract from the
TreeSummarizedExperiment object (default
an optional logical value. If
A data frame containing the continuity corrected logarithms of the
average fitted values for each row of the matrix of counts in the
column, and the estimated probability to observe a zero in the
# Generate some random counts counts = matrix(rnbinom(n = 60, size = 3, prob = 0.5), nrow = 10, ncol = 6) # Fit model on the counts matrix ZIG <- fitZIG(counts) head(ZIG)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.