This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DESeq2 results comparison

I am going through the manual of DESeq2, there are some places I still don't understand:

dds <- DESeqDataSetFromMatrix(countData =countdata, colData = coldata, design = ~ cell +dex)

in the design the cell are cell type and dex contains untreat and treat. in the results extraction part, the example in the manual is results(dds, contrast=c("cell", "N061011", "N61311) So this is the results from differential expression between 2 cell types: N061011 and N61311, no matter treat or untreat? if I want to compare the differential expresssion between treat and untreat in a certain cell type, say, N061011 what should I do? With this, I was guided to the page of results function in detail:here. and it raised me more questions : 1. what is the difference between design = ~ genotype + condition + genotype:condition and design = ~ genotype + condition 2. take 2 genotypes and 2 conditons for example; in the manual it is:

dds <- makeExampleDESeqDataSet(n=100,m=12)
dds$genotype <- factor(rep(rep(c("I","II"),each=3),2))

design(dds) <- ~ genotype + condition + genotype:condition
dds <- DESeq(dds) 
resultsNames(dds)

# the condition effect for genotype I (the main effect)
results(dds, contrast=c("condition","B","A"))

# the condition effect for genotype II
# this is, by definition, the main effect *plus* the interaction term
# (the extra condition effect in genotype II compared to genotype I).
results(dds, list( c("condition_B_vs_A","genotypeII.conditionB") ))

# the interaction term, answering: is the condition effect *different* across genotypes?
results(dds, name="genotypeII.conditionB")

How can I know in my real experiments which is the main effects? why the command for condition effect on genotype I and genotype II is so different? the last command results(dds, name="genotypeII.conditionB), in only mentions genotypeII and conditionB, when I see the comments explained it as the condition effect "differences"across genotypes, I am a little bit confused. So my next question is:how to interpret the meaning of resultsNames(dds)?

Too many questions during this process. Appreciate a lot for your patience and any advice during my learning process.

r rna-seq deseq2

Cross-posted: https://support.bioconductor.org/p/133274/

Note: for such an extensive question, you really should have mentioned all other web-sites where you posted. Otherwise, different users could have used up a substantial amount of their time in answering.

Hi Kevin:

How can I delete this one? It looks that there are too many questions and too long. Maybe delete it is an option?

0 answers

No answers yet.

Log in to answer this question.