Hard coding in DiffBind DBA.R line 968
1
0
Entering edit mode
Ou, Jianhong ★ 1.3k
@ou-jianhong-4539
Last seen 13 hours ago
United States

Hi Rory and Gord,

A student report a potential bug in DiffBind3. At DBA.R line 967-968,

pv.DBA(DBA, method,bTagwise=bTagwise,
           minMembers=3, bParallel=bParallel)

will it be better to be replaced by

pv.DBA(DBA, method,bTagwise=bTagwise,
           minMembers=minMembers, bParallel=bParallel) # minMembers = parameter from pv.contrast

or keep minMembers parameter?

Jianhong.

DiffBind • 665 views
ADD COMMENT
0
Entering edit mode
Rory Stark ★ 5.1k
@rory-stark-5741
Last seen 7 weeks ago
Cambridge, UK

Thank you for noticing and reporting this potential issue relating to the minimum number of sample replicates required to form a contrast.

In this case, the code is correct, as the line appears in dba.analyze(), which does not take minMembers as a parameter. If this function is called on an object with no established contrasts, the defined behavior is to set up the default contrasts (with the default being minMembers=3). If a user needs to automatically establish contrasts using a different value for minMembers, they should call dba.contrast() explicitly, providing the desired value for minMembers, prior to calling dba.analyze().

ADD COMMENT

Login before adding your answer.

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