* for deseq2, use `design = ~ batch + condition` * for limma-voom, may need to do batch correction on the count matrix before running limma