I am currently using deseq2 for differential abundance analysis. Consider a control-treatment type study where animals are housed in cages/pens and all animals in a given pen receive the same treatment. The only way to account for this using deseq2 is to calculate a median for each OTU and then use cage as the experimental unit. This works but some power is lost. How is this type of analysis handled with gneiss or ancom ?