Question: Interaction with DESeq2 via DiffBind
gravatar for rbronste
23 days ago by
rbronste20 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 16 days ago by Rory Stark2.1k • written 23 days ago by rbronste20
gravatar for Rory Stark
16 days ago by
Rory Stark2.1k
CRUK, Cambridge, UK
Rory Stark2.1k 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 16 days ago by Rory Stark2.1k

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 15 days ago by rbronste20
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: 355 users visited in the last hour