MGLMdtmFit: Fit the Dirichlet-tree multinomial (DTM) distribution.

In liudoubletian/phyloMDA: An R package for phylogeny-aware microbiome data analysis

Fit the Dirichlet-tree multinomial (DTM) distribution.

Description

Fit the Dirichlet-tree multinomial (DTM) distribution.

Usage

MGLMdtmFit(otu.tab, tree)

Arguments

otu.tab	a data frame or matrix containing the count data. Rows of the matrix represent observations and columns are the taxa.
tree	the phylogenetic tree.

Value

Returns a list of MGLM objects ('MGLMfit') corresponding to internal nodes of the tree.

References

Wang, T., Zhao, H.: A Dirichlet-tree multinomial regression model for associating dietary nutrients with gut microorganisms. Biometrics 73(3), 792–801 (2017).

Examples

library(phyloseq)
otu.tab <- t(otu_table(combo.phyloseq.obj))
tree <- phy_tree(combo.phyloseq.obj)
fit <- MGLMdtmFit(otu.tab, tree)
Extract_logL(fit)

liudoubletian/phyloMDA documentation built on April 28, 2022, 6:51 p.m.