Question: Interaction with DESeq2 via DiffBind
gravatar for rbronste
11 months ago by
rbronste60 wrote:

Hi wondering how I can use the following DESeq2 logic within the framework of DiffBind for calling differential peaks and looking at multiple interactions:

to look at interaction--
dds<-DESeqDataSetFromMatrix(countData = countData,
                           colData = colData,
                           design = ~ sex + treatment + sex:treatment)
dds$sex<-relevel(dds$sex, ref="female")
dds$age<-relevel(dds$treatment, ref="CO")

dds <- DESeq(dds)
results(dds, name="sexmale.treatment")
ADD COMMENTlink modified 11 months ago by Rory Stark2.6k • written 11 months ago by rbronste60
gravatar for Rory Stark
11 months ago by
Rory Stark2.6k
CRUK, Cambridge, UK
Rory Stark2.6k wrote:

Currently DiffBind only handles fairly simple designs. I'm looking to upgrade this in future, but the original idea was that if you wanted to do more complex modelling, you should use the underlying DE package (DESeq2) directly.



ADD COMMENTlink written 11 months ago by Rory Stark2.6k

Thanks Rory! Would it be fair to use the dba.count output matrix and feed that into DESeq2 via DEseqDataSetFromMatrix?

So as to get around having to manually create a read depth adjusted matrix for all of your peaks?

Thanks again.

ADD REPLYlink written 11 months ago by rbronste60
Please log in to add an answer.


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