Search
Question: Multi-factor blocking in DiffBind
0
gravatar for jared.andrews07
6 months ago by
jared.andrews070 wrote:

I was wondering if it's possible to account for multiple factors when analyzing samples in DiffBind. I want to block for potential batch effects as well as treatment conditions, but I'm unsure how to do so (or if it's even possible). The docs make it clear how to block for one factor, but don't make any mention of multiple. Hoping Rory will pop in here and clarify.

The relevant portions of my samplesheet (I'd like to block for Treatment and Tissue while using Condition for the actual contrast):

                                     
                                     
                                     
                                     
                                     
                                     

 

ADD COMMENTlink modified 6 months ago by Rory Stark2.4k • written 6 months ago by jared.andrews070
0
gravatar for Rory Stark
6 months ago by
Rory Stark2.4k
CRUK, Cambridge, UK
Rory Stark2.4k wrote:

You can specify multiple factors for blocking:

> h3k27ac <- dba.contrast(h3k27ac, categories=DBA_CONDITION, 
                          block=c(DBA_TISSUE,DBA_TREATMENT))

While this may be useful to get a quick look, it probably doesn't do what you really want however. It will create a factor for each unique TISSUE-TREATMENT combination and fit a model for each. It would be more correct to fit a TISSUE model and a TREATMENT model. I'm looking at updating DiffBind to handle factor designs more appropriately, but in the meantime you may want to model this directly in DESeq2 or edgeR.

-Rory

ADD COMMENTlink written 6 months ago by Rory Stark2.4k
Please log in to add an answer.

Help
Access

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