Hi there, I noticed that the lfcshrink() function in DESeq2 doesn't seem to work with more complex contrast designs, eg: those with interaction terms. Is it possible for this to be implemented? Thanks! :)

This is implemented in the latest release:

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#extended-section-on-shrinkage-estimators

You can use type="apeglm" or type="ashr" with designs with interaction terms.