Search
Question: Multi-factor blocking in DiffBind
0
gravatar for jared.andrews07
6 weeks 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 5 weeks ago by Rory Stark2.1k • written 6 weeks ago by jared.andrews070
0
gravatar for Rory Stark
5 weeks ago by
Rory Stark2.1k
CRUK, Cambridge, UK
Rory Stark2.1k 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 5 weeks ago by Rory Stark2.1k
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: 158 users visited in the last hour