Question: Interaction with DESeq2 via DiffBind
gravatar for rbronste
4 months ago by
rbronste50 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 4 months ago by Rory Stark2.3k • written 4 months ago by rbronste50
gravatar for Rory Stark
4 months ago by
Rory Stark2.3k
CRUK, Cambridge, UK
Rory Stark2.3k 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 4 months ago by Rory Stark2.3k

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 4 months ago by rbronste50
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: 345 users visited in the last hour