Question: Multi-factor blocking in DiffBind
gravatar for jared.andrews07
12 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 12 months ago by Rory Stark2.6k • written 12 months ago by jared.andrews070
gravatar for Rory Stark
12 months ago by
Rory Stark2.6k
CRUK, Cambridge, UK
Rory Stark2.6k wrote:

You can specify multiple factors for blocking:

> h3k27ac <- dba.contrast(h3k27ac, categories=DBA_CONDITION, 

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.


ADD COMMENTlink written 12 months ago by Rory Stark2.6k
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: 240 users visited in the last hour