I have a performed a typical differential expression analysis using limma voom and I want to extract the log cpm values to draw ROC curves. But I want to extract the log cpm values produced specifically by limma because they are normalised and free of batch effect that was added to the design matrix. Is there a function to do that?

Thanks in advance

The log-cpm values from voom() are not free of batch effect. If you want batch-corrected log-cpm values, we recommend cpm() followed by removeBatchEffect().