Question: Bioconductor - DESeq
0
gravatar for soheilazareie
21 months ago by
soheilazareie10 wrote:

Hi We have 20 Acinar, 20 Duct, 4 control samples for RNA seq analysis, the aim is to see how Acinar are different from Duct, I have a couple of questions:

  1. Is it necessary to include control samples if just the comparison of Acinar vs Duct (res_A_D <- results(ddsMF, contrast=c("condition","A","D")) is important for us? and also, how important is to have control in our DESeq analysis?
  2. When I do contrast, my codes are

    res_A_WT <- results(ddsMF, contrast=c("condition","A","WT"))

    res_D_WT <- results(ddsMF, contrast=c("condition","D","WT"))

    res_A_D <- results(ddsMF, contrast=c("condition","A","D"))

My question is: why I get the same list of genes when I do the above contrast?, is it correct? if so, what is the point of doing contrast?

  1. For making a heatmap of Differentially expressed genes, which one is more important, fold change or p-adjusted value?

Thanks, Soheila

rna-seq • 547 views
ADD COMMENTlink modified 21 months ago by JC9.5k • written 21 months ago by soheilazareie10

Hi Soheila, I just formatted your post and tagged as Question, because this is not a "Tutorial".

ADD REPLYlink written 21 months ago by JC9.5k
1
gravatar for JC
21 months ago by
JC9.5k
Mexico
JC9.5k wrote:

Answering your 3 questions:

  1. Yes, you can compare Acinar (A), Duct (D) and controls (C) but I would contrast (A / C) vs (D / C) and then get DEG between groups, that will mean to get genes differentially expressed to the control and differential expressed between A / D. But if you only care about the differences between A / D is fine just to compare both without consider C, just keep in mind that both are similar questions but not necessary are the same.

  2. Are you sure the lists are the same? I would expect different genes in each class. I guess you are observing the same genes when you compare ( A/C vs D/C ) and (A/D), that can be explained if the C doesn't have an impact in the DEG.

  3. Depends, if you want to show how strong are the changes, use fold-change, if you want to show how statistically significant are, use the p-value.

ADD COMMENTlink written 21 months ago by JC9.5k
Please log in to add an answer.

Help
Access

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