Hi! I hope you can help me to solve this issue.
I'm using DESeq to look for differentially expressed genes betweem male and female specimes, comparing my samples (trated) with the control (untreated). If I understood correctly the procedure, it should be done like this:
cds <- DESeqDataSetFromMatrix(dataset, sample_metadata, ~ Treatment + Sex + Treatment:Sex)
But this gives me the results for all the samples that I put in the analysis. Since I have 500 samples, I'm looking for a way to get a list (or at least the number) of DE genes for each one of my sample. Is there a way to do this?