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

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 10 weeks 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: 171 users visited in the last hour