Condition shown as 'Untreated vs. Treated' in output of results() in DESeq2
Entering edit mode
treebig44 ▴ 30
Last seen 19 months ago
United States of America

I am trying to find differentially expressed genes using DESeq2 on some RNA-seq data. In the pheno data, there is a column named 'condition' with factored values of 'Treated' and 'Untreated'. I want to use just this column for my experiment.

However, when I use the results() function in DESeq2 with my data, I get the first line as: 'log2 fold change (MLE): condition Untreated vs Treated'

Now, I am new to this, so I am not sure if this is the same as 'Treated vs Untreated' or not. And I don't know why it decided to show 'Untreated vs Treated'.

I would love some clarification. Are the 2 things (A vs B and B vs A) the same? If not, what does each one mean? and how to control what is compared?

exp_1 = DESeqDataSetFromMatrix(
  countData = Counts,
  colData = Phenotypic[match,],
  design = ~ condition

my_deg = DESeq(exp_1)

fin_results = results(my_deg, alpha=0.05)
RNASeqData DESeq2 • 970 views
Entering edit mode

I have understood that the difference is in the swapping of the numerator and the denominator during the calculation of log2FC. But would still like to know how this affects my results, and how to control this.

Entering edit mode
Last seen 17 hours ago
United States

This is covered in the vignette (note on factor levels).


Login before adding your answer.

Traffic: 711 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6