For a dataset I'm working on, we would like to use the limma package for differential expression analysis. We are using `limma::removeBatchEffect`

(doing it for downstream analysis) before using limma for differential expression analysis. Would it be okay to still include the batch covariate in the model to account for degrees of freedom, even though it's been already corrected for?

Thank you in advance