Question: RNA-seq with DESeq2 and different conditions
gravatar for martapatriciaribeiroferreira
14 months ago by

Hi everyone!

I am doing a deferentially expressed analysis with DESeq2, but now i have a doubt about using the contrast.

I have my info table of all samples that for example condition A, I have x, y,w,z different variables but I am only interested in analyse the deferentially expressed genes in y vs z, and what I do was using all samples and all condition, but I define the contrast=c("A","y",z") , but when i subset the info table and counts table for only samples and condition that i am interested and again i use the same contrast but this time the only variables is y and z, the results are different, that is in one analysis i have more genes than the other analysis. Someone can tell why this may happen and what I have to use, because I do not understand the differences.

Thanks for your help!

rna-seq R • 325 views
ADD COMMENTlink written 14 months ago by martapatriciaribeiroferreira0
gravatar for russhh
14 months ago by
UK, U. Glasgow
russhh5.4k wrote:

Why did you subset the counts table? Although the additional samples might not contribute to the contrast that you are interested in, they help estimate things like the average expression level and the dispersion-relationship. These things influence the model fitting and any downstream things like contrast calculation.

ADD COMMENTlink written 14 months ago by russhh5.4k

Thanks for your reply! But I thought that aren't important use the samples that are not in my contrast and because I am doing a PCA and heatmap and I only can show in these graphs the samples of these two variables. So you thing that for a set of samples I should use always the complete matrix of counts but defining the contrast with only the variables that I want??

Thanks a lot!

ADD REPLYlink written 14 months ago by martapatriciaribeiroferreira0

Provided all the rest of your samples are of good quality.

ADD REPLYlink written 14 months ago by russhh5.4k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1235 users visited in the last hour