I have a question related to DESeq2, just to make sure I am not missing anything. We have RNA-seq data and we generated the count matrix using featureCounts. Then, during our analysis and research, we found interesting data available online and wanted to use it and __compare our data to the data found__ (perform differential expression analysis).
I was checking the DESeq2 vignette and my understanding is that I should include in the __colData__ an additional column, and then include it in the formula, for instance,
~source + condition. This will allow me to compare the conditions while also accounting for the source. Am I right?
Thank you in advance.