Question: Bioconductor - DESeq
0
gravatar for soheilazareie
10 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 • 341 views
ADD COMMENTlink modified 10 months ago by JC7.7k • written 10 months ago by soheilazareie10

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

ADD REPLYlink written 10 months ago by JC7.7k
1
gravatar for JC
10 months ago by
JC7.7k
Mexico
JC7.7k 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 10 months ago by JC7.7k
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: 1348 users visited in the last hour