This minimal example on ANCOM (Mandal et al. 205) can be run with main.R (in R). This generates the summary in ancom.md