Using block and design together in Diffbind 3.0?
1
0
Entering edit mode
fleur_p90 • 0
@fleur_p90-10888
Last seen 8 days ago
Netherlands

Hi,

I am performing DiffBind on ATACseq data and I have a question regarding the use of block and design in setting a dba.contrast. From a previous question a couple of years ago I learned how to use block (DiffBind paired samples) but I think it might be different now in the new version as the design setting has changed.

I would like to block for replicate (or individual patient in my case) whilst also using a design that includes Tissue as confounder: ~Tissue + Condition. However, when I do that it says that block is ignored when design is used. Should I just include replicate in my design instead of as block or is there another way to analyze my data in a paired way while also including confounders in the design?

Data <- dba.contrast(Data, block = DBA_REPLICATE, design = "~Tissue + Condition", reorderMeta = list(Condition = "Isolated")) 

Output:

Computing results names...
Warning message:
block ignored when design is used.

Many thanks for any help!

DiffBind R ATACSeq • 89 views
ADD COMMENT
1
Entering edit mode
Rory Stark ★ 4.0k
@rory-stark-5741
Last seen 1 day ago
CRUK, Cambridge, UK

The original block parameter was there to get around DiffBind not supporting arbitrary design formulae. It has been retained for backward compatibility.

The best way to do this now is to put everything (including Replicate) explicitly in the design. So don't use block, and set design="~Replicate + Tissue + Condition".

ADD COMMENT
0
Entering edit mode

I see, okay thanks for the advice!

ADD REPLY
0
Entering edit mode

Great, I was about to ask pretty much the same, so thanks for both question + answer :-)

ADD REPLY

Login before adding your answer.

Traffic: 421 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6