Diffbind batch effect blocking error
1
0
Entering edit mode
rm1238 • 0
@rm1238-19204
Last seen 5.3 years ago

Hello Rory,

I am experiencing trouble trying to perform a dba with blocking factor (replicate in my case). My current dba object looks like this, and to my eyes, it has the correct format, indicating the Replicate metadata in the corresponding column

Escapeall ID Condition Treatment Replicate Caller Intervals FRiP 1 CtrlD81 Ctrl 0d 1 counts 127618 0.30 2 CtrlD82 Ctrl 0d 2 counts 127618 0.38 7 TreatD81 Treat 8d 1 counts 127618 0.36 8 TreatD8_2 Treat 8d 2 counts 127618 0.36

However, once I try to perform the analysis using the blocking factor, the following error appears:

Escapeallblock <- dba.analyze(Escapeall,block=DBAREPLICATE) Error in dba.analyze(Escapeall, block = DBAREPLICATE) : unused argument (block = DBA_REPLICATE)

I am failing to understand where the error is coming from, as the Replicate column is present in the sample sheet.

I am using R studio 1.2.1335 on a Mac OS 10.14.6

Thank you very much for your help,

Ricardo

Diffbind ATAC batch correction blocking factor • 1.1k views
ADD COMMENT
0
Entering edit mode
Rory Stark ★ 5.2k
@rory-stark-5741
Last seen 6 weeks ago
Cambridge, UK

block is not a parameter to dba.analyze(). If you want to do a blocking contrast, you need to add it explicitly using dba.contrast().

ADD COMMENT
0
Entering edit mode

Thanks a lot for the reply.

Cheers,

Ricardo

ADD REPLY
0
Entering edit mode

Thanks a lot for the reply.

Cheers,

Ricardo

ADD REPLY

Login before adding your answer.

Traffic: 559 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