Search
Question: DiffBind "no lines available in input" error
0
gravatar for igor
16 months ago by
igor20
United States
igor20 wrote:

Sometimes my ChIP-seq samples do not get any significant peaks by MACS, so the output peaks file is empty. When I use that file in DiffBind sample sheet, I get an error:

Error in read.table(fn, blank.lines.skip = T, header = T) : 
  no lines available in input
Calls: dba ... pv.peakset -> pv.readPeaks -> pv.macs -> read.table

Is there a way to suppress that error?

I understand if all the samples have no peaks, there should definitely be an error, but sometimes you have one sample that just doesn't have strong peaks that you still want to include in the analysis.

 

 

ADD COMMENTlink modified 16 months ago by Rory Stark2.1k • written 16 months ago by igor20
2
gravatar for Rory Stark
16 months ago by
Rory Stark2.1k
CRUK, Cambridge, UK
Rory Stark2.1k wrote:

Actually there is a workaround. If you leave the Peaks file empty (blank, no value) in the samplesheet for that sample, it should work! dba.count() will still count reads for this sample in the consensus peaks, and these reads can be used for a differential analysis.

We'll have a look at the error message so we can fix this in a future release.

-Rory

ADD COMMENTlink written 16 months ago by Rory Stark2.1k
1
gravatar for Gord Brown
16 months ago by
Gord Brown540
United Kingdom
Gord Brown540 wrote:

To the best of my knowledge there's no way to suppress that error.  I can see how it might be useful, though, if (for example) one condition completely eliminated binding.  I'll add it to the list of feature requests, but for now I can't see how to bypass it without hacking the source code.  It wouldn't surprise me if, even if you could read the file, other problems would arise downstream, since DiffBind hasn't been tested with empty peak files (that I know of).  Sorry I can't be more helpful.

Wait... just reading the manual page for "dba.peakset", it looks like you can add a peak set that is explicitly empty by setting "peaks=NA".  No idea what might happen downstream, though.

ADD COMMENTlink modified 16 months ago • written 16 months ago by Gord Brown540
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: 131 users visited in the last hour