Question: DESeq2: interaction term
gravatar for salamandra
4 months ago by
salamandra0 wrote:
If we look at the example of Deseq2 vignette on interaction terms:
   sampleName         fileName genotype condition group
1       AI_s1 A1_ATCACG_counts        I         A    AI
2       AI_s2 A2_TGACCA_counts        I         A    AI
3      AII_s1 A3_CAGATC_counts       II         A   AII
4      AII_s2 A4_TAGCTT_counts       II         A   AII
5     AIII_s1 B1_CGATGT_counts      III         A  AIII
6     AIII_s2 B2_ACAGTG_counts      III         A  AIII
7       BI_s1 B3_ACTTGA_counts        I         B    BI
8       BI_s2 B4_GGCTAC_counts        I         B    BI
9      BII_s1 C1_TTAGGC_counts       II         B   BII
10     BII_s2 C2_GCCAAT_counts       II         B   BII
11    BIII_s1 C3_GATCAG_counts      III         B  BIII
12    BIII_s2 C4_CTTGTA_counts      III         B  BIII

and we do:

dds3 <- DESeqDataSetFromHTSeqCount(sampleTable = Table, design= ~ genotype + condition + genotype:condition)
dds3 <- dds3[rowSums(counts(dds3)) > 1, ]
ddsHTSeq3 <- DESeq(dds3)

[1] "Intercept"              "genotype_II_vs_I"       "genotype_III_vs_I"      "condition_B_vs_A"       "genotypeII.conditionB"  "genotypeIII.conditionB"
results(ddsHTSeq3, name="genotypeII.conditionB")

"genotypeII.conditionB" gives us the difference between the condition effect for a genotype II and the condition effect for the reference genotype I.

But what if we want to know the difference between the genotype effect for conditionB and the genotype effect for the reference condition A? What should we do?

P.S: this is a repost of

deseq2 • 153 views
ADD COMMENTlink modified 4 months ago • written 4 months ago by salamandra0

I’m out of the office for winter break, but will reply when I’m back.

ADD REPLYlink written 4 months ago by Michael Love23k
Answer: C: DESeq2: interaction term
gravatar for Kevin Blighe
4 months ago by
Kevin Blighe90
Guy's Hospital, London
Kevin Blighe90 wrote:

Hey Michael, I had been assisting the user on Biostars:

On more than 1 occasion, I had advised to read the DESeq2 vignette.

salamandra, if you do post on both Biostars and Bioconductor, can you try to indicate this when you re-post? Many users are on both web-sites.

ADD COMMENTlink modified 4 months ago • written 4 months ago by Kevin Blighe90
Please log in to add an answer.


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